code for plot below

> hor:=seq(plot([[-1,j/10],[1, j/10]]),j=-10 ..10):
ver:=seq(plot([[j/10,-1],[ j/10,1]]),j=-10..10):
xaxis:=plot([[0,-1], [0,1]], thickness=2):
yaxis:=plot([[-1,0], [1,0]], thickness=2):
gp:=plots[display]([hor,ver,xaxis,yaxis]):
grph:=plot([x^3,x, x=-1..1],thickness=3 ):
plots[display]({ grph,gp},title=`y = f(x)`,titlefont=[TIMES,BOLD,16], axesfont=[TIMES,BOLD,12], xtickmarks=[-1,-.5,0,.5,1],ytickmarks=[-1,-.5,0,.5,1], scaling=constrained);