二階阮奇庫塔法微分方程計算

Second order Runge-Kutta Method

程式編寫日期: 2007年5月4日 修改日期: 2008年8月3日

注意: πr 是按 Shift EXP Shift Ans 2

程式 (70 bytes,不包括綠色的函數方程 )

Mem clear: ?→X: ?→Y: ?→A: ?→M:

M-1(A - X→B: Lbl 0: Y + . 5 BC: 1 + 2X + Ans→C:

X + . 5B→X: sin Mπr => Y + BC→Y: . 5M-:

M => Goto 0: Y

注意: 若果想計算其它微分方程,只要修改綠色的部份即可(變數x使用記憶X,變數y使用記憶Ans)。

 

例題: 利用二階阮奇庫塔法計算微分方程 dy/dx = 1 + 2x + y, y(1)=2的解在 x =3 的近似值,其中區間數目為12。

按 Prog 1再按 1 EXE 2 EXE 3 EXE 12 EXE (顯示 y(3)的近似值為42.74692888)

計算完成後,按 RCL Y 顯示近似值答案。

 

返回 fx-3650P及SC185程式集

Free Web Hosting