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"

SageDays/Kreyszig-2.8-13 (last edited 2010-08-11 13:26:29 by zubinmehta)