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