Вычислить интеграл по формуле прямоугольников - Turbo Pascal

uses crt;
function f(x:real):real;
begin
f:=x+1;
end;
var a,b,p,h,x,s:real;
    t,k,n,i:integer;
begin
clrscr;
t:=2;{t=1*2}
a:=t*t*t-sin(t);
if a>t then b:=a else b:=t;
p:=b+t*t;
write('Число разбиений n=');{чем больше n, тем выше точнсть}
readln(n);
h:=(p-a)/n;
x:=a+h/2;
s:=0.0;
for i:=1 to n do
 begin
  s:=s+f(x);
  x:=x+h;
 end;
s:=s*h;
write('Интеграл = ',s:0:4);
readln;
end.