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

program integral;
var
 n,i : integer;
 a,b,shag,sum,itog : real;
 
function F(x:real):real;
 begin
   F:=ввести интегрируемую функцию 
 end;
{======================================}
 
begin
 write('Начало интегрирования a = '); readln(a);
 write('Конец  интегрирования b = '); readln(b);
 write('Количество разбиений интервала n = '); readln(n);
 
 shag:=(b-a)/n;
 sum:=0;
 for i:=1 to n-1 do
   sum := sum + F(shag*i+a);
 sum := sum + (F(a)+F(b))/2;
 
 itog:=(b-a)/n * sum;
 writeLn('Интеграл = ', itog:0:5);
 readln
end.

program integral;
var
 n,i : integer;
 a,b,shag,sum,itog,c,d : real;
 
function F(x:real):real;
 begin
   F:=(1/c*x+d);
 end;
{======================================}
 
begin
 write('Начало интегрирования a = '); readln(a);
 write('Конец  интегрирования b = '); readln(b);
 write('Количество разбиений интервала n = '); readln(n);
 write('Введите коэффициент при x, c= '); readln(c);
 write('Введите свободный коэффициент d= '); readln(d);
 shag:=(b-a)/n;
 sum:=0;
 for i:=1 to n-1 do
   sum := sum + F(shag*i+a);
 sum := sum + (F(a)+F(b))/2;
 
 itog:=(b-a)/n * sum;
 writeLn('Интеграл = ', itog:5:5);
 readln
end.