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

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