private static void FillFreeVariables(LinearEquations equations)
{
int numberFreedomVariables = equations.Matrix.Cols - equations.Matrix.Rows;
for (int i = 1; i < numberFreedomVariables + 1; i++)
{
equations.Variables[^ i] = 1.0;
public static Vector FindAnswer(double[,] leftSide, double[] rightSide)
{
LinearEquations equations = new LinearEquations(leftSide, rightSide);
for (int i = 0; i < equations.Matrix.Rows; i++)
{
(int row, int col) = equations.Matrix.FindIndexMaxAbsItem(i, i);
equations.SwapRows(i, row);
equations.SwapColumns(i, col);
equations = ExcludeMaximumElement(equations, i);
}
if (equations.Matrix.Rows < equations.Matrix.Cols)
{
FillFreeVariables(equations);
}
return(FindVariables(equations));
}
private static LinearEquations ExcludeMaximumElement(LinearEquations equations, int startSubMatrix)
{
int j = startSubMatrix;
for (int i = j + 1; i < equations.Matrix.Rows; i++)
{
double multiplier = equations.Matrix[i, j] / equations.Matrix[j, j];
equations.Matrix[i, j] = 0;
for (int k = j + 1; k < equations.Matrix.Cols; k++)
{
equations.Matrix[i, k] -= multiplier * equations.Matrix[j, k];
}
equations.Vector[i] -= multiplier * equations.Vector[j];
}
return(equations);
}
private void Show(LinearEquations equations, string text)
{
Console.Out.WriteLine(text);
Console.Out.WriteLine("---------------");
for (int i = 0; i < equations.Matrix.Rows; i++)
{
for (int j = 0; j < equations.Matrix.Cols; j++)
{
Console.Out.Write("{0:0.###} ", equations.Matrix[i, j]);
}
Console.Out.Write(" | " + equations.Vector[i]);
Console.Out.WriteLine("");
}
Console.Out.WriteLine("_______________________");
for (int j = 0; j < equations.Matrix.Cols; j++)
{
Console.Out.Write("{0:0.###} ", equations.Variables[j]);
}
}