using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using L = Science.Mathematics.LinearAlgebra;

namespace Strang3Ed.Chapter06.Section7
{
    public class ExampleA
    {
        public ExampleA()
        {
        }
        private string result;
        public string Result
        {
            get { return result; }
        }
        public void Compute()
        {
            double[,] x = { {1.0, 1.0, 0.0}, 
                            {0.0, 1.0, 1.0} };

            L.Matrix A = new L.Matrix(x);

            L.FactorizationAeqUSVt obj
                = new L.FactorizationAeqUSVt(A);
            obj.Compute();

            L.Matrix B = obj.ColumnOrthogonal * obj.Diagonal * obj.Orthogonal.Transpose;

            result += B.ToString() + "\r\n";
        }
    }
}
/*
1     1     1.11022302462516E-16     
1.66533453693773E-16     1     1     
*/