% 正方形領域上の関数の表示(Neumann境界条件用)
% 使用例
%  [v,d]=eigs(eigp2nsp(nx,ny),10,0); 小さい方から１０個の固有値、固有関数を求める
%  plotn(v(:,10),nx,ny);             最小の固有値に属する固有関数を表示
function plot_nzb(v,nx,ny)
  % vvは ny+1,nx+1
  vvv=zeros(nx+1,ny+1);
  vvv(:)=v;
  vv=vvv';
  vv(1,:)=vv(1,:)*sqrt(2);
  vv(ny+1,:)=vv(ny+1,:)*sqrt(2);
  vv(:,1)=vv(:,1)*sqrt(2);
  vv(:,nx+1)=vv(:,nx+1)*sqrt(2);
  x=0:1/nx:1;
  y=0:1/ny:1;
% 右側に等高線
  contour(x,y,vv,[0 0],'-k');
  end
