>
sum('k^4','k=1..n');factor(%);
>
simplify(%-(1/30)*(n*(n+1)*(2*n+1)*(3*n^2+3*n-1)));
>