假设有耦合PDE的一个系统,诸如在Python中使用sympy可以解决一个PDE系统吗?
第一PDE在F(A,B)
第二PDE在F(A,B)
以下代码能够分别解决每个PDE:
import numpy as np
import sympy as sp
# definition of variables
a, b = sp.symbols('a b')
f = sp.Function('f')
F = f(a, b)
Fda = F.diff(a)
Fdb = F.diff(b)
# definition of PDEs
eq1 = Fda - 2
eq2 = Fda + Fdb + 2
# solution of separated PDEs
sp.pprint(sp.pdsolve(eq1))
sp.pprint(sp.pdsolve(eq2))
是否有可能解决PDEs系统?语法可能类似于sp.pprint(sp.pdsolve([eq1, eq2]))
。我试过[eq1, eq2]
,{eq1, eq2}
,np.array([eq1, eq2])
等我看看help(sp.pdsolve)
和help(sp.pde)
,但还没有找到解决办法。