x = var('x')    # define a variable x
y = function('y',x)   # define x to be a function of that variable
DE=diff(y, x) +y
f=desolve(DE, y,ics=(0,1/3)); f
plot(f)

SageDays/Kreyszig-1.2-3 (last edited 2010-08-09 14:40:48 by dassarun)