Book
Advanced Engineering Mathematics
Author
Erwin Kreyszig
Edition
8th Edition
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"
Solution by:
- Zubin Mehta, Student, IIT Bombay