Помогите написать программу на Паскале, пожалуйста, Метод Симпсона на языке pascal Опции

Здравствуйте, помогите пожалуйста написать программу на языке pascal

Написать подпрограмму численного интегрирования функции f(x) на отрезке [a,b] с достижением заданной точности методом Симпсона. Вычисление функции оформить в виде подпрограммы с входным параметром sx - значением аргумета x и выходным параметром sf - значением функции f(x).
Алгоритм достижения точности:
- применим составную формулу Симпсона для количества подотрезков на отрезке [a,b], равном n=1. Получим приблежение к интегралу R1;
- удвоим количество подотрезков n и получим очередное приблежение к интегралу R2;
- используя значения R1 и R2 проверяем достижение точности. Если достигнута требуемая точность, то в качестве ответа выбираем R2. Если точность не достигнута, повторяем процесс удвоения количества подотрезков n.
Входные параметры:
- sa,sb - отрезок [a,b], на котором вычисляется интеграл;
- sd,se - абсолютная и относительная точности интегрирования;
- sl - максимальное разрешенное количество обращений к функции f(x).
Выходные параметры:
- sr - ответ
- sn - реально потребовавшееся количество обращений к f(x);
- si - индикатор ошибки: =0 - точность достигнута, =1 - использованы все разрешенные sl вычислений функции, но точность не достигнута, =2 - sd и se оюа =<0, либо sl<5.