MathGroup Archive 1999

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

Search the Archive

ExpIntegralEi

  • To: mathgroup at smc.vnet.net
  • Subject: [mg18463] ExpIntegralEi
  • From: Lionel ARNAUD <arnaud at lmt.ens-cachan.fr>
  • Date: Wed, 7 Jul 1999 00:11:42 -0400
  • Organization: ENS Cachan, France
  • Sender: owner-wri-mathgroup at wolfram.com

Hello,

	I am from LMT-Cachan FRANCE, working with MATHEMATICA V.3, I made this
calculation:

c2 = -0.05018627683354541 - 0.153047656745338 I;
c3 = -0.7828709924214918 + 0.2780791279205129 I;
c5 = -0.6758555487562639 - 0.04753624179417532 I;

Integrate[Exp[beta*c2+s*(c3+beta*c5)],	{s,0,1},{beta,0,1}]

NIntegrate[Exp[beta*c2+s*(c3+beta*c5)],	{s,0,1},{beta,0,1}]

The results given are:
-0.228103 + 10.5644 I
 0.587252 +  0.0191685 I

Not the same !

Options of accuracy, algorithm,... don't change much the result, even if
you change c2, c3 or c5 a little. In fact the Exp[beta...] is very regular
and you can plot it and observe that the good result is given by the Nintegral.
It seems that the formal integration given by Mathematica 
is not totally correct.

If you can tell me more, I am interested..


  • Prev by Date: Absence of Compile warning in Mathematica 4
  • Next by Date: Kronecker product of matrices
  • Previous by thread: Absence of Compile warning in Mathematica 4
  • Next by thread: RE: ExpIntegralEi