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

Gauss-Legendre 2-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 : D + A + D √ 3 -1 ( 1 + 4 C – 4 Ans :
Goto 0 : Lbl 2 : Ans + X → X :
C + . 5 → C : Rnd : Ans ≠ C => Goto 1 :
A + 2 D → A : M > C => Goto 1 :
Norm1 : D X → X
>112 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 : D + A + D √ 3 -1 ( 1 + 4 C – 4 Ans :
Goto 0 : Lbl 2 : Ans + X → X :
C + . 5 → C : Rnd : Ans ≠ C => Goto 1 :
A + 2 D → A : M > C => Goto 1 :
Norm1 : D X → X
118 Bytes
啟動程式 按 Prog 1
輸入上下限
a EXE b EXE		 [A?] 0 EXE
[B?] π EXE
輸入區間數
n EXE		 [M?] 12 EXE
顯示答案
Answer		 [DX→X] 1.821529591
(正確答案 = 1.8225763879)

記憶體

A b
B b
C n
D h/2
X Answer
Y  
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 3-Point Quadrature / 高斯-勒讓德三點求積法
Hosted by www.Geocities.ws

1