disp("Ef-Ec = K*T*((log(n/Nc))+((1/sqrt(8))*(n/Nc)))");
a=4.4;  //say n/Nc=a
Nc=2.78*10^19;
n=a*Nc;
a1=4.51;
n1=a*a1*10^17;
printf('\n For GaAs, the density for degeneracy is %f',n1);