| 查看: 708 | 回复: 2 | ||
[求助]
1stopt，求积分，为什么每次求的值不一样 已有1人参与
|
|
f3和f4每次计算值不一样，有时0有时正有时负，或者帮我看看是我模型的问题，还是程序的问题？求大佬帮帮我。 Constant q=1.6e-19, mox=0.42*9.11e-31, h=6.63e-34, d=7.4e-7, mm=0.1*9.11e-31, F1=0.2, ss=5.2, s0=8.85E-14, Tem=298, K=1.38E-23, pi=3.1415926, Fb0=2.5,Ft=1.48,Nt=1.5e19, E=1.5e6, //t=2.2, A=8*pi*sqrt(2*q*mox)/(3*h), Fb=Fb0-sqrt((q*E)/(4*pi*ss*s0)), Ct=(mm/mox)^(5/2)*(16*pi*q*F1^(3/2))/(3*h*sqrt(Ft-F1)), C=(q*Ct*Nt)/E; conststr P1=exp((-A/E)*(t^(3/2)-Ft^(3/2))), //P2_tri=exp((-A/E)*Ft^(3/2)), //P2_tra=exp((-A/E)*(Ft^(3/2)-(t-E*d)^(3/2))), fFD=1/(1+exp(q*(Fb-t)/(K*Tem))); C1=1E160,C2=1E-150; FUNCTION //f0=A, //f1=C2, //f2=int(exp(x^(3/2)),x=1,2), f3=(8*pi*sqrt(2*q*mox)/(3*h))*int(1/(1/(1/(1+exp(q*(Fb-t)/(K*Tem)))*exp((-A/E)*(t^(3/2)-Ft^(3/2))))+1/exp((-A/E)*(Ft^(3/2)-(t-E*d)^(3/2)))),t=Ft,2.5), //f=1e200*Int((1/(fFD*P1)+1/P2_tra)^(-1),t=Ft,Ft+E*d); f4=Int((1/((1/(1+exp(q*(Fb-t)/(K*Tem))))*(exp((-A/E)*(t^(3/2)-Ft^(3/2)))))+1/(exp((-A/E)*(Ft^(3/2)-(t-E*d)^(3/2)))))^(-1),t=Ft,Ft+E*d); |
回复此楼» 猜你喜欢
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有87人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
-
求助NH4V4O10晶体的CIF文件
已经有0人回复
2楼2018-07-18 19:54:04
3楼2018-07-18 20:58:15