x = var('x')    
y = function('y',x)   
DE= diff(y,x,2)+(9*y)
desolve(DE,[y,x],ics=[0,4,-2])

SageDays/Kreyszig-2.1-17 (last edited 2010-08-10 13:09:14 by manojmanu)