i1 : R = ZZ/7[x,y,u,v]/ideal(x*y-u*v) o1 = R o1 : QuotientRing
i2 : D = divisor(x) o2 = 1*Div(v, x) + 1*Div(u, x) of R o2 : WDiv
i3 : getGBList(D) o3 = {{v, x}, {u, x}} o3 : List