Finding Coefficients of Cross Terms

[Umut-Can-Physics]

When I try to run the following code, the print function always returns the integer 0 (zero). However, you can easily see that the coefficient of the z[1] * z[2] term equals -2. The Symbolics.jl module is robust for just single variables, such as coefficients of z[1]^2 and z[2]^2.

using Symbolics

@variables z[1:2]

Ψ = z[1]^2 - 2 * z[1] * z[2] + z[2]^2

# Extract coefficient
coeff_value = Symbolics.coeff(Ψ, z[1]*z[2])

println(coeff_value)

But interestingly Symbolics.coeff(Ψ,z[2]^0*z[1]^2) coefficient is true. The module doesn’t work for Symbolics.coeff(Ψ, z[1]*z[2]) Why?

[Umut-Can-Physics]

Anyone who solve this problem?

[zengmao]

Duplicate of #1041

[ChrisRackauckas]

Fixed in #1240