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=(1,2)); f
plot(f)

SageDays/Kreyszig-1.2-S (last edited 2010-08-09 14:35:03 by dassarun)