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]; 