% çàäàâàìå äàäåíèòå óñëîâèÿ çà çàäà÷àòà, êîèòî ñà 
% ó''' + 7ó'' + 2ó' + ó = 0
% ó(0) = -2
% ó'(0) = 12
% ó''(0) = -3

sol = dsolve('D3y + 7*D2y + 2*Dy + y = 0,y(0)= -2,Dy(0) = 12 , D2y(0) = -3');
x = -2:0.1:0;             % -2 < õ < 0 ñúñ ñòúïêà 0.1
y = subs(sol, 't' , x);   % ðåøåíèå ñïðÿìî t
plot(x,y)