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

Hardy’s Rule

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

<COMP>
Goto 3 : Lbl 0 : f(A) → Y : Goto 2 :
Lbl 3 : MemClear : ? → A : ? → B : ? → M :
M -1 ( B – A → D :
Lbl 1 : 6 nCr C = 15 => Goto 4 : Goto 0 : Lbl 2 : 6 nCr C :
X – Y ( 793 Ans ² – 23373 Ans + 3960 → X :
Lbl 4 : A + D → A : C + 1 → C : 6 ≥ C => Goto 1 :
0 → C : A – D → A : B > A + D => Goto 1 :
D X ÷ 66500 → X
>136 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 : 6 nCr C = 15 => Goto 4 : Goto 0 : Lbl 2 : 6 nCr C :
X – Y ( 793 Ans ² – 23373 Ans + 3960 → X :
Lbl 4 : A + D → A : C + 1 → C : 6 ≧ C => Goto 1 :
0 → C : A – D → A : B > A + D => Goto 1 :
D X ÷ 66500 → X
142 Bytes
啟動程式 Prog 1
輸入上下限
a EXE b EXE
[A?] 0 EXE
[B?] π EXE
輸入區間數
n EXE
[M?] 18 EXE
顯示答案
Answer
[DX÷66500→X] 1.857697813
(正確答案 = 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. Weddle’s Rule
  7. Gauss-Legendre 2-Point Quadrature / 高斯-勒讓德兩點求積法
  8. Gauss-Legendre 3-Point Quadrature / 高斯-勒讓德三點求積
Hosted by www.Geocities.ws

1