public void USolve()
{
di = new double[] { 1, 2, 3 };
al = new double[] { 1, 2, 3 };
au = new double[] { 3, 2, 1 };
ia = new int[] { 1, 1, 2, 4 };
skylineMatrix = new SkylineMatrix(di, ia, al, au);
Vector vector = new Vector(new double[] { 4, 4, 1 });
var result = skylineMatrix.USolve(vector, false);
Vector resultActual = new Vector(new double[] { -7, 3, 1 });
for (int i = 0; i < result.Size; i++)
{
Assert.Equal(result[i], resultActual[i], 8);
}
}