Разработать программу для вычисления значения определенного интеграла от функции f(X) по заданному

uses crt;
const
  a = 1;{нижний предел}
  b = 5;{верхний предел}
function F(x:real):real;
begin
F:=x*sqrt(2+sin(pi*x));
end;
 
var n,i:longint;
    e,h,k,s: real;
begin
clrscr;
repeat
write('Введите точность от 0.0001 до 0.01 e=');
readln(e);
until (e>0.00099)and(e<=0.0101);
n:=round((b-a)/e);
s:=0;
k:=a;
for i:=1 to n-1 do
 begin
  k:=k+e;
  if i mod 2 = 0 then s:=s+2*f(k)
  else s:=s+4*f(k);
 end;
s:=(s+f(a)+f(b))*e/3;
write('Интеграл=',s:0:6);
readln
end.

var a,b:real;
.....
write('Начало интервала a=');
readln(a);
repeat
write('Конец интервала >',a:0:2,' b=');
readln(b);
until b>a;

uses crt;
const
  a = 1;{нижний предел}
  b = 5;{верхний предел}
function F(x:real):real;
begin
F:=x*sqrt(2+sin(pi*x));
end;
 
var n,i:longint;
    e,h,k,s: real;
 
procedure vvod;
begin
repeat
write('Введите точность от 0.0001 до 0.01 e=');
readln(e);
until (e>0.00099)and(e<=0.0101);
n:=round((b-a)/e);
s:=0;
k:=a;
for i:=1 to n-1 do
 begin
  k:=k+e;
  if i mod 2 = 0 then s:=s+2*f(k)
  else s:=s+4*f(k);
 end;
s:=(s+f(a)+f(b))*e/3;
write('Интеграл=',s:0:6);
end;
 
begin
clrscr;
vvod;
readln
end.