y = var('y')   
u = function('u',y)   
sol=desolve(diff(u,y,y) - u, [u,y])
print sol
#Here, a PDE is asked to be solved as ODE, treating other variables as parameters
#So, we get the constant 'k1' and 'k2' in the solution as functions of other variables

SageDays/Kreyszig-11.1-20 (last edited 2010-08-12 10:19:43 by Aashita)