| 查看: 755 | 回复: 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); |
回复此楼» 猜你喜欢
［调剂信息］211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有236人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有2人回复
-
070200求调剂，一志愿某211，288分
已经有18人回复
-
法国博士后职位
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
广州大学光电信息工程专业调剂，招收物理学专业学生
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
2楼2018-07-18 19:54:04
3楼2018-07-18 20:58:15
5