proc choosesubsetsv1(int n,int k)
{

  if(k<0 or k>n)
  {
    return(L);
  }

  // the list of all non empty subsets out of 1,...,n
  list l = Realrad::subset(n); 
  list L,li;

  int i,j;

  for(i=1;i<=size(l);i++)
  {
   li = l[i];
   if (size(li)==k)       // select those subsets with k elements
   {
    j++;
    L[j]= intvec(li[1..k]);   // convert this list to intvec
   }
  }
  return(L);
}

LIB "realrad.lib";
ring r0 = 0,x,dp;
> choosesubsetsv1(5,3);
[1]:
   1,2,3
[2]:
   1,2,4
[3]:
   1,3,4
[4]:
   2,3,4
[5]:
   1,2,5
[6]:
   1,3,5
[7]:
   2,3,5
[8]:
   1,4,5
[9]:
   2,4,5
[10]:
   3,4,5

proc choosesubsets(int n, int k)
"USAGE: choosesubsets(n,k); n, k int 
RETURN: list,  all subsets of 1,...,n with k elements
EXAMPLE: example choosesubsets; shows an example
"
{
  list L;
  int i,j,m;
  intvec v,w;

  int flipped = (n-k) < k;
  if (flipped) { k = n-k;}

  ring r= 0,(x(1..n)),dp;
  ideal J;

  for (i=1;i<=n;i++) {J[i] = x(i)^2;}  // J is a std-basis
  attrib(J,"isSB",1);
  ideal K = reduce(maxideal(k),J)+0;   //

  for(j=1;j<=size(K);j++)     // transform the exponents to the indices
  {
     v = leadexp(K[j][1]);
     m = 0;
     for(i=1;i<=n;i++)
     {
       if(v[i]==!flipped)
       {
         m++;
         w[m]=i;
       }   
     }
     L[j] = w;
  }
return(L);
}