Метод Эйлера

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

Метод Эйлера основан на разложении функции y в ряд Тейлора в окрестности x₀:

Если h мало, то члены, содержащие производные второго и более высоких порядков, можно отбросить. Тогда .
y’(x₀) находим из дифференциального уравнения, подставив в него начальное условие. Таким образом, можно получить приближённое значение y при малом смещении h от начальной точки x₀. Этот процесс можно продолжить, используя следующую рекуррентную формулу:

Рис.17.2. Принцип метода Эйлера

y_i+1=y_i + hf(x_i,y_i), i=1,2,…

Графически метод Эйлера показан на рис.17.2. Ошибка метода имеет порядок h².

Пример 1. Составим программу для решения дифференциального уравнения y¢=2x²+2y при начальном условии y(0)=1; 0£ x£1 и h=0,1.

Program Euler;

Uses Crt;

Var

xn,xk,yn,h,x,y:real;

i:integer;

Function f(x,y:real):real;

begin

{ Здесь приводим выражение для вычисления функции f(x,y) }

f:=2*x*x+2*y;

end;

Begin

ClrScr;

Writeln(' Решение дифференциального уравнения ');

Writeln(' dy/dx=2x^2+2y методом Эйлера ');

{ Ввод исходных данных }

xn:=0; yn:=1; xk:=1; h:=0.1;

{ Выводим шапку таблицы и первую точку }

Writeln('--------------------');

Writeln('| № | x | y |');

Writeln('--------------------');

{ Начинаем расчет }

x:=xn; y:=yn; i:=1;

Writeln('|', i:2, ' |', x:5:2, ' |', y:7:4, ' |');

Repeat

y:=y+h*f(x,y);

Writeln('|', i:2, ' |', x:5:2, ' |', y:7:4, ' |');

x:=x+h;

i:=i+1;

Until x>xk;

Writeln('--------------------');

Readln;

End.