blob: 60c4aca7fd06855f0b54be697246b56b509f6f78 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
|
// EXAMPLE (PG 512)
A = [1 2 1;2 2 3;-1 -3 0] // Coefficient matrix
b = [0 3 2]' // Right hand matrix
[l,u] = lu(A)
// l is lower triangular matrix & u is upper triangular matrix
l*u
if(A==l*u)
disp('A = LU is verified')
end
det(A)
det(u)
if(det(A)==det(u))
disp('Determinant of A is equal to that of its upper triangular matrix')
// Product rule of determinants is verified
|
