MathGroup Archive 2006

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: does anybody know how to find the inverse Laplace transform of this wierd thing?

  • To: mathgroup at smc.vnet.net
  • Subject: [mg64127] Re: [mg64100] does anybody know how to find the inverse Laplace transform of this wierd thing?
  • From: Daniel Lichtblau <danl at wolfram.com>
  • Date: Thu, 2 Feb 2006 00:05:37 -0500 (EST)
  • References: <200602010934.EAA23043@smc.vnet.net>
  • Sender: owner-wri-mathgroup at wolfram.com

gino wrote:
> Want to find the inverse Laplace transform of the following term:
> 
> H(s)=1/s^2*exp(s^2*a^2/2)*integrate(exp(-u^2/2), u from s*a to +infinity)
> 
> How to do that?
> 
> ------------------------------
> 
> Making relaxation to the problem, if I have to find only certain sampled 
> values of the inverse Laplace transform of H(s), let's denote it as h(t),
> 
> I just need to find h(1), h(2), h(3), etc.
> 
> Is there a short cut for it?
> 
> Thanks a lot! 
> 

We start with your input.

In[9]:= InputForm[hH[s_] =
         1/s^2*Exp[s^2*a^2/2]*Integrate[Exp[-u^2/2],
         {u,s*a,Infinity}, Assumptions->Element[a,Reals]]]
Out[9]//InputForm= (E^((a^2*s^2)/2)*Sqrt[Pi/2]*Erfc[(a*s)/Sqrt[2]])/s^2

Then the inverse Laplace transform is computed as below.

In[11]:= InputForm[InverseLaplaceTransform[hH[s], s, t,
          Assumptions->Element[a,Reals]]]
Out[11]//InputForm=
Sqrt[Pi/2]*(a*(-1 + E^(-t^2/(2*a^2)))*Sqrt[2/Pi] + t*Erf[t/(Sqrt[2]*a)])

You will need to decide for yourself whether this is what you had in mind.

Daniel Lichtblau
Wolfram Research


  • Prev by Date: New MathCode Functionality Expands Options for Prototyping C++ and Fortran 90 Applications in Mathematica
  • Next by Date: Re: Utilizing the Result From Solve[]
  • Previous by thread: does anybody know how to find the inverse Laplace transform of this wierd thing?
  • Next by thread: Re: does anybody know how to find the inverse Laplace transform of this wierd thing?