sum 1/(sqrt((z+i*d)^2+r^2)^3), i=0 to n