MathCad 2000. Шейкер Т.Д. - 63 стр.

                    ⎛    −8 ⋅ Y0 + 8 ⋅ Y1      ⎞
                    ⎜                          ⎟
      D ( t , Y) :=
                    ⎜ 30 ⋅ Y0 + Y1 − Y0 ⋅ Y2   ⎟
                    ⎜              8           ⎟
                    ⎜ Y0 ⋅ Y1 − ⋅ Y2           ⎟
                    ⎝              3           ⎠
         Определим параметры, необходимые для получения
решения:
         t0, t1 – начальное и конечное значения независимой
переменной;
         Y0 – вектор начальных условий;
         N – количество шагов (при необходимости) на интервале
[t0,t1].
            ⎛⎜ −1 ⎟⎞
      Y0 := ⎜ 0 ⎟          t0 := 0         t1 := 10    N := 1000
             ⎜ 1 ⎟
             ⎝    ⎠
      Получим решение, вызвав одну из функций: rkfixed,
Bulstoer, Rkadapt:
               S := Rkadapt ( Y0 , t0 , t1 , N , D )
      Выделим из матрицы решения S независимую переменную и
функции y0, y1, y2, учитывая, что аргумент расположен в
начальном столбце (номер 0), а функции в последующих столбцах.

          〈〉               〈〉                     〈〉           〈〉
   t := S 0        y0 := S 1              y1 := S 2    y2 := S 3

       По найденным t, y0, y1, y2 можно построить графики.
Однако формирование этих векторов для построения графиков не
является обязательным. Для примера создадим график по двум
функциям решения, не выделяя их предварительно из матрицы
(рис. 7.).
       Следует понимать, что при реализации рассмотренного
метода решения системы дифференциальных уравнений, нет
необходимости записывать систему и начальные условия в
документе Mathcad, так как в вычислениях участвуют только
параметры, заданные в используемой функции.
                                     63