Example program for dgelsy
This example solves the linear least squares problem
min
x
b
-
A
x
2
for the solution,
x
, of minimum norm, where
A
=
-
0.09
0.14
-
0.46
0.68
1.29
-
1.56
0.20
0.29
1.09
0.51
-
1.48
-
0.43
0.89
-
0.71
-
0.96
-
1.09
0.84
0.77
2.11
-
1.27
0.08
0.55
-
1.13
0.14
1.74
-
1.59
-
0.72
1.06
1.24
0.34
and
b
=
7.4
4.2
-
8.3
1.8
8.6
2.1
.
A tolerance of
0.01
is used to determine the effective rank of
A
.
Note that the block size (NB) of
64
assumed in this example is not realistic for such a small problem, but should be suitable for large problems.
f90 source code
Input data file
Expected results