Графическое представление решений
Для отображения численного решения задачи Коши для дифференциального уравнения используется команда odeplot(). Её первым параметром является построенная командой dsolve() с опцией type=numeric процедура численного решения задачи Коши, а вторым списковым параметром задаются отображаемые по горизонтальной и вертикальной осям величины. Например, чтобы отобразить искомую функцию решения y(x), следует задать список [x,y(x)].
Нарисуем график точного решения:
Пример 7.
>
>
>
>
> assign(y60):y:=evalf(y(x));
> py60:=plot(y,x=-1..1,legend=`zadacha Cauchy`,color=blue,thickness=3):
plots[display]({py60});
Здесь решению присвоено имя у 60, а соответствующему графику – ру 60.
Теперь построим график решения краевой задачи.
Пример 8.
>
>
>
>
>
>
> assign(y71):y:=evalf(y(x));
> py71:=plot(y,x=-1..1,legend=`krajevaya zadacha`,color=brown,thickness=3):
plots[display]({py71});
Графическое представление приближённых решений осуществляется с помощью функции odeplot.
1 | 2 | 3 | 4 | 5 | Поиск по сайту:
|