I created a benchmark based on this problem in https://github.com/mattpap/sympy/commit/ac818efbbe69100e31691daf21de93ec7fb7df40, of course using sparse polynomials. It takes 30 seconds on my machine to solve this system, but `prun` shows that it can be improved.