[Date Index]
[Thread Index]
[Author Index]
Re: need help with integration
*To*: mathgroup at smc.vnet.net
*Subject*: [mg52019] Re: need help with integration
*From*: "David W. Cantrell" <DWCantrell at sigmaxi.org>
*Date*: Sun, 7 Nov 2004 01:03:50 -0500 (EST)
*References*: <cmhv6l$pqv$1@smc.vnet.net>
*Sender*: owner-wri-mathgroup at wolfram.com
"symbio" <symbio at sha.com> wrote:
> Below I define two functions x[t] and h[t], then in eq1 I integrate
> the integrand ( x[tau] h[t-tau] ) from -inf to +inf and get one set of
> results and plots, then I integrate the same integrand as before but this
> time in two steps, once from -inf to 0 and once from 0 to +inf, but the
> results and plots from the integration performed in two steps are NOT the
> same as the results and plots from integration in one step. Is this a
> bug and if so what's the work around?
This looks essentially like what you asked here about two weeks ago in
"Integration of UnitStep has bugs!? help!"
<http://groups.google.com/groups?threadm=clneor%24obq%241%40smc.vnet.net>
I would have answered your question in the previous thread had you not
already gotten two good answers. Look at them carefully. I suspect that the
trouble you're having now is the same as you were having previously; if so,
then your question has already been answered, twice.
David
> Try this in the input cell
>
> In[19]:=
> x[t_] = (3Exp[-.5t]UnitStep[t]) + DiracDelta[t + 3]
> h[t_] = UnitStep[t] - UnitStep[t - 2]
> eq1 = Integrate[x[\[Tau]]h[t - \[Tau]], {\[Tau], -\[Infinity],
> \[Infinity]}] Plot[eq1, {t, -10, 10}, PlotRange -> All]
> eq2 = Integrate[x[\[Tau]]h[t - \[Tau]], {\[Tau], 0, \[Infinity]}]
> Plot[eq2, {t, -10, 10}, PlotRange -> All]
> eq3 = Integrate[x[\[Tau]]h[t - \[Tau]], {\[Tau], -\[Infinity], 0}]
> Plot[eq3, {t, -10, 0}, PlotRange -> All]
>
> You will get the following output expressions (I couldn't paste the plots
> here, but if you run the above input cells you should get the plots too,
> then it will be more obvious what the problem is)
>
> Out[21]=
> \!\(1\/2\ \((1 + \(6.`\ \@\((\(-2\) + t)\)\^2\)\/\(\(\(2.`\)\(\
> \[InvisibleSpace]\)\) - 1.`\ t\) + \(6.`\ \@t\^2\)\/t + \@\((3 +
> t)\)\^2\/\(3 \
> + t\) - \(1 + t + \@\((1 + t)\)\^2\)\/\(1 + t\))\)\)
>
> Out[23]=
> \!\(\((\(-6.`\) +
> 16.30969097075427`\ \[ExponentialE]\^\(\(-0.5`\)\ t\))\)\ UnitStep[\
> \(-2\) + t] + \((\(\(6.`\)\(\[InvisibleSpace]\)\) -
> 6.`\ \[ExponentialE]\^\(\(-0.5`\)\ t\))\)\ UnitStep[t]\)
>
> Out[25]=
> \!\(1\/2\ \((\(-\(\(1 +
> t\)\/\@\((1 + t)\)\^2\)\) + \(3 + t\)\/\@\((3 + t)\)\^2)\)\)
Prev by Date:
**Re: MathGroup /: Descriptive headings**
Next by Date:
**Re: need help with integration**
Previous by thread:
**Re: need help with integration**
Next by thread:
**Re: need help with integration**
| |