| 查看: 1084 | 回复: 3 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
wbjhpp金虫 (初入文坛)
|
[交流]
【求助】用mathematic编个二分法的程序 已有2人参与
|
||
| 请用mathematic编个二分法的程序，谢谢喽！如果方程类似反函数的话，应该怎么排除这种情况？ |
回复此楼» 猜你喜欢
遇见不省心的家人很难过
已经有10人回复
-
博士延得我，科研能力直往上蹿
已经有3人回复
-
退学或坚持读
已经有24人回复
-
免疫学博士有名额，速联系
已经有14人回复
-
面上基金申报没有其他的参与者成吗
已经有4人回复
-
多组分精馏求助
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 编辑意见回复
已经有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人回复
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
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