fixed make multipole method, centre was being incorrectly calculated
fixed multipole calculation and tested for rotational invariance