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