x = var('x')
y = function('y',x)
deq1 = y*diff(y,x,4)-y ==0
soln = desolve_laplace(deq1 ,  [y,x] , ics=[-1,7,-1,7])
soln
plot(y)