高斯-勒讓德七點定積分(I)
程式編寫日期: 2008年3月24日
以下程式使用了高斯-勒讓德法七點方法(Gauss-Legendre Quadrature)計算一個定積分的近似值，現時不少圖像計算機亦是使用相似的方法計算定積分的近似值。
程式 (183 bytes，不包括綠色的函數方程)
Mem clear: ?→A: ?→B: 2-1(B - A→B:
Lbl 0: A + B + BX: ln Ans: Ans ( . 4179591837(Y=0)
+ . 1294849662(Y2+2=3Y) + . 2797053915(Y>2)(5>Y)
+ . 3818300505(Y>4M+: 1+Y→Y: -X→X:
Y=1 => . 9491079123→X: Y=3 => . 7415311856→X:
Y=5 => . 4058451514→X: Y-7 => Goto 0: MB
註: 綠色的 ln Ans是函數方程(變數是Ans)，若果想計算其它函數的積分，只要修改綠色的部份。
例題: 利用高斯-勒讓德七點方法(Gauss-Legendre Quadrature)，計算
按 Prog 1 再按 1 EXE 2 EXE (顯示答案為0.386294361)
參考資料:
http://mathworld.wolfram.com/Legendre-GaussQuadrature.html
http://en.wikipedia.org/wiki/Gaussian_quadrature
