# HP Forums

Full Version: New variant for the Romberg Integration Method
You're currently viewing a stripped down version of our content. View the full version with proper formatting.

Hi All,

I just posted, on my web site, the following article for a new variant for the Romberg method. The article actually looks at several variants and selects the best one.

Enjoy!!

Namir

Edited: 18 Apr 2012, 12:49 a.m.

Thanks. This should be intriguing

Using both the HP-15c and Wolfram Alpha I get values that differ from yours for these examples:

ln(x)/x integrate 1 to 100 = 10.60378

sin(x) integrate 1e-10 to pi/4 = 0.2928932

Nick

Edited: 18 Apr 2012, 5:27 a.m.

Nick,

Thanks for the corrections. I n the case of the sin(x), I meant sin(x)/x. I posted the article with the corrected results.

Namir

Edited: 18 Apr 2012, 7:29 a.m. after one or more responses were posted

Quote:
Nick,

Thanks for the corrections. I n the case of the sin(x), I meant sin(x)/x. I should posted a corrected article very shortly.

Namir

The standard name for the sin(x)/x function is sinc(x), an abbreviation of sinus cardinalis (i.e.: cardinal sine).

Regards from V.

Right you are! And I learned a new function name. Alpha Worlfram recognized the sinc(x) funcion!!

:=)

So does the 34S :-)

- Pauli

This brings back memories as the sinc function was one of the first things I plotted as a teenager on my newly acquired Sinclair ZX81 computer.

Nick

Quote:
Right you are! And I learned a new function name. Alpha Worlfram recognized the sinc(x) funcion!!

:=)

I'm glad you did, I also learn new things each and every day.

About the sinc(x) function, it has many interesting properties and quirks but the one
I find most uncanny is this: a little computation or theoretical work will quickly stablish the following results:

• I1 = Integral( 0, Infinity, sinc(x) dx) = Pi/2

• I2 = Integral( 0, Infinity, sinc(x)*sinc(x/3) dx) = Pi/2

• I3 = Integral( 0, Infinity, sinc(x)*sinc(x/3)*sinc(x/5) dx) = Pi/2

• I4 = Integral( 0, Infinity, sinc(x)*sinc(x/3)*sinc(x/5)*sinc(x/7) dx) = Pi/2

• I5 = Integral( 0, Infinity, sinc(x)*sinc(x/3)*sinc(x/5)*sinc(x/7)*sinc(x/9) = Pi/2

• I6 = Integral( 0, Infinity, sinc(x)*sinc(x/3)*sinc(x/5)*sinc(x/7)*sinc(x/9)*sinc(x/11) = Pi/2

• I7 = Integral( 0, Infinity, sinc(x)*sinc(x/3)*sinc(x/5)*sinc(x/7)*sinc(x/9)*sinc(x/11)*sinc(x/13) = Pi/2

but lo and behold, we unexpectedly find that

• I8 = Integral( 0, Infinity, sinc(x)*sinc(x/3)*sinc(x/5)*sinc(x/7)*sinc(x/9)*sinc(x/11)*sinc(x/13)*sinc(x/15) = Pi/2.0000000000294+ !!

You might want to check this amazing fact by trying and computing said integrals I1, I2, ..., I8 using the 34S' extreme precision capabilities, it would be a fine test for any numerical integration procedure such as yours ! ... XD

Best regards from V.