proc fitting(matrix M, int n)
{
  n=nrows(M)-n;
  if (n<=0) {return(ideal(1));}
  if ((n>nrows(M))||(n>ncols(M))) {return(ideal(0));}
  return(std(minor(M,n)));
}


proc flatteningStrat(matrix M)
{
   list L;
   int k,r;
   ideal F;
   while(1)
   {
      F=interred(fitting(M,r));
      if (F[1]==1) {return(L);}
      if (size(F)!=0) {k++; L[k]=F;}
      r++;
   }
}