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

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