help in Numerical Intergation
- To: mathgroup at smc.vnet.net
- Subject: [mg36518] help in Numerical Intergation
- From: "Salman Durrani" <dsalman at itee.uq.edu.au>
- Date: Tue, 10 Sep 2002 06:24:52 -0400 (EDT)
- Organization: The University of Queensland, Australia
- Sender: owner-wri-mathgroup at wolfram.com
Hi
I am having some trouble with numerical integration.
\!\(\(\(Timing[
ansOUT = {}; \[IndentingNewLine]For[EbNo = 0, EbNo <= 2,
EbNo += 10\/10, \[IndentingNewLine]linEbNo =
10\^\(EbNo\/10\); \[IndentingNewLine]M = 64; \[IndentingNewLine]K =
Log[2, M]; \[IndentingNewLine]p = \(M/2\)\/\(M - 1\)\ Sum[
Binomial[M - 1,
n] \((\(-1\))\)\^\(n + 1\)\/\(1 + n\)\ Exp[\(-\
linEbNo\)\ K\ n\/\(1 + n\)], {n, 1,
M - 1}]; \[IndentingNewLine]Pm = \(1\/\@\(2\ \[Pi]\)\)
NIntegrate[
1 - \(\((\(1\/\@\(2\ \[Pi]\)\)
NIntegrate[
Exp[\(-x\^2\)\/2], {x, \(-\[Infinity]\), y},
WorkingPrecision -> 32\ ])\)\^\(M - 1\)\)
Exp[\(\(-1\)\/2\) \((y - \@\(2\ K\ linEbNo\))\)\^2], {y,
\(-\
\[Infinity]\), \[Infinity]},
WorkingPrecision ->
32]; \[IndentingNewLine]Pb = \(\(M/2\)\/\(M - 1\)\)
Pm; \[IndentingNewLine]AppendTo[
ansOUT, {N[a, 4], N[p, 8],
N[Pb, 8]}];\[IndentingNewLine]];]\)\(\[IndentingNewLine]\)\)\)
TableForm[ansOUT]
The code can also be downloaded from
http://www.itee.uq.edu.au/~dsalman/_downloads/Proakis.nb
Col. 2(symblic formula) and 3 (result of numerical integration) of the
variable ansOUT should be identical.
However the numerical integration is giving a 0 answer.
Can anyone please suggest how a solution to this problem.
Thanks
Salman