% 正方形領域上の関数の表示
function plot_d3(v,nx,ny)
  vv=zeros(ny+1,nx+1);
  vvv=zeros(nx-1,ny-1);
  vvv(:)=v;
  vv(2:ny,2:nx)=vvv';
  x=0:1/nx:1;
  y=0:1/ny:1;
% 左側にグラフの鳥瞰図
  subplot(2,2,[1 3]);
  colormap hsv;
  mesh(x,y,vv);
% 右側に等高線
  subplot(2,2,2);
  contour(x,y,vv);
  subplot(2,2,4);
  contour(x,y,vv,[0 0]);
  end
