Efficient Laplace Transform Solution for Sin Square Root of T

[mechprog]

Laplace Transform of...

Few months back my friend called me to help him (and his whole class!) finding the Laplace Transform of \sin \sqrt{t}.
Since then I found myself wasting hours and days (even before non-mathematics exams!) trying to solve the problem. I also checked the provided answer by online mathematica (wolfram alpha) and found it to be correct.
One of the most apparent way is to use Euler's formula, it does not work, however. I don't even remember the different ways I tried.
When my friend asked his instructor (who gave him the problem), he was quick to suggest the use of series expansion, the method I personally dislike.
Is there no other better way?

 

[Ben Niehoff]

You don't need a series expansion. It can be done using contour integration. I won't just give you the answer, but here's a description of how to do it. You should follow along with your own work so you can learn how to do it. Then check against your Mathematica answer:

We want to integrate

\int_0^{\infty} \sin{\sqrt{t}} \; e^{-st} \; dt

First do a substitution u = \sqrt{t}. This will put an e^{-su^2} inside your integral, but don't fret. There will also be some additional stuff in there; keep track of it.

Next, you'll find that the integrand in terms of u is an even function. So, instead of integrating over u from 0 to \infty, you can integrate from -\infty to \infty and then divide by 2. This is good, because we know how to integrate e^{-su^2} from -\infty to \infty.

Now, express \sin u as \Im \{ e^{iu} \}. You can pull the \Im all the way outside your integral, since everything else in the integrand is a real function.

Now you can combine your exponentials into a single exponential of a quadratic polynomial in u. Complete the square in this polynomial, so you have

e^{-A(u - \lambda)^2 + B}[/itex]<br /> <br /> where A, B, and \lambda are functions of s only.<br /> <br /> Now do a second substitution v = u - \lambda. You should end up with a (finite) sum of terms of the form<br /> <br /> \Im \left\{ C_k(s) \int_{-\infty + i\alpha}^{\infty + i\alpha} v^k e^{-A(s) v^2} \; dv \right\}<br /> <br /> where \alpha is the imaginary part of \lambda. This integral can be done by contour integration: we can simply shift the contour downward to the real axis, since the integrand has no poles. Effectively, we can set \alpha = 0. Then we&#039;re left with a standard Gaussian integral, which you should look up if you don&#039;t know how to evaluate it.

 

[mechprog]

Thanks for solution! It worked elegantly!
One interesting thing about this is that the transform is e^{\frac{-1}{4s}} times the transform of \sqrt{t} .
Is there any link/hint to some general formula/theorem?