Index: tol_tests/tol/LinearAlgebra/PolMatrix/test_0001/check.tol
===================================================================
--- tol_tests/tol/LinearAlgebra/PolMatrix/test_0001/check.tol	(revision 0f7ecd14b1b7f646776844e785a34eb7bf35c3ab)
+++ tol_tests/tol/LinearAlgebra/PolMatrix/test_0001/check.tol	(revision 747dbd9cb8b55afbc9d4eec2fdebff2f09d7e9d5)
@@ -42,3 +42,12 @@
 }));
 
+Matrix x = GetNumeric([[X]]);
+Matrix xP = GetNumeric([[XP]]);
 
+Real quality.dif.xP = 1-MatMax(Abs(x-xP));
+
+Set X2 = DeepCopy([[X]],x*2);
+Matrix x2 = GetNumeric([[X2]]);
+Real quality.dif.x2 = 1-MatMax(Abs(x2-x*2));
+
+
Index: tol_tests/tol/LinearAlgebra/PolMatrix/test_0002/check.tol
===================================================================
--- tol_tests/tol/LinearAlgebra/PolMatrix/test_0002/check.tol	(revision 0f7ecd14b1b7f646776844e785a34eb7bf35c3ab)
+++ tol_tests/tol/LinearAlgebra/PolMatrix/test_0002/check.tol	(revision 747dbd9cb8b55afbc9d4eec2fdebff2f09d7e9d5)
@@ -1,16 +1,47 @@
 #Require VariableChange;
+#Require AlgLib;
+#Require MPM;
 
 Real dim = 2;
-Real deg = 3;
+Real deg = 2;
+Real T = 100;
 
-Matrix x = VariableChange::StableMatrixPolyn::
+Matrix x = VariableChange::@StableMatrixPolyn::
   SquareSetGen.Rand(dim, deg, -1, 1, 0.001);
 
-Set PHI = VariableChange::StableMatrixPolyn::
+Set PHI = VariableChange::@StableMatrixPolyn::
   DuffinSchurBinFact(dim, deg, (x));
 Polyn determinant = VecSer::PolynMatrixDeterminant(PHI);
-
 
 PolMatrix P = TabPol2PolMat(PHI);
 
 Polyn determinant_ = Determinant(P);
+Set PHI.det.roots = MPM::@ModuleVAR::PolynRootsReport(determinant_);
+Real root.abs.min = PHI.det.roots[1][3];
+Real P.isStationary = PolMatIsStationary(P);
+
+Real inequality_constraint_ = -Log(StationaryValue(determinant));
+
+Matrix X = Rand(T,dim,-1,1);
+
+Matrix Y = P:X;
+
+
+VMatrix C = CompanionMatrix(P);
+Matrix eigen = AlgLib::CppTools::rmatrixevd.eigen(VMat2Mat(C));
+Matrix eigen.module = Sqrt(eigen^2 * Col(1,1));
+Real eigen.module.max = MatMax(eigen.module);
+Real eigen.module.max.inv = 1/eigen.module.max;
+
+
+/*
+Real m = VRows(C);
+Real S=10000*m;
+VMatrix U = Gaussian(m,S,0,1);
+VMatrix V = C*U;
+VMatrix u2 = (U^2)*Constant(S,1,1/S);
+VMatrix v2 = (V^2)*Constant(S,1,1/S);
+VMatrix r2 = v2 $/ u2;
+//VMatrix R = MtMSqr(V) $/ MtMSqr(U);
+//Real VMatMax(R)
+/* */