#Kreyszig Problem set 2.9 Question 16
var('x')
y=function('y',x)
print desolve(diff(y,x,2)-6*diff(y,x)+13*y==4*exp(3*x),[y,x],ics=[0,2,4])