scalar m=10[ kg ];
scalar g=9.81[ m/ s2 ];
scalar v0=3[ m/ s ];
scalar b=5[ N/ s ];
scalar T=2[ s ];
scalar mu=1 [ N s/ m ];
point point1=point( 0 [ m ], 0 [ m ], 0 [ m ], visible = hide: );
point point2=point( 0 [ m ], -0.5 [ m ], 0 [ m ], visible = hide: );
point point3=point( 0 [ m ], 0[ m ], 0 [ m ], visible = hide: );
point point4=point( 0 [ m ], -0.25 [ m ], 0 [ m ], visible = hide: );
solid solid3=cylinder( point2, point1, 0.1 [ m ], mass = m );
color color1=RGB( 0, 0, 0 );
body body1=body( color = color1, visible = show: );
body body1 < ( solid3, point1, point2, point4 );
color color2=index( 7 );
body body2=body( color = color2, visible = hide: );
set ground = body2;
body body2 < ( point3 );
condition condition1=transVelocity( body2, reverse( projectY ), body1, point2, v0 );
function function1(v[m/s])=-mu*#v;
function function2(t[ s ])=b*#t;
sensor v_theor=(v0+m*b/(mu*mu))*exp(-mu*T/m)+b*(T-m/mu)/mu;
sensor v=transVelocity( body2, reverse( projectY ), body1, point2 );
sensor v_delta=abs((v_theor-v)/v_theor);
gravity gravity1=parallel( reverse( projectY ), g = g );
force force1=force( body1, point4, reverse( projectY ), function1, list( v, v, v ), work = on:, visible = show: );
force force2=force( body1, point4, projectY, m*g, list( time ), work = on:, visible = hide: );
force force3=force( body1, point4, reverse( projectY ), function2, list( time ), visible = hide: );
event event1=reformsBySensor( list( stop(  ) ), time, T );

/\///////////////////////////////////////////////////////////////////////////////////
/\   Единицы измерения;
set units = SI;