Решение задачи Коши для обыкновенных дифференциальных уравнений, страница 7

Вычисляя значения yi приближенными методами, занесем в таблицу 1 максимальные погрешности каждого метода при различных значениях M (10, 20, 40).

Таблица 1

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

Метод

М

Max погрешность |y(xi) - yi|

Порядок точности

Метод Эйлера

10

0,187377344

1

20

0,10798817

40

0,0581612

Модифицированный метод Эйлера (предиктор-корректор)

10

0,015824958

2

20

0,004282

40

0,00111125

Метод Рунге-Кутты

10

0,00055952

3

20

7,3E-05

40

9,27699E-06

    Построим график на отрезке [1,3] функции точного решения (y(x) = 4/x) и вынесем точки (xi, yi), соответствующие массивам приближенных решений, полученных по приближенным методам для M = 10:

       Вывод: из таблицы 1 видно, что в сравнении с точным решением самым точным является метод Рунге-Кутты; при увеличении M погрешность данного метода уменьшается в 9 раз, что говорит о том, что приближенный метод имеет третий порядок точности. Погрешность модифицированного метода Эйлера (предиктор-корректор) при увеличении M уменьшается в 4 раза, из этого следует, что метод имеет второй порядок точности. Метод Эйлера из перечисленных приближенных методов является менее точным, его погрешность с увеличением M уменьшается в 2 раза, метод имеет первый порядок точности. Сделанные выводы подтверждает и график, который показывает, что ломанные, соответствующие методу Рунге-Кутты (третий порядок аппроксимации) и методу предиктор-корректора (второй порядок аппроксимации), проходят ближе к графику точного решения, чем ломанная, полученная по методу Эйлера первого порядка аппроксимации. Наиболее ближе к графику точного решения проходит ломанная метода Рунге-Кутты, что еще раз свидетельствует об его третьем порядке аппроксимации.