function theat2d3_loop(m,dt,loop) % 温度分布を2次元グラフで表示 n=m^2; u = zeros(n,1); for k=1:loop u = theat2d3(u,m,dt); 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) axis([1 m 1 m 0 40]) pause(dt) end