Формула Ньютона (правило трех восьмых):

где

 

Остаточный член имеет вид

 

Заметим, что в формуле Ньютона число узлов обязательно равно , т.е. .

 

 

Рис. 4.

 

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

 

 

где под подразумевается арифметическое среднее значение разностей соответствующего порядка.

Система MathCad вычисляет определенные интегралы методом Ромберга. Не описывая его подробно, отметим лишь, что он является вариантом метода трапеций с делением на два интервала интегрирования с итерационным уточнением решения до достижения заданной точности (она определяется значением системной переменной TOL).

Если за заданное число итераций точность не достигнута, используется более точный метод Ромберга с открытыми концами. При нем число интервалов утраивается на каждом шаге интегрирования. Этот метод увеличивает число шагов интегрирования там, где подынтегральная функция меняется более резко (например, если она имеет разрыв).

К достоинствам метода можно отнести то, что он делает все возможное, чтобы вычислить интеграл даже сложной функции. Но для простых функций это ведет к увеличению времени вычислений. При наличии у подынтегральной функции особенностей время вычисления может резко возрастать из-за перехода от одной реализации метода Ромберга к другой. Поэтому нередко бывает оправданным применение достаточно точных формул интегрирования, например формул Ньютона-Котесса с легко предсказуемыми узлами, которые можно выбрать вдали от особых точек подынтегральной функции.

Для вычисления определенных интегралов в Mathcad нажмите & или вставьте знак интегрирования с соответствующей панели.

Функция возвращает определенный интеграл f(x) от а до b, где F – любая скалярная функция, определенная в замкнутом интервале [a, b]; x – аргумент этой скалярной функции. а и b должны быть действительными скалярными величинами, но f(х) может быть комплексной величиной. а и b должны иметь одинаковые размерности, если они есть. может быть функцией любого числа переменных.

Подобно всем численным методам, алгоритм интегрирования Mathcad может иметь трудности с неправильными подынтегральными выражениями. Если подынтегральное выражение имеет особенности, разрывы, или большую и быструю флуктуацию, решение Mathcad может быть неточно.

Чтобы вычислить двойные или кратные интегралы, нажмите & несколько раз.

В следующем примере вычисляется определенный интеграл действительной функции f(x,y) в некоторой плоской области. Вводятся границы области интегрирования, принимаются а < x < b и с(x)<y<d(x) для всех x:

Вводится подынтегральная функция, описывающая плотность треугольника:

Тогда его масса вычисляется следующим образом:

3. ПРИБЛИЖЕННЫЕ МЕТОДЫ РЕШЕНИЯ
ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ

Задача Коши для дифференциального уравнения n- го порядка

,

заключается в отыскании функции удовлетворяющей этому уравнению и начальным условиям:

 

где - заданные числа.

Найти общее решение задачи Коши удается в редких случаях, поэтому чаще всего приходится решать задачу Коши приближенно. Приближенные методы в зависимости от формы, в которой они представляют решение, можно разделить на две группы.

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

2) Численные методы, дающие приближенное решение в виде таблицы.

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

К наиболее часто используемым численным методам решения дифференциальных уравнений относятся методы Эйлера и Рунге-Кутта, имеющие, в свою очередь, несколько модификаций.

Метод Эйлера относится к численным методам, дающим решение в виде таблицы приближенных значений искомой функции y(x).

Пусть дано дифференциальное уравнение

(1)

с начальным условием . (2)

Выбрав достаточно малый шаг h, рассмотрим систему равноотстоящих точек

В методе Эйлера приближенные значения вычисляются последовательно по формулам

.

При этом искомая интегральная кривая , проходящая через точку , заменяется ломаной с вершинами ; каждое звено этой ломаной, называемой ломаной Эйлера, имеет направление, совпадающее с направлением той интегральной кривой уравнения (1), которая проходит через точку

Если правая часть уравнения (1) в некотором прямоугольнике удовлетворяет условиям

(3)

то имеет место следующая оценка погрешности:

где - значение точного решения уравнения при а - приближенное значение, полученное на n- м шаге.

Формула (3) имеет лишь теоретическое применение. На практике иногда оказывается более удобным двойной просчет: расчет повторяют с шагом и погрешность более точного значения (при шаге ) оценивают приближенно так: . Блок-схема алгоритма решения дифференциального уравнения методом Эйлера приведена на рис. 5.

Метод Эйлера легко распространяется на системы дифференциальных уравнений и на дифференциальные уравнения высших порядков. Последние должны быть предварительно приведены к системе дифференциальных уравнений первого порядка.

Рассмотрим систему двух уравнений первого порядка

с начальными условиями

 

 

Рис. 5.

 

Приближенные значения и вычисляются последовательно по формулам

Существует несколько модификаций метода Эйлера.

Первый улучшенный метод Эйлера для решения задачи (1), (2) состоит в том, что сначала вычисляют промежуточные значения, а затем полагают

Второй улучшенный метод – метод Эйлера–Коши -заключается в том, что сначала определяют "грубое приближение"

затем вычисляют и приближенно полагают

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

где - точное решение дифференциального уравнения.

Метод Эйлера-Коши решения задачи (1), (2) можно ещё более уточнить, применяя итерационную обработку каждого значения y₁ (метод Эйлера с последующей итерационной обработкой). Исходя из грубого приближения

рассмотрим итерационный процесс

Итерации продолжаем до тех пор, пока в двух последовательных приближениях не совпадут соответствующие десятичные знаки. После этого полагаем

Как правило, при достаточно малом h итерации быстро сходятся. Если после трёх–четырёх итераций не произошло совпадения нужного числа десятичных знаков, то следует уменьшить шаг расчета h.

Блок-схемы модифицированных методов Эйлера легко получить самостоятельно, по аналогии с рис. 5.

При решении методом Рунге-Кутта дифференциального уравнения (1) с начальными условиями (2) через обозначают приближенное значение искомого решения в точке x₁ и вычисление приближенного значения в следующей точке производится по формулам

где

Схема метода Рунге-Кутта приведена в таблице 1.

Таблица 1.

Порядок заполнения таблицы (выполнения вычислений по методу Рунге-Кутта):

1. Выбираются x и y.

2. Вычисляются . Определяются

3. Определяются

4. Вычисляются и .

5. Принимаются

6. Вычисляются

7. Определяются

8. Вычисляются и .

9. Суммируются , делим на 6 и получаем таким образом

10. Вычисляются

Затем все вычисления продолжаются в том же порядке, принимая за начальную точку .

Заметим, что шаг расчета можно менять при переходе от одной точки к другой. Для контроля правильности выбора шага h рекомендуется вычислять дробь

.

Величина не должна превышать нескольких сотых. В противном случае шаг h следует уменьшить.

Метод Рунге-Кутта имеет порядок точности на всем отрезке .Оценка погрешности метода очень затруднительна. Грубую оценку погрешности можно получить с помощью двойного просчета по формуле

_,

где - значение точного решения уравнения в точке а - приближенные значения, полученные с шагом и .

При реализации на ЭВМ метода Рунге-Кутта на ЭВМ с автоматическим выбором шага обычно в каждой точке делают двойной просчет - сначала с шагом h, затем с шагом .Если полученные при этом значения различаются в пределах допустимой точности, то шаг для следующей точки удваивают, в противном случае берут половинный шаг.

Блок-схему решения обыкновенного дифференциального уравнения методом Рунге-кутта с автоматическим выбором шага можно получить самостоятельно.