помогите с решением интеграла в программе Delphi7

[Imax2]

Условие задания:
Программа вычисляет определенный интгерал функции: е(в степени Х) соs(в квадрате) X, на промежутке интегрирования [0;П].тремя способами: прямоугольников, трапеций и параболы(Симпсона) при заданном числе разбиений отрезка интегрирования.Вычесляет точное значение интеграла по формуле НьютонаоЛейбница.Программа строит график подинтегральной функции и на одном рисунке 3 графика зависимости относительной погрешности вычесления интегралов для каждого способа от числа отрезков разбиения.

Сложность:
Не могу(не умею) написать программный код для этой задачи.
Пожалуйста помогите.Буду рад любой помощи.Заранее спасибо