Computes the kernel of a matrix M numerically using singular value decomposition.
i1 : M = matrix {{1., 1, 1}}
o1 = | 1 1 1 |
1 3
o1 : Matrix RR <--- RR
53 53
|
i2 : numericalKernel(M, 0.01)
o2 = | -.57735 -.57735 |
| .788675 -.211325 |
| -.211325 .788675 |
3 2
o2 : Matrix RR <--- RR
53 53
|
Singular values less than the tolerance are treated as zero.
i3 : M = matrix {{1., 1}, {1.001, 1}}
o3 = | 1 1 |
| 1.001 1 |
2 2
o3 : Matrix RR <--- RR
53 53
|
i4 : numericalKernel(M, 0.01)
o4 = | .70693 |
| -.707284 |
2 1
o4 : Matrix RR <--- RR
53 53
|