1.5 Исследование динамической системы с использованием пакета Mathcad

Характерной особенностью пакета является использование привычных стандартных математических обозначений, то есть документ на экране выглядит точно так же обычный математический расчет.

Определяем параметры системы, время и количество точек, в которых будет искаться решение: а:=0Л      6:= 0.2      n:=800     i:=0..n Р:=0.05    е:=0.15    у:-0.03    t:=164

1

Решаем систему с помощью функции rkfixed, реализующую метод Рунге-Кутты. Первый параметр - начальные условия, второй и третий -диапазон изменения времени, последний — функция, задающая систему уравнений.

z := rkfixed(y, 0, t, n, D)              zl := rkfixed(y 1, 0, t, n, D)