For comparison with other CAS choose from: Axiom Derive Macsyma Maple Mathematica
decompose(x**6+9x**5+52x**4+177x**3+435x**2+630x+593);
3 2 2 {u + 25*u + 210*u + 593,u=x + 3*x}
decompose(x**4+2x**3*y + 3x**2*y**2 + 2x*y**3 + y**4 + 2x**2*y +2x*y**2 + 2y**3 + 5 x**2 + 5*x*y + 6*y**2 + 5y + 9);
2 2 2 {u + 5*u + 9,u=x + x*y + y + y}