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

[nicky729 (御風)]

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)
推 of:評翰太威了阿阿阿阿阿阿!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
→ nicky729:還好啦...CJ應該也寫得出來，他不想寫而已
推 of:我寫不出來阿 哭哭
推 icarusming:大感謝~~!!
推 nicky729:不客氣
推 llewxam:評翰太厲害了 <(_ _)>
推 cjcat2266:酷喔~        (不過你還有閒情逸致做花俏的按鈕呀...XD)
推 nicky729:回樓上～我對自己的作品有底限要求 XD
推 jerrypotato:有用有推XD
推 sobluesk:推~~評翰太威了!!!
→ nicky729:不敢...學長說..微方課本后面的CD有法寶，不知道是真的嗎