11. Решение обыкновенных дифференциальных уравнений.

Функция для решения дифференциальных уравнений имеет вид:

         Z: =rkfixed (y, x1, x2, npoints, D), где у – вектор начальных значений решений

                                                        X1 – начало отрезка интегрирования

                                                        X2 – конец отрезка интегрирования

                                                        Npoints – число точек интегрирования

                                                        D – вектор правых частей

Число строк вектора D равно порядку уравнения. При этом дифференциальное уравнение записывается как бы в виде системы. На рис. 18  y1 - значение первой производной. Следовательно, первая строка вектора D – это первая производная, вторая строка – вторая производная.

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

Для построения графика найденной функции необходимо по оси абсцисс отложить первый столбец таблицы, по оси ординат – второй. Столбец указывается верхним индексом (горячие клавиши Ctrl + 6) (см. рис. 18).