automatic_names('True')
y=function('y',x)
d=diff(y,x,4)-diff(y,x,3)+100*diff(y,x,1)-100*y
desolve_laplace(d,[y,x])  

sol:
        

-1/30300*((100^(2/3) - 100^(1/3) + 1)*D[0, 0, 0](y)(0) - (100^(2/3) +
100*100^(1/3) - 100)*y(0) + 101*100^(1/3)*D[0](y)(0) - 101*D[0,
0](y)(0))*100^(1/3)*e^(-100^(1/3)*x) + 1/101*(100*y(0) + D[0, 0,
0](y)(0))*e^x + 1/9090000*(((2*(100^(2/3) - 50*100^(1/3) + 50)*y(0) -
(2*100^(2/3) + 100^(1/3) - 1)*D[0, 0, 0](y)(0) +
101*100^(1/3)*D[0](y)(0) - 101*D[0, 0](y)(0))*100^(2/3) - 2*((100^(2/3)
+ 2*100^(1/3) - 100)*D[0, 0, 0](y)(0) + 100*(100^(2/3) + 2*100^(1/3) +
1)*y(0) - 101*100^(2/3)*D[0](y)(0) - 202*100^(1/3)*D[0,
0](y)(0))*100^(1/3))*sqrt(3)*100^(2/3)*sin(1/2*sqrt(3)*100^(1/3)*x) +
300*(2*(100^(2/3) - 50*100^(1/3) + 50)*y(0) - (2*100^(2/3) + 100^(1/3) -
1)*D[0, 0, 0](y)(0) + 101*100^(1/3)*D[0](y)(0) - 101*D[0,
0](y)(0))*100^(1/3)*cos(1/2*sqrt(3)*100^(1/3)*x))*e^(1/2*100^(1/3)*x)

Kreyszig-2.14_35 (last edited 2010-12-18 10:51:26 by santosh)