Численные методы

помощь делать задачу численных методов по MATLAB
Используя квадратурную формулу трапеций Эйлера-Маклорена, вычислите интеграл intlimits_{-1}^{1} dfrac{dx}{x^2 + 1} с точностью 1%, используя сетку с четырьмя интервалами (N = 4). елси бы для вычисления интеграла применяли формулу Рунге-Ромберга не однократно, а рекуррентно. Т.е. после очередного применения этой формулы, в очередной раз повышая порядок точности схемы.
программа уже есть надо поменять.