I've tried all sorts of initial conditions but this just won't converge. It is an easy Lagrange multipliers problem.
var("y z L")
F(x, y, z, L) = (x-3)^2 (y-1)^2 (z 1)^2 L*(x^2 y^2 z^2 - 4)
minimize(F, [4, 3, 2, 1], algorithm='ncg' )
Any thoughts? Are there ways of forcing it to use more iterations, or something??