Вычисление определенных интегралов

program integral;uses crt;
var i,j,m,n:integer;
h,x,sum,a,b,e,it1,it2,it3:real;
function f(x:real): real;
begin
f:=sin(x);end;
procedure trap(n,i:integer;var it:real);
begin
h:=(b-a)/n;x:=0;
sum:=0;
for i:=1 to n-1 do
begin
x:=a+h*i;
sum:=sum+f(x);
end;
it:=(sum+(f(a)+f(b))/2)*h;
end;
procedure sim(m,i:integer; var sim:real);
begin
x:=a;
h:=(b-a)/n;sum:=0;
for i:=1 to 2*n-1 do
begin
x:=x+(h/2);
if (i/2=i div 2) then sum:=sum+f(x)/3;
if (i/2<>i div 2) then sum:=sum +2*f(x)/3;
end;
sim:=(sum+(1/6)*(f(a)+f(b)))*h;
end;
procedure pram(n,i:integer;var ip:real);
begin
h:=(b-a)/n;x:=a-h/2;sum:=0;
for i :=0 to n-1 do
begin
x:=x+h;sum:=sum+h*f(x);
end;
ip:=sum;
end;
begin
clrscr;
{write('введите начало отрезка');
readln(a);
write('введите конец отрезка');
readln(b);}
a:=-2*arctan(1);
b:=2*arctan(1);
write('Количество разбиений = ');
readln(m);
write('Точность =');
readln(e);
n:=m;
repeat
trap(n,i,it1);trap(2*n,i,it2);
n:=2*n;
until ((abs(it1-it2))/3)<e;
write(' Интеграл по методу трапеции: ');
writeln(it1);
repeat
sim(n,i,it1);sim(2*n,i,it2);
n:=2*n;
until ((abs(it1-it2))/15)<e;
write(' Интеграл по методу Симпсона: ');
writeln(it1);
repeat
pram(n,i,it1);
pram(2*n,i,it2);
n:=2*n;
until ((abs(it1-it2))/3)<e;
write(' Интеграл по методу прямоугольника: ');
writeln(it1);readln;
end.

program z8;
const a=1; Eps=0.01;
var n,I,k:integer; D,h,S,S1,b:real;
   begin
      begin
b:=exp();
n:=10;
S:=0;
  repeat
h:=(b-a)/n;
S1:=S;
for i:=1 to n do 
S:=S+h*((1+ln((i-0.5)*h+a))/
((i-0.5)*h+a));
n:=n*2;
  until abs(S1-S)<Eps;
writeln (‘S=’, S:8:3);
      end;
    begin
k:=10;
h:=(b-a)/k;
D:=0;
for i:=1 to k do 
D:=D+sqrt(sqr(a+i*h-(a+(i-1)*h))+
sqr((1+ln(a+i*h)/(a+i*h))-
((1+ln(a+(i-1)*h))/(a+(i-1)*h)));
writeln (‘D=’, D:6:3);
    end;
end.