sympy.solveについて

Question

sympy.solveにて、解けない場合は？

現在sumpy.solveにて方程式を解こうとしているのですが、永遠に解が出てこず(エラーも無)、
解がないorコードミスor他、なぜ解けないのかが分からない状態です。どのようにデバッグすればいいのかアドバイスいただきたいです。

簡易的にコードを示します

def fx(x,y)
return sy.log(x)-(y**2)/2

def gx(x,y):
return x-w*y/(1+0.01)

p1 = 0.00000001
p2 = 0.000001
p3 = 0.1
def Ux(x,y):
p = p1x+p2y+p3
return fx(x,y)+β((1-p)*fx(x,y))

f=fx(x,y)
g=gx(x,y)
L = f-g*l

diffx=diff(L,x)
diffy=diff(L,y)
diffl=diff(L,l)
A = solve([diffx,diffy,diffl)

---

要するに、関数Uxを制約gxのもと最大になるようなx,yの組み合わせを、ラグランジュの未定乗数法により解こうとしています。
処理的に大きな計算量ではないので解けるとしたら数秒で解は出ると思うのですが...

Answer 1

質問欄にあるのは実際に動かしたコードとは別ですか？　そこら中にある文法ミス？コピペミス？を修正したところ、普通に動きました。

import sympy as sy
(a,b,c,d,e,f,g,h,i,j,k,l,m,n,o,p,q,r,s,t,u,v,w,x,y,z) = sy.symbols('a b c d e f g h i j k l m n o p q r s t u v w x y z')

def fx(x,y):
    return sy.log(x)-(y**2)/2

def gx(x,y):
    return c-w*x*y/(1+0.01)

p1 = 0.00000001
p2 = 0.000001
p3 = 0.1
def Ux(x,y):
    p = p1*x+p2*y+p3
    return fx(x,y)+β((1-p)*fx(x,y))

f=fx(x,y)
g=gx(x,y)
L = f-g*l

diffx=sy.diff(L,x)
diffy=sy.diff(L,y)
diffl=sy.diff(L,l)
A = sy.solve([diffx,diffy,diffl])
print(A)
# [{c: -0.99009900990099*I*w*x, l: -1.01*I/(w*x), y: -1.0*I},
#  {c: 0.99009900990099*I*w*x, l: 1.01*I/(w*x), y: 1.0*I}]

Comments

[junjunn]

ご回答ありがとうございます。実際のコードより雑に書いてしまっています、申し訳ありません。
また、gxに少し誤りがありましたので修正しました。解中のIというのはどこから出てきたものか分かりますでしょうか。

[jeanbiego]

Iは虚数単位と思います