Вычисление интеграла функции методом трапеций

program Integral_trap;

{$APPTYPE CONSOLE}

var a,b,h,s,x,y,eps: real;

i,n: integer;

function f(c: real): real;     // вычисляемая функция

begin

f:=1/(sqrt(2*sqr(x)+0.3)); 

end;

begin

writeln('a=');      // ввод значений

readln(a);

writeln('b=');

readln(b);

writeln('eps=');

readln(eps);

n:= Round((b-a)/eps);

s:=0;

h := (b-a)/n;

for i:=1 to n-1 do       // расчет и вывод результатов

begin

x:= a + i* h;

s:= s + f(x);

writeln ('x=',x: 5 : 3, '  s=',s: 5: 3);

end;

y:= h*((f(a)+f(b))/2+s);

writeln ('n=',n, ' y=',y: 8: 3);

readln;

end.