t = var('t')    # define a variable x
r = function('r',t)   # define y to be a function of that variable
DE = diff(r,t) == -2*t*r
print 'Differential Equation:'
show(DE)
soln = desolve(DE, [r,t])
print 'Solution:'
show(soln)

c = solve(soln.subs(t=0)==2.5, soln.variables()[0])
show(c)

print 'Explicit solution:'
soln.subs(c=c[0].rhs())

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