MATLAB

clear % Define thermal diffusivity for aluminum alpha = 97.1 * 10^ (-6); Define the total time, in seconds totalTime 300; Define delta t, in seconds dt -0.1; Define number of nodes in the bar nodes-5; % Define the length of the bar, in meters len = 0.2; Calculate the distance between nodes, delta x dx = (len / (nodes - 1)); Define the number of time steps % The plus 1 is necessary because the time starts at t-0 sec timesteps = totalTime/dt + 1; % Define the Temperature Matrix, TM = zeros (timesteps, nodes); Suppose that the initial temperature of the bar is 20C, but starting at time t-0 seconds, the ends of the bar are heated to 200C and kept at a constant 200C for the duration of the experiment. Initialize TM to reflect the initial temperature conditions. % Initialize the temperature matrix % The first row is the initial 20 degrees c. The first and last columns are 200 degrees c. TM(1, :) = 20; TM(s, 1) = 200; TM (, nodes) 200; % initial temperature % constant temperature on left side % constant temperature on right side The commands above initialize the TM matrix with the values shown below 200 20 20 20 200