x = var('x')
y = function('y',x)
eq1 = diff(y,x,4)==0
soln = desolve_laplace(eq1 ,  [y,x] , ics=[1,16,-4,24])
soln
plot(y)

Kreyszig-2.14-11-U (last edited 2010-12-17 11:37:47 by harsha)