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

Kreyszig-2.14-17-U (last edited 2010-12-18 07:50:12 by harsha)