% TempProfile_v2
% Coded by Nigel F. Reuel on 9.13.2016
% This code demonstrates contourf, nested loops, animation
% also has an anonymous function!
H = 100; %mm
L = 200; %mm
kappa = 23; % mm^2/s
Ufunc = @(x,y,t) 1600/pi^2*sin(pi*x/L)*sin(pi*y/(2*H))*exp(-(1/(4*H^2)+1/L^2)*kappa*pi^2*t);
%
Temp = zeros(100,200);
%
for k = 1:100  %<-- Outer loop is for the time
    t = k;  
    %
    for i = 1:100
        for j = 1:200
            x = j;
            y = i;
            Temp(i,j) = Ufunc(x,y,t);
        end
    end
    contourf(Temp)
    colorbar
    caxis([0 150])
    title(['Time =',num2str(k),'s'])
    pause(.1)
end