梯形及辛卜生(森遜)法則定積分計算

程式新版 (舊版保留在網頁下方)

程式可以計算梯形法則或辛卜生法則定積分問題。第二個程式雖然長一點,但函數的修改部份接近程式的開始位置,因此較容易及快速修改函數,使整體速度加快。 注意新版方程變數使用X,與舊版使用Ans不同。

程式編寫日期: 2009年2月16日

注意: πr 是按 SHIFT EXP SHIFT Ans 2

第一個程式 (69 bytes,不包括綠色函數方程)

ClrMemory: ?→B: ?→X: ?→Y:

?→A: A-1(Y-X→Y: For 0→C To A: ln( X:

Ans(3 - cos( πrBC ) - (C2=CAM+:

X+Y→X: Next: YM÷( 2+B

 

第二個程式 (83 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: ClrMemory: ?→B: ?→X: ?→Y: ?→A:

A-1(Y-X→Y: Goto 0

 

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

註2: 進行三角函數的積分計算時請將度單位設定為Rad(弧度),否則答案將會錯誤。

註3: 若果想保留一些記憶用作儲存臨時數據,第一個可將程式中"ClrMemory"改為"MM-", 第二個程式將程式中"ClrMemory"改為"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 終止程式。

 

注意: 梯形法則的區間數目必為正整數,而辛卜生法則的區間數目必為正雙數

 

返回 CASIO fx-50FH、fx-3650P II、fx-50FH II及fx-50F PLUS 程式集

 

 



程式舊版  (另一個版本)

程式編寫日期: 2008年2月12日

注意: πr 是按 SHIFT EXP SHIFT Ans 2

程式 (68 bytes,不包括綠色函數方程)

ClrMemory: ?→X: ?→A: ?→B: ?→C: C-1(B-A→D:

For 0→B To C: A + BD: ln( Ans: Ans(3 - cos( πrBX ) - (B2=BCM+:

Next: DM÷( 2+X

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

註2: 若果想保留一些記憶用作儲存臨時數據,可將程式中"ClrMemory"改為"MM-",而程式所使用的記憶為A,B,C,D,X及M。