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

Weddle’s Rule

本式可以 Weddle 法則求出定積分的值

<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 : 6 nCr C :
X + Y ( 121 Ans ³ – 3726 Ans ² + 30455 Ans – 14880 → X :
A + D → A : C + 1 → C : 6 ≥ C => Goto 1 :
0 → C : A – D → A : B > A + D => Goto 1 :
D X ÷ 39900 → X
>129 Bytes

操作方法

一般操作 例子


int(f(x),a,b)


整體程式為
<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 : 6 nCr C :
X + Y ( 121 Ans ³ – 3726 Ans ² + 30455 Ans – 14880 → X :
A + D → A : C + 1 → C : 6 ≧ C => Goto 1 :
0 → C : A – D → A : B > A + D => Goto 1 :
D X ÷ 39900 → X
137 Bytes
啟動程式 Prog 1
輸入上下限
a EXE b EXE
[A?] 0 EXE
[B?] π EXE
輸入區間數
n EXE
[M?] 18 EXE
顯示答案
Answer
[DX÷39900→X] 1.843608248
(正確答案 = 1.8225763879)

記憶體

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

其他資料

n 必為六之倍數

相關程式

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

1