elif (cx - u) * (u - ulim) > 0.0:
            fu = fct(u)
            if fu < fc:
                bx, cx = (cx, u)
                u = cx + GOLD * (cx - bx)