program prog17
      implicit none
      double precision Xa, Ya, Xb, Yb
      double precision vcos
      
      read *, Xa, Ya
      read *, Xb, Yb
      print *, 'Косинус угла между векторами:', vcos(Xa,Ya,Xb,Yb)
      end
      
      
      double precision function vcos(a1, a2, b1, b2)
      implicit none
      double precision a1, a2, b1, b2
           
      vcos = (a1*b1 + a2*b2)/sqrt(a1**2 + a2**2)/sqrt(b1**2 + b2**2)
      end