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

SageDays/Kreyszig-2.8-4 (last edited 2010-08-11 12:03:46 by zubinmehta)