This function approaches the Dirac delta-function (hence its name) in the limit of large multipoles: lim l-->+o o dl(z) = d(z). Indeed the source looks more and more like a point mass as we increase the multipolar order l.