| 查看: 812 | 回复: 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); |
回复此楼» 猜你喜欢
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有673人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
物理学I论文润色/翻译怎么收费?
已经有62人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有4人回复
2楼2018-07-18 19:54:04
3楼2018-07-18 20:58:15
