файлы матлаб / parameters

ft=kB*T/q;

Eg0=1.1695;
Alfa=4.73e-4;
Betta=636;
Eg=Eg0-Alfa*T^2/(Betta+T);

mnt=0.19;
mnl=0.98;
m1n=0.009;
mpl=0.16;
mph=0.49;
m1p=0.6;
m2p=-0.1;
m0n=(mnt^2*mnl)^(1/3);
m0p=(mpl^1.5+mph^1.5)^(2/3);
mn=m0n+m1n*(T/300);
mp=m0p+m1p*(T/300)+m2p*(T/300)^2;

Nc=2*(2*pi*mn*me*kB*T/h^2)^1.5;
Nv=2*(2*pi*mp*me*kB*T/h^2)^1.5;
ni=sqrt(Nc*Nv)*exp(-Eg/2/ft);

muL300n=1430*1e-4;
muL300p=460*1e-4;
g0n=-2;
g0p=-2.18;
g1n=-0.45;
g1p=-0.45;
g3n=3.2;
g3p=3.2;
g4n=0.065;
g4p=0.065;
muMin300n=80*1e-4;
muMin300p=45*1e-4;
Cref300n=1.12e17*1e6;
Cref300p=2.23e17*1e6;
a300n=0.72;
a300p=0.72;
muLn=muL300n*(T/300)^g0n;
muLp=muL300p*(T/300)^g0p;
muMinn=muMin300n*(T/300)^g1n;
muMinp=muMin300p*(T/300)^g1p;
Crefn=Cref300n*(T/300)^g3n;
Crefp=Cref300p*(T/300)^g3p;
an=a300n*(T/300)^g4n;
ap=a300p*(T/300)^g4p;
MUn=muMinn+(muLn-muMinn)./(1+((ND+NA)./Crefn).^an);
MUp=muMinp+(muLp-muMinp)./(1+((ND+NA)./Crefp).^ap);