Project 2 Solving Linear Systems
This project has two parts .
1. A Matlab function LowerSolve for solving lower triangular systems by
column oriented forward substitution
and the Matlab script TestLowSolve.m for testing LowerSolve are provided on the course
web page. Modify LowerSolve to produce a Matlab function UpperSolve that solves upper triangular
systems using row oriented backward substitution . Modify TestLowSolve.m to give a Matlab script
TestUppSolve.m which compares the result of using both your function UpperSolve and the Matlab
backslash operator on randomly generated upper triangular systems of sizes 4, 8 and 16.
You must submit printouts of TestUppSolve.m and UpperSolve.m, and the output
from running the
script TestUppSolve.m. Also, provide marked up listings of TestLowSolve.m and LowerSolve showing
the changes that you made .
2. You are to compute solutions to systems of linear equations , to calculate
the errors in the computed
solutions, and to calculate the residuals for the computed solutions. Each linear system has n equations
in the n unknowns x 1, x2, . . . , xn.
that is, in matrix vector form, where
The matrix A is to be set to the Hilbert matrix of order n . Use the Matlab command hilb to build A.
One way to form
is first to form
using the above values, then to compute
using the simple
b = A * x
In a single program, for linear systems of size n = 6, 8, 10 and 12 in turn.
(a) Calculate the matrix A, the true solution
and right hand side
(b) Find an approximate solution of the linear system using the Matlab backslash operator
to solve the linear systems, and print in a table the maximum magnitude of the components of
(c) Compute and print the maximum magnitude of the components of the residual
The Matlab script vande on the course web page may be of assistance as a template for this project.
To complete this part of the assignment, you must submit.
•A printout of the program that you used to compute the approximate
solutions, errors and
•The output containing the approximate solutions, errors and residuals.
•A brief report outlining the problem and what you did to solve it. (Either give a list of the
changes you made to vande or provide a marked up copy of vande showing your changes.) Your
report should discuss the magnitudes of the errors and residuals, their behavior as the value of n
increases , and the probable cause of their sizes.