Ask Your Question

Solving Lagrangians in Sage

asked 2020-05-03 08:55:31 +0200

EconJohn gravatar image

updated 2020-07-16 18:06:59 +0200

I have been running the following Sage code on to solve a standard utility maximization problem that is in an intermediate microeconomics course. However I'm having difficulty solving these problems in the context of Sage.

Below is the code that I have run to solve the problem and so far I believe it's correct.

x, y, l = var('x, y, l')
U = x^0.7 * y^0.3; U
m = 2*x+2*y; m
solve(m == 4000, y)
L = U - l * (m - 4000); L
dLdx = L.diff(x); dLdx
dLdy = L.diff(y); dLdy
dLdl = L.diff(l); dLdl
solve([dLdx == 0, dLdy == 0, dLdl == 0], x, y, l)

All my code seems to be running fine, however the last line of:

solve([dLdx == 0, dLdy == 0, dLdl == 0], x, y, l)

is not outputting the final demand equations. Is there any reason why this is happening?

edit retag flag offensive close merge delete


Regarding your note for mods: visit your profile page and find all your questions and answers.

To visit your profile page, click on your name near this question. or at the very top left from any page.

slelievre gravatar imageslelievre ( 2020-05-03 20:42:27 +0200 )edit

1 Answer

Sort by » oldest newest most voted

answered 2020-05-03 15:56:21 +0200

tmonteil gravatar image

Sage symbolics handles fractional powers better than numerical ones, so if you replace




it seems to work much better.

edit flag offensive delete link more

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

1 follower


Asked: 2020-05-03 08:55:31 +0200

Seen: 208 times

Last updated: Jul 16 '20