[Matlab,매트랩] 2계 비제차 미분방정식 라플라스 변환 소스

>> clc; clearvars;
syms  t s Y y(t) Dy(t)
Dy = diff(y, t);
D2y = diff(Dy, t);
LS = 18*exp(-t)*sin(3*t);
EQN=D2y-Dy-2*y-LS;
LEQN=laplace(EQN,t,s);
LT_Y=subs(LEQN,laplace(y,t,s),Y);
LT_Y=subs(LT_Y, {'y(0)', 'D(y)(0)'}, {0, 3});
ys=solve(LT_Y,Y);
y=ilaplace(ys,s,t)
ezplot(y); shg