code for diagram below

> plot(x^3+x,x=-1.5..1.5,xtickmarks=[0],ytickmarks=[0],labelfont=[TIMES,BOLD,16],thickness=3,axesfont=[TIMES,BOLD,12]);