n=200;
A=1;
y=10;
for i=1:n
   A = log(y) - log(A);
   nothing = sprintf('%d  %e  %e', i, A, A - Ans);
   disp(nothing)
end
Ans = A;
