%RUN_INERFC
%
global ab0
f0='%8.1f %4.0f %12.4e %12.4e %4.0f\n';
digits(32); dig=32;
ab=loadvpa('ab_hrhermite',dig,200,2);
ab0=subs(ab);
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:-17.5
  for n=0:25:150
%  for n=0:10:50
    if x>27-.084*n|x<-17.9+.024*n, break; end
    [y,N]=inerfc(n,x);
    N1=N+5; 
    if N1>200, break; end
    y1=quad_inerfc(N1,n,x);
    err=abs((y-y1)/y1);
    fprintf(f0,x,n,y,err,N)
    if err>errmax, errmax=err; end
  end
  fprintf('\n')
end
errmax
