Что считает калькулятор
Определённый интеграл ∫f(x)dx от a до b равен знаковой площади между графиком функции и осью абсцисс. Калькулятор берёт функцию f(x), границы a и b и число разбиений n, а на экран выдаёт значение интеграла, сравнение с методом трапеций, оценку погрешности и график функции с заштрихованной площадью.
Как заполнить поля. В строку «Функция f(x)» впишите выражение на JavaScript: x**2, sin(x), sqrt(x), exp(x), ln(x). Степень записывается как x**2 или x^2. В полях «Нижний предел a» и «Верхний предел b» задайте границы, в списке «Число разбиений n» выберите число отрезков. Вертикальные маркеры на графике двигаются мышью и стрелками ←/→, после чего пределы и результат пересчитываются.
Откуда берётся число
Отрезок [a; b] делится на n равных частей, шаг h = (b − a) / n. Метод Симпсона заменяет функцию на каждой паре соседних отрезков квадратичной параболой и складывает площади под ними: S(n) = h/3 · [f(a) + f(b) + 4·Σf(нечётные) + 2·Σf(чётные)]. Для метода трапеций рядом считается T(n) = h · [ (f(a) + f(b))/2 + Σf(внутренние) ].
Погрешность оценивается правилом Ричардсона: то же значение считается с удвоенным числом разбиений и берётся err = |S(n) − S(2n)| / 15. Множитель 15 связан с тем, что ошибка Симпсона убывает как O(h⁴): при уменьшении шага вдвое она падает примерно в шестнадцать раз.
Пример, который можно повторить
f(x) = x², a = 0, b = 1, n = 1000. Шаг равен 0,001. Калькулятор показывает 0,33333333. Это ровно 1/3 с точностью отображения. Метод трапеций даёт 0,3333335, а оценка Ричардсона на этом примере близка к нулю: функция гладкая, и удвоение разбиений почти не меняет результат. Подынтегральная функция здесь неотрицательна, поэтому строка под результатом сообщает, что знак площади не меняется.
Второй пример: f(x) = x³, a = −1, b = 1. Интеграл равен нулю, но площадь не равна нулю, потому что часть графика лежит выше оси X, а часть ниже. Калькулятор разделяет эти области. Красная заливка показывает вклад выше оси X, голубая — ниже, а строка под результатом сообщает, что функция меняет знак.
Где модель не работает
Функция должна быть определена и конечна на всём отрезке [a; b]. Разрыв, деление на ноль или логарифм неположительного числа дают ошибку вычисления, а не молчаливый ноль. Число разбиений берётся чётным, нечётное значение округляется вверх, как требует формула Симпсона. Численное интегрирование приближённое: результат зависит от n и от гладкости функции, а не только от её значения в отдельных точках.
Составная формула Симпсона и правило Ричардсона — стандартные приёмы численного анализа; их вывод есть в курсах вычислительной математики, например A. Ralston, P. Rabinowitz, A First Course in Numerical Analysis.