Z pomocą idzie sagemath
http://doc.sagemath.org/html/en/reference/calculus/sage/symbolic/expression.html

Są również inne programy np. Derive, dostępne tutaj [1] [2]
Wystarczy zdefiniować poszczególne równania dla x i y
Zaznaczyć #1 #2
Wykonać SOLVE dla równań #1 i #2 czyli Solve->System (Ctrl+Shift+Y), wpisując w oknie ‘Solve 2 equation(s)’ w polu 1 – #1, w polu 2 – #2
Na równaniu #3 wykonać Simplify (Ctrl+B)
Po otrzymaniu wyniku na #4 wykonać Approximate (Ctrl+G)