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

Running the code costs a lot of time, I need some suggestions to make faster and more accurate. Thanks!

sonkok.mw

Hi! I'm trying to find the way to plot the solution with series representation. I need some help to find the easiest way.

Note: I realized some typing errors, which do not change the question a lot ,and I corrected them.

plot.mw

I would like to plot the following singular double integral, but I cannot due to singularities...

where x>0, t=0.2 and m=0.2.

I defined f(y) function as f:=y->exp(-(y-4.68)^2/0.4):

I attached my file:
1st_try.mw

Thank you !

I would like integrate and plot the following double integral:

I enclosed tryingdf.mw