Даны действительные числа a,h, целое n. Вычислить:f(a) + f(a+h)+f(a+2*h)+…+ f(a+n*h), где f(x)=(x^2+1)/x — Pascal(Паскаль)

program p1;
uses crt;
var
a,h,sum:real;
n:integer;
function F(x:real):real;
begin
F := (sqr(x)+1)/x;
end;
function Summa(a,h:real;n:integer):real;
var i:integer;
    s:real;
begin
    s:=0.0;
    for i:=0 to n do 
    begin
     s := s + F(a+i*h);
    end;
    Summa := s;
end;
var 
begin
clrscr;
writeln('Vvedite dannie');
write('a=');readln(a);
write('h=');readln(h);
write('n=');readln(n);
sum := Summa(a,h,n);
writeln('Summa S =',sum:8:2);
readln;
end.