%SNPLUS_POS  Symbolic analogue of Nplus_pos.m
%
global ab ab0
digits(32); dig=32;
ab=loadvpa('ab_hrhermite',dig,200,2);
dN=1; eps0=.5e-30;
n0=0;
%n0=25:
%n0=50;
%n0=75;
%n0=100;
%n0=125;
i=0; N0=10;
for x=0:3:27
  if x<=14.4
    tic
    i=i+1; xp(i)=x;
    Nm=0;
    for n=n0:5:n0+25
      N=sN_inerfc(dig,n,x,eps0,N0,dN);
      [N (n-n0)/5+1]
      if N>Nm, Nm=N; end
      N0=N-1;
    end
    Np(i)=Nm;
    Nm
    toc
  else
    nmax=floor((27-x)/.084);
    if nmax>=n0
      tic
      i=i+1; xp(i)=x;
      Nm=0;
      for n=n0:5:min(n0+25,nmax)
        [N,y]=sN_inerfc(dig,n,x,eps0,N0,dN);
        [N (n-n0)/5+1]
        if N>Nm, Nm=N; end
        N0=N-1;
      end
      Np(i)=Nm;
      Nm
      toc
    end
  end
end  
[xp' Np']
plot(xp,Np);set(gca,'FontSize',14)
