x = var('x')
y = function('y',x)
deq1 = diff(y,x,3)+9*diff(y,x,2)+27*diff(y,x,1)+27*y==0
soln = desolve_laplace(deq1 ,  [y,x])
soln