with(plots):
g:=x->x/2+2/x^2:
q[1]:=plot({x,g(x)},x=1.45..1.7,y=1.45..1.7,color=blue):
x[1]:=1.5:
x[2]:=g(x[1]):
p[2]:=[[x[1],0],[x[1],x[2]]]:
for i from 2 to 10 do
  p[2*i-1]:=[[x[i-1],x[i]],[x[i],x[i]]]:
  x[i+1]:=g(x[i]):
  p[2*i]:=[[x[i],x[i]],[x[i],x[i+1]]]
od:
for i from 2 to 10 do
  q[i]:=display({q[1],plot({seq(p[j],j=2..2*(i-1))},color=red)})
od:
display([seq(q[i],i=1..10)],insequence=true,view=[1.45..1.7,1.45..1.7]);
p:='p': q:='q': i:='i':