function yh = odezva3(t,y0) global M1 B K1 x=y0(1:9); v=y0(10:18); yh = [v;-inv(M1)*B*v-inv(M1)*K1*x];