This preview shows page 1. Sign up to view the full content.
Unformatted text preview: rithm, with no pivoting. It must compute L row-by-row and U column-by-column. Start with the following 2-by-2 block matrix form: " L 11 L 21 L 22 # " U 11 U 12 U 22 # = " A 11 A 12 A 21 A 22 # where L 22 , U 22 , and A 22 , are each scalars. The L 11 , U 11 , and A 11 , are each ( n-1)-by-( n-1) in size. Derive the algorithm (show your work) and then write a MATLAB function. Write a test code that tests your LU factorization function. Check the norm of L*U-A , divided by the norm of A , using a range of square random test matrices (from rand ). 1...
View Full Document
This note was uploaded on 01/18/2010 for the course COT 4501 taught by Professor Davis during the Fall '08 term at University of Florida.
- Fall '08