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 [email protected].
To post to this group, send email to [email protected].
Visit this group at http://groups.google.com/group/sage-support.
For more options, visit https://groups.google.com/groups/opt_out.