Интегрирование с помощью встроенных функций MathCad

Результат численного интегрирования с помощью встроенных функций MathCad зависит от точности вычислений TOL. По умолчанию TOL=0.001. Эту константу можно изменять программно, объявив, например TOL:=0.0001, а можно на соответствующей панели (Математика – Опции – TOL) изменить ее численное значение. Интегрирование, как и многие другие математические действия в MathCad устроено по принципу «как пишется, так и вводится». Знак интеграла (определенный и неопределенный) можно найти на панели исчислений. Останется только заполнить несколько местозаменителей – ввести подынтегральную функцию, пределы интегрирования и переменную интегрирования. Чтобы получить результат интегрирования, следует ввести знак численного равенства. Символьный процессор MathCad может для небольшого круга несложных подынтегральных функций провести точное интегрирование. Для этого вместо знака = надо ввести знак символьных вычислений . Несколько примеров интегрирования приведены в программе 19.

Программа 19

Изменение точности TOL вычислений не приводит к изменению количества значащих цифр в результате. Чтобы изменить формат вывода результата надо обратиться к панели Форматирование – результат… – number of desimal places.

Разработчиками MathCad запрограммированы 4 численных метода интегрирования: Ромберга (Romberg), Адаптивный (Adaptive), Бесконечный предел (Infinite limit), C неопределенностью (Singular Endpoint). Если щелкнуть правой кнопкой мыши по знаку интеграла, то в появившемся окне можно найти отметку какой метод использован по умолчанию. По желанию можно изменить метод расчета. Последнее лучше не делать, оставить выбор за MathCad, т. е. оставить флаг на AutoSelect.

Основные идеи, положенные в основу итерационного алгоритма Ромберга, который применяется для большинства случаев, следующий.

1. Сначала строится несколько интерполирующих полиномов, которые заменяют на интервале интегрирования подынтегральную функцию f(x): прямую, параболу, и т. д.

2. Находят ряд значений интеграла I₁,I₂,I₃… по правилу трапеций, Симпсона, и т. д.

3. Экстраполяцией этой зависимости к нулевому шагу находят J‑тое значение вычисляемого интеграла.

4. Осуществляют переход к новой итерации с помощью более частого разбиения интервала интегрирования и добавлением нового члена интерполирующего полинома (J_N – приближение Ромберга).

5. Выход из итерационного процесса происходит тогда, когда разница между двумя последующими приближениями по абсолютной величине будет меньше TOL.