7.22 Решение дифференциальных уравнений с помощью степенных рядов

Пример:

Решить дифференциальное уравнение y' = cosx+e^y с начальным условием y(0) = 1.

Решение:

Обращение к dsolveвозвращает сообщение о том, что решение не найдено:

>> dsolve('Dy=cos(x)+exp(y)','y(0)=1','x')

ans =

[ empty sym ]

Команда dsolve не нашла аналитического решения в MATLAB. Известно, что решения этого дифференциального уравнения в аналитическом виде не существует. Найти разложение решения в степенной ряд (до 6-й степени по умолчанию) можно с помощью команды dsolve системы Maple.

>> maple('dsolve({diff(y(x),x)=cos(x)+exp(y(x)),y(0)=1},y(x),series)')

ans =

y(x) = series(1+(1+exp(1))*x+(1/2*exp(1)*(1+exp(1)))*x^2+(-1/6+1/3*exp(1)*(3/2*exp(1)+exp(1)^2+1/2))*x^3+(1/2*exp(1)^3+1/4*exp(1)^4+7/24*exp(1)^2)*x^4+(1/2*exp(1)^4+1/5*exp(1)^5+5/12*exp(1)^3+1/120-1/40*exp(1)+1/12*exp(1)^2)*x^5+O(x^6),x,6)

Имеется возможность управлять порядком разложения. Найдем разложение решения в степенной ряд до 3 -й степени:

>> maple('Order:=3;dsolve({diff(y(x),x)=cos(x)+exp(y(x)),y(0)=1},y(x),series)')

ans =

Order := 3y(x) = series(1+(1+exp(1))*x+(1/2*exp(1)*(1+exp(1)))*x^2+O(x^3),x,3)