program colon

 implicit none
 real::a,b,h,m
 real(8),parameter::pi=3.1415926535897932384626433D0
 real(8)::rad_h1,rad_h2,rad_m,cos_h,cos_m,sita1,sita2,sita,x
 read*,a,b,h,m

 rad_h1=30.0D0*pi/180.0D0
 rad_h2=0.50D0*pi/180.0D0
 rad_m=6.0D0*pi/180.0D0

 cos_h=rad_h1*h+rad_h2*m
 cos_m=rad_m*m

 sita1=abs(cos_h-cos_m)
 sita2=abs(360.0D0-abs(cos_h-cos_m))
 sita=min(sita1,sita2)
 
 x=sqrt(a**2+b**2-2*a*b*cos(sita))
 print*,x

end program colon
 
