Курсовая работа: Программа вычисления значения и вывода графика определенного интеграла

dx:=(b-a)/DivNo;

for i:=1 to DivNo-1 do

begin

x1:=a+dx*i;

if MainFunction(x1,y1) then sum:=sum+y1 else Raise EAbort.Create('');

end;

sum:=sum*2;

if MainFunction(a,y1) and MainFunction(b,y2) then

sum:=sum+y1+y2 else Raise EAbort.Create('');

Res:=sum*dx/2;

except

PTrap:=False;

end;

end;

{-----}

var

T: integer;

i1, i2: real;

begin

countTrap:=true;

try

if PTrap(MainFunction,a,b,2,i1) and PTrap(MainFunction,a,b,4,i2) then

begin

T:=4;

while abs(i1-i2)>eps do

begin

if PTrap(MainFunction,a,b,t,i1) and PTrap(MainFunction,a,b,t*2,i2) then t:=t*2

else EAbort.Create('');

T:=T*2;