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