| 查看: 1619 | 回复: 7 | ||
| 【悬赏金币】回答本帖问题，作者问天卖饼2将赠送您 100 个金币 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
振荡的波函数如何求积分 已有1人参与
|
||
我用薛定谔方程分别求出199个phi和199个psi，现在的问题是phi*psi后求积分求不出来，代码如图所示。我尝试了很多种method都不行。求助广大强大的网友，该怎么求积分呢？ @月只蓝 @月只蓝 @beefly 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
散金币祈福
已经有186人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有182人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有114人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有70人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有39人回复
4楼2020-11-12 10:17:43
2楼2020-11-12 10:09:41
3楼2020-11-12 10:13:41
|
L = 1000*10^-10; F = 10^7; m = 0.0665*9.1*10^-31; mhole = 0.34*9.1*10^-31; e = 1.6*10^-19; \[HBar] = 1.05*10^-34; Ne = 199;eqn1 = (-(\[HBar]^2/(2 m))* (\[Phi]^\[Prime]\[Prime])[x] + e*F*x*\[Phi][x]);{phivals, phifuns} = NDEigensystem[{eqn1, DirichletCondition[\[Phi][x] == 0, x \[GreaterSlantEqual] L/2 || x <= -L/2]}, \[Phi][x], {x, -L/2, L/2}, Ne, Method -> {"SpatialDiscretization" -> {"FiniteElement", \ {"MeshOptions" -> {MaxCellMeasure -> 10^-9}}}}];eqn2 = (-(\[HBar]^2/(2 mhole))* (\[Psi]^\[Prime]\[Prime])[x] - e*F*x*\[Psi][x]);{psivals, psifuns} = NDEigensystem[{eqn2, DirichletCondition[[Psi][x] == 0, x [GreaterSlantEqual] L/2 || x <= -L/2]}, [Psi][x], {x, -L/2, L/2}, Ne, Method -> {"SpatialDiscretization" -> {"FiniteElement", \ {"MeshOptions" -> {MaxCellMeasure -> 10^-9}}}}];Etable = Table[(psivals[] + phivals[[j]])/(1.6*10^-22), {i, Ne}, {j, Ne}];Itable = Table[ Abs[NIntegrate[psifuns[]*phifuns[[j]], {x, -L/2, L/2}, Method -> {"ExtrapolatingOscillatory"}]], {i, Ne}, {j, Ne}]; |
5楼2020-11-12 11:06:59