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