Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
2004
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*November
*December
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 2004

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

Search the Archive

Re: Strange PV results from Integrate

  • To: mathgroup at smc.vnet.net
  • Subject: [mg51288] Re: Strange PV results from Integrate
  • From: carlos at colorado.edu (Carlos Felippa)
  • Date: Tue, 12 Oct 2004 01:57:54 -0400 (EDT)
  • References: <ckd61e$522$1@smc.vnet.net>
  • Sender: owner-wri-mathgroup at wolfram.com

carlos at colorado.edu (Carlos Felippa) wrote in message news:<ckd61e$522$1 at smc.vnet.net>...
> Could somebody pls check if these results persist in 5.0? Thx.
> 
> $Version
> "4.2 for Power Macintosh (August 27, 2002)"
> 
> Integrate[Log[  2+2*Cos[x]  ], {x,-Pi,Pi}] 
> 0   (* correct *)
> 
> Integrate[Log[ 2*(1+Cos[x]) ], {x,-Pi,Pi}] 
> identical integrand returns unevaluated
> 
> Integrate[Log[ -2-2*Cos[x]  ], {x,-Pi,Pi}]//InputForm; 
> (2*I)*Pi^2   (* correct *)
> 
> Integrate[Log[ -2*(1+Cos[x]) ], {x,-Pi,Pi}] 
> identical integrand returns unevaluated
> 
> Integrate[Log[ 1+ Cos[x]  ], {x,-Pi,Pi}]
> (I/2)*Pi^2 - 2*Pi*Log[1 + I] - Pi*Log[2]    (* wrong *)
> 
> FullSimplify[Integrate[Log[ 1+ Cos[x]  ], {x,-Pi,Pi}]]
> -(Pi*Log[4])     (* wrong *)
> 
> These integrals arise on applying the first Szego theorem 
> to some benchmark infinite Toeplitz matrices.  Results 
> labelled "wrong" contradict the theorem.

An update: a student with 5.0.1 did the tests and the unevaluated
integral cases return the values 0 and 2*I*Pi^2.

The last value is fine and agrees with Szego theorem.  I had the wrong
sign in the real part of his Toeplitz asymptotic expansion.  The PV
for Log[ a*(2+2*Cos[x])] in [-Pi,Pi] should be 2*Pi*Log[a], 
so the - sign is correct if a=1/2.


  • Prev by Date: Re: cross-product in cylindrical problem
  • Next by Date: Re: Outer product in mathematica
  • Previous by thread: Strange PV results from Integrate
  • Next by thread: Re: Strange PV results from Integrate