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

Gauss-Legendre 3-Point Quadrature
高斯-勒讓德三點求積法

本式可以高斯-勒讓德求積法 (使用三點) 求出定積分的值

<COMP>
Goto 3 : Lbl 0 : f(Ans) : Goto 2 :
Lbl 3 : MemClear : ? → A : ? → B : ? → M :
. 5 M -1 ( B – A → D : Fix0 :
Lbl 1 : C : Rnd : √ 5.4 ( C – Ans → Y : D Y + A + D : Goto 0 :
Lbl 2 : X + 9 -1 Ans ( 8 – 5 Y ² → X :
C + 3 -1 → C : Rnd : Ans ≠ C => Goto 1 :
A + 2 D → A : M > C => Goto 1 :
Norm1 : D X → X
>121 Bytes

操作方法

一般操作 例子

求
int(f(x),a,b)

求
整體程式為
<COMP>
Goto 3 : Lbl 0 : ( cos Ans ² ) ² : Goto 2 :
Lbl 3 : MemClear : ? → A : ? → B : ? → M :
. 5 M -1 ( B – A → D : Fix0 :
Lbl 1 : C : Rnd : √ 5.4 ( C – Ans → Y : D Y + A + D : Goto 0 :
Lbl 2 : X + 9 -1 Ans ( 8 – 5 Y ² → X :
C + 3 -1 → C : Rnd : Ans ≠ C => Goto 1 :
A + 2 D → A : M > C => Goto 1 :
Norm1 : D X → X
127 Bytes
啟動程式 按 Prog 1
輸入上下限
a EXE b EXE		 [A?] 0 EXE
[B?] π EXE
輸入區間數
n EXE		 [M?] 8 EXE
顯示答案
Answer		 [DX→X] 1.822629468
(正確答案 = 1.8225763879)

記憶體

A b
B b
C n
D h/2
X Answer
Y -√0.6
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. Hardy’s Rule
  8. Gauss-Legendre 2-Point Quadrature / 高斯-勒讓德兩點求積法
Hosted by www.Geocities.ws

1