scalar R=1 [ m ];
scalar h0=0.9 [ m ];
scalar alpha0=45 [ deg ];
scalar m=1[ kg ];
scalar v0=0.5[ m/ s ];
scalar h=0.6[ m ];
scalar g=9.81 [ m/ s2 ];
solid solid1=sphere( point( 0 [ m ], 0 [ m ], 0 [ m ] ), R );
point M=pointCyl( alpha0, sqrt(R*R-h0*h0), h0 );
solid solid2=sphere( M, 0.05 [ m ], mass = m );
body body1=body( color = RGB( 229, 229, 229 ) );
set ground = body1;
body body1 < ( solid1 );
body body2=body( color = RGB( 255, 153, 153 ) );
body body2 < ( M, solid2 );
joint joint1=contactPointSurface( M, solid1 );
gravity gravity1=parallel( reverse( projectZ ) );
vector vector1=rotate( rotate( rotate( vectorPP( point( 0 [ m ], 0 [ m ], 0 [ m ] ), M ), projectZ, -alpha0 ), projectY, 90 [ deg ] ),
	projectZ, alpha0 );
condition condition1=transVelocity( body1, vector1, body2, M, v0 );
sensor sensor1=bodyDisplacement( body1, point( 0 [ m ], 0 [ m ], 0 [ m ] ), projectZ, body2, M );
event event1=reformsBySensor( list( stop(  ) ), sensor1, h );
sensor teor=m*g/R*(3*h-2*h0-v0*v0/g);
sensor sensor2=jointForce( force:, body2, M, projectX, joint1 );
sensor sensor3=jointForce( force:, body2, M, projectY, joint1 );
sensor sensor4=jointForce( force:, body2, M, projectZ, joint1 );
sensor N=-sqrt(sensor2*sensor2+sensor3*sensor3+sensor4*sensor4);
sensor delta=abs(teor-N);

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