50 setxunit
7.5 setborder
90 setangle_repere
0 4 setxrange
0 4 setyrange
1 1 settkstep
.5 .5 setsubtkstep
1 1 setmkstep
ticks
subticks
ymarks
traceaxes
axesarrow
unites
/f {setxvar
#rpn#   x*(4-x)
} def
gsave
   bleu
   .75 setlinewidth
   pointilles
   {f} courbe
grestore
/u0 .5 def           %% 1er terme
/n 500 def             %% nb de termes
/l1 [u0 n {dup f} repeat] def  %% l1 = le tableau [u0 u1 ... u(n-1)]
/l2 l1 {f} apply def           %% l2 = le tableau [u1 u2 ... un]
/L1 l1 dup fuz def             %% L1 = (l1, l1) soit [u0 u0 u1 u1 ... un un]
/L2 l1 l2 fuz def              %% L2 = (l1, l2) soit [u0 u1 u1 u2 ... u(n-1) un]
L2 {point} plot                %% on represente L2
{} courbe                      %% trace de la ligne de rappel
gsave
   .25 setlinewidth
   pointilles 
   gris 
   [                                      
   u0 0                           %% le point (u0 0)        
   /i 0 def                               
   L2 length 2 idiv {            
      L1 i getp                   %% le point ( u(i-1) , ui)
      L2 i getp
      /i i 1 add store                    
   } repeat                               
   ] ligne                        %% on relie par des pointilles   
grestore

u0 xtick
setTimesItalic 
   (u)   u0 0 dtext
setTimesRoman 
   (0)   indice

0 1 10 1 sub {
   /i exch def
   setTimesItalic 
   (u)
   l2 i get 
   0 dctext
   setTimesRoman 
   i 1 add (   ) cvs
   indice
} for