i1 : gndR = QQ[x,y,z];
i2 : subR = subring({x^2, y^2, z^2});
i3 : gens subR o3 = | x2 y2 z2 | 1 3 o3 : Matrix gndR <--- gndR