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

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