Julia 如何从 SymPy 求解集中获取系数

Question

我在 Julia 中使用符号数学。当我进行微分时，结果非常好，但我无法计算出系数

using SymPy
@vars x y
z = x*(10 + 5*x + 4*y)
dz = diff(z,x)
x_s = solveset(dz,x)

如何从 x_s 中获取系数？

Answer 1

您可以使用elements 将有限集的元素作为数组获取:

julia> elements(s)
1-element Array{Sym,1}:
 -2*y/5 - 1

可以通过不同的方式获取系数，但这里我们将值转换为 Polynomial 类型，然后使用其 coeffs 方法:

julia> p = sympy.Poly(elements(s)[1], y)
Poly(-2/5*y - 1, y, domain='QQ')

julia> p.coeffs()
2-element Array{Sym,1}:
 -2/5
   -1