plot temperature distribution of a plane wall undergoing constant convection via a fluid on the right side and large heat flux from the left side for 0.2 seconds every 1 second in matlab

main.m
% Define parameters
k = 100; % thermal conductivity of the wall (W/mK)
rho = 7800; % density of the wall (kg/m^3)
cp = 450; % specific heat capacity of the wall (J/kgK)
h = 50; % convective heat transfer coefficient with fluid (W/m^2K)
Tinf = 300; % temperature of the fluid (K)
q = 50000; % heat flux from the left side (W/m^2)

L = 0.1; % length of the wall (m)
n = 50; % number of segments to discretize the wall
dx = L/n; % length of each segment (m)
A = dx*dx; % area of each segment (m^2)
t = 0:0.1:10; % time steps to iterate over

% Initialize temperature distribution
T = ones(1,n)*300; % initial temperature distribution of the wall (K)

% Solve heat equation using finite difference method
for i=1:length(t)
    dt = t(i) - t(i-1); % time step
    for j=1:n
        if j==1 % left boundary
            T(j) = T(j) + q/k*A*dt/(rho*cp*dx);
        elseif j==n % right boundary
            T(j) = T(j-1) - h*A*dt/k*(T(j-1)-Tinf)/dx;
        else % interior nodes
            T(j) = T(j) + k*dt/(rho*cp*dx*dx)*(T(j-1)-2*T(j)+T(j+1));
        end
    end
    
    % Plot temperature distribution of the wall
    figure(1); clf; hold on
    plot(linspace(0, L, n), T, 'r-', 'LineWidth', 2)
    xlabel('Distance (m)')
    ylabel('Temperature (K)')
    title(sprintf('Temperature distribution at t=%0.1f sec', t(i)))
    ylim([200 500])
    hold off
end
1357 chars
40 lines