Linear algebraic equations can arise in the solution of differential equations. In this problem you will develop a MATLAB code to solve the heat transfer problem in a 1-dimensional bar. The following equation derives from a heat balance for a long, thin rod (see the figure): dx where Tis the temperature at location x along the rod, his the heat transfer coefficient between the rod and ambient air, and Ta is the ambient temperature. This equation can be transformed to a set of linear algebraic equations by using a finite divided difference approximation for the second derivative T T-27+ T1 where Ti designate the temperature at node i. By substituting this equation into the first differential equation, we can have This equation can be writen for each node i and is linear in terms of "Ts. The temperature at the first and last node is fixed and is given by the boundary conditions. 1. First formulate the problem as Ax = b where x is the unknown temperatures at the nodes. Assume we have n nodes along the rod 2. use MATLAB to automatically create the A and b matrices based on the number of desired nodes, the boundary conditions, which are given as the input to your function. The rod has the length of 10 [m], heat transfer rate of h = 0.01 [m-2], and ambient temperature is T-20C 3. Plot the temperature distribution vs. x. Note that the figure below is an example with "n=6" nodes, and temperature at node zero is "TO 40" and at the end node "Tend 200". a 20 Your Function C Reset EE MATLAB Documentation function TmyFunction (n, Te, Tend) % n is the number of nodes along the rod % Te is the temperature at node zero at x-e % Tend is the temperature at the last node at x=L % is a vector of Temperature starting at node zero and going to the last node with Tend % T should have n elements, as the temperature at each node L 10; % [m] % Write your code here: end Code to call your function C Reset 1 my Function(x)