x = var('x')
y = function('y', x)

DE = 4*diff(y,x,x) - 15*diff(y,x) - 4*y

desolve(DE, [y,x])