Calling Mathematica into Julia to symbolically solve a system of non-linear equations

Seems to me you are simply solving a system of linear equations here. You can use JuMP for that.