定積分計算 ( 梯形法則及森遜法則 )
( Numerical Integration : Trapezoidal
Rule and Simpson's Rule )
注意：此程式只能在 HP-32SII 上使用。
|
LBL A
|
INPUT A
|
INPUT B
|
STO H
|
RCL A
|
STO –H
|
|
LBL B
|
INPUT N
|
x ≦
0?
|
GTO B
|
RCL N
|
IP
|
|
x ≠ y?
|
GTO B
|
0
|
STO J
|
RCL N
|
STO ÷H
|
|
STO C
|
XEQ C
|
RCL J
|
R/S
|
STO I
|
RCL N
|
|
2
|
÷
|
IP
|
RCL N
|
2
|
÷
|
|
x ≠ y?
|
GTO B
|
2
|
STO ÷N
|
STO ×H
|
RCL N
|
|
STO C
|
0
|
STO J
|
XEQ C
|
RCL I
|
4
|
|
×
|
RCL –J
|
3
|
÷
|
STO J
|
GTO E
|
|
LBL C
|
RCL A
|
STO X
|
LBL D
|
XEQ F
|
STO +J
|
|
RCL H
|
STO +X
|
XEQ F
|
STO +J
|
1
|
STO –C
|
|
RCL C
|
x >
0?
|
GTO D
|
RCL H
|
2
|
÷
|
|
STO ×J
|
RTN
|
LBL F
|
RAD
|
RCL
X
|
COS
|
|
RCL
X
|
ex
|
÷
|
RTN
|
LBL E
|
PRGM
|
|
最少
106.5 bytes
|
|||||
( 程式連同數字記憶共使用最少
162.5 bytes。)
( 如果日後想繼續輸入其他程式，則請在離開 PRGM
模式前按 RTN。)
| 例：利用梯形法則 ( Trapezoidal Rule ) 及森遜法則 ( Simpson's Rule )，用 10 個區間 ( subintervals ) 計算積分 |
 |
先按 XEQ A
再按 0 R/S ( 積分下限 (
Lower Limit ) )
再按 1 R/S ( 積分上限 (
Upper Limit ) )
再按 10 R/S ( 區間數量 (
Number of subintervals )，注意此數必須是正整數 )
顯示 0.555807
( 使用梯形法則計算的積分近似值 )
再按 R/S 顯示 0.555396
( 使用森遜法則計算的積分近似值 )
( 如果輸入的區間數量 ( Number of subintervals
) 不是雙數，則此時會顯示 N?，並要求重新輸入區間數量。)
| 亦即是說使用梯形法則 |  |
| 而使用森遜法則 |  |
注意：如想計算其他積分，只需修改褐色部份。( X 是函數變數 )
程式執行完成後，按 RCL I 會顯示使用梯形法則計算的積分近似值，而按 RCL J 會顯示使用森遜法則計算的積分近似值。