me = 9.109*10^-31; mh = me; me2 = me; k = 1.38*10^-23; T = 300; 
    a = 5.43*10^-10; em = 1.602*10^-19; Eg = 1*em; Ea = 0.2*em; 
    h = 6.626*10^-34; NC = 2*(2*π*me2*k*T/(h^2))^(3/2); 
    NV = 2*(2*π*mh*k*T/(h^2))^(3/2); NA := 10^24; 
    mup = 10^(-2);

ni = Sqrt[NC*NV]*Exp[-Eg/2/k/T]

p = 2*NA*1/(1 + Sqrt[1 + 4*NA/NC*Exp[Ea/k/T]])

p/NA*100

mup*p*em