системы дифференциальных уравнений (СДУ) в нормальной форме и коэффициенты СДУв виде векторов а и b, которые можно менять в интерактивном режиме и всякий раз получать решение для анализа и принятия рекомендаций по оптимизации его поиска:
0<FRAME<50
Запишем правую часть СДУ под именем D(t,X) в следующем виде:
задав начальные условия:
Запишем процедуру нахождения решения:
Приведем матрицу решений Z с числовыми данными
В полученной матрице решенияZ<1> aналог времени , заданный в виде вектора.
Прочие векторы-столбцы Z<j> представляют решения СДУ (1).
Между введенными обозначениями установлены следующие соответствия:
которые позволяют построить графики найденных решений (см. рис.1-2 и др.)
Нахождение решения с помощью более мощной процедуры rkfixed
Проведем решение СДУ (1) с помощью опции Radau(y0,t0,t1,N,D), записав её в виде:
Полученное решение СДУ (1) методом RADAUS в данном случае ничем не отличается от полученных ранее решений с использованием других опций. Это своего рода проверка, которая подтверждает достоверность полученных решений.