function [y,flag]=newton2(f,fprime,y0,tol)
%
%  function [y,flag]=newton2(f,fprime,y0,tol)
%
count=0;
quit = 0;
flag=0;
y=y0
x=zeros(1,20)
while quit ==0 
  oldy=y;
  fy=eval([f,'(y)'])
  fprimey=eval([fprime,'(y)'])
  y=y-fy/fprimey
  m=tol;
  if abs(y)>m
    m=abs(y) 
   end;
  if abs(y-oldy) < m*tol
    quit=1
  end;
  count=count+1
  x(count)=y;
  if count >19 
    quit=1;
    flag=1;
  end;
end;
x(1:count)