辛卜生第二法則定積分計算

辛卜生第二法則定積分計算

程式編寫日期: 2006年6月21日

程式 (最少81 bytes)

Mem clear: ?→A: ?→B: ?→C: C^-1(B - A→B:

Fix 0 : Lbl 0: A + BM: ln Ans→ X: M ÷ 3: Rnd:

M=3Ans: D + X(3 - Ans - (M²=CM→D: 1M+:

C≧M => Goto 0: 3DB ÷ 8: Norm 1

注意: 綠色的 ln Ans是函數方程(變數是Ans)，若果想計算其它函數的積分，只要修改綠色的部份。

例題: 利用辛卜生第二法則，區間數目為12，計算

按 Prog 1 再按 1 EXE 2 EXE 12 EXE (顯示答案為0.386293328)

注意: 辛卜生第二法則的區間數目必須為3的倍數，否則計算將會錯誤。