sympyで関数の二回微分と代入ができない。

前提・実現したいこと

sympyを使って、関数の二回微分の行列を作りたい
変数がa1,a2,s1,s1,s2,s3,s4,s5
ベクトルｘの要素にそれぞれ上の変数の値を入れている。

発生している問題・エラーメッセージ

sequence[i] = None if _aresame(*s) else tuple(s)

TypeError: _aresame() takes 2 positional arguments but 5 were given

該当のソースコード

x = [3000, 3000, 500, 500, 500, 500, 500]
def d2l(x):
a1, a2, s1, s2, s3, s4, s5 = sym.symbols("a1, a2, s1, s2, s3, s4, s5")
L = a1 * length[0] + a2 * length[1] - lamda * ((-(P1 + P2) / a1 + sigma - s1) + (-(P2) / a2 + sigma - s2) + (-(P1 + P2)length[0] / a1E - P2 * length[1] / a2 * E + U - s3) + (a1 - s4) + (a2 - s5))
dLa1a1 = sym.diff(L, a1, 2).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLa1a2 = sym.diff(L, a1, a2).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLa1s1 = sym.diff(L, a1, s1).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLa1s2 = sym.diff(L, a1, s2).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLa1s3 = sym.diff(L, a1, s3).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLa1s4 = sym.diff(L, a1, s4).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLa1s5 = sym.diff(L, a1, s5).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLa2a1 = sym.diff(L, a2, a1).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLa2a2 = sym.diff(L, a2, 2).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLa2s1 = sym.diff(L, a2, s1).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLa2s2 = sym.diff(L, a2, s2).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLa2s3 = sym.diff(L, a2, s3).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLa2s4 = sym.diff(L, a2, s4).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLa2s5 = sym.diff(L, a2, s5).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs1a1 = sym.diff(L, s1, a1).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs1a2 = sym.diff(L, s1, a2).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs1s1 = sym.diff(L, s1, 2).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs1s2 = sym.diff(L, s1, s2).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs1s3 = sym.diff(L, s1, s3).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs1s4 = sym.diff(L, s1, s4).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs1s5 = sym.diff(L, s1, s5).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs2a1 = sym.diff(L, s2, a1).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs2a2 = sym.diff(L, s2, a2).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs2s1 = sym.diff(L, s2, s1).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs2s2 = sym.diff(L, s2, 2).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs2s3 = sym.diff(L, s2, s3).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs2s4 = sym.diff(L, s2, s4).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs2s5 = sym.diff(L, s2, s5).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs3a1 = sym.diff(L, s3, a1).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs3a2 = sym.diff(L, s3, a2).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs3s1 = sym.diff(L, s3, s1).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs3s2 = sym.diff(L, s3, s2).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs3s3 = sym.diff(L, s3, 2).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs3s4 = sym.diff(L, s3, s4).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs3s5 = sym.diff(L, s3, s5).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs4a1 = sym.diff(L, s4, a1).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs4a2 = sym.diff(L, s4, a2).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs4s1 = sym.diff(L, s4, s1).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs4s2 = sym.diff(L, s4, s2).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs4s3 = sym.diff(L, s4, s3).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs4s4 = sym.diff(L, s4, 2).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs4s5 = sym.diff(L, s4, s5).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs5a1 = sym.diff(L, s5, a1).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs5a2 = sym.diff(L, s5, a2).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs5s1 = sym.diff(L, s5, s1).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs5s2 = sym.diff(L, s5, s2).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs5s3 = sym.diff(L, s5, s3).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs5s4 = sym.diff(L, s5, s4).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
dLs5s5 = sym.diff(L, s5, 2).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
d2L = np.array([[dLa1a1, dLa1a2, dLa1s1, dLa1s2, dLa1s3, dLa1s4, dLa1s5], [dLa2a1, dLa2a2, dLa2s1, dLa2s2, dLa2s3, dLa2s4, dLa2s5],
[dLs1a1, dLs1a2, dLs1s1, dLs1s2, dLs1s3, dLs1s4, dLs1s5], [dLs2a1, dLs2a2, dLs2s1, dLs2s2, dLs2s3, dLs2s4, dLs2s5],
[dLs3a1, dLs3a2, dLs3s1, dLs3s2, dLs3s3, dLs3s4, dLs3s5], [dLs4a1, dLs4a2, dLs4s1, dLs4s2, dLs4s3, dLs4s4, dLs4s5],
[dLs5a1, dLs5a2, dLs5s1, dLs5s2, dLs5s3, dLs5s4, dLs5s5]])
d2L = d2L.astype(np.float64)
return d2L

print(d2l(x))

試したこと

補足情報（FW/ツールのバージョンなど）

ここにより詳細な情報を記載してください。

toast-uz

2020/10/06 11:57 編集

コードは```に挟まれた中に記載ください。インデントが消えて、分からなくなってしまっています。 length、sigma、E、Uといった変数が定義されているようですが、示されたコードに定義がありません。必ず、貼り付けたコード単体で動作するようにしてください。またimport文も省略しないでください。エラーが発生した場所を明確にしてください。

行動規範の内容に同意します

回答2件

自己解決

コードミスです。すいません。

投稿2020/10/06 15:37

gingingin

総合スコア2

今出ているエラーは、カッコの位置が間違っているために発生しています。

誤: dLa1a1 = sym.diff(L, a1, 2).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])

正: dLa1a1 = sym.diff(L, a1, 2).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3]), (s3, x[4]), (s4, x[5]), (s5, x[6])])

x[3]の後ろの閉じカッコが無く、その分最後のx[6]の閉じカッコが多いです。たくさんの行で同様の間違いがあります。

L=の行は定義が不足して再現不可能ですので、L=1として、再現を試みたところ、以下のエラーが出ました。

  File "main.py", line 8, in d2l
    dLa1a2 = sym.diff(L, a1, a2).subs([(a1, x[0]), (a2, x[1]), (s1, x[2]), (s2, x[3], (s3, x[4]), (s4, x[5]), (s5, x[6]))])
  File "/Users/（ユーザ名）/python/.teratail/lib/python3.8/site-packages/sympy/core/basic.py", line 939, in subs
    sequence[i] = None if _aresame(*s) else tuple(s)
TypeError: _aresame() takes 2 positional arguments but 5 were given

質問ではエラーの下の2行だけを投稿されていますが、上の行があれば、dLa1a1 = の行のsubsというメソッドでエラーが発生していることが明確になり、TypeError: _aresame() takes 2 positional arguments but 5 were givenから、エラーの内容はsubsの引数の数が間違っていることが、すぐにわかります。エラーの内容を削って投稿されると、切り分けが困難になりますので、ご注意ください。

投稿2020/10/06 12:08

toast-uz

総合スコア3266