x=var('x')
y=function('y',x)
desolve(diff(y,x)+4*y-17*sin(x),y)