Hardy's rule定積分計算

程式編寫日期: 2007年5月2日

程式 (最少104 bytes)

注意: 10x是按shift log。

Mem clear: ?→M: ?→Y: ?→A: (Y - M) ÷ 6A→D:

Lbl 0: ln M: B + Ans( 28(C2=6C) + 162(C2+5=6C) + 220(C=3→B:

DM+: 1+C→C: C=7 => DM- => 1+X→X => 0→C:

A>X => Goto 0: BD÷10x2

 

註1: 程式要輸入的區間數目n是指將定積分計算的範圍平分為n個區間,而每一個區間都會應用Weddle's rule計算,即是將該區間再分六分計算七個對應函數值,然後再應用法則計算定積分。

註2: 綠色的 ln M是函數方程(變數是X),若果想計算其它函數的積分,只要修改綠色的部份。

 

例題: 利用Hardy's rule,區間數目為2,計算

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

 

有關Weddle's rule 可以參考以下網址:

http://mathworld.wolfram.com/HardysRule.html

 

返回 fx-3650P及SC185程式集

Free Web Hosting