ОГЛАВЛЕНИЕ
ВВЕДЕНИЕ ……………………………………………………………….…………3
ГЛАВА I. РЕШЕНИЕ ЗАДАЧИ КОШИ ДЛЯ СИСТЕМ ДИФФЕРНЕЦИАЛЬНЫХ УРАВНЕНИЙ МЕТОДОМ РУНГЕ-КУТТА-МЕРСОНА …………………………………………………………………..………4
1.1. Постановка задачи Коши для систем обыкновенных дифференциальных уравнений………….……………..................…...4
1.2. Метод Рунге-Кутта-Мерсона………………………………………….5
1.3. Алгоритм решения задачи Коши для систем обыкновенных дифференциальных уравнений методом Рунге-Кутта-Мерсона……………….........................................................………..….7
ГЛАВА II. ВЫЧИСЛИТЕЛЬНЫЙ ЭКСПЕРЕМЕНТ ………….………...…....8
2.1. Постановка задачи……………………………………………………..8
2.2. Анализ результатов……………………………………………….....…9
ЗАКЛЮЧЕНИЕ ………………………………………………………………..….10 СПИСОК ЛИТЕРАТУРЫ ………...……………..................................................11
ПРИЛОЖЕНИЕ.......................................................................................................12
Приложение 1........................................................................................12
Приложение 2........................................................................................16
ВВЕДЕНИЕ
Обыкновенные дифференциальные уравнения широко используются для математического моделирования процессов и явлений в различных областях науки и техники. Множество переходных процессов в радиотехнике, кинетика химических реакций, динамика биологических популяций, движение космических объектов, модели экономического развития исследуются с помощью обыкновенных дифференциальных уравнений.
Актуальность темы курсовой работы состоит в том, что ОДУ имеют аналитически сложное решение и составление программы, реализируещей численное решение облегчило бы эту задачу.
В курсовой работе решается задача разработки программы поиска решения системы дифференциальных уравнений методам Рунге-Кутта-Мерсона.
Выбор метода решения системы дифференциальных уравнений объясняется тем, что метод Кутта-Мерсона сочетает хорошую точность и высокую скорость.
Цель работы: составить программу для решения задачи Коши для системы дифференциальных уравнений методом Рунге-Кутта-Мерсона на примере, проверить полученное решение в MathCad и проанализировать результаты.
ГЛАВА I. РЕШЕНИЕ ЗАДАЧИ КОШИ ДЛЯ СИСТЕМ ДИФФЕРНЕЦИАЛЬНЫХ УРАВНЕНИЙ МЕТОДОМ РУНГЕ-КУТТА-МЕРСОНА
Постановка задачи Коши для систем обыкновенных дифференциальных уравнений
Задача Коши заключается в решении систем обыкновенных дифференциальных уравнений первого порядка, представляемых в виде:
(1)
Где j=1 
N-номер каждой зависимой переменной yj, x-независимая переменная.
Решение системы (1) при заданных начальных условиях
x=x0,
y1(x0)=y10,…, (2)
y2(x0)=y20,
yN(x0)=yN0
сводиться к нахождению зависимостей (интегральных кривых ) y1(x),…,y2(x), yN(x), проходящих через точки (x0,y10), (x0,y20),…, (x0,yN0). Задача Коши сводиться к интегрированию дифференциальных уравнений. Порядок метода численного интегрирования при этом определяется и порядок метода решения.
Метод Рунге-Кутта-Мерсона
Автоматическое изменение шага в ходе решения систем дифференциальных уравнений необходимо, если решение требуется получить с заданной точностью. При высокой точности (погрешность 
) и решении в виде кривых с сильно различающейся крутизной автоматическое изменение шага обеспечивает уменьшение общего числа шагов в несколько раз, резко уменьшается вероятность числовой неустойчивости, даёт более равномерное расположение точек графика кривых (решений) при их выводе на печать. Данный метод обеспечивает приближённую оценку погрешностей на каждом шаге интегрирования. Погрешность интегрирования имеет порядок h5. Этот метод реализуется следующим алгоритмом: Задаём число уравнений N, погрешность ε=E, начальный шаг интегрирования h=H и начальное значение y10,…,yN0. С помощью пяти циклов с управляющей переменной J=1,2,..,N вычисляем коэффициенты:
(3)
(4)
(5)
(6)
(7)
Находим (в последнем цикле) значение (8):
(8)
И погрешность
(9)
Проверяем выполнения условий
(10)
(11)
Если условие (10) не выполняется, то делим шаг h на 2 и повторяем вычисления. Если это условие выполняется и выполняется условие (11), значение xi+1=xi+h и Yj(i+1), то считаем, что решение системы дифференциальных уравнений найдено с заданной точностью. Если условие (11) не выполняется, шаг h увеличивается вдвое и вычисления повторяются.