Численное интегрирование (Pascal)

Листинг программы
{$R-}
program Lab6;
uses crt;
var
a,b:real; { Границы отрезка }
r:real;
n:integer;

function f(x:real):real;
begin
f:= cos(sqr(x));
end;
...

Цель: Изучить формулы приближенного вычисления интегралов и погрешностей интегрирования.

Задание: Вариант3.
Разработать алгоритм и составить программу для приближенного вычисления интеграла при n=4, заданным методом для функции из лабораторной работы 3.
Функция cos(x^2). Интервал: (0; 0,5).
Метод вычисления: метод Симпсона(парабол).

p.s. В работе также содежится исходный код, блок схема, результаты выполнения и теоретические cведения.