scalar l=0.3[m];
scalar r=0.1[m];
scalar m=5[kg];
scalar alfa=45[deg];
scalar w=17.453292[rad/s];
scalar h_rel=0.01;
scalar h=1[ m ];
point point1=point( 0[m], -1 [ m ], 0[m] );
point point2=point( 0[m], 1 [ m ], 0[m] );
solid solid1=cylinder( point1, point2, 0.01[m] );
point point3=point( 0[m], 0[m], 0[m] );
vector vector1=rotate( projectY, projectZ, alfa );
point point4=move( point3, vector1, l );
point point5=move( point3, vector1, -l );
solid solid2=cylinder( point4, point5, 0.1[ m ], mass = m );
color color1=index( 54 );
body body1=body( color = color1 );
set ground = body1;
color color2=index( 64 );
body body2=body( color = color2 );
body body2 < ( solid1, solid2 );
joint joint1=spherical( body1, body2, point1 );
node node1=nodePoint( point2 );
joint joint2=user( body1, body2, node1, noX:, yesY:, noZ:, yesFiX:, yesFiY:, yesFiZ: );
condition condition1=rotVelocity( body1, projectY, body2, w );
sensor F1x=jointForce( force:, body1, point1, projectX, joint1 );
sensor F1z=jointForce( force:, body1, point1, projectZ, joint1 );
sensor F1=sqrt(F1x*F1x+F1z*F1z);
sensor F2x=jointForce( force:, body1, point2, projectX, joint2 );
sensor F2z=jointForce( force:, body1, point2, projectZ, joint2 );
sensor F2=sqrt(F2x*F2x+F2z*F2z);
sensor Q_theoretical=m*w*w*sin(2*alfa)/(4*h)*(l*l/3-r*r/4)/1[rad2];
sensor Q_delta_rel=abs((F1-Q_theoretical)/Q_theoretical);

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