I have an expression like
e_p = simplify(e_p) print(type(e_p)) print(e_p) -(-(b[omega] - 1)*k[i, omega]*r[i]/(b[omega]*w[