...(2x2)) to resolve system (33x33)
mt=[0 0 30 0;0,0,0,50;30,0,0,0;0,50,0,80];
kt=[-40,0,0,0;0,-50,0,0;0,0,35000,-25000;0,0,-25000, 4000];
Z=inv(mt)*kt;
[V,D]=eig(Z);
disp('Eigenvalues')
DS=[D(1,1),D(2,2),D(3,3),D(4,4)]
disp('Eigenvectors')
V
x0=[0;0;0.01;0];
S=inv(V)*x0...