| 查看: 741 | 回复: 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); |
回复此楼» 猜你喜欢
深圳大学2026年秋博士招生-物理学-活性胶体方向-高永祥课题组
已经有18人回复
-
论物质与能量的统一模型及物理现象解释
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有228人回复
-
基于基元I统一理论的数学相关应用推导
已经有0人回复
-
基元I统一理论：宇宙本质、层级演化与修炼文明的本源规律
已经有1人回复
-
基元I理论下三大核心空间现象精准推导与细节解析
已经有0人回复
-
基于基元 I 统一理论的反重力理论推导
已经有0人回复
-
基于基元I统一理论的量子力学本源推导
已经有1人回复
-
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，AI润色，全免费
已经有20人回复
-
【EI|Scopus 双检索】第六届智能机器人系统国际会议（ISoIRS 2026）
已经有0人回复
-
2026年第四届电动车与车辆工程国际会议（CEVVE 2026）
已经有0人回复
3楼2018-07-18 20:58:15
2楼2018-07-18 19:54:04