Задание необходимо выполнить в Паскале используя оператор цикла с параметром (for).Даны

Задание необходимо выполнить в Паскале используя оператор цикла с параметром (for).
Даны действительные a, b и естественное n. Вывести таблицу значений функции f(x). Таблица имеет две колонки, в i-ю строчку таблицы заносятся значения xi и yi, где xi=a+ih, h=(b-a)/n, yi=f(xi)(i=0,1...,n).
f(x)=(x2+1)cos(x)
Вычислить f(a)+f(a+h)+...+f(a+nh)
Делал сам, но вышло криво:

var a, b, h, x, f, xi, yi: real;
n, i: integer;
Begin
write('a=');
read(a);
write('b=');
read(b);
write('n=');
read(n);
for var i := 1 to 10 do
writeln(i);

h:=(b-a)/n;
writeln('h=', h);
f:=(sqr(x)+1)*cos(x);
xi:=a+i*h;
writeln('xi=',xi);
yi:=xi;
writeln('yi=',yi);
writeln('f=', f);
End.