[Date Index]
[Thread Index]
[Author Index]
Re: Calculate a numerical integral with enough precision
*To*: mathgroup at smc.vnet.net
*Subject*: [mg114736] Re: Calculate a numerical integral with enough precision
*From*: Bob Hanlon <hanlonr at cox.net>
*Date*: Tue, 14 Dec 2010 06:57:56 -0500 (EST)
You want WorkingPrecision to be at least as large as PrecisionGoal. However, the default is to use a WorkingPrecision that is twice the PrecisionGoal. For this particular integral, the oscillations near zero cause problems eventhough the automatic method used appears to be OscillatorySelection (identical results).
(tab = Table[{n, NIntegrate[Exp[Sin[1/x]], {x, 0, \[Pi]/2},
WorkingPrecision -> n, PrecisionGoal -> 15,
Exclusions -> (x == 0)] // Quiet}, {n, 15, 50, 5}]) // Grid
Bob Hanlon
---- Alexei Boulbitch <alexei.boulbitch at iee.lu> wrote:
=============
NIntegrate[Exp[Sin[1/x]], {x, 0, \[Pi]/2}, WorkingPrecision -> 12,
Exclusions -> (x == 0), PrecisionGoal -> 15] // Quiet
3.09892833696
* /Subject/: [mg114668] Calculate a numerical integral with enough
precision
* /From/: alphatest <iliurarfwpuap at mailinator.com>
* /Date/: Sun, 12 Dec 2010 05:45:17 -0500 (EST)
------------------------------------------------------------------------
How can we calculate the following integral up to 10-12 decimal places?
integrate exp(sin(1/x)) , x=0..pi/2
It's as if no integration method or precision option lets you calculate more than 5-6 decimal places. Is it possible?
--
Alexei Boulbitch, Dr. habil.
Senior Scientist
Material Development
IEE S.A.
ZAE Weiergewan
11, rue Edmond Reuter
L-5326 CONTERN
Luxembourg
Tel: +352 2454 2566
Fax: +352 2454 3566
Mobile: +49 (0) 151 52 40 66 44
e-mail: alexei.boulbitch at iee.lu
www.iee.lu
Prev by Date:
**Re: LessEqual vs Inequality, was ..Re: Replacement Rule with Sqrt in denominator**
Next by Date:
**Re: Integral no longer evaluated in Version 7, 8**
Previous by thread:
**Re: Calculate a numerical integral with enough precision**
Next by thread:
**Re: Calculate a numerical integral with enough precision**
| |