function mygraf2d(m,u)
% 温度分布を２次元グラフで表示

w = zeros(m);
for jj=1:m
  for ii=1:m
    i = ii + (jj - 1)*m;
    w(ii,jj) = u(i);
  end
end
surf(w)