微分方程 ( 歐拉計算法 ) ( Differential Equations : Euler Method )

注意:此程式只能在 HP-32SII 上使用。

這個程式使用歐拉計算法 ( Euler Method ) 計算微分方程 dy/dx = f(x, y) , y(a) = b 的解的近似值,其中 ab 是常數。

LBL A
R/S
STO A
+/-
STO C
R/S
STO B
R/S
STO +C
R/S
STO D
STO ÷C
1
STO E
LBL B
XEQ F
RCL C
STO +A
×
STO +B
1
STO +E
RCL D
ENTER
RCL E
x y?
GTO B
GTO E
LBL F
3
RCL ×A
RCL +B
2
+
RTN
LBL E
RCL B
PRGM
最少 48 bytes

( 程式連同數字記憶共使用最少 88 bytes。)
( 如果日後想繼續輸入其他程式,則請在離開 PRGM 模式前按 RTN。)

例:使用歐拉計算法計算微分方程 dy/dx = 2 + 3x + y, y(0) = 1 的解在 x = 1 的近似值,其中區間數量 ( Number of steps )10

XEQ A,再按

0  R/S  ( x 的最初數值  Initial Value )
1  R/S  ( y 的最初數值  Initial Value )
1  R/S  ( x 的最後數值  Final Value )
10  R/S  ( 區間數量  Number of steps )

顯示   7.56245  ( y(1) 的近似值 )

亦即是說 y(1) 的近似值是 7.56245。( y(1) 的準確數值是 8.30969 )

注意:如想解其他微分方程,只需修改褐色部份。( AxBy )

程式執行完成後,按 RCL B 會顯示微分方程的解在指定點的近似值。

微分方程 ( 改良歐拉計算法 ) ( Differential Equations : Improved Euler Method )

微分方程 ( Runge-Kutta 計算法 ) ( Differential Equations : Runge-Kutta Method )

返回選擇程式頁

返回頁首

Hosted by www.Geocities.ws

1