[問題] findroot後繼續運算

看板Mathematica作者reforevering (nick)時間9年前 (2015/07/17 00:01)推噓1(1推 0噓 1→)

留言2則, 2人參與討論串1/1

各位前輩大家好,我現在的問題是這樣子我需要解類似的方程式好幾次首先會假設一個初始條件的數字F0 然後去找到一個解X1,使得下列的2個方程式成立 (1)F1=F0-N1 (2) PC0*J0= J1 --> 第一組 N1是X1的函數, PC0是F0的函數, J0(or J1)是跟F0(or F1)有關的積分式(NIntegrate) 找到X1之後就繼續做同樣的事情解 F2=F1-N2 & PC1*J0=J2 找到解X2 --> 第二組到Fn=Fn-1 - Nn & PCn*J0=Jn 找到解Xn --> 第n組解了n組之後看 F0+F1+...+Fn有沒有等於我想的值,沒有的話就重新假設F0,然後再重複繼續解直到找到一個正確的F0使得F0+F1+...+Fn等於我想要的值 -------------------------------------------------------------------- 我現在可以用Findroot找到X1 但我不知道怎麼讓程式知道要解第2組的時候的F1已經等於F0-N1(帶入X1得到的結果) 就是把上一組的解(X1)帶入得到相關參數的值(F1 / PC1)之後再進行下一組的計算有什麼相關指令可以做到這件事情嗎? (假如沒爬到文請見諒....) 另外我這整組解包含最後F0+F1+...+Fn到重新輸入F0應該是可以用迴圈相關指令來完成? 不過我對這方面不太熟悉...... 不知後有沒有前輩可以給個提示麻煩大家了謝謝 <__ __> -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.113.55.150 ※ 文章網址: https://www.ptt.cc/bbs/Mathematica/M.1437062502.A.D4B.html