PROGRAM c168c

IMPLICIT NONE

REAL*16 :: xa, xb, ya, yb, sa, sb, z
INTEGER :: A, B, H, M
REAL*16, PARAMETER :: p = 3.1415926535897932384626433832795

READ *, A, B, H, M

sa = 2.0q0 * p / 7.20q2 * (H * 6.0q1 + M * 1.0q0 -1.5027q-5) 
sb = 2.0q0 * p / 6.0q1 * M * 1.0q0
xa = A * 1.0q0 * COS(-sa)
xb = B * 1.0q0 * COS(-sb)
ya = A * 1.0q0 * SIN(-sa)
yb = B * 1.0q0 * SIN(-sb)
z = ((xa - xb) ** 2.0q0 + (ya - yb) ** 2.0q0) ** 5.0q-1

PRINT *, z

END PROGRAM c168c