Re: [請益] 微方3.3第16題

看板NTUEE110HW作者nicky729 (御風)時間18年前 (2007/10/31 01:32)推噓9(9推 0噓 2→)

留言11則, 7人參與討論串4/5 (看更多)

http://www2.ee.ntu.edu.tw/~b95901189/DE/ 如果大家不介意就用我寫的小程式算吧... 要留意的事項如下： 1. 沒有檢查下限，所以s,i,r,n是否為負數都沒錯 2. r的值 = n - s[t] - i[t] 3. dt 為　每次代值取的間隔 4. 由於懶惰問題，當n>10時，描點將畫出界限外 5. 此程式的算法為 s[t+dt] = s[t] + s'[t]*dt, i如前所述，r如2所計算 6. 如果有bug，請別怪我.... <在林清富的作業海的同時，完成電網導的PA2...還要弄這個....waler...不要太計較了> <本來不想公開的...希望大家開開心心地用唄，不再用手描....> ※ 引述《icarusming (icarus)》之銘言： : 試了很久還是試不出來>"< : NDSolve[{s'[x]==-0.2i[x]s[x],i'[x]==-0.7i[x]+0.2s[x]i[x],i[0]==8,s[0]==2} : ,{i[x],s[x]},x] : 答案還是無法算出來，會變成把後面兩個initial condition 代入前面再印出來>"< : 請問該如何改呢? : ※ 引述《thewindjuei (well...)》之銘言： : ※ 引述《shinchen (starrydawn星晨)》之銘言： : : b小題到底要我們幹麼啊@@ : : 不可能要手動解啊 : : 但是我手邊的程式又不吃這種聯立方程式 : : (應該說我不知道要怎麼說服程式s跟i是互相獨立的= =) : : 解答卻碰碰碰的就跑出函式圖形了 : : 怎麼做出來的 : : 有人知道嗎@@ : 參考下面網站 : http://www.math.sjtu.edu.cn/Mathematica%E6%95%99%E7%A8%8B/6.1.1.htm : DSolve[{s'[t] == -0.2×s[t]×i[t], i'[t] == -0.7×i[t] + 0.2×s[t]×i[t], : i[0] == 8, s[0] == 2}, {i[t], s[t]}, t] : 可是mathmatica說找不到符合此條件的解， : 這是我輸入錯誤嗎？還是那邊出問題了呢？ : 歡迎多多討論~ -- -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.112.249.18 ※ 編輯: nicky729 來自: 140.112.249.18 (10/31 01:33) ※ 編輯: nicky729 來自: 140.112.249.18 (10/31 01:34)

推

10/31 01:40, , 1^F

10/31 01:40, 1^F

→

nicky729

10/31 01:47, , 2^F

10/31 01:47, 2^F

推

10/31 02:49, , 3^F

10/31 02:49, 3^F

推

icarusming