%CONDENU
%
function [cond,lcond]=condEnu(N,nu)
f0='%8.0f %12.4e %12.4e\n';
dig0=48; dig1=64; 
cond=zeros(N,1); lcond=zeros(N,1);
tic
ab0=sr_Enu(dig0,N,nu); 
ab1=sr_Enu(dig1,N,nu);
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))
