plot[y=x^3+3x+1,x=-1,x=0,x=1,y=-3,y=5, {x,-1,1}]