proc Substitute(int m, ideal I) 
{
  int i;
  int n = nvars(basering);
  int k;
  ideal P = I;
  if (m==0)
  {
    k = n/2;
  }
  else
  {
    k = n/(m+1);
  }
  for (i=0; i<k; i++)
    {
      if (m==0)
      {
	P = subst(P,var(i*2+1),1-var(2*i+2));
      }
      if (m==1)
	{
	  P = subst(P,var(i*(m+1)+1),1-3*var(i*(m+1)+2));
	}
      if (m==2)
	{
	  P = subst(P,var(i*(m+1)+1),1-2*var(i*(m+1)+2)-var(i*(m+1)+3));
	}
      if (m==3)
	{
	  P = subst(P,var(i*(m+1)+1),1-var(i*(m+1)+2)-var(i*(m+1)+3)-var(i*(m+1)+4));
	}
    }
  return(P);
}

proc Suma(ideal P)
{
  int i; 
  poly suma=P[1];
  for(i=2; i<=size(P); i++)
    { 
      suma=suma+P[i];
    }
  return(suma);
}