curr_h = (D[ind] / 2) - (((1 + i * cmath.sqrt(3)) * cube_root(2 * cmath.sqrt(6) * cmath.sqrt(6 * V[ind] ** 2 - cmath.pi * D[ind] ** 3 * V[ind]) + cmath.pi * D[ind] ** 3 - 12 * V[ind])) / (4 * cube_root(cmath.pi))) - (((1 - i * cmath.sqrt(3)) * cube_root(cmath.pi) * D[ind] ** 2) / (4 * cube_root(2 * cmath.sqrt(6) * cmath.sqrt(6 * V[ind] ** 2 - cmath.pi * D[ind] ** 3 * V[ind]) + cmath.pi * D[ind] ** 3 - 12 * V[ind])))