f0='%8.0f %5.0f %14.4e %12.4e\n';
disp('      N     dig      erra         errb')
for N=[10 25 50:5:100]
  abp=zeros(N,2);
  i=0;
  dig1=8; dig2=80; nd=1+(dig2-dig1)/4;
  err=zeros(nd,3);
  for dig=dig1:4:dig2
    i=i+1;
    abp0=abp;
    digits(dig);
    ab0=sr_jacobi01(dig,N,-1/2,0);
    gam=vpa(1/2);
    c=vpa(sqrt(1+(gam/2)^2));
    r0=vpa((2/gam)*sqrt(1+c));
    rho0=vpa(2*r0*atan(r0));
    ab=schri4(dig,N,ab0,c,rho0);
    abp=subs(ab);
    err(i,2:3)=max(abs(abp-abp0));
  end
  d=(dig1:4:dig2)';
  err(:,1)=d;
  for n=1:nd
    fprintf(f0,N,err(n,1),err(n,2),err(n,3))
  end
  fprintf('\n')
end
