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

Simpson’s Rule
森遜法則

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

<COMP>
Goto 3 : Lbl 0 : f(A) → Y : Goto 2 :
Lbl 3 : MemClear : ? → A : ? → B : ? → M :
M -1 ( B – A → D :
Goto 0 : Lbl 2 :
X + Y (2 nCr C) ² → X : A + D → A : C + 1 → C : 2 ≥ C => Goto 0 :
0 → C : A – D → A : B ≥ A + D => Goto 0 :
D X ÷ 3 → X
>100 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 :
Goto 0 : Lbl 2 :
X + Y (2 nCr C) ² → X : A + D → A : C + 1 → C : 2 ≥ C => Goto 0 :
0 → C : A – D → A : B ≥ A + D => Goto 0 :
D X ÷ 3 → X
106 Bytes
啟動程式 按 Prog 1
輸入上下限
a EXE b EXE		 0 EXE
π EXE
輸入區間數
n EXE		 18 EXE
顯示答案
Answer		 1.828054552
(正確答案 = 1.8225763879)

記憶體

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

相關程式

  1. Simpson’s 3/8 Rule / 森遜第二法則
  2. Trapezoidal Rule / 梯形法則
  3. Mid-Point Rule / 中點法則
  4. Boole’s Rule / 保爾法則
  5. Weddle’s Rule
  6. Hardy’s Rule
  7. Gauss-Legendre 2-Point Quadrature / 高斯-勒讓德兩點求積法
  8. Gauss-Legendre 3-Point Quadrature / 高斯-勒讓德三點求積法
Hosted by www.Geocities.ws

1