function [r]=iter(f,x0,N)
// 
//  Methode du point fixe selon la formule d appel
//  r=iter(f,x0,N)
//  
xp = x0;
r = [xp];
for i=1:1:N
  xp = f(xp);
  r = [r, xp];
end
endfunction