Casio⁵⁰ » fx-3x50P » Programs 程式集 » Numerical Integration 數值法積分 »

Simpson’s 3/8 Rule
森遜第二法則

本式可以森遜第二法則求定積分

<COMP>
Goto 3 : Lbl 0 : f(A) → Y : Goto 2 :
Lbl 3 : MemClear : ? → A : ? → B : ? → M :
M ^-1 ( B – A → D : Lbl 1 : Goto 0 : Lbl 2 :
X + 3 nCr C × Y → X : A + D → A : C + 1 → C : 3 ≥ C => Goto 1 :
0 → C : A – D → A : B ≥ A + D => Goto 1 :
3 D X ÷ 8 → X
>102 Bytes

操作方法

一般操作

例子

求

int(f(x),a,b)

求

int(cos^2 x^2,0,pi)

整體程式為

<COMP>
Goto 3 : Lbl 0 : (cos A ² ) ² → Y : Goto 2 :
Lbl 3 : MemClear : ? → A : ? → B : ? → M :
M ^-1 ( B – A → D : Lbl 1 : Goto 0 : Lbl 2 :
X + 3 C C × Y → X : A + D → A : C + 1 → C : 3 ≧ C => Goto 1 :
0 → C : A – D → A : B ≧ A + D => Goto 1 :
3 D X ÷ 8 → X
108 Bytes

啟動程式

按 Prog 1

輸入上下限
a EXE b EXE

[A?] 0 EXE
[B?] π EXE

輸入區間數
n EXE

[M?] 18 EXE

顯示答案
Answer

[3DX÷8→X] 1.808612432
(正確答案 = 1.8225763879)

記憶體

A	b
B	b
C	0
D	h
X	Answer
Y	f(b)
M	n

相關程式