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