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