% 正方形領域上の関数の表示
function plot_d(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(1,2,1);
  colormap hsv;
  mesh(x,y,vv);
% 右側に等高線
  subplot(1,2,2);
  contour(x,y,vv);
  end
