%RUN_SINERFC
%
global ab
digits(32); dig=32;
ab=loadvpa('ab_hrhermite',dig,200,2);
disp('      x    n        y           err       N')
errmax=0;
for x=0:.2:1
%for x=1.5:.5:5
%for x=5.5:.5:8
%for x=8.5:.5:12
%for x=13:19
%for x=20:27
%for x=-.5:-.5:-5
  for n=0:25:150
%  for n=0:10:50
    [y,N]=sinerfc(dig,n,x);
    y1=squad_inerfc(dig,N+5,n,x);
    if N+5>200, break; end
    err=subs(abs((y-y1)/y1));
    fprintf('%8.1f %4.0f %12.4e %12.4e %4.0f\n',x,n,subs(y),err,N)
    if err>errmax, errmax=err; end
  end
  fprintf('\n')
end
errmax
