%
% TESTCONJ
%
f0='%6.2f %12.8f %6.0f\n';
n=50; db=.02;
aa=zeros(25,1); bb=zeros(25,1);
for ia=1:25
  a=.02+(ia-1)*.02;
  [b0,diff]=jacconj(n,a,db);
  fprintf(f0,a,b0,diff)
  aa(ia)=a; bb(ia)=b0; 
end