问题出在哪？循环后数值越来越大

A=[0,1;-4,0];
K=inv(A);
C=[1,1;1,1];
B=[0;1];
syms dt p;
y0=[0;0];
for T=0:0.01:0.2
y=exp(Adt)y0+K[exp(Adt)-C]Bp;
m=subs(y,{dt,p},{0.01,sin(4*T)});
y0=m;
vpa(m,5)
end