import sympy as sp
from IPython.display import display

x = sp.symbols('x',real=True)
eq1=x*sp.sqrt(x+1)
display(sp.Eq(sp.Integral(eq1, x), sp.factor(sp.Integral(eq1, x).doit())))