{ reset() x=var('x') y=function('y',x) DE=diff(y,x,2) + 2*diff(y,x) +101*y -10.4*exp(x) q=desolve(DE,[y,x]) }