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

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

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

            double[,] y = { {-1.0, 0.0, 0.0, 0.0}, 
                            {0.0, 1.0, 0.0, 0.0}, 
                            {0.0, 0.0, -1.0, 0.0}, 
                            {0.0, 0.0, 0.0, 1.0} };

            L.Matrix M = new L.Matrix(y);

            L.Matrix B = A.Similar(M);

            result += B.ToString() + "\r\n";

            result += A.Inverse.ToString() + "\r\n";
        }
    }
}
/*
1     0     0     0     
-1     1     0     0     
1     -2     1     0     
-1     3     -3     1     


1     0     0     0     
-1     1     0     0     
1     -2     1     0     
-1     3     -3     1    
*/