umfpack_simple.m.out
1.03 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
umfpack_simple
UMFPACK_SIMPLE a simple demo
Example:
umfpack_simple
Copyright 1995-2007 by Timothy A. Davis.
UMFPACK License:
Your use or distribution of UMFPACK or any modified version of
UMFPACK implies that you agree to this License. UMFPACK is
is free software; you can redistribute it and/or
modify it under the terms of the GNU Lesser General Public
License as published by the Free Software Foundation; either
version 2.1 of the License, or (at your option) any later version.
A =
2 3 0 0 0
3 0 4 0 6
0 -1 -3 2 0
0 0 1 0 0
0 4 2 0 1
b =
8
45
-3
3
19
Solution to Ax=b via UMFPACK:
x1 = umfpack2 (A, '\', b)
x1 =
1.0000
2.0000
3.0000
4.0000
5.0000
Solution to Ax=b via MATLAB:
x2 = A\b
x2 =
1.0000
2.0000
3.0000
4.0000
5.0000
norm (x1-x2) should be small: 1.28037e-15
Type 'umfpack_demo' for a full demo of UMFPACK
diary off