Re: Complex integration
- To: mathgroup at smc.vnet.net
- Subject: [mg55816] Re: Complex integration
- From: dh <dh at metrohm.ch>
- Date: Thu, 7 Apr 2005 05:09:55 -0400 (EDT)
- References: <d305d3$rv7$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
Hi Berhard
considering your two integrands:
1) 2^(-(-12500 + w)^2/31250 - (-12500 + w + W)^2/
31250)
and
2) 2^(-(-12500 + w)^2/31250 - (-12500 + w + W)^2/
31250)*E^((-I)*Pi*Cos[(Pi*w)/64] +
I*Pi*Cos[(Pi*(w + W))/64])
it is not surprising that the second integral is much more difficult.
HOWEVER, if I try to do the second integral, Mathematica simply dies!!!
Here my versio of Mathematica:
$Version
"5.1 for Microsoft Windows (October 25, 2004)"
If somebody else could very this we will have a case for WRI
Sincerely, Daniel
Bernhard von Vacano wrote:
> Dear Group,
>
> I am experiencing problems with a Mathematica simulation I am trying to
> perform.
>
> The first part of my notebook evaluates without any problems:
>
>
> ETL[w_] := Exp[-(w - w0)^2/width^2*2 Log[2]]
>
> EPer[w_] := ETL[w]* Exp[I Pi Cos[2 Pi w/Wm]]
> EPerCC[w_] := ETL[w]* Exp[-I Pi Cos[2 Pi w/Wm]]
>
> parameter := {w0 -> 12500, width -> 250, Wm -> 128}
>
> AmpTL[W_] = Integrate[ETL[w]ETL[w + W] /. parameter, {w, 10000, 20000}]
>
>
> This is the first expression I need for plotting (in the Range of W =
> 0..1000), and gives the correct results.
> With my so far limited mathematica skills, I could not get the following
> modified expression to be evaluated:
>
> AmpPer[W_] = Integrate[EPerCC[w]EPer[w + W] /. parameter, {w, 10000, 20000}]
>
> Like this, it only produces "$RecursionLimit::reclim: "Recursion depth of
> 256 exceeded." und "$IterationLimit::itlim: Iteration limit of 4096
> exceeded."
>
> I have tried other limits for the integration, indefinite integration and
> NIntegrate, but so far it did not work. (At the moment, I am using
> Mathematica 4).
>
> Could anybody help me with this? I don't really know which parameters to use
> or how to restate the problem to get it working. It doesn't seem too hard to
> compute, does it?
>
>
> Thank you very much in advance,
>
> Yours
> Bernhard
>
>