Курсовая работа: Дослідження чисельних методів інтегрування

while (x<b-0.000001) do

begin

j:=j+1;

if ((j mod 2) = 0) then

I:=I+2*f(x)

else

I:=I+4*f(x);

x:=x+h

end;

Simpson:=I*h/3;

end;

{------------------------------------------------}

function difChebushev (a, b: extended): extended;

const

t4:array [1..4] of extended = (-0.794654, -0.187592, 0.187592, 0.794654);

n=4;

var

x, I: extended;

j: integer;

begin

I:=0;

for j:=1 to 4 do

begin

x:=(a+b)/2+t4[j]*(b-a)/2;

I:=I+f(x);

end;

difChebushev:=I*(b-a)/n;

end;

{------------------------------------------------}