using System; using System.Collections.Generic; using System.Linq; using System.Text; using L = Science.Mathematics.LinearAlgebra; namespace Strang3Ed.Chapter02.Section2 { public class Example02 { public Example02() { } private string result; public string Result { get { return result; } } public void Compute() { double[,] x = {{1.0,-2.0000001}, {3.0,-6.0}}; L.Matrix A = new L.Matrix(x); L.FactorizationPAeqLU obj = new L.FactorizationPAeqLU(A); obj.Compute(); L.Matrix LU = obj.LowerTriangular * obj.UpperTriangular; L.Matrix PA = obj.RowPermutation * A; result += PA.ToString() + "\r\n"; result += LU.ToString() + "\r\n"; result += obj.LowerTriangular.ToString() + "\r\n"; result += obj.UpperTriangular.ToString() + "\r\n"; result += A.Determinant.ToString(); } } } /* 3 -6 1 -2.0000001 3 -6 1 -2.0000001 1 0 0.333333333333333 1 3 -6 0 -9.99999998363421E-08 2.99999999509026E-07 */