x = var('x')    # define a variable x
y = function('y',x)   # define y to be a function of that variable
DE = diff(y,x) == 1 + 4*y^2
print 'Differential Equation:'
show(DE)
soln = desolve(DE, [y,x], ics=(0,0))
print 'Solution:'
show(soln)

print 'Explicit solution:'
solve(soln,y)

SageDays/Kreyszig-1.3-16 (last edited 2010-08-09 14:04:05 by PankajPandey)