# part a
y = var('y')    # define a variable y
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=(0,1)); f
plot(f)

# part b
y = var('y')    # define a variable y
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=(0,2)); f
plot(f)