Графическое представление решений дифференциальных уравнений
Применение функции odeplot пакета plots
Для обычного графического представления результатов решения дифференциальных уравнений может использоваться функция odeplot из описанного выше пакета plots. Эта функция используется в следующем виде:
odeplot(s,vars.r,o)
Где s – запись (в выходной фирме) дифференциального уравнения или системы дифференциальных уравнений, решаемых численно функцией dsolve, vars – переменные, г – параметр, задающий пределы решения (например, а..b), и о – необязательные дополнительные опции.
На рис. 13.5 представлен пример решения одиночного дифференциального уравнения с выводом решения у(х) с помощью функции odeplot.
Рис. 13.5. Пример решения одиночного дифференциального уравнения
В этом примере решается дифференциальное уравнение: y'(x)=cos(x2 y(x)).
При у(0) = 2 и x, меняющемся от -5 до 5. Левая часть уравнения записана с помощью функции вычисления производной diff. Результатом построения является график решения у(х).