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)
Need more control? Try SymPy Live.