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

SageDays/Kreyszig-1.1-7 (last edited 2010-08-09 14:26:26 by GauravLuthra)