x = var('x')
y = function('y',x)
deq1 = 16*diff(y,x,4)-40*diff(y,x,2)+9==0
soln = desolve_laplace(deq1 ,[y,x])
soln