code for diagram above

> f:=x->x^5-5*x^4+5*x^3:
p1:=plot(f(x),x= -1..3.7,thickness=1):

> p2:=plot({[.63,.57],[1,1],[2.37,-16.3],[0,0],[.63,.57]},
style=POINT, symbol=circle,thickness=3,color=blue):

> plots[display]({p1,p2},ytickmarks=[5,0,-10,-20,-30],
axesfont=[TIMES,BOLD,14]);