Maple Code

> fixedpoint := proc(x0,max)
  local k,p0,p1;
  k := 0;
  p0 := evalf(x0);
  printf(" P%g  = %g \n",k,p0);
  p1 := p0;
  for k from 1 to max do
    p0 := p1;
    p1 := g(p0);
    printf(" P%g  = %g \n",k,p1);  
  end do;
  print(`g(x) = `,g(x));
  printf("  P  = %g \n",p1);
  printf("g(P) = %g \n",g(p1));
  RETURN(p1);
  end;

> g := proc(x)
  1+x-x^2/3
  end proc;

(c) John H. Mathews 2004