Pro polar, coords1, coords2
 
magnitude = sqrt((coords_begin[0]-coords_end[0])^2+(coords_begin[1]-coords_end[1])^2)
slope = (coords_end[1]-coords_begin[1])/(coords_end[0]-coords_begin[0])
angle = 180/!DPI*atan(slope)
 
print, 'magnitude of vector is: ', magnitude
print, 'angle of vector is: ', angle
end