specGu

 computes and plots neig eigenvalues nearest to origin in figure p.ifig

0001 function muv=specGu(p,varargin)
0002 % computes and plots neig eigenvalues nearest to origin in figure p.ifig
0003 r=resi(p,p.u,p.lam); Gu=getGu(p,p.u,p.lam,r);
0004 if (nargin==2)
0005     dum=p.neig; p.neig=varargin{1}; dum2=p.eigmeth; p.eigmeth='eigs';
0006     [ineg,muv]=spcalc(Gu,p); p.neig=dum; p.eigmeth=dum2;
0007 elseif (nargin==3)
0008     dum=p.eigint; p.eigint=[varargin{1},varargin{2}]; 
0009     dum2=p.eigmeth; p.eigmeth='sarn';
0010     [ineg,muv]=spcalc(Gu,p); p.neig=dum; p.eigmeth=dum2;
0011 end
0012 figure(p.ifig); plot(real(muv),imag(muv),'*'); axis tight; 
0013 fprintf('#-EV = %i\n',ineg);
0014 end