y = var('y')    # define a variable x
x = function('x',y)   # define x to be a function of that variable
DE = diff(x, y) - x*y
f=desolve(DE, [x,y],ics=(1,2)); f
plot(f)