qgaus.txt

来自「1.梯形求积法 2.辛普森求积法 3.龙贝格求积法 4.反常积分 5.高」· 文本 代码 · 共 21 行

procedure QGAUS(A, B:real;var SS:real);
var
    X, W:array[0..5] of real;
    XM,XR,DX:real;  J:integer;
begin
    X[1]:=0.1488743389;    X[2]:=0.4333953941;
    X[3]:=0.6794095682;    X[4]:=0.8650633666;
    X[5]:=0.9739065285;
    W[1]:=0.2955242247;    W[2]:=0.2692667193;
    W[3]:=0.2190863625;    W[4]:=0.1494513491;
    W[5]:=0.0666713443;
    XM:=0.5 * (A + B);
    XR:=0.5 * (B - A);
    SS:=0;
    For J:=1 To 5 do
    begin
        DX:=XR * X[J];
        SS:=SS + W[J] * (FUNC(XM + DX) + FUNC(XM - DX));
    end;
    SS:=XR * SS;
end;