A = 2*b*(a**2*c**2 + b**2*s**2)/(a*d)
A*A_M reduced = 8*a**4*b**4/(a**2 + b**2)**2
free symbols of the reduced product: {b, a}
difference at two members (should be 0): 0
c2 (x^2 coeff) - c2 (y^2 coeff) = 0  ; xy coefficient = 0
A*c0 reduced = 8*a**4*b**4/(a**2 + b**2)**2   free symbols: {b, a}
A*c1x reduced = 0   free symbols: set()
A*c1y reduced = 0   free symbols: set()
A*c2 reduced = 0   free symbols: set()
