%Fermi Function Calculation, f(ΔE, T) % Constant k=8.617e-5; %Computation proper for ii=l: 4; T= l00*ii; kT=k*T; dE(ii,1)=-5*kT; for ij=1: 101; f(ii,jj)= 1/(1 +exp(dE(ii,jj)/kT)); dE(ii,jj+1)=dE(ii,jj)+(0.1 *kT) end end dE=dE(1:4,1 :101); %This step strips the extra dE value %plotting results close plot(dE,f);grid; xlabel(E-EF(eV)); text(0.05,0.2,'T= 400K'); text(-.03,0.1,'T=100K');
Best Answer