% Sitzung S. 210 unten

d = 1+ceil(2*sum(log10(diag(A))))
alpha = 2^30; K = ceil((d+7)/log10(alpha));
tic, [X1,fail] = wan(A,b,M,alpha,K); toc

% Store the list of xi values to the file 'X1.txt'

file= fopen('X1.txt','w');
fprintf(file,'%20.16g\n',n);
for k=1:length(X1);
    fprintf(file,'%20.16g\n',X1(k));
end
fclose(file);