%SNPLUS_NEG  Symbolic analogue of Nplus_neg.m
%
global ab ab0
digits(32); dig=32;
ab=loadvpa('ab_hrhermite',dig,200,2);
dN=1; eps0=.5e-30; Np=zeros(6,1); np=(0:5)';
N0=10;
for x=-1:-1:-5
  i=0; tic
  for n=0:10:50
    i=i+1;
    Np(i)=sN_inerfc(dig,n,x,eps0,N0,dN);
  end
  toc
  N0=Np(1)-1;
  plot(np,Np);set(gca,'FontSize',14)
  axis([0 5 10 180])
  hold on
end
