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

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

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

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

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

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

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

[an error occurred while processing this directive]

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

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

 y

 M₂

 M₁ M₃

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

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

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

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

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

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

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

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

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

  [an error occurred while processing this directive]

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

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

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

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

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