| 查看: 1096 | 回复: 3 | |||
wbjhpp金虫 (初入文坛)
|
[交流]
【求助】用mathematic编个二分法的程序 已有2人参与
|
| 请用mathematic编个二分法的程序，谢谢喽！如果方程类似反函数的话，应该怎么排除这种情况？ |
回复此楼» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 编辑意见回复
已经有6人回复
- 【求助】给Journal of computer 发了3个评审员的联系方式，主编给回复
已经有7人回复
- 由于语法问题，编辑要求回稿修改，并要告诉他怎么改的，请问回信该放coverletter里吗
已经有30人回复
- 【软件】Ginger 世界一流的英语语法和拼写检查程序
已经有1213人回复
- 在科研过程中如何验证自己编译程序的正确性
已经有3人回复
- 虫友们，请问如何使用endnote编写aip期刊的参考文献？
已经有6人回复
- 求助一段代码解读，有奖励
已经有3人回复
- 用c语言编程怎来表达时滞！谢谢！
已经有6人回复
- 招聘石油化学，有机化学方面的兼职编辑
已经有21人回复
- matlab 积分计算 算不出结果，求帮助！！！
已经有10人回复
- 如何编写一个应用程序
已经有6人回复
- 一个从来没有编程经验的人，你们觉得VB作为入门语言怎么样？
已经有23人回复
- 【求助】帮我编个程啊~C语言，C++，Matlab，VB...都行~
已经有14人回复
- 【求助】帮忙编个用于VASP小程序
已经有22人回复
mshwangg
至尊木虫 (正式写手)
- 程序强帖: 5
- 应助: 206 (大学生)
- 金币: 10702.8
- 散金: 100
- 红花: 19
- 帖子: 597
- 在线: 195.4小时
- 虫号: 576702
- 注册: 2008-06-21
- 专业: 物理学I
2楼2010-11-05 13:21:15
wbjhpp
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 935.4
- 散金: 16
- 帖子: 28
- 在线: 7小时
- 虫号: 1114459
- 注册: 2010-10-05
- 专业: 凝聚态物性 II ：电子结构
具体问题
nono2009:建议通使用“引用回复该帖”，不然别人收不到你的回复提示。 2010-11-08 12:44:39
|
h[x]是薛定谔方程的波函数（不知道具体形式），x表示能量，下面是我编的找出h[x]==0的解的二分法程序，v1、v2是已知的势能；但如果此波函数是类似反函数的情况，也满足h[p1]*h[p3]<0，但它与x轴没交点，也就是说并没有解，那怎么用程序排除这种情况！！谢谢 total = 1000; y1 = v1;y2 = v2; y0 = (y2 - y1)/total; rt[x_, y_] := Module[{p1 = x, p2 = y, p3 = (x + y)/2}, While[Abs[p1 - p2] > 0.0001, If[h[p1]*h[p2] == 0, If[h[p1] == 0, p3 = p1, p3 = p2]; Break[]]; If[Negative[h[p1]*h[p3]], p2 = p3; p3 = (p1 + p2)/2, p1 = p3; p3 = (p1 + p2)/2]]; p3]; last = Table[ If[h[y1 + w*y0]*h[y1 + (w + 1) y0] < 0, rt[y1 + w*y0, y1 + (w + 1) y0], 0], {w, 1, total - 2}]; sz = Module[{i = 0}, Do[If[Part[last, j] != 0, i = i + 1, i], {j, 1, Length[last]}]; i]; netlast = Table[0, {i, 1, sz}]; Module[{j = 1}, Do[If[Part[last, i] != 0, Part[netlast, j] = Part[last, i]; j = j + 1], {i, 1, Length[last]}]]; netlast |
3楼2010-11-08 11:16:20
wbjhpp
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 935.4
- 散金: 16
- 帖子: 28
- 在线: 7小时
- 虫号: 1114459
- 注册: 2010-10-05
- 专业: 凝聚态物性 II ：电子结构
具体问题
|
h[x]是薛定谔方程的波函数（不知道具体形式），x表示能量，下面是我编的找出h[x]==0的解的二分法程序，v1、v2是已知的势能；但如果此波函数是类似反函数的情况，也满足h[p1]*h[p3]<0，但它与x轴没交点，也就是说并没有解，那怎么用程序排除这种情况！！谢谢 total = 1000; y1 = v1;y2 = v2; y0 = (y2 - y1)/total; rt[x_, y_] := Module[{p1 = x, p2 = y, p3 = (x + y)/2}, While[Abs[p1 - p2] > 0.0001, If[h[p1]*h[p2] == 0, If[h[p1] == 0, p3 = p1, p3 = p2]; Break[]]; If[Negative[h[p1]*h[p3]], p2 = p3; p3 = (p1 + p2)/2, p1 = p3; p3 = (p1 + p2)/2]]; p3]; last = Table[ If[h[y1 + w*y0]*h[y1 + (w + 1) y0] < 0, rt[y1 + w*y0, y1 + (w + 1) y0], 0], {w, 1, total - 2}]; sz = Module[{i = 0}, Do[If[Part[last, j] != 0, i = i + 1, i], {j, 1, Length[last]}]; i]; netlast = Table[0, {i, 1, sz}]; Module[{j = 1}, Do[If[Part[last, i] != 0, Part[netlast, j] = Part[last, i]; j = j + 1], {i, 1, Length[last]}]]; netlast |
4楼2010-11-09 08:44:58
