Даны действительные числа a, h, натуральное число n. Вычислить f(a) + 2f(a + h) + 2f(a + 2h) + … + 2f(а + (n-1)h) + f(a + nh), где f (х) = (х2 + 1) cos2 x- Pascal(Паскаль)

program d;
var a,s,h:real;
    n,i: integer;
function f(x:real): real;
begin
f:=(sqr(x)+1)*sqr(cos(x))
end;
begin
writeln('Enter a, h, n:');
readln(a,h,n);
s:=0;
for i:=0 to n do s:=s+f(a+i*h);
writeln('S = ',s:2:2);
end.