%CONDENUFIN
%
function [cond,lcond]=condEnufin(N,nu,c)
f0='%8.0f %12.4e\n';
dig0=64; dig1=120;
cond=zeros(N,1); lcond=zeros(N,1);
tic
ab0=sr_Enufin(dig0,N,nu,c);
ab1=sr_Enufin(dig1,N,nu,c);
toc
err=subs(abs((ab0-ab1)./ab1));
eps0=.5*10^(-dig0); c=err/eps0;
for n=1:N
  cond(n)=max(c(n,1),c(n,2));
  lcond(n)=log10(cond(n));
%  fprintf(f0,n,cond(n),lcond(n))
end
plot((2:N)',lcond(2:N));set(gca,'FontSize',14)
axis([0 42 -10 60])
