+- HP Forums (https://www.hpmuseum.org/forum)
+-- Forum: HP Calculators (and very old HP Computers) (/forum-3.html)
+--- Forum: General Forum (/forum-4.html)
+--- Thread: Fractional Part - A Difficult Integral (/thread-1739.html)
Fractional Part - A Difficult Integral - Gerald H - 06-29-2014 11:36 AM
Here some results for integrating FP from zero to 5.8:
Maple V Release 3.0
evalf(int(frac(x), x=0..5.8)) 2.820001887
Maxima 5.23.2
quad_qags(x-floor(x), x, 0, 5.8) [2.720000001263673,1.5265947617137954*10^-
8,3171,0]
PARI 2.4.2
intnum(x=0,5.8,frac(x)) 2.770657207613055496607250297
WolframAlpha 2011-5-28
integrate(frac(x),x,0,5.8) 2.82
CASIO fx-9860GII (& fx-5800P)
∫_0^5.8▒〖Frac X〗 dx 2.72
HP 42S
Gauss-Lagrange 16 point:
Divisions: 1 2.32
2 2.71441362769
3 2.95110930932
4 2.84122268726
5 2.72066778956
6 2.8405122545
HP 50G
Inbuilt integration programme: ∫_0^5.8▒〖FP(X)〗 dX
FIX 6 2.820116 in 1,232.15 seconds
Gauss-Lobatto 4 point formula with 7 and 13 point Kronrod extensions:
FIX 2 2.80
3 2.819
4 2.8199
5 2.82000
6 2.820000 in 12.34 seconds
SHARP EL-9650
∫_0^5.8▒〖fpart X〗 dxX 2.419454779
TI-84 Plus 2.53MP
fnInt(fPart(X),X,0,5.8) 2.720014333
TI-86 (tol = 1E-5)
fnInt(fPart x,x,0,5.8) 2.71992436738
TI-89 (& voyage 200 & 92(plus))
∫_0^5.8▒〖fPart (x)〗 dx 2.72000000036
RE: Fractional Part - A Difficult Integral - Massimo Gnerucci - 06-29-2014 11:55 AM
(06-29-2014 11:36 AM)Gerald H Wrote: Here some results for integrating FP from zero to 5.8:
Another example (apart from WolframAlpha) where a pencil and paper approach is better.
RE: Fractional Part - A Difficult Integral - pito - 06-29-2014 05:18 PM
(06-29-2014 11:36 AM)Gerald H Wrote: Here some results for integrating FP from zero to 5.8:WP34s:
2.819612562 after several minutes
RE: Fractional Part - A Difficult Integral - kakima - 06-29-2014 06:17 PM
Change the upper limit to 6.4 and try evaluating with the built-in integrator on any HP calculator that has one.
RE: Fractional Part - A Difficult Integral - Thomas Klemm - 06-29-2014 06:47 PM
(06-29-2014 06:17 PM)kakima Wrote: Change the upper limit to 6.4 and try evaluating with the built-in integrator on any HP calculator that has one.
Link for the lazy: Numerical Integration on the 35S.
This is another example where most HP-calculators fail.
Cheers
Thomas
RE: Fractional Part - A Difficult Integral - Gerald H - 06-29-2014 07:57 PM
(06-29-2014 06:17 PM)kakima Wrote: Change the upper limit to 6.4 and try evaluating with the built-in integrator on any HP calculator that has one.
HP 50G Gauss-Lobatto 4 point formula with 7 and 13 point Kronrod extensions:
3.0489 in 20.8 sec
RE: Fractional Part - A Difficult Integral - pito - 06-29-2014 09:40 PM
(06-29-2014 06:47 PM)Thomas Klemm Wrote:(06-29-2014 06:17 PM)kakima Wrote: Change the upper limit to 6.4 and try evaluating with the built-in integrator on any HP calculator that has one.
Link for the lazy: Numerical Integration on the 35S.
This is another example where most HP-calculators fail.
Cheers
Thomas
What about to use a simple Monte Carlo method as a "pre-processor" to get a "rough estimate", ie.:
Code:
NIntegrate[x - Floor[x], {x, 0, 6.4},
Method -> {"MonteCarloRule", "Points" -> 1000}]returns always nice results around 3.08. So we can see whether the specific "precise" integration method with some "difficult" integrands has failed or not..
PS: another example from above posts:
Code:
NIntegrate[
x*x*(x*x - 47*47)*(x*x - 88*88)*(x*x - 117*117), {x, -128, 128},
Method -> {"MonteCarloRule", "Points" -> 1000}]RE: Fractional Part - A Difficult Integral - ttw - 06-30-2014 01:17 PM
Numerical integration of things like Frac(x)dx do not gain by using high-order rules (Simpson's, Gauss, etc.) because the first derivative is not continuous.
Monte Carlo (and quasi Monte Carlo) does no worse than its usual performance (which isn't that good anyway) as the variance (or variation) is small.
One QMC example is to use the points (2j-1)/2N for j=1,N as an N-point integration formula. Another is to use Frac(N*Sqrt(2)) as a sequence to integrate these types of functions.
The last sequence is easy to compute; just add Frac(Sqrt(2)) at each step and reduce below 1.