SolomonTeraoIdeal(A: LIST, f: RINGELEM): IDEAL
/**/ use QQ[x,y]; /**/ A := [x, y, x-y]; /**/ f := x^2+y^2; /**/ SolomonTeraoIdeal(A, f); ideal(2*x^2 +2*y^2, 2*x*y^2 -2*y^3)