Casio50 » fx-3x50P » Programs 程式集 » Differential Equations 微分方程 »

Second Order Runge-Kutta Method
二階 Runge-Kutta 計算法

本式可以二階 Runge-Kutta 計算法計算微分方程

<COMP>
Goto 3 : Lbl 0 : f(X,Ans) : Goto 2 :
Lbl 3 : MemClear : ? → X : ? → Y : ? → B : ? → M :
M -1 ( B – X → D :
Lbl 1 : Y + C : Goto 0 : Lbl 2 :
. 5 D Ans → C : X + . 5 D → X :
A => Y + 2 C → Y : A = 0 → A :
B > X => Goto 1 : Y
>91 Bytes

操作方法

一般操作 例子
解
y(xn)
其中:
dy/dx = f(x,y); y(x0) = y0 		解
y(5)
其中:
dy/dx = ln|y^2-x|; y(2)=1
整體程式為:
<COMP>
Goto 3 : Lbl 0 : ln √ ( Ans ² – X ) ² : Goto 2 :
Lbl 3 : MemClear : ? → X : ? → Y : ? → B : ? → M :
M -1 ( B – X → D :
Lbl 1 : Y + C : Goto 0 : Lbl 2 :
. 5 D Ans → C : X + . 5 D → X :
A => Y + 2 C → Y : A = 0 → A :
B > X => Goto 1 : Y
100 Bytes

註:
Ans 是代表 y(x) 的變量而 X 則代表 x.
√x² 是指取 x 的絕對值.
啟動程式 按 Prog 1
輸入數據
x0 EXE y0 EXE xn EXE n EXE 		2 EXE 1 EXE 5 EXE 36 EXE (分成36份)
顯示答案
y(xn)		 1.91980773
(正確答案 = 1.9182575899)

記憶體

A 1
B xn
C K2/2
D h
X xn
Y y(xn)
M n

相關程式

  1. Euler’s Method / 歐拉計算法
  2. Fourth Order Runge-Kutta Method / 四階 Runge-Kutta 計算法