automatic_names('true')
y=function('y',x)
d=diff(y,x,3)+3*diff(y,x,2)+3*diff(y,x,1)+y-30*(e^(-x))==0
desolve_laplace(d,[y,x],ics=(0,3,-3,-47))


sol:
        

5*x^3*e^(-x) - 25*x^2*e^(-x) + 3*e^(-x)

Kreyszig-2.15_1 (last edited 2010-12-18 09:26:33 by eshwara)