plot f(x)=(-x^3+x^2+3x+1)/(x+1)