syms x a
>> factor(x^4+x^2-2*a*x-a^2+1)
ans =
-(a+x+1+x^2)*(a+x-1-x^2)
MATLAB 编程解得：
整理后(x^2+x+a+1)(x^2-x-a+1)