微分方程 ( 歐拉計算法 ) ( Differential Equations : Euler Method )
注意:此程式只能在 HP-32SII 上使用。
這個程式使用歐拉計算法 ( Euler Method ) 計算微分方程 dy/dx = f(x, y) , y(a) = b 的解的近似值,其中 a 和 b 是常數。
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 )
注意:如想解其他微分方程,只需修改褐色部份。( A 是 x,B 是 y )
程式執行完成後,按 RCL B 會顯示微分方程的解在指定點的近似值。
微分方程 ( 改良歐拉計算法 ) ( Differential Equations : Improved Euler Method )
微分方程 ( Runge-Kutta 計算法 ) ( Differential Equations : Runge-Kutta Method )