Pascal/Turbo Pascal - integracion de funciones

HOLA COMPAÑEROOOS!!

Ante todo estado buscando inforamcion por todas partes sobre como realizar este programa, pero resulta que estoy encasquildo y ni si quiera se empezar, pero me dijeron como empezar y no se seguir, y me gustaria saber si podrian ayudarme a solucionar este programilla, soy un principiante sobre este tema en pascal, y lei este problemilla en una hoja de problemas de pascal de internet, pero la verdad que no tengo ni idea de como hacerlo, si me pudieran ayudar se lo agradeceria.
El problema se dictaba mas o menos de esta manera:


definid una funcion Pascal como la del siguiente ejemplo:

function f(x: real) : real ;
begin
f := sqr(x) ; {o lo que se quiera}
end;

e implementad la funcion:

function integral(a, b :real ; n : integer) : real;

que integra la funcion f en el intervalo [a, b] usando la aproximacion

(integral de ''a'' a ''b'' de f(X) d(x) = sumatorio (n-1) y (i=0) de f(a+ibeta) donde beta= (beta-a) /n. Observando la convergencia de este metodo de integracion variando n.

Cualquier ayuda que sea aportada sera de muchisima utilidad muchisisa gracias !!

Un saludo!!