四階阮奇庫塔法微分方程計算
Fourth order Runge-Kutta Method
程式編寫日期: 2008年3月2日 修改日期: 2008年8月2日
程式 (114 bytes，不包括綠色的函數方程)
Mem clear: ?→X: ?→A: ?→B: ?→Y:
Y-1(B - X→B: Lbl 0: A→C: Lbl 1: 1 + 2X + Ans→D:
A + BD ÷ ( 3 + 3(M2=3M→A: M2=2M => X + B ÷ 2→X:
1M+: M=4 => Goto 2: C + DB ÷ (2 - (M=3: Goto 1: Lbl 2:
MM-: Y - 1→Y => Goto 0: A
注意: 若果想計算其它微分方程，只要修改綠色的部份即可(變數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.72281362)
計算完成後，按 RCL A 顯示近似值答案。