tic

itmax=100;

k=0
if x ~= 0
  xold=0;
else 
  xold=1;
end
[x abs(x-xold)]

S=size(A);
n=S(2)-1;

while abs(x-xold)>eps*max(abs(x),1.0)
     if k+1>itmax
       break
     end 
     xold=x;
     k=k+1
     b=A(1);
     c=A(1);
     for i=1:n-1
	b=A(i+1)+x*b;
        c=b+x*c;
     end
     b=A(n+1)+x*b;
     x=x-b/c;
     [x abs(x-xold)]
end  

toc