# interleaf following page 134 of Chapter 1
# see # lines in il1

with(plots):

f := 16*x^3-7*x^2*y+55*x^2-104*x*y^2+45*x*y+59*x-65*y^3+130*y^2+65*y-130:
l := [[-5, -1], [-3, 1], [-13/16, 51/16], [-60/31, 26/31], [0, 2], 
[0, 1], [0, -1]]:
a := [[-5, -1.2, `P`], [-3.17, 1.2, `O'`], [-.3, .8, O], [-.7, -1.1, `P+Q`], 
[-.4, 2, R], [-1.25, 3.2, `-P`], [-2, .6, Q], [-4, -2, `Plot by MapleV3`]]:

p1:=implicitplot({f},x=-6..4,y=-2..4,grid=[120,105]\
         ,scaling=CONSTRAINED):
p2:=implicitplot({y=1},x=-6..4,y=-2..4,grid=[120,105]\
         ,scaling=CONSTRAINED):
p3:=implicitplot({y=x+4},x=-6..4,y=-2..4,grid=[120,105]\
         ,scaling=CONSTRAINED):
p4:=implicitplot({y=(3/5)*x+2},x=-6..4,y=-2..4,grid=[120,105]\
         ,scaling=CONSTRAINED):
p5:=plot(l,style=point,symbol=CIRCLE,scaling=CONSTRAINED):
p6:=textplot(a,align = RIGHT, font = [TIMES, ITALIC, 12]):

display({p.(1..6)});