using System; using System.Collections.Generic; using System.Linq; using System.Text; using L = Science.Mathematics.LinearAlgebra; namespace Strang3Ed.Chapter02.Section6 { public class Example02 { public Example02() { } private string result; public string Result { get { return result; } } public void Compute() { double[,] x = {{1.0,1.0,0.0,0.0}, {1.0,2.0,1.0,0.0}, {0.0,1.0,2.0,1.0}, {0.0,0.0,1.0,2.0}}; L.Matrix A = new L.Matrix(x); L.FactorizationPAeqLU obj = new L.FactorizationPAeqLU(A); obj.Compute(); L.Matrix L = obj.LowerTriangular; result += L.ToString(); L.Matrix U = obj.UpperTriangular; result += U.ToString(); L.Matrix P = L * U; result += P.ToString(); } } } /* 1 0 0 0 1 1 0 0 0 1 1 0 0 0 1 1 1 1 0 0 0 1 1 0 0 0 1 1 0 0 0 1 1 1 0 0 1 2 1 0 0 1 2 1 0 0 1 2 */