Re: numerical inversion of laplace transform

• To: mathgroup at smc.vnet.net
• Subject: [mg74861] Re: numerical inversion of laplace transform
• From: "Roman" <rschmied at gmail.com>
• Date: Mon, 9 Apr 2007 06:09:11 -0400 (EDT)
• References: <ev7j3g\$l2j\$1@smc.vnet.net>

```Dan,
Since your p[t] is a step function, your problem can be solved
exactly.
G[t] = Integrate[p[x]*F[t-x], {x,-Infinity,Infinity}]
Given your choice of p this becomes
G[t] = Integrate[F[t-x], {x,-Infinity,a}]
where a=30. Substitute z=t-x, you get
G[t] = Integrate[F[z], {z,t-a,Infinity}]
The first derivative of this is
G'[t] = -F[t-a]
So you can find F simply from
F[t] = -G'[t+a] = -G'[t+30]

Roman.

On Apr 7, 10:00 am, "dantimatter" <dantimat... at gmail.com> wrote:
> hello all,
>
> i have a function G(t) which is a convolution of F(t) and p(t):
> G=F*p .
> now p(t) is a step function (p(t)=UnitStep[30-t]) and i know what G(t)
> is, but i'd like to get at F.  what i've been trying to do is take the
> laplace transform of G and dividing by the laplace transform of p,
> then inverting that to get F. i'm having a lot of trouble.   the built
> in inverseLaplaceTransform function is unable to do it, and i haven't
> found any numerical routines that are reliable.  is there some reason
> why i shouldn't be able to do this?  perhaps there are better ways to