homotopy: use ((x-x0)**2 + (y-y0)**2 + (z-z0)**2 - r0**2)**2 for homotopy objective instead
homotopy uses collocation points instead of equal spacing