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

SageDays/Kreyszig-1.2-5 (last edited 2010-08-10 11:42:20 by dassarun)