Составить подпрограмму-функцию вычисления значений интеграла методом Симпсон.

Не знаю как кого, а меня лично просто поражает вся абсурдность такого вот задания. И дело тут не в Фортране как таковом, я всё равно на нём ничего здесь писать не буду, а в отсутствии всякого смысла в таком расчёте: очень простой полином на очень простом сегменте интегрировать по Симпсону с точностью ε=1е-4 или ε=1е-5 ! Да какой даун это придумал ? А какой олигофрен в формулу для расчёта интегральной суммы вписал вот эту чушь: ⅓•N•(...), где N -чётное число? Он что -вообще, что-ли, не понимает что означает интеграл и как он приблизительно (a иногда хоть и приближёнными методами, но совершенно точно) вычисляется ?
Запустим на плюсах код:
#include <iostream>
using namespace std;
double Simpson(double (*Y)(double),
double a, double b, int n) { double i,
step = (b - a) / n, halfstep = step / 2,
u = Y(b - halfstep), v = 0., x; for (i = 1; i < n;
i++) { x = a + i * step; u += Y(x - halfstep);
v += Y(x); } return (Y(a) + Y(b) + 4. * u + 2. * v)
* step / 6.; }
double F(double x)
{ return -3. + x * (2. + x * x); }
int main()
{ int N; double s, S = 166.25;
cout.precision(16); while (true)
{ cout << "N = "; cin >> N;
s = Simpson(F, 0., 5., N); cout << s
<< ", error = " << s - S << endl; } }
Уже при количестве интервалов разбиения N=1 отрезка [0;5] ошибка получается нулевая -вот насколько этот метод точный ! Всё понятно ?

" а в отсутствии всякого смысла в таком расчёте: очень простой полином на очень простом сегменте интегрировать по Симпсону с точностью ε=1е-4 или ε=1е-5 ! Да какой даун это придумал ?"
Тупость задания соответствует тупости студентов. Посмотрите внимательно на вопросы в разделе программирования.

тут что-то в условиях напутано. или опечатка