Оценка погрешности по правилу Рунге.

Функция, реализующая явный метод Эйлера; возвращает вектор решения:

Входные параметры:

f - функция правой части;

y0 - начальное значение;

t0 - начальная точка отрезка;

h - шаг сетки;

N - число узлов сетки.

Правая часть уравнения:
Начальное значение: Отрезок:   
Шаг сетки: Число узлов сетки:   

ТОЧНОЕ РЕШЕНИЕ

Точное решение в узлах сетки:      

ВЫЧИСЛЕНИЕ РЕШЕНИЯ ПО МЕТОДУ ЭЙЛЕРА:

Функция, реализующая явный метод Эйлера; возвращает вектор решения:

      

ВЫЧИСЛЕНИЕ ПОГРЕШНОСТИ ПО ПРАВИЛУ РУНГЕ:

Вычисление приближенных решений с шагом h/2:

       

Вычисление погрешностей:

     

Графики приближенных и точного решений: