function theat1d2_loop(n,dt,loop)

u = zeros(n,1);

for k=1:loop

    u = theat1d2(u,n,dt);

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

end