I'm trying to use Sage for optimization using the Lagrange multipliers method, but it does not solve this apparently simple problem:
var('x y l') f = x + 2*y g = 2*x^(1/4)*y^(1/2) L = f - l * (16 - g) dfdx = L.diff(x) dfdy = L.diff(y) dfdl = L.diff(l) solve([dfdx == 0, dfdy == 0, dfdl == 0], x,y,l) First Sage takes a very long time and then it just returns the equations, not the solutions. Is this a bug or am I doing something wrong? -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-support+unsubscr...@googlegroups.com. To post to this group, send email to sage-support@googlegroups.com. Visit this group at http://groups.google.com/group/sage-support. For more options, visit https://groups.google.com/groups/opt_out.