function [ TIA, TIS, stFull, stHalf, stEmpty, stAll ] = SpinStatistic( EkT,x )
% [ TIA TIS stFull stHalf stEmpty stAll ] = SpinStatistic( EkT,x ); 
% x is spin polarization  EkT  is electron energy /kT
% TIA -number of spin-polarized electrons
% TIS -number of spin-unpolarized electrons
% stFull - number of electrons in full-filled states
% stHalf -number of electrons in half-filled states (spin-polarized +
% spin-unpolarized)
%stEmpty - number of states, which are not filled
% stAll - number of all electrons (The Fermi-Dirac distribution)
% use for gragh EkT=-5:0.05:5; EkT=EkT';
% 
%   

stAll=zeros(length(EkT),1);
stEmpty=zeros(length(EkT),1);
TIA=zeros(length(EkT),length(x));
TIS=zeros(length(EkT),length(x));
stFull=zeros(length(EkT),length(x));
stHalf=zeros(length(EkT),length(x));


for jE=1:length(EkT)
    EkT1=EkT(jE);
    stAll(jE)=1/(1+exp(EkT1));
    stEmpty(jE)=1-stAll(jE);
    for jx=1:length(x)
        SP=x(jx);
temp=(1-sqrt(1-2*(1-0.3832*(1-SP)^2)/(1+cosh(EkT1))))/2/(1-0.3832*(1-SP)^2);
stHalf(jE,jx)=temp;
TIA(jE,jx)=SP*temp;
TIS(jE,jx)=(1-SP)*temp;
stFull(jE,jx)=stAll(jE)-temp;

    end
end
end