function [x,e,it]=odwracanie(a,x0=.75,tol=1e-14)
#odwracanie bez dzielenia
#m Newtona dla 1/x-a=0
#dla a\in[1,2) x0=.75 ok

 it=0;x=x0;
 e=1-x*a;
 while ((abs(e)>tol)&&(it<7))
  x=x*(1+e);
  it++;
  e=1-x*a;
 endwhile
endfunction
