Question: how to evaluate numerical inverse Laplace transform and plot it?

I would like to apply inverse Laplace transform to U(x,p), which is defined by

For simplicity with my calculations, I assumed p:=i*beta^2. That is why I have the following equation after applying Laplace transform

(beta=0 is not a pole, that is why I removed the last term in my calculations later. Because there is no contribution) where

Here p and beta are complex values, we can write Re(p)=-2*Re(beta)*Im(beta), Im(p)=(Re(beta))^2-(Im(beta))^2 due to p:=i*beta^2. I numerically compute the roots of h(beta), you can find the numerical values of beta (I assumed digits are 50 due to accuracy ) betap.mw

Finally, I would like to plot U(x,t) with the values t=0.8, lambda=1, L=10, k=1. For checking the figure give t=0 and observe that U(x,0)=0.

I am expecting the plot is more or less like the following figure

PS: I already tried to solve and plot the problem, but I could not find where I make a mistake. I  share the worksheet below. Thank you!

complexplot.mw

Please Wait...