Oh yeah, and it can tell the PMI VECTOR3 from any other VECTOR3. Compilers are so smart nowadays :P
It will look like this (assuming the axes are similarly aligned for both vessels):
#define sqr(x) (x*x)
double m1, r1, m2, r2; // Mass and coordinates of the vessels
VECTOR3 PMI1, PMI2...