Book
Advanced Engineering Mathematics
Author
Erwin Kreyszig
Edition
8th Edition
x = var('x') y=function('y',x) DE = x^2*diff(y,x,2)-3*x*diff(y,x)+3*y-3*log(x)+4 print "GS:", desolve(DE, [y,x], ics=[1,0,1]) #check for particular solution yp = log(x) V = x^2*diff(yp,x,2)-3*x*diff(yp,x)+3*y-3*log(x)+4 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