m1=0.999
m2=1.01
m3=1
m4=1.02
m5=0.998
k11=1
k22=1
k33=0.9
k44=1
k55=1
b11=10.2
b22=10.3
b33=10.1
b44=10.1
b55=10.5
k12=10.6
k23=10.5
k34=11.7
k45=10.1
k51=10
b12=0.1
b23=0.2
b34=0.5
b45=0.0
b51=0.3
freq=0.4
t=0:0.05:100/freq
u=(1/5)*sin(freq*t)
F=[0 1 0 0 0 0 0 0 0 0;
-(1/m1)*(k12+k51+k11) (-1/m1)*(b12+b51+b11) k12/m1 b12/m1 0 0 0 0 k51/m1 b51/m1;
0 0 0 1 0 0 0 0 0 0; 
k12/m2 b12/m2 (-1/m2)*(k12+k23+k22) (-1/m2)*(b12+b23+b22) k23/m2 b23/m2 0 0 0 0; 
0 0 0 0 0 1 0 0 0 0;
0 0 k23/m3 b23/m3 (-1/m3)*(k23+k34+k33) (-1/m3)*(b23+b34+b33) k34/m3 b34/m3 0 0;
0 0 0 0 0 0 0 1 0 0;
0 0 0 0 k34/m4 b34/m4 (-1/m4)*(k34+k45+k44) (-1/m4)*(b34+b45+b44) k45/m4 b45/m4;
0 0 0 0 0 0 0 0 0 1;
k51/m5 b51/m5 0 0 0 0 k45/m5 b45/m5 (-1/m5)*(k45+k51+k55) (-1/m5)*(b45+b51+b55)]
G=[1/m1; 0; 1/m2; 0; 1/m3; 0; 1/m4; 0; 1/m5; 0]
H=[1,0,1,0,1,0,1,0,1,0]
H1=[1 0 0 0 0 0 0 0 0 0]
H2=[0 0 1 0 0 0 0 0 0 0]
H3=[0 0 0 0 1 0 0 0 0 0]
H4=[0 0 0 0 0 0 1 0 0 0]
H5=[0 0 0 0 0 0 0 0 1 0]
J =0 
ss_blades=ss(F,G,H,J); [F,G,H,J]=ssdata(ss_blades);
[blades_out,blades_ss]=lsim(F,G,H1,J,u,t);
lsim(F,G,H1,J,u,t);
X1=lsim(F,G,H1,J,u,t);
X2=lsim(F,G,H2,J,u,t);
X3=lsim(F,G,H3,J,u,t);
X4=lsim(F,G,H4,J,u,t);
X5=lsim(F,G,H5,J,u,t);