To compute the grobner basis of three polynomials the command is well known: R.ideal(x^2+y^2-1, (x+y)^2-1).groebner_basis(), where the ring is defined as R.<x,y,z> = PolynomialRing(QQ). I want to compute this for n polynomials $f_1, f_2,\ldots,f_n$ where $ n$ is as my choice, so how to use the loop for to solve this problem