#Kreyszig Problem set 2.9 Question 1
var('x')
y=function('y',x)
print desolve(diff(y,x,2)+4*y==sin(3*x),[y,x])