x = var('x')    # define a variable x
y = function('y',x)   # define y to be a function of that variable
DE = diff(y,x) == - x/y
print 'Differential Equation:'
show(DE)
soln = desolve(DE, [y,x], ics=(sqrt(3),1))
print 'Solution:'
show(soln)
print 'Explicit solution:'
solve(soln,y)

SageDays/Kreyszig-1.3-12 (last edited 2010-08-09 14:00:55 by PankajPandey)