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

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

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

  Рассмотрим некоторые из них.

 

Метод Эйлера.

(Леонард Эйлер (1707 – 1783) швейцарский математик )

Предел монотонной функции. Математика лекции и задачи

  Известно, что уравнение  задает в некоторой области поле направлений. Решение этого уравнения с некоторыми начальными условиями дает кривую, которая касается поля направлений в любой точке.

  Если взять последовательность точек х₀, х₁, х₂, …. и заменить на получившихся отрезках интегральную кривую на отрезки касательных к ней, то получим ломаную линию. 

 y

 

 M₂

 M₁ M₃

 

 При подстановке заданных начальных условий (х₀, у₀) в дифференциальное уравнение получаем угловой коэффициент касательной к интегральной кривой в начальной точке

  Заменив на отрезке [x₀, x₁] интегральную кривую на касательную к ней, получаем значение

  Производя аналогичную операцию для отрезка [x₁, x₂], получаем:

  Продолжая подобные действия далее, получаем ломаную кривую, которая называется ломаной Эйлера.

  Можно записать общую формулу вычислений:

 

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

 

Следует отметить, что точность метода Эйлера относительно невысока. Увеличить точность можно, конечно, уменьшив шаг вычислений, однако, это приведет к усложнению расчетов. Поэтому на практике применяется так называемый уточненный метод Эйлера или формула пересчета.

Суть метода состоит в том, что в формуле  вместо значения

 берется среднее арифметическое значений f(x₀, y₀) и f(x₁, y₁). Тогда уточненное значение:

  Пределы Интеграл Типовые задачи С помощью дифференциала функции вычислить приближённо   при x = 7,76. Справочный материал и примеры к выполнению контрольной работы по математике

 Затем находится значение производной в точке . Заменяя f(x₀, y₀) средним арифметическим значений f(x₀, y₀) и , находят второе уточненное значение у₁.

 Затем третье:

и т.д. пока два последовательных уточненных значения не совпадут в пределах заданной степени точности. Тогда это значение принимается за ординату точки М₁ ломаной Эйлера.

  Аналогичная операция производится для остальных значений у.

 

Подобное уточнение позволяет существенно повысить точность результата.