定積分計算 ( 森遜法則及保爾法則 )
( Numerical Integration : Simpson's
Rule and Boole's Rule )
( 程式版本：1.0，最後更新日期 14 AUG 2003。)
|
LBL 0
|
R/S
|
STO 1
|
R/S
|
STO 2
|
STO 4
|
|
RCL 1
|
STO –4
|
R/S
|
STO ÷4
|
÷
|
2
|
|
=
|
STO 3
|
STO 5
|
n!
|
XEQ 1
|
RCL 6
|
|
STO 7
|
R/S
|
STO 8
|
1
|
6
|
STO ×8
|
|
2
|
STO ÷3
|
STO ×4
|
RCL 3
|
STO 5
|
n!
|
|
XEQ 1
|
RCL 6
|
STO –8
|
1
|
5
|
STO ÷8
|
|
RCL 8
|
GTO E
|
LBL 1
|
RCL 1
|
STO 0
|
XEQ F
|
|
STO 6
|
XEQ C
|
RCL 5
|
x = 0?
|
GTO 3
|
LBL 2
|
|
RCL 4
|
STO +0
|
XEQ F
|
STO +6
|
STO +6
|
XEQ C
|
|
RCL 5
|
x = 0?
|
GTO 3
|
GTO 2
|
LBL 3
|
RCL 2
|
|
STO 0
|
XEQ F
|
STO +6
|
RCL 4
|
STO ×6
|
3
|
|
STO ÷6
|
RTN
|
LBL C
|
RCL 4
|
STO +0
|
XEQ F
|
|
×
|
4
|
=
|
STO +6
|
1
|
STO –5
|
|
RTN
|
LBL F
|
RAD
|
RCL
0
|
COS
|
÷
|
|
RCL
0
|
ex
|
=
|
RTN
|
LBL E
|
PRGM
|
|
最少
83 steps
|
|||||
| 例：利用森遜法則 ( Simpson's Rule ) 及保爾法則 ( Boole's Rule )，用 12 個區間 ( subintervals ) 計算積分 |
 |
先按 XEQ 0
再按 0 R/S ( 積分下限 (
Lower Limit ) )
再按 1 R/S ( 積分上限 (
Upper Limit ) )
再按 12 R/S ( 區間數目 (
Number of subintervals )，注意此數必須是正整數及雙數，否則會顯示錯誤 Error-Func。)
顯示 0.555396 ( 使用森遜法則計算的積分近似值 )
再按 R/S 顯示 0.555397
( 使用保爾法則計算的積分近似值 )
( 如果輸入的區間數量 ( Number of subintervals
) 不是 4 的倍數，則此時會顯示錯誤
Error-Func。)
| 亦即是說使用森遜法則 |  |
| 而使用保爾法則 |  |
注意：如想計算其他積分，只需修改褐色部份。( R0 是函數變數 )
程式執行完成後，按 RCL 7 會顯示使用森遜法則計算的積分近似值，而按 RCL 8 會顯示使用保爾法則計算的積分近似值。