heat transfer by convection for wall in matlab

main.m
% Input parameters
T_room = 25;              % [C] Room temperature
T_wall = 100;             % [C] Wall temperature
h = 10;                   % [W/m^2-C] Convective heat transfer coefficient
k = 0.5;                  % [W/m-C] Thermal conductivity of the wall
L = 0.1;                  % [m] Thickness of the wall
A = 1;                    % [m^2] Area of the wall
t_final = 3600;           % [s] Final time
dt = 1;                   % [s] Time step

% Initialization
t = (0:dt:t_final)';
N = length(t);
T = zeros(N,1);
T(1) = T_wall;

% Main loop
for i = 2:N
  q_conv = h*A*(T(i-1)-T_room);   % [W/m^2] Heat transfer by convection
  q_cond = k*A/L*(T(i-1)-T_room); % [W/m^2] Heat transfer by conduction
  q = q_conv - q_cond;
  dTdt = q/(A*L);                 % [C/s] Rate of temperature change
  T(i) = T(i-1) + dTdt*dt;        % [C] Temperature
end

% Plotting
plot(t,T)
xlabel('Time [s]')
ylabel('Temperature [C]')
title('Temperature Evolution of Wall with Convection Heat Transfer')
989 chars
31 lines