function yh = odezva4(t,y0) global M2 B K2 x=y0(1:9); v=y0(10:18); yh = [v;-inv(M2)*B*v-inv(M2)*K2*x];