scalar m=1 [kg];
scalar F=20 [N];
scalar k=1[kg/s];
scalar mu=1[kg/s];
scalar v0=20[m/s];
scalar g=9.81[m/ s2];
function Q(v_x[m/s])=k*#v_x;
function R_x(v_x[m/s])=-mu*#v_x;
function R_y(v_y[m/s])=-mu*#v_y;
function Ft()=m*g;
function Fy_sum(v_y[m/s],v_x[m/s])=Q(#v_x)+R_y(#v_y)-Ft();
function Fx_sum(v_x[m/s])=F+R_x(#v_x);
color color1=RGB( 0, 0, 128 );
point point1=point( 0 [ m ], 0 [ m ], 0 [ m ], pointStyle = thickdot:, color = color1 );
color color2=RGB( 0, 0, 128 );
body bodyM=body( color = color2 );
point pointM=point1;
body bodyM < ( pointM );
MIP MIP1=massPoint( pointM, m );
body bodyM < ( MIP1 );
body body1=body( color = index( 1 ) );
set ground = body1;
body body1 < ( point1 );
sensor v_x=componentVelocity( pointM, projectX );
sensor S=-time*mu/m;
sensor v_y=componentVelocity( pointM, projectY );
sensor v=velocity( pointM );
sensor y_m=displacement( point1, pointM, projectY );
sensor y_theoretical=((m/mu)*(k*F/(mu*m)-g+(k*(F-mu*v0)/(mu*m))*exp(S))*time-(m/mu)*(m/mu)*(k*(2*F-mu*v0)/(mu*m)-g)*(1-exp(S)));
sensor y_delta=abs(y_m-y_theoretical);
condition cond1=transVelocity( body1, projectX, bodyM, pointM, v0 );
force forcex=force( bodyM, pointM, projectX, Fx_sum, list( v_x ) );
force forcey=force( bodyM, pointM, projectY, Fy_sum, list( v_y, v_x ) );

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