梯形及辛卜生(森遜)法則定積分計算
程式可以計算梯形法則或辛卜生定積分問題。第二個程式雖然長一點，但函數的修改部份接近程式的開始位置，因此較容易及快速修改函數，使整體速度加快。
程式新版
程式編寫日期: 2009年2月11日
注意: πr 是按 SHIFT EXP SHIFT Ans 2
第一個程式 (73 bytes，不包括綠色函數方程)
Mem clear: ?→B: ?→X: ?→Y:
?→A: A-1(Y-X→Y: Lbl 0: ln X:
C + Ans(3 - cos πrBM - (M2=MA→C:
X+Y→X: A>MM+ => Goto 0: YC÷( 2+B
第二個程式 (82 bytes，不包括綠色函數方程)
Goto 1: Lbl 0: ln X: C + Ans(3 - cos πrBM - (M2=MA→C:
X+Y→X: A>MM+ => Goto 0: YC÷( 2+B◢
Lbl 1: Mem clear: ?→B: ?→X: ?→Y: ?→A:
A-1(Y-X→Y: Goto 0
註1: 綠色的 ln X 是函數方程(變數是X)，若果想計算其它函數的積分，只要修改綠色的部份。
註2: 進行三角函數的積分計算時請將度單位設定為Rad(弧度)，否則答案將會錯誤。
註3: 若果想保留一些記憶用作儲存臨時數據，第一個可將程式中"Mem clear"改為"MM-"， 第二個程式將程式中"Mem clear"改為"MM-: 0→B"，而程式所使用的記憶為A，B，C，X，Y及M。
例題1: 利用梯形法則，區間數目為10，計算
按 Prog 1 再按 0 EXE (輸入0代表計算梯形法則)
1 EXE 2 EXE 10 EXE (顯示答案為0.385877936)
若果是第二個程式計算完結後，請按 AC 終止程式。
例題2: 利用辛卜生法則，區間數目為10，計算
按 Prog 1 再按 1 EXE (1代表計辛卜生法則)
1 EXE 2 EXE 10 EXE (顯示答案為0.386293403)
若果是第二個程式計算完結後，請按 AC 終止程式。
注意: 梯形法則的區間數目必為正整數，而辛卜生法則的區間數目必為正雙數
