function y=testnewton(f,df,x,sol=0)
#testujemy met Newtona
#f - wsk do funkcji y=f(x)
#df wsk do pochodnej funkcji y=ff(x)
#x- przyblizenie startowe
#sol - rozwiazanie f(sol)=0

  tol=1e-8;
  y=f(x);
  e=x-sol;it=0;
 while(abs(y)>tol)
  dy=df(x);it++;
  x=x-y/dy;
  ep=e;e=x-sol;
  y=f(x);
  printf("[%d] e=%5.3g e/ep=%g e/(ep)^2=%g \n",it,e,e/ep,e/(ep*ep));
 endwhile
endfunction
