New variant for the Romberg Integration Method



#10

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.


#11

Thanks. This should be intriguing

#12

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

x in radians

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

Nick

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


#13

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


#14

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.


#15

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

:=)


#16

So does the 34S :-)


- Pauli

#17

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.

#18

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


Possibly Related Threads...
Thread Author Replies Views Last Post
  Integration question and "RPN" mode comment Craig Thomas 16 2,022 12-05-2013, 02:32 AM
Last Post: Nick_S
  A fast Bernoulli Number method for the HP Prime Namir 16 1,652 11-22-2013, 04:46 PM
Last Post: Namir
  WP34s integration trapped in infinite loop Bernd Grubert 25 2,130 10-17-2013, 08:50 AM
Last Post: Dieter
  HP Prime integration Richard Berler 1 401 10-06-2013, 10:52 PM
Last Post: Helge Gabert
  A brand new calculator benchmark: "middle square method seed test" Pier Aiello 25 2,194 09-13-2013, 01:58 PM
Last Post: Pier Aiello
  Downhill Simplex Method (Nelder Mead) for WP-34S Marcel Samek 0 272 07-07-2013, 03:05 AM
Last Post: Marcel Samek
  integration on 39gII emulator Wes Loewer 29 2,381 06-07-2013, 05:58 PM
Last Post: Chris Smith
  a new variant of the 200LX? robert rozee 1 412 06-03-2013, 10:33 AM
Last Post: BShoring
  WP-34S Integration Richard Berler 15 1,356 03-08-2013, 02:29 AM
Last Post: Walter B
  HP 34S integration Richard Berler 16 1,263 02-18-2013, 04:42 PM
Last Post: Marcus von Cube, Germany

Forum Jump: