x = var('x')
y = function('y',x)
DE = diff(y, x) + y - x^2+2
a=desolve(DE, y,x)
print a
print
# validating the output
s= e^(-x)+x^2-2*x
diff(s,x) + s - x^2+2

SageDays/Kreyszig-1.1-5 (last edited 2010-08-09 14:19:31 by GauravLuthra)