Code: Select all
import sympy as sp
a,b,c,d = sp.symbols('a, b, c, d', positive = True)
expr = (a*b*c+2*d*a*c+b*c*d)/(2*a+b)
correct_simpl = c*(d + a*b/(2*a+b))
display(expr)
display(expr.simplify())
display(correct_simpl)
correct_simpl.equals(expr)