plot y=x^3+2x+1