The convergence behaviour of the algorithm when fitting multipole based electrostatic models is also analysed against the seven test problems of this study. The steps followed by the algorithm are essentially identical to that followed for the point-charge model except that lattice minimisations and second derivative calculations are performed with the inclusion of all multipole interactions up to l 1 + l 2 = 4.