| 查看: 772 | 回复: 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); |
回复此楼» 猜你喜欢
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有250人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有6人回复
-
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有22人回复
-
考博自荐
已经有4人回复
-
山东大学第二批博士研究生招生
已经有0人回复
-
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
2楼2018-07-18 19:54:04
3楼2018-07-18 20:58:15