Метод Симпсона и трапеций

uses crt;
const
  a=0;{нижний предел}
  b=1;{верхний предел}
  e=0.001;{точность}
function f(x:real):real;
begin
  f:=sqrt(x)+cos(x*x);
end;
 
var s,k,y1,y2:real;
    it,is:integer;
begin
clrscr;
writeln('Метод трапеций:');
s:=(F(a)+F(b))/2;
it:=1;
k:=a-e;
while k<=b+e do
 begin
  k:=k+e;
  s:=s+F(k);
  it:=it+1;
 end;
s:=s*e;
writeln('S=',s:0:3,' количество итераций=',it);
writeln('Метод Симпсона:');
y1:=0;
k:=a+e;
is:=0;
while k<=b-e do
 begin
  y1:=y1+F(k);
  k:=k+2*e;
  is:=is+1;
 end;
y2:=0;
k:=a+2*e;
is:=0;
while k<=b-2*e do
 begin
  y2:=y2+F(k);
  k:=k+2*e;
  is:=is+1;
 end;
s:=((b-a)/(3*is))*(F(a)+F(b)+2*y1+y2);
writeln('S=',s:0:3,' количество итераций=',is);
if it>is then write('Количество итераций по методу трапеций больше на ', it-is)
else write('Количество итераций по методу Симпсона больше на ', is-it);
readln
end.

uses crt;
function F(x:real):real;
begin
F:=exp(x)-exp(-x)-2;//функция
end;
function F1(x:real):real;
begin
F1:=exp(x)-exp(-x);//первая производная
end;
 
var x1,x2,eps,b:real;
begin
clrscr;
repeat
write('Введите начальное приближение [0;1] a=');
readln(x1);
until(x1>=0)and(x1<=1);
eps:=0.001;
x2:=x1;
repeat
b:=x2;
x2:=b-F(b)/F1(b);
until abs(x2-b)<eps;
write('X=',x2:0:3);
readln
end.