Note please that this is really not a Julia question per se, since all SymPy.jl does is call out to Python.

I let it run in the background and forgot all about it, but the solve did finish. **Edit:** But it took twenty minutes.

```
julia> @time solve([eq1,eq2,eq3],[x,y,z], domain=SymPy.sympy.S.Reals)
1193.533067 seconds (29.31 M allocations: 333.714 GiB, 0.06% gc time)
10-element Vector{Tuple{Sym, Sym, Sym}}:
(0, 0, 0)
(0, 0, 1)
(0, 1, 0)
(0, 1, 1)
(1, 0, 0)
(1, 0, 1)
(1, 1, 0)
(1, 1, 1)
(-(c₃ - g₁)/(g₁ + s₁), (c₁ - e₁)/(c₂ - e₂ + g₁ + s₁), 0)
(-(c₃ - g₁)/(g₁ + s₁), (c₁ - c₂ - e₁ + e₂ - g₂)/(g₁ + s₁), 1)
```

Your `fx,...`

are called `eq1,...`

here, and I asked for real solutions. Although that probably won’t make a difference, because no assumptions about the other variables are made.