Index: tol_tests/tol/LinearAlgebra/PolMatrix/test_0002/check.tol
===================================================================
--- tol_tests/tol/LinearAlgebra/PolMatrix/test_0002/check.tol	(revision 747dbd9cb8b55afbc9d4eec2fdebff2f09d7e9d5)
+++ tol_tests/tol/LinearAlgebra/PolMatrix/test_0002/check.tol	(revision e1969fb93af47c5fa0d238f608965d3e9ac9963d)
@@ -3,7 +3,9 @@
 #Require MPM;
 
+Real PutRandomSeed(23453);
+
 Real dim = 2;
 Real deg = 2;
-Real T = 100;
+Real T = 1000;
 
 Matrix x = VariableChange::@StableMatrixPolyn::
@@ -23,7 +25,24 @@
 Real inequality_constraint_ = -Log(StationaryValue(determinant));
 
-Matrix X = Rand(T,dim,-1,1);
+Matrix X0 = Rand(deg,dim,-1,1);
+Matrix X  = Rand(T,  dim,-1,1);
 
-Matrix Y = P:X;
+Matrix Y = P:(X0<<X);
+
+Matrix Y_ = Group("ConcatColumns",For(1, dim, Matrix(Real i)
+{
+  SetSum(For(1, dim, Matrix(Real j)
+  { 
+    DifEq(PHI[i][j]/1,SubCol(X,[[j]]),SubCol(X0,[[j]]))
+  }))
+}));
+
+Matrix Y_t = Group("ConcatColumns",For(1, dim, Matrix(Real i)
+{
+  SetSum(For(1, dim, Matrix(Real j)
+  { 
+    DifEq(PHI[j][i]/1,SubCol(X,[[j]]),SubCol(X0,[[j]]))
+  }))
+}));
 
 
@@ -33,4 +52,5 @@
 Real eigen.module.max = MatMax(eigen.module);
 Real eigen.module.max.inv = 1/eigen.module.max;
+