Помогите пожалуйста с программированием определенного интеграла в Паскале(методом прямоугольников)

Вот программный код:
Program Integral;
type func=function(x:real):real;
Var a,b:real;n:integer;
{$F+}
function f1(x:real):real;
Begin
f1:=x*(exp(ln(1+x)/3 ));
End;
{$F-}
Function integ(a,b:real;n:integer;f:func):real;
Var s,h,x:real;i:integer;
Begin
h:=(b-a)/n;
s:=0;
For i:=1 to n do
s:=s+h*f(a+i*h);
integ:=S;
End;
Begin {main}
readln(a,b,n);
writeln(itnteg(a,b,n,f1):7:4);
readln;
End.
Программа работает но считает интеграл неверно
пределы интегрирования - от 8 до 1, использую метод прямоугольников)
В проге a,b -пределы интегрирования, n - число разбиений, h-шаг.
И да еще если не сложно, можете помочь посчитать этот интеграл с заданной точностью( я видел готовую блок-схему там вроде нужно использовать цикл for, но для чего все это я так и не понял).
Заранее спасибо!!!