y = var('y')   
u = function('u',y)   
sol=desolve(diff(u,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 'c' in the solution as a function of other variables

SageDays/Kreyszig-11.1-15 (last edited 2010-08-12 10:16:55 by Aashita)