using System;
using L = Science.Mathematics.LinearAlgebraBasic;

namespace ScienceTest.MathematicsTest.LinearAlgebraTest
{
	/// <summary>
	/// SingularLinearEquationTest
	/// </summary>
	public class SingularLinearEquationTest
	{
		public SingularLinearEquationTest()
		{
		}
		private string result;
		public string Result
		{
			get{return result;}
		}
		public void Compute()
		{
			double[,] a = new Double[4,3];
			a[0,0]=2.0;
			a[0,1]=4.0;
			a[0,2]=3.0;
			a[1,0]=1.0;
			a[1,1]=4.0;
			a[1,2]=3.0;
			a[2,0]=1.0;
			a[2,1]=4.0;
			a[2,2]=8.0;
			a[3,0]=2.0;
			a[3,1]=8.0;
			a[3,2]=11.0;
			L.Matrix mm = new L.Matrix(a);

			double[] c = new Double[4];
			c[0] = 5.3;
			c[1] = 6.3;
			c[2] = 7.3;
			c[3] = 13.6;
			L.Vector vv = new L.Vector(c);

			L.SingularLinearEquation eq 
				= new L.SingularLinearEquation(mm,vv);
			eq.Solve();
			foreach(double k in eq.Solution)
				result += Convert.ToString(k)+"\r\n";		
		}
	}
}