integrate (2*x+y,y)

y*(4*x + y)/2

integrate(2*x + y, y)

plot(2*x*y + y**2/2)

solve(2*x*y + y**2/2, y)

diff(2*x*y + y**2/2, y)

series(2*x*y + y**2/2, y, 0, 10)

Experiment with SymPy at SymPy Live.