Ask Your Question
1

Solving Lagrangians in Sage

asked 2020-05-03 01:55:31 -0500

EconJohn gravatar image

updated 2020-07-16 11:06:59 -0500

I have been running the following Sage code on cocalc.com 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

Comments

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 13:42:27 -0500 )edit

1 answer

Sort by » oldest newest most voted
3

answered 2020-05-03 08:56:21 -0500

tmonteil gravatar image

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

U=x^0.7*y^0.3;U

with

U=x^(7/10)*y^(3/10);U

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

Stats

Asked: 2020-05-03 01:55:31 -0500

Seen: 181 times

Last updated: Jul 16