added another way of computation for thick multipoles