Question: Why does simplify produce wrong results for this polynomial?

When defining a polynomial as follows

p:=x->x^(r+1)-(r+1)*Sum((-1)^(r+k)*(r+k)!*(r+1)!*r!/((k!)^2*(r-k+1)!*(2*r+1)!)*x^k,k=0..r);

the result of

simplify((eval(p(1),r=1)));

is given as 1 whereas hand evaluation and also plotting the polynomial clearly shows that the result is 0. How can this behaviour be explained?

Please Wait...