Reference no: EM13309453
Include all your programs and their input and output in ONE .doc or .pdf file. Also, include all your comments on the same file. Submit the file using turn-it-in to be found in the assessment menu on your L@G account. Please be reminded that turn-it-in allows submission only once. The marks for the two questions are 10 + 10.
1. Engineers often use thermodynamics extensively in their work. The following polynomial can be used in relating the zero-pressure specific heat of dry air cp kJ/(kg K) to temperature (K):
cp = 0.99403 + 1.671 * 10-4T + 9.7215 * 10-8T2 - 9.5838* 10-11T3 + 1.9520 * 10-14T4.
Write a MATLAB program that does the following things.
a) Convert the equation into f(T) = 0 when cp of 1.1kJ/kgK. In our MATLAB program, develop a plot of f(T) versus a range of T = 0 to 1200 K. MATLAB. On the axes, put proper labelling.
b) Use your plot to graphically select two points that bracket a solution. In your program write a bisection algorithm to find the value of T such that f(T) = 0 using the two selected points as your initial values.
c) Your algorithm will find the value of T up to 10 decimal places. Determine the level of the absolute approximate relative error to be achieved in this case. Use the error level as the stopping condition in your program.
d) In each iteration, your program should display the current root and the current value of relative approximate error.
e) Show your selected points and the root in the plot as well with their x values.
2. Consider the following system of linear equation:
2x1 - 6x2 - x3 = -38
-3x1 -x2 + 7x3 = -34
-8x1 + x2 - 2x3 = -20
a) Check whether system has no solution, infinite number of solutions or ill-conditions. Perform this checking from the coefficient and constant matrices.
b) Manually solve the system of linear equation using Gauss elimination method with partial pivoting.
Consider 3 decimal places. Show all steps of your computation in detail.
c) Once you obtain the solution, go to MATLAB's command window to compute the solutions by using Cramer's Rule and verify the correctness of the solution.
d) Considering the values obtained from MATLAB as true value, compute the absolute relative errors for the values that you manually calculated for each of variables.
e) Include all your computation details, and MATLAB command window activities in your file to be submited.