dz=[z(2); 2*z(4)*wc/wn+wc^2*z(1)/(wn)^2-ks*z(1)/(ms*wn^2)-cs*z(2)/(ms*wn)+ks1*a*sin(-20)/(ms*wn^ 2)+ks2*b*sin(100)/(ms*wn^2)+ks3*c*sin(220)/(ms*wn^2)+cs1*z(14)*sin(-20)/(ms*wn)+cs2*z( 16)*sin(100)/(ms*wn)+cs3*z(18)*sin(220)/(ms*wn)-qs1*u*ks1*a*cos(-20)/(ms*wn^2)-qs2*u*ks 2*b*cos(100)/(ms*wn^2)-qs3*u*ks3*c*cos(220)/(ms*wn^2)-qs1*u*cs1*z(14)*cos(-20)/(ms*wn) -qs2*u*cs2*z(16)*cos(100)/(ms*wn)-qs3*u*cs3*z(18)*cos(220)/(ms*wn); z(4); -2*wc*z(2)/wn+wc^2*z(3)/wn^2-ks*z(3)/(ms*wn^2)-cs*z(4)/(ms*wn)-ks1*a*cos(-20)/(ms*wn^ 2)-ks2*b*cos(100)/(ms*wn^2)-ks3*c*cos(220)/(ms*wn^2)-cs1*z(14)*cos(-20)/(ms*wn)-cs2*z(1 6)*cos(100)/(ms*wn)-cs3*z(18)*cos(220)/(ms*wn)-qs1*u*ks1*a*sin(-20)/(ms*wn^2)-qs2*u*ks2 *b*sin(100)/(ms*wn^2)-qs3*u*ks3*c*sin(220)/(ms*wn^2)-qs1*u*cs1*z(14)*sin(-20)/(ms*wn)-q s2*u*cs2*z(16)*sin(100)/(ms*wn)-qs3*u*cs3*z(18)*sin(220)/(ms*wn); z(6); 2*wc*z(8)/wn+(wc)^2*z(5)/(wn)^2-kc*z(5)/(mce*wn^2)-cc*z(6)/(mce*wn)-ks1*a*sin(-20)/(mce *wn^2)-ks2*b*sin(100)/(mce*wn^2)-ks3*c*sin(220)/(mce*wn^2)-cs1*z(14)*sin(-20)/(mce*wn)- cs2*z(16)*sin(100)/(mce*wn)-cs3*z(18)*sin(220)/(mce*wn)+kr1*d*sin(20)/(mce*wn^2)+kr2*e* sin(140)/(mce*wn^2)+kr3*f*sin(260)/(mce*wn^2)+cr1*z(20)*sin(20)/(mce*wn)+cr2*z(22)*sin( 140)/(mce*wn)+cr3*z(24)*sin(260)/(mce*wn)-qr1*u*kr1*d*cos(20)/(mce*wn^2)-qr2*u*kr2*e*c os(140)/(mce*wn^2)-qr3*u*kr3*f*cos(260)/(mce*wn^2)-qr1*u*cr1*z(20)*cos(20)/(mce*wn)-qr2 *u*cr2*z(22)*cos(140)/(mce*wn)-qr3*u*cr3*z(24)*cos(260)/(mce*wn)-qs1*u*ks1*a*cos(-20)/( mce*wn^2)-qs2*u*ks2*b*cos(100)/(mce*wn^2)-qs3*u*ks3*c*cos(220)/(mce*wn^2)-qs1*u*cs1* z(14)*cos(-20)/(mce*wn)-qs2*u*cs2*z(16)*cos(100)/(mce*wn)-qs3*u*cs3*z(18)*cos(220)/(mce *wn); z(8); -2*wc*z(6)/wn+(wc)^2*z(7)/(wn)^2-kc*z(7)/(mce*wn^2)-cc*z(8)/(mce*wn)-ks1*a*cos(-20)/(mc e*wn^2)-ks2*b*cos(100)/(mce*wn^2)-ks3*c*cos(220)/(mce*wn^2)-cs1*z(14)*cos(-20)/(mce*w n)-cs2*z(16)*cos(100)/(mce*wn)-cs3*z(18)*cos(220)/(mce*wn)+kr1*d*cos(20)/(mce*wn^2)+kr 2*e*cos(140)/(mce*wn^2)+kr3*f*cos(260)/(mce*wn^2)+cr1*z(20)*cos(20)/(mce*wn)+cr2*z(22) *cos(140)/(mce*wn)+cr3*z(24)*cos(260)/(mce*wn)-qs1*u*ks1*a*sin(-20)/(mce*wn^2)-qs2*u*k s2*b*sin(100)/(mce*wn^2)-qs3*u*ks3*c*sin(220)/(mce*wn^2)-qs1*u*cs1*z(14)*sin(-20)/(mce* wn)-qs2*u*cs2*z(16)*sin(100)/(mce*wn)-qs3*u*cs3*z(18)*sin(220)/(mce*wn)-qr1*u*kr1*d*sin (20)/(mce*wn^2)-qr2*u*kr2*e*sin(140)/(mce*wn^2)-qr3*u*kr3*f*sin(260)/(mce*wn^2)-qr1*u* cr1*z(20)*sin(20)/(mce*wn)-qr2*u*cr2*z(22)*sin(140)/(mce*wn)-qr3*u*cr3*z(24)*sin(260)/(mc e*wn); z(10); 2*wc*z(12)/wn+wc^2*z(9)/(wn)^2-kr*z(9)/(mr*wn^2)-cr*z(10)/(mr*wn)+kr1*d*sin(20)/(mr*wn ^2)+kr2*e*sin(140)/(mr*wn^2)+kr3*f*sin(260)/(mr*wn^2)+cr1*z(20)*sin(20)/(mr*wn)+cr2*z(2 2)*sin(140)/(mr*wn)+cr3*z(24)*sin(260)/(mr*wn)-qr1*u*kr1*d*cos(20)/(mr*wn^2)-qr2*u*kr2* e*cos(140)/(mr*wn^2)-qr3*u*kr3*f*cos(260)/(mr*wn^2)-qr1*u*cr1*z(20)*cos(20)/(ms*wn)-qr2 

*u*cr2*z(22)*cos(140)/(mr*wn)-qr3*u*cr3*z(24)*cos(260)/(mr*wn); z(12); -2*wc*z(10)/wn+wc^2*z(11)/(wn)^2-kr*z(11)/(mr*wn^2)-cr*z(12)/(mr*wn)-kr1*d*cos(20)/(mr* wn^2)-kr2*e*cos(140)/(mr*wn^2)-kr3*f*cos(260)/(mr*wn^2)-cr1*z(20)*cos(20)/(mr*wn)-cr2*z( 22)*cos(140)/(mr*wn)-cr3*z(24)*cos(260)/(mr*wn)-qr1*u*kr1*d*sin(20)/(mr*wn^2)-qr2*u*kr2 *e*sin(140)/(mr*wn^2)-qr3*u*kr3*f*sin(260)/(mr*wn^2)-qr1*u*cr1*z(20)*sin(20)/(ms*wn)-qr2 *u*cr2*z(22)*sin(140)/(mr*wn)-qr3*u*cr3*z(24)*sin(260)/(mr*wn); z(14); -2*wc*(z(4)-z(8))*sin(-20)/wn-wc^2*(z(1)-z(5))*sin(-20)/wn^2-kc*z(5)*sin(-20)/(mce*wn^2)+k c*z(7)*cos(-20)/(mce*wn^2)+ks*z(1)*sin(-20)/(ms*wn^2)-ks*z(3)*cos(-20)/(ms*wn^2)-2*wc*(z (2)-z(6))*cos(-20)/wn-wc^2*(z(7)-z(3))*cos(-20)/wn^2-(rbs^2/js+rbc^2/jce)*(ks1*a/wn^2+cs1*z( 14)/wn)-(rbs^2/js+rbc^2/jce)*(ks2*b/wn^2+cs2*z(16)/wn)-(rbs^2/js+rbc^2/jce)*(ks3*c/wn^2+cs3 *z(18)/wn)+rbc^2/jce*(kr1*d/wn^2+cr1*z(20)/wn)+rbc^2/jce*(kr2*e/wn^2+cr2*z(22)/wn)+rbc^2 /jce*(kr3*f/wn^2+cr3*z(24)/wn)-(mce+ms)*sin(-20)*sin(-20)/(mce*ms)*(ks1*a/wn^2+cs1*z(14) /wn)-(mce+ms)*sin(-20)*sin(100)/(mce*ms)*(ks2*b/wn^2+cs2*z(16)/wn)-(mce+ms)*sin(-20)*si n(220)/(mce*ms)*(ks3*c/wn^2+cs3*z(18)/wn)+sin(-20)*sin(20)/mce*(kr1*d/wn^2+cr1*z(20)/wn) +sin(-20)*sin(140)/mce*(kr2*e/wn^2+cr2*z(22)/wn)+sin(-20)*sin(260)/mce*(kr3*f/wn^2+cr3*z( 24)/wn)+(mce+ms)*cos(-20)*cos(-20)/(mce*ms)*(ks1*a/wn^2+cs1*z(14)/wn)+(mce+ms)*cos(-2 0)*cos(100)/(mce*ms)*(ks2*b/wn^2+cs2*z(16)/wn)+(mce+ms)*cos(-20)*cos(220)/(mce*ms)*(k s3*c/wn^2+cs3*z(18)/wn)-cos(-20)*cos(20)/mce*(kr1*d/wn^2+cr1*z(20)/wn)-cos(-20)*cos(140) /mce*(kr2*e/wn^2+cr2*z(22)/wn)-cos(-20)*cos(260)/mce*(kr3*f/wn^2+cr3*z(24)/wn)-rb1^2/j1* (ks1*a/wn^2+cs1*z(14)/wn)-rb1^2/j1*(kr1*d/wn^2+cr1*z(20)/wn)-rbc*tc/(jce*bc*wn^2)+rbs*ts/ (js*bc*wn^2)+we^2*0.74*cos(we*t-1.05)/(bc*wn^2)-sin(-20)*cc*z(6)/(mce*wn)+cos(-20)*cc*z( 8)/(mce*wn)+sin(-20)*cs*z(2)/(ms*wn)-cos(-20)*cs*z(4)/(ms*wn)+rbs/js*qs1*u*ls1*(ks1*a/wn^ 2+cs1*z(14)/wn)+rbs/js*qs2*u*ls2*(ks2*b/wn^2+cs2*z(16)/wn)+rbs/js*qs3*u*ls3*(ks3*c/wn^2 +cs3*z(18)/wn)+rb1/j1*qs1*u*ls1*(ks1*a/wn^2+cs1*z(14)/wn)+rb1/j1*qs2*u*ls2*(ks2*b/wn^2 +cs2*z(16)/wn)+rb1/j1*qs3*u*ls3*(ks3*c/wn^2+cs3*z(18)/wn)-rb1/j1*qr1*u*lr1*(kr1*d/wn^2+ cr1*z(20)/wn)-rb1/j1*qr2*u*lr2*(kr2*e/wn^2+cr2*z(22)/wn)-rb1/j1*qr3*u*lr3*(kr3*f/wn^2+cr3 *z(24)/wn)-rbc^2/jce*qs1*u*ls1*(ks1*a/wn^2+cs1*z(14)/wn)-rbc^2/jce*qs2*u*ls2*(ks2*b/wn^2 +cs2*z(16)/wn)-rbc^2/jce*qs3*u*ls3*(ks3*c/wn^2+cs3*z(18)/wn); z(16); -2*wc*(z(4)-z(8))*sin(100)/wn-wc^2*(z(1)-z(5))*sin(100)/wn^2-kc*z(5)*sin(100)/(mce*wn^2)+ kc*z(7)*cos(100)/(mce*wn^2)+ks*z(1)*sin(100)/(ms*wn^2)-ks*z(3)*cos(100)/(ms*wn^2)-2*wc *(z(2)-z(6))*cos(100)/wn-wc^2*(z(7)-z(3))*cos(100)/wn^2-(rbs^2/js+rbc^2/jce)*(ks1*a/wn^2+cs 1*z(14)/wn)-(rbs^2/js+rbc^2/jce)*(ks2*b/wn^2+cs2*z(16)/wn)-(rbs^2/js+rbc^2/jce)*(ks3*c/wn^2 +cs3*z(18)/wn)+rbc^2/jce*(kr1*d/wn^2+cr1*z(20)/wn)+rbc^2/jce*(kr2*e/wn^2+cr2*z(22)/wn)+r bc^2/jce*(kr3*f/wn^2+cr3*z(24)/wn)-(mce+ms)*sin(100)*sin(-20)/(mce*ms)*(ks1*a/wn^2+cs1* z(14)/wn)-(mce+ms)*sin(100)*sin(100)/(mce*ms)*(ks2*b/wn^2+cs2*z(16)/wn)-(mce+ms)*sin(1 00)*sin(220)/(mce*ms)*(ks3*c/wn^2+cs3*z(18)/wn)+sin(100)*sin(20)/mce*(kr1*d/wn^2+cr1*z( 20)/wn)+sin(100)*sin(140)/mce*(kr2*e/wn^2+cr2*z(22)/wn)+sin(100)*sin(260)/mce*(kr3*f/wn^ 2+cr3*z(24)/wn)+(mce+ms)*cos(100)*cos(-20)/(mce*ms)*(ks1*a/wn^2+cs1*z(14)/wn)+(mce+m 

s)*cos(100)*cos(100)/(mce*ms)*(ks2*b/wn^2+cs2*z(16)/wn)+(mce+ms)*cos(100)*cos(220)/(mc e*ms)*(ks3*c/wn^2+cs3*z(18)/wn)-cos(100)*cos(20)/mce*(kr1*d/wn^2+cr1*z(20)/wn)-cos(100) *cos(140)/mce*(kr2*e/wn^2+cr2*z(22)/wn)-cos(100)*cos(260)/mce*(kr3*f/wn^2+cr3*z(24)/wn) -rb2^2/j2*(ks2*b/wn^2+cs2*z(16)/wn)-rb2^2/j2*(kr2*e/wn^2+cr2*z(22)/wn)-rbc*tc/(jce*bc*wn^ 2)+rbs*ts/(js*bc*wn^2)+we^2*0.74*cos(we*t-1.05)/(bc*wn^2)-sin(100)*cc*z(6)/(mce*wn)+cos( 100)*cc*z(8)/(mce*wn)+sin(100)*cs*z(2)/(ms*wn)-cos(100)*cs*z(4)/(ms*wn)+rbs/js*qs1*u*ls1 *(ks1*a/wn^2+cs1*z(14)/wn)+rbs/js*qs2*u*ls2*(ks2*b/wn^2+cs2*z(16)/wn)+rbs/js*qs3*u*ls3*( ks3*c/wn^2+cs3*z(18)/wn)+rb2/j2*qs1*u*ls1*(ks1*a/wn^2+cs1*z(14)/wn)+rb2/j2*qs2*u*ls2*(k s2*b/wn^2+cs2*z(16)/wn)+rb2/j2*qs3*u*ls3*(ks3*c/wn^2+cs3*z(18)/wn)-rb2/j2*qr1*u*lr1*(kr1 *d/wn^2+cr1*z(20)/wn)-rb2/j2*qr2*u*lr2*(kr2*e/wn^2+cr2*z(22)/wn)-rb2/j2*qr3*u*lr3*(kr3*f/ wn^2+cr3*z(24)/wn)-rbc^2/jce*qs1*u*ls1*(ks1*a/wn^2+cs1*z(14)/wn)-rbc^2/jce*qs2*u*ls2*(ks 2*b/wn^2+cs2*z(16)/wn)-rbc^2/jce*qs3*u*ls3*(ks3*c/wn^2+cs3*z(18)/wn); z(18); -2*wc*(z(4)-z(8))*sin(220)/wn-wc^2*(z(1)-z(5))*sin(220)/wn^2-kc*z(5)*sin(220)/(mce*wn^2)+ kc*z(7)*cos(220)/(mce*wn^2)+ks*z(1)*sin(220)/(ms*wn^2)-ks*z(3)*cos(220)/(ms*wn^2)-2*wc *(z(2)-z(6))*cos(220)/wn-wc^2*(z(7)-z(3))*cos(220)/wn^2-(rbs^2/js+rbc^2/jce)*(ks1*a/wn^2+cs 1*z(14)/wn)-(rbs^2/js+rbc^2/jce)*(ks2*b/wn^2+cs2*z(16)/wn)-(rbs^2/js+rbc^2/jce)*(ks3*c/wn^2 +cs3*z(18)/wn)+rbc^2/jce*(kr1*d/wn^2+cr1*z(20)/wn)+rbc^2/jce*(kr2*e/wn^2+cr2*z(22)/wn)+r bc^2/jce*(kr3*f/wn^2+cr3*z(24)/wn)-(mce+ms)*sin(220)*sin(-20)/(mce*ms)*(ks1*a/wn^2+cs1* z(14)/wn)-(mce+ms)*sin(220)*sin(100)/(mce*ms)*(ks2*b/wn^2+cs2*z(16)/wn)-(mce+ms)*sin(2 20)*sin(220)/(mce*ms)*(ks3*c/wn^2+cs3*z(18)/wn)+sin(220)*sin(20)/mce*(kr1*d/wn^2+cr1*z( 20)/wn)+sin(220)*sin(140)/mce*(kr2*e/wn^2+cr2*z(22)/wn)+sin(220)*sin(260)/mce*(kr3*f/wn^ 2+cr3*z(24)/wn)+(mce+ms)*cos(220)*cos(-20)/(mce*ms)*(ks1*a/wn^2+cs1*z(14)/wn)+(mce+m s)*cos(220)*cos(100)/(mce*ms)*(ks2*b/wn^2+cs2*z(16)/wn)+(mce+ms)*cos(220)*cos(220)/(mc e*ms)*(ks3*c/wn^2+cs3*z(18)/wn)-cos(220)*cos(20)/mce*(kr1*d/wn^2+cr1*z(20)/wn)-cos(220) *cos(140)/mce*(kr2*e/wn^2+cr2*z(22)/wn)-cos(220)*cos(260)/mce*(kr3*f/wn^2+cr3*z(24)/wn) -rb3^2/j3*(ks3*c/wn^2+cs3*z(18)/wn)-rb3^2/j3*(kr3*f/wn^2+cr3*z(24)/wn)-rbc*tc/(jce*bc*wn^ 2)+rbs*ts/(js*bc*wn^2)+we^2*0.74*cos(we*t-1.05)/(bc*wn^2)-sin(220)*cc*z(6)/(mce*wn)+cos( 220)*cc*z(8)/(mce*wn)+sin(220)*cs*z(2)/(ms*wn)-cos(220)*cs*z(4)/(ms*wn)+rbs/js*qs1*u*ls1 *(ks1*a/wn^2+cs1*z(14)/wn)+rbs/js*qs2*u*ls2*(ks2*b/wn^2+cs2*z(16)/wn)+rbs/js*qs3*u*ls3*( ks3*c/wn^2+cs3*z(18)/wn)+rb3/j3*qs1*u*ls1*(ks1*a/wn^2+cs1*z(14)/wn)+rb3/j3*qs2*u*ls2*(k s2*b/wn^2+cs2*z(16)/wn)+rb3/j3*qs3*u*ls3*(ks3*c/wn^2+cs3*z(18)/wn)-rb3/j3*qr1*u*lr1*(kr1 *d/wn^2+cr1*z(20)/wn)-rb3/j3*qr2*u*lr2*(kr2*e/wn^2+cr2*z(22)/wn)-rb3/j3*qr3*u*lr3*(kr3*f/ wn^2+cr3*z(24)/wn)-rbc^2/jce*qs1*u*ls1*(ks1*a/wn^2+cs1*z(14)/wn)-rbc^2/jce*qs2*u*ls2*(ks 2*b/wn^2+cs2*z(16)/wn)-rbc^2/jce*qs3*u*ls3*(ks3*c/wn^2+cs3*z(18)/wn); z(20); 2*wc*(z(8)-z(12))*sin(20)/wn-wc^2*(z(5)-z(9))*sin(20)/wn^2+sin(20)*kr*z(9)/(mr*wn^2)-cos(2 0)*kr*z(11)/(mr*wn^2)+2*wc*(z(6)-z(10))*cos(20)/wn+wc^2*(z(11)-z(7))*cos(20)/wn^2-kc*z(5) *sin(20)/(mce*wn^2)+kc*z(7)*cos(20)/(mce*wn^2)-cos(20)*cos(-20)/mce*(ks1*a/wn^2+cs1*z(1 4)/wn)-cos(20)*cos(100)/mce*(ks2*b/wn^2+cs2*z(16)/wn)-cos(20)*cos(220)/mce*(ks3*c/wn^2+ cs3*z(18)/wn)+rbc^2/jce*(ks1*a/wn^2+cs1*z(14)/wn)+rbc^2/jce*(ks2*b/wn^2+cs2*z(16)/wn)+r bc^2/jce*(ks3*c/wn^2+cs3*z(18)/wn)-(rbr^2/jr+rbc^2/jce)*(kr1*d/wn^2+cr1*z(20)/wn)-(rbr^2/jr 

+rbc^2/jce)*(kr2*e/wn^2+cr2*z(22)/wn)-(rbr^2/jr+rbc^2/jce)*(kr3*f/wn^2+cr3*z(24)/wn)-rb1^2/ j1*(kr1*d/wn^2+cr1*z(20)/wn)-rb1^2/j1*(ks1*a/wn^2+cs1*z(14)/wn)-rbc*tc/(jce*bc*wn^2)+we ^2*2.33*cos(we*t-0.85)/(bc*wn^2)+(sin(20)*cr*z(10)/(mr*wn)-cos(20)*cr*z(12)/(mr*wn)-sin(20) *cc*z(6)/(mce*wn)+cos(20)*cc*z(8)/(mce*wn))*1-sin(20)*sin(-20)/mce*(ks1*a/wn^2+cs1*z(14) /wn)-sin(20)*sin(100)/mce*(ks2*b/wn^2+cs2*z(16)/wn)-sin(20)*sin(220)/mce*(ks3*c/wn^2+cs3 *z(18)/wn)+(mce+mr)*sin(20)*sin(20)/(mce*mr)*(kr1*d/wn^2+cr1*z(20)/wn)+(mce+mr)*sin(20) *sin(140)/(mce*mr)*(kr2*e/wn^2+cr2*z(22)/wn)+(mce+mr)*sin(20)*sin(260)/(mce*mr)*(kr3*f/ wn^2+cr3*z(24)/wn)+(mce+mr)*cos(20)*cos(20)/(mce*mr)*(kr1*d/wn^2+cr1*z(20)/wn)+(mce+ mr)*cos(20)*cos(140)/(mce*mr)*(kr2*e/wn^2+cr2*z(22)/wn)+(mce+mr)*cos(20)*cos(260)/(mce *mr)*(kr3*f/wn^2+cr3*z(24)/wn)+2*(mce+mr)/(mce*mr)*qr1*u*(kr1*d/wn^2+cr1*z(20)/wn)*si n(20)*cos(20)+2*(mce+mr)/(mce*mr)*qr2*u*(kr2*e/wn^2+cr2*z(22)/wn)*sin(140)*cos(140)+2* (mce+mr)/(mce*mr)*qr3*u*(kr3*f/wn^2+cr3*z(24)/wn)*sin(260)*cos(260)+rbr/jr*qr1*u*lr1*(kr 1*d/wn^2+cr1*z(20)/wn)+rbr/jr*qr2*u*lr2*(kr2*e/wn^2+cr2*z(22)/wn)+rbr/jr*qr3*u*lr3*(kr3*f/ wn^2+cr3*z(24)/wn)-rb1/j1*qs1*u*ls1*(ks1*a/wn^2+cs1*z(14)/wn)-rb1/j1*qs2*u*ls2*(ks2*b/w n^2+cs2*z(16)/wn)-rb1/j1*qs3*u*ls3*(ks3*c/wn^2+cs3*z(18)/wn)-rb1/j1*qr1*u*lr1*(kr1*d/wn^ 2+cr1*z(20)/wn)-rb1/j1*qr2*u*lr2*(kr2*e/wn^2+cr2*z(22)/wn)-rb1/j1*qr3*u*lr3*(kr3*f/wn^2+c r3*z(24)/wn)+rbc^2/jce*qr1*u*lr1*(kr1*d/wn^2+cr1*z(20)/wn)+rbc^2/jce*qr2*u*lr2*(kr2*e/wn ^2+cr2*z(22)/wn)+rbc^2/jce*qr3*u*lr3*(kr3*f/wn^2+cr3*z(24)/wn); z(22); 2*wc*(z(8)-z(12))*sin(140)/wn-wc^2*(z(5)-z(9))*sin(140)/wn^2+sin(140)*kr*z(9)/(mr*wn^2)-c os(140)*kr*z(11)/(mr*wn^2)+2*wc*(z(6)-z(10))*cos(140)/wn+wc^2*(z(11)-z(7))*cos(140)/wn^2 -kc*z(5)*sin(140)/(mce*wn^2)+kc*z(7)*cos(140)/(mce*wn^2)-cos(140)*cos(-20)/mce*(ks1*a/w n^2+cs1*z(14)/wn)-cos(140)*cos(100)/mce*(ks2*b/wn^2+cs2*z(16)/wn)-cos(140)*cos(220)/mce *(ks3*c/wn^2+cs3*z(18)/wn)+rbc^2/jce*(ks1*a/wn^2+cs1*z(14)/wn)+rbc^2/jce*(ks2*b/wn^2+cs 2*z(16)/wn)+rbc^2/jce*(ks3*c/wn^2+cs3*z(18)/wn)-(rbr^2/jr+rbc^2/jce)*(kr1*d/wn^2+cr1*z(20) /wn)-(rbr^2/jr+rbc^2/jce)*(kr2*e/wn^2+cr2*z(22)/wn)-(rbr^2/jr+rbc^2/jce)*(kr3*f/wn^2+cr3*z(2 4)/wn)-rb2^2/j2*(kr2*e/wn^2+cr2*z(22)/wn)-rb2^2/j2*(ks2*b/wn^2+cs2*z(16)/wn)-rbc*tc/(jce*b c*wn^2)+we^2*2.33*cos(we*t-0.85)/(bc*wn^2)+(sin(140)*cr*z(10)/(mr*wn)-cos(140)*cr*z(12)/ (mr*wn)-sin(140)*cc*z(6)/(mce*wn)+cos(140)*cc*z(8)/(mce*wn))*1-sin(140)*sin(-20)/mce*(ks 1*a/wn^2+cs1*z(14)/wn)-sin(140)*sin(100)/mce*(ks2*b/wn^2+cs2*z(16)/wn)-sin(140)*sin(220)/ mce*(ks3*c/wn^2+cs3*z(18)/wn)+(mce+mr)*sin(140)*sin(20)/(mce*mr)*(kr1*d/wn^2+cr1*z(20) /wn)+(mce+mr)*sin(140)*sin(140)/(mce*mr)*(kr2*e/wn^2+cr2*z(22)/wn)+(mce+mr)*sin(140)*s in(260)/(mce*mr)*(kr3*f/wn^2+cr3*z(24)/wn)+(mce+mr)*cos(140)*cos(20)/(mce*mr)*(kr1*d/w n^2+cr1*z(20)/wn)+(mce+mr)*cos(140)*cos(140)/(mce*mr)*(kr2*e/wn^2+cr2*z(22)/wn)+(mce+ mr)*cos(140)*cos(260)/(mce*mr)*(kr3*f/wn^2+cr3*z(24)/wn)+2*(mce+mr)/(mce*mr)*qr1*u*(k r1*d/wn^2+cr1*z(20)/wn)*sin(20)*cos(20)+2*(mce+mr)/(mce*mr)*qr2*u*(kr2*e/wn^2+cr2*z(2 2)/wn)*sin(140)*cos(140)+2*(mce+mr)/(mce*mr)*qr3*u*(kr3*f/wn^2+cr3*z(24)/wn)*sin(260)* cos(260)+rbr/jr*qr1*u*lr1*(kr1*d/wn^2+cr1*z(20)/wn)+rbr/jr*qr2*u*lr2*(kr2*e/wn^2+cr2*z(22) /wn)+rbr/jr*qr3*u*lr3*(kr3*f/wn^2+cr3*z(24)/wn)-rb2/j2*qs1*u*ls1*(ks1*a/wn^2+cs1*z(14)/w n)-rb2/j2*qs2*u*ls2*(ks2*b/wn^2+cs2*z(16)/wn)-rb2/j2*qs3*u*ls3*(ks3*c/wn^2+cs3*z(18)/wn) -rb2/j2*qr1*u*lr1*(kr1*d/wn^2+cr1*z(20)/wn)-rb2/j2*qr2*u*lr2*(kr2*e/wn^2+cr2*z(22)/wn)-rb 2/j2*qr3*u*lr3*(kr3*f/wn^2+cr3*z(24)/wn)+rbc^2/jce*qr1*u*lr1*(kr1*d/wn^2+cr1*z(20)/wn)+r bc^2/jce*qr2*u*lr2*(kr2*e/wn^2+cr2*z(22)/wn)+rbc^2/jce*qr3*u*lr3*(kr3*f/wn^2+cr3*z(24)/w 

n); z(24); 2*wc*(z(8)-z(12))*sin(260)/wn-wc^2*(z(5)-z(9))*sin(260)/wn^2+sin(260)*kr*z(9)/(mr*wn^2)-c os(260)*kr*z(11)/(mr*wn^2)+2*wc*(z(6)-z(10))*cos(260)/wn+wc^2*(z(11)-z(7))*cos(260)/wn^2 -kc*z(5)*sin(260)/(mce*wn^2)+kc*z(7)*cos(260)/(mce*wn^2)-cos(260)*cos(-20)/mce*(ks1*a/w n^2+cs1*z(14)/wn)-cos(260)*cos(100)/mce*(ks2*b/wn^2+cs2*z(16)/wn)-cos(260)*cos(220)/mce *(ks3*c/wn^2+cs3*z(18)/wn)+rbc^2/jce*(ks1*a/wn^2+cs1*z(14)/wn)+rbc^2/jce*(ks2*b/wn^2+cs 2*z(16)/wn)+rbc^2/jce*(ks3*c/wn^2+cs3*z(18)/wn)-(rbr^2/jr+rbc^2/jce)*(kr1*d/wn^2+cr1*z(20) /wn)-(rbr^2/jr+rbc^2/jce)*(kr2*e/wn^2+cr2*z(22)/wn)-(rbr^2/jr+rbc^2/jce)*(kr3*f/wn^2+cr3*z(2 4)/wn)-rb3^2/j3*(kr3*f/wn^2+cr3*z(24)/wn)-rb3^2/j3*(ks3*c/wn^2+cs3*z(18)/wn)-rbc*tc/(jce*b c*wn^2)+we^2*2.33*cos(we*t-0.85)/(bc*wn^2)+(sin(260)*cr*z(10)/(mr*wn)-cos(260)*cr*z(12)/ (mr*wn)-sin(260)*cc*z(6)/(mce*wn)+cos(260)*cc*z(8)/(mce*wn))*1-sin(260)*sin(-20)/mce*(ks 1*a/wn^2+cs1*z(14)/wn)-sin(260)*sin(100)/mce*(ks2*b/wn^2+cs2*z(16)/wn)-sin(260)*sin(220)/ mce*(ks3*c/wn^2+cs3*z(18)/wn)+(mce+mr)*sin(260)*sin(20)/(mce*mr)*(kr1*d/wn^2+cr1*z(20) /wn)+(mce+mr)*sin(260)*sin(140)/(mce*mr)*(kr2*e/wn^2+cr2*z(22)/wn)+(mce+mr)*sin(260)*s in(260)/(mce*mr)*(kr3*f/wn^2+cr3*z(24)/wn)+(mce+mr)*cos(260)*cos(20)/(mce*mr)*(kr1*d/w n^2+cr1*z(20)/wn)+(mce+mr)*cos(260)*cos(140)/(mce*mr)*(kr2*e/wn^2+cr2*z(22)/wn)+(mce+ mr)*cos(260)*cos(260)/(mce*mr)*(kr3*f/wn^2+cr3*z(24)/wn)+2*(mce+mr)/(mce*mr)*qr1*u*(k r1*d/wn^2+cr1*z(20)/wn)*sin(20)*cos(20)+2*(mce+mr)/(mce*mr)*qr2*u*(kr2*e/wn^2+cr2*z(2 2)/wn)*sin(140)*cos(140)+2*(mce+mr)/(mce*mr)*qr3*u*(kr3*f/wn^2+cr3*z(24)/wn)*sin(260)* cos(260)+rbr/jr*qr1*u*lr1*(kr1*d/wn^2+cr1*z(20)/wn)+rbr/jr*qr2*u*lr2*(kr2*e/wn^2+cr2*z(22) /wn)+rbr/jr*qr3*u*lr3*(kr3*f/wn^2+cr3*z(24)/wn)-rb3/j3*qs1*u*ls1*(ks1*a/wn^2+cs1*z(14)/w n)-rb3/j3*qs2*u*ls2*(ks2*b/wn^2+cs2*z(16)/wn)-rb3/j3*qs3*u*ls3*(ks3*c/wn^2+cs3*z(18)/wn) -rb3/j3*qr1*u*lr1*(kr1*d/wn^2+cr1*z(20)/wn)-rb3/j3*qr2*u*lr2*(kr2*e/wn^2+cr2*z(22)/wn)-rb 3/j3*qr3*u*lr3*(kr3*f/wn^2+cr3*z(24)/wn)+rbc^2/jce*qr1*u*lr1*(kr1*d/wn^2+cr1*z(20)/wn)+r bc^2/jce*qr2*u*lr2*(kr2*e/wn^2+cr2*z(22)/wn)+rbc^2/jce*qr3*u*lr3*(kr3*f/wn^2+cr3*z(24)/w n)];