CoCoA » CoCoA-5

Sometimes LinearSimplify could produce a simpler result:

f := x^4+4*x^3+5*x^2+2*x+1;
LinearSimplify(f);
record[LinearChange := x, SimplePoly := x^4 +4*x^3 +5*x^2 +2*x +1]

subst(f,x,x-1);
x^4 -x^2 +1

Improve it (without making it measurably slower).