y=x(x^2+1)