function theat1d3_loop(n,dt,loop)

u = zeros(n,1);

for k=1:loop

    u = theat1d3(u,n,dt);

    plot(u)
    axis([1 n 0 80])
    pause(dt)
end

end