The idealproduct command computes a basis for the product of ideals A and B. An optional argument P is the modulus, which should be a machine size prime.
A = [x-1, y^2]; B = [y]; idealproduct(A,B); [-y+x*y, y^3]