x = var('x')
y = function('y',x)
DE = diff(diff(y, x),x) - x^(-4)
desolve(DE, y)

SageDays/Kreyszig-1.1-3 (last edited 2010-08-09 14:13:07 by GauravLuthra)