point point1=point( 0 [ m ], 0 [ m ], 0 [ m ] );
scalar l=0.5[ m ];
scalar m=2[ kg ];
scalar omega0=1000[ deg/ s ];
point point2=point( 3*l, 0 [ m ], 0 [ m ] );
solid solid1=cylinder( point1, point2, 0.01 [ m ], mass = m );
scalar d=0*l;
point point3=point( d, 0 [ m ], 0 [ m ] );
point point4=point( d+l, 0 [ m ], 0 [ m ] );
point point5=point( 0 [ m ], 0 [ m ], 0 [ m ] );
body body1=body( color = RGB( 229, 229, 229 ) );
set ground = body1;
body body1 < ( point1 );
body body2=body( color = RGB( 255, 153, 153 ) );
body body2 < ( point2, solid1, point5 );
solid solid2=sphere( point3, 0.03 [ m ], mass = m );
solid solid3=sphere( point4, 0.03 [ m ], mass = m );
joint joint1=rotational( body1, body2, point1, projectZ );
condition condition1=rotVelocity( body1, projectZ, body2, omega0 );
body body3=body( color = RGB( 229, 229, 229 ) );
body body3 < ( solid2, point3 );
body body4=body( color = RGB( 255, 153, 153 ) );
body body4 < ( solid3, point4 );
sensor T_teor=m*l*omega0*omega0/32;
sensor v_teor=sqrt(3/2)*omega0*l/1[ rad ];
vector vector1=vectorPP( point1, point2 );
joint joint2=translational( body2, body3, point3, vector1 );
joint joint3=translational( body2, body4, point4, vector1 );
sensor sensor1=displacement( point4, point2, vector1 );
joint joint4=sphericalSpherical( body3, point3, body4, point4 );
event event1=reformsBySensor( list( stop(  ) ), sensor1, 0 [ m ] );
sensor sensorT=jointForce( force:, body4, point4, vector1, joint4, fixing = unlock: );
sensor T=abs(sensorT)*1[ rad2 ];
sensor delta_T=abs((T_teor-T)/T_teor);
sensor v=transVelocity( body2, vector1, body4, point4 );
sensor delta_v=abs((v_teor-v)/v_teor);

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