Updated version of the multipole calculation where the matrix is computed on the fly. Algorithm is now O(N).
