x, a, b = symbols('x a b')
fx = '3*x - 3*x**3'
gx = 'x'
m = 1 / 3
int_fx = integrate(fx, (x, a, b))
int_gx = integrate(gx, (x, a, b))
solve([Eq(int_fx, m), Eq(int_gx, m)], (a, b))
решение для пределов интегрирования численно?