救助：请大家帮忙看看程序哪里出错

412

这是根据单自由度振动系统写的，但是程序运行老出错。请大家帮忙看看
m1=3;k1=3;c1=0.6;y0=[0,1]'
M=diag([m1]);
K=[k1];
C=[c1];
[m,n]=size(M);
A=[zeros(m),eye(m);-K/M,-C/M];
T=50;dt=0.01;
for i=1:round(T/dt)+1
    t(i)=dt*(i-1);
    y(:,i)=expm(A*t(i))*y0;
end
subplot(2,1,1),plot(t,y(1,:))
title('波形图')
subplot(2,1,2),plot(y(1,:),y(2,:))% 相图
title('相图')
[p,d]=eig(A)
for i=1:4
    f=p(i,:)
end

ChaChing

报错不是告诉LZ了吗?!
Index exceeds matrix dimensions.
p为2*2矩阵, for i=1:4, f=p(i,:); end有误!