function y = phi(x)
  y = x + 1 / (x + 1);
endfunction


getf('mpf.sci');

[raiz, iter] = mpf(phi, 1)
