scalar m=2 [ kg ];
scalar c=100 [ N/ m ];
scalar lambda=0.2 [ m ];
scalar l=2[ m ];
scalar R=0.1[ m ];
point C=point( 0 [ m ], 0 [ m ], 0 [ m ] );
point B=point( l+lambda, 0 [ m ], 0 [ m ] );
MIP MIP1=massPoint( B, m );
line line1=polyLine( list( move( C, reverse( projectY ), R ), move( B, reverse( projectY ), R ) ) );
color color1=index( 23 );
line line2=circle( C, projectZ, R, color = color1 );
solid solid1=extrude( line2, move( C, projectZ, 0.05 [ m ] ), move( C, reverse( projectZ ), 0.05 [ m ] ), mass = 2*m );
point A=point( -3.1000e-001 [ m ], 2.4000e-001 [ m ], 0.0000e+000 [ m ] );
body body1=body( color = RGB( 229, 229, 229 ) );
body body1 < ( line2, solid1, A );
body body2=body( color = RGB( 255, 153, 153 ) );
body body2 < ( C );
body body3=body( color = RGB( 255, 153, 255 ) );
body body3 < ( B, MIP1 );
body body4=body( color = RGB( 255, 255, 153 ) );
set ground = body4;
body body4 < ( line1 );
joint joint1=contactLineLine( line2, line1, lines_slip = linesNoSlip: );
force force1=spring( body1, C, body3, B, c, s0 = l );
joint joint2=rotational( body2, body1, C, projectZ );
sensor sensor1=bodyDisplacement( body2, C, projectX, body3, B, fixingVector = unlockVector: );
sensor sensor2=abs(sensor1);
event event1=reformsBySensor( list( stop(  ) ), sensor2, l );
sensor sensor3=transVelocity( body4, projectX, body2, C, fixingVector = unlockVector: );
sensor sensor4=(lambda/2)*sqrt(c/(3*m));
sensor pogr_otn=abs(sensor3-sensor4)/sensor4;

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