Implements multipoles. Multipole integrals are correct at least up to quadrupole

test update included.