24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

wbjhpp

金虫 (初入文坛)

应助: 0 (幼儿园)
金币: 935.4
散金: 16
帖子: 28
在线: 7小时
虫号: 1114459
注册: 2010-10-05
专业: 凝聚态物性 II ：电子结构

[交流] 【求助】用mathematic编个二分法的程序已有2人参与

请用mathematic编个二分法的程序，谢谢喽！如果方程类似反函数的话，应该怎么排除这种情况？

回复此楼

» 猜你喜欢

一个化合物的合成路线：CAS:367929-02-0 名称：8β-乙烯基雌二醇已经有4人回复
太白金星有点烦已经有3人回复
中国地质大学（北京）博士招生补录，数理学院材料科学与工程专业和材料与化工专业已经有3人回复
河北省自然基金已经有8人回复
收到国自然专家邀请后几年才会有本子送过来评已经有3人回复
考博已经有3人回复
有没有快的中文核心比较快录用的，纳米材料光催化已经有4人回复
本人42，博士刚毕业，现在找不到工作，怎么办？:( 已经有21人回复
有人投过CCC中国控制会议吗？已经有3人回复
3,4-二羟基苯乙酮如何纯化？已经有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人回复

1楼 2010-11-04 20:24:41

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

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