Курсовая работа: Інтеграли зі змінними границями

for i:=1 to n do

GridOfInt[i]:=integralSimpsona(a,Xgrid[i]);

end;

//--------------------------------------------------------------

function integralSimpsona(aSim,bSim:real):real;

var X:one_array; i:integer; sum1, sum2, h:real;

begin

if(aSim<>bSim)then

begin

h:=(bSim-aSim)/nSim;

SetLength(X,nSim+1);

for i:=0 to nSim do

X[i]:=aSim+h*i;

sum1:=0;

sum2:=0;

for i:=1 to (nSim div 2) do

sum1:=sum1+y(X[2*i-1]);

for i:=1 to ((nSim div 2)-1) do

sum2:=sum2+y(X[2*i]);

Result:=(bSim-aSim)*(y(X[0])+y(X[nSim])+4*sum1+2*sum2)/(3*nSim);

end

else

Result:=0;

end;

//--------------------------------------------------------------

function integralTrapeciay(aTrap,bTrap:real):real;

var i:integer; sum, h :real; X:one_array;

begin

h:=(bTrap-aTrap)/nTrap;