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

Kreyszig-2.14-7-U (last edited 2010-12-18 07:46:29 by shankar)