請問各位大神，下面兩個matlab程式，哪里出問題了，運行好久，沒結果
程式1
T=0:1:100;
fplot(T,Ce)
function C=Ce(T)
syms T E
h=6.6261*10^-34; c=3*10^8; me=9.10956*10^-31; ne=8.5*10^28; kb=1.38065*10^-23;
f=c/800*10^-9;
w=2*pi*f;
Ef=(h*c)/(8*me*c^2)*(3/pi)^(2/3)*ne^(2/3);
U=Ef*[1-(pi^2/12)*(kb*T/Ef).^2+(pi^2/80)*(kb*T/Ef).^4];
B=1./(kb*T);
np=1./(exp(B.*(E-U))+1);
P=(8*sqrt(2)*pi*me^(3/2)*sqrt(E))/(h^3);
f1=np.*P.*E;
f2=np.*P;
C=ne*diff(Ep,T);
end
程式2
T=0:1:100;
fplot(T,F1)
function F=F1(T)
h=6.6261*10^-34; c=3*10^8; me=9.10956*10^-31; ne=8.5*10^28; kb=1.38065*10^-23;
f=c/800*10^-9;e=1.6*10^-19; Z=1;
syms T t
Ef=(h*c)/(8*me*c^2)*(3/pi)^(2/3)*ne^(2/3);
U=Ef*[1-(pi^2/12)*(kb*T/Ef).^2+(pi^2/80)*(kb*T/Ef).^4];
F=int(t^0.5/(1+exp(t-U/(kb*T))),t,0,inf);
end
轉載請註明出處，本文鏈接：https://www.uj5u.com/houduan/58139.html
標籤：C++ 語言