x=var('x')
y=function('y',x)
deq1=diff(y,x)+2*y-e^x*(3*sin(2*x)+2*cos(2*x))==0
r=desolve(deq1,[y,x],y)
show(r.simplify_full())

Kreyszig-1.6-36-ex-3 (last edited 2010-12-17 04:37:30 by mlec)