Casio⁵⁰ » fx-3x50P » Programs 程式集 » Differential Equations 微分方程 »

Euler’s Method
歐拉計算法

本式可以歐拉計算法計算微分方程

<COMP>
Goto 3 : Lbl 0 : f(X,Y) : Goto 2 :
Lbl 3 : ? → X : ? → Y : ? → B : ? → M :
M ^-1 ( B – X → D :
Lbl 1 : Goto 0 : Lbl 2 :
Y + D Ans → Y : X + D → X :
B > X => Goto 1 : Y
>68 Bytes

操作方法

一般操作

例子

解
y(x_n)
其中:
dy/dx = f(x,y); y(x0) = y0

dy/dx = f(x,y); y(x0) = y0

解
y(5)
其中:
dy/dx = ln|y^2-x|; y(2)=1

dy/dx = ln|y^2-x|; y(2)=1

整體程式為:

<COMP>
Goto 3 : Lbl 0 : ln √ ( Y ² – X ) ² : Goto 2 :
Lbl 3 : ? → X : ? → Y : ? → B : ? → M :
M ^-1 ( B – X → D :
Lbl 1 : Goto 0 : Lbl 2 :
Y + D Ans → Y : X + D → X :
B > X => Goto 1 : Y
77 Bytes

註:
Y 是代表 y(x) 的變量而 X 則代表 x.
√x² 是指取 x 的絕對值.

啟動程式

按 Prog 1

輸入數據
x₀ EXE y₀ EXE x_n EXE n EXE

2 EXE 1 EXE 5 EXE 36 EXE (分成36份)

顯示答案
y(x_n)

1.919091783
(正確答案 = 1.9182575899)

記憶體

A
B	x_n
C
D	h
X	x_n
Y	y(x_n)
M	n

相關程式