x = var('x')
y=function('y',x)
DE = diff(y,x,2)-y-8*exp(-3*x)
print "GS:", desolve(DE, [y,x])
#check for particular solution
yp = exp(-3*x)
V = diff(yp,x,2)-yp-8*exp(-3*x)
if V == 0:
    print "yp is a particular solution"
else:
    print "yp is not a particular solution"

SageDays/Kreyszig-2.8-1 (last edited 2010-08-11 11:35:34 by zubinmehta)