Learn to use Symbolics.jl

General Usage
,
2

If you do the same thing as SymPy it seems to work fine:

using Symbolics, Groebner, Nemo, LinearAlgebra
@variables x[1:2] y[1:2] b[1:2]
eqs = collect([one.(x) x] * b .~ y)
symbolic_solve(eqs, collect(b))

1-element Vector{Any}:
 Dict{Num, Any}(b[2] => (y[1] - y[2]) / (x[1] - x[2]), b[1] => (x[1]*y[2] - x[2]*y[1]) / (x[1] - x[2]))
1 Like