create a for loop that accounts for cyclic transient heat conduction in a plane wall in matlab

main.m
% Define inputs
Tm = 20;
Ta = 5;
k = 0.5;
h = 10;
L = 0.2;
t = 0:1:3600;
alpha = k/(p*c);

% Preallocate memory for temperature matrix
T = NaN(numel(t), 2);

% Set initial conditions
T(1,:) = [Tm, Tm];

% Calculate temperature for each time step using for loop
for i = 2:numel(t)
    % Calculate temperature difference between inner and outer wall
    dT = T(i-1,1) - T(i-1,2);
    
    % Use cyclic boundary condition to update temperature matrix
    T(i,1) = T(i-1,1) - ((h*L)/(p*c))*(T(i-1,1) - Ta)*dt + ((k*L)/(p*c))*(dT/L)*dt;
    T(i,2) = T(i-1,2) + ((h*L)/(p*c))*(T(i-1,1) - Ta)*dt - ((k*L)/(p*c))*(dT/L)*dt;
end

% Plot results
plot(t, T(:,1), 'r', t, T(:,2), 'b');
xlabel('Time (s)');
ylabel('Temperature (°C)');
legend('Inner wall', 'Outer wall');
758 chars
31 lines