# 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)

SageDays/Kreyszig-1.2-1 (last edited 2010-08-09 14:38:22 by dassarun)