Thanks to visit codestin.com
Credit goes to www.scribd.com

Scribd
Upload
16 views2 pages

LU Decomposition

Uploaded by

gotikow387
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
16 views2 pages

LU Decomposition

Uploaded by

gotikow387
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 2

wxMaxima document 1 / 2

Q1

(%i6) $load(linearalgebra)$
'a=a:matrix( [5,2,1] , [3,7,4] , [1,1,9] );
b=b:matrix( [10],[21],[12]);
'A=A:lu_factor(a);
'x=lu_backsub(A,b);
C=C:get_lu_factors(A);
5 2 1

(%o2) a = 3 7 4
1 1 9

10

(%o3) b = 21
12

5 2 1
3 29 17
(%o4) A = [ 5 5 5 , [ 1 , 2 , 3 ] , generalring ]
1 3 245
5 29 29

(%o5) x = 2
1

1 0 0 5 2 1
1 0 0
3 29 17
1 0 0
(%o6) C = [ 0 1 0 , 5 , 5 5 ]
0 0 1 1 3 245
1 0 0
5 29 29

(%i7) L:second(C);
1 0 0
3
1 0
(%o7) 5

1 3
1
5 29

(%i8) U:third(C);
5 2 1
29 17
0
(%o8) 5 5

245
0 0
29
wxMaxima document 2 / 2

(%i9) P:first(C);
1 0 0

(%o9) 0 1 0
0 0 1

(%i10) inl:invert(L);
1 0 0
3
− 1 0
(%o10) 5

4 3
− − 1
29 29

(%i11) inu:invert(U);
1 2 1
−
5 29 245

5 17
(%o11) 0 −
29 245

29
0 0
245

(%i12) z:inl.P.b;
10
15
(%o12)
245
29

(%i13) x:inu.z;
1
(%o13) 2
1

You might also like