x = var('x')   
y = function('y',x)   
sol=desolve(diff(y,x) + 4*y - cos(x), [y,x])
sol