    function xdot = vozik(t,x)
    global u m kt kp
    xdot = [x(2); (u/m)-kt*(x(1)/m)-kp*(x(2)/m)];
