x = var('x')
y = function('y',x)
deq1 = diff(y,x,3)+3*diff(y,x,2)+3*diff(y,x)+y-30*e**(-x)==0
soln = desolve_laplace(deq1 ,  [y,x] , ics=[3,-3,-47])
soln

Kreyszig-2.15-6 (last edited 2010-12-17 10:20:54 by harsha)