| 4 | 1/1 | 返回列表 |
| 查看: 1106 | 回復(fù): 3 | |||
wbjhpp金蟲 (初入文壇)
|
[交流]
【求助】用mathematic編個二分法的程序 已有2人參與
|
| 請用mathematic編個二分法的程序,謝謝嘍!如果方程類似反函數(shù)的話,應(yīng)該怎么排除這種情況? |
至尊木蟲 (正式寫手)
金蟲 (初入文壇)
|
h[x]是薛定諤方程的波函數(shù)(不知道具體形式),x表示能量,下面是我編的找出h[x]==0的解的二分法程序,v1、v2是已知的勢能;但如果此波函數(shù)是類似反函數(shù)的情況,也滿足h[p1]*h[p3]<0,但它與x軸沒交點,也就是說并沒有解,那怎么用程序排除這種情況。≈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 |
金蟲 (初入文壇)
|
h[x]是薛定諤方程的波函數(shù)(不知道具體形式),x表示能量,下面是我編的找出h[x]==0的解的二分法程序,v1、v2是已知的勢能;但如果此波函數(shù)是類似反函數(shù)的情況,也滿足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 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +13 | Mr. Z 2026-03-25 | 13/650 |
|
|---|---|---|---|---|
|
[考研] 291求調(diào)劑 +5 | Y-cap 2026-03-29 | 6/300 |
|
|
[考研] 一志愿鄭州大學(xué),080500學(xué)碩,總分317分求調(diào)劑 +8 | 舉個栗子oi 2026-03-24 | 9/450 |
|
|
[基金申請] 面上5B能上會嗎? +3 | redcom 2026-03-29 | 3/150 |
|
|
[考研] 343求調(diào)劑 +6 | 愛羈絆 2026-03-29 | 6/300 |
|
|
[考研] 311求調(diào)劑 +5 | 冬十三 2026-03-24 | 5/250 |
|
|
[碩博家園] 求調(diào)劑 有機化學(xué)考研356分 +3 | Nadiums 2026-03-25 | 4/200 |
|
|
[考研] 085600,材料與化工321分求調(diào)劑 +9 | 大饞小子 2026-03-28 | 9/450 |
|
|
[考研] 085701環(huán)境工程,267求調(diào)劑 +16 | minht 2026-03-26 | 16/800 |
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
[考研] 266求調(diào)劑 +11 | 陽陽哇塞 2026-03-27 | 12/600 |
|
|
[考研] 考研化學(xué)308分求調(diào)劑 +10 | 你好明天你好 2026-03-23 | 12/600 |
|
|
[考研] 0703化學(xué)338求調(diào)劑! +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
[考研] 中國科學(xué)院深圳先進技術(shù)研究院-光纖傳感課題組招生-中國科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|