Numerical Integration « Next Oldest | Next Newest »

 ▼ Namir Posting Freak Posts: 2,247 Threads: 200 Joined: Jun 2005 11-21-2011, 02:38 PM Hi All, I posted an article (as a pdf file) on my web site. You can download it if you click here. The article is a result of some basic investigation to enhance the Trapezoidal method by using additional points inside each integration interval. My new approach produces some improvements over the basic Trapezoidal rule using weights with points inside the interval. Interestingly, the new method does very well when these weights match those of Simpson's rule, Simpson's 3/8 rule, and other similar methods. The article explains how you can mix the Trapezoidal method with better methods when the scheme of the latter are applied to each integration interval in an isolated or compartmentalized manner. The article shows the source code for VBA functions that implement the various algorithms covered by the article. Namir ▼ Allen Senior Member Posts: 562 Threads: 27 Joined: Oct 2006 11-22-2011, 07:56 AM Thank you for the interesting article! You are the quintessential quadrature expert. ▼ Namir Posting Freak Posts: 2,247 Threads: 200 Joined: Jun 2005 11-22-2011, 08:07 AM Glad you enjoyed the article. I was just asking some basic questions about numerical integration. Namir Namir Posting Freak Posts: 2,247 Threads: 200 Joined: Jun 2005 11-22-2011, 09:27 AM I updated the article to correct a mistake and clarify a few things. Namir Mike Morrow Posting Freak Posts: 758 Threads: 9 Joined: Jul 2007 11-22-2011, 07:37 PM How does it compare to Gauss-Kronrod quadrature? The implementation of G-K quadrature on the fx-115/991ES seems to give much better and faster results than the quadrature methods of the various HP calculators, including the HP 50g and my beloved HP 42S. ▼ Namir Posting Freak Posts: 2,247 Threads: 200 Joined: Jun 2005 11-22-2011, 10:59 PM The Gauss-Kronrod quadrature is a better method. Allen Senior Member Posts: 562 Threads: 27 Joined: Oct 2006 11-22-2011, 09:57 PM I wonder if you have considered using something other than the monomial basis for your objective function for your numeric integration?? For example, whether using uniform weighting or individual weights, you can use the method of undetermined coefficients to solve the abscissa locations. For example to find a three point, equal-weighted solution on the following range: With 4 unknowns, you can just use this equation: to find +/-inv(sqrt(2)), and 0 are the nodes, with the weight w=2/3. Any number of weights or points can be calculated this way, but this only uses the monomial basis. With a different (perhaps higher-order?) basis function, Could it be more accurate?? Edited: 22 Nov 2011, 10:07 p.m. ▼ Namir Posting Freak Posts: 2,247 Threads: 200 Joined: Jun 2005 11-23-2011, 09:14 AM I really have not seen any numerical integration algorithm with a single weight. Being a pragmatic person, if you can find a method with a constant weight and it works well, that is fine. Namir ▼ Allen Senior Member Posts: 562 Threads: 27 Joined: Oct 2006 11-23-2011, 08:37 PM I was just reading about this method last week. The key is using the monomial basis to actually move the nodes around to compensate for the equal weight. I first heard about Chebyshev Quadrature in the book I am reading. It's simpler numerically than Simpson and a few other methods because you save n-multiplications by factoring out the weighting function from the summation. It's interesting to solve, too. I hope you all have a Merry Thanksgiving!!

 Possibly Related Threads... Thread Author Replies Views Last Post Integration question and "RPN" mode comment Craig Thomas 16 2,863 12-05-2013, 02:32 AM Last Post: Nick_S HP Prime numerical restrictions? Alasdair McAndrew 4 806 11-16-2013, 05:32 PM Last Post: Alasdair McAndrew HP Prime numerical precision in CAS and HOME Javier Goizueta 5 1,032 11-16-2013, 03:51 AM Last Post: Paul Dale WP34s integration trapped in infinite loop Bernd Grubert 25 3,007 10-17-2013, 08:50 AM Last Post: Dieter HP Prime integration Richard Berler 1 605 10-06-2013, 10:52 PM Last Post: Helge Gabert [HP-Prime] AMBIGUITY between Numerical Calculation (HOME) and Numerical/Symbolic Calculation (CAS mode) CompSystems 2 671 08-18-2013, 07:06 PM Last Post: CompSystems OT: My brain is failing me again. Help with numerical / mechanical problem required. Harald 4 866 07-01-2013, 10:31 AM Last Post: Harald integration on 39gII emulator Wes Loewer 29 3,297 06-07-2013, 05:58 PM Last Post: Chris Smith WP-34S Integration Richard Berler 15 1,839 03-08-2013, 02:29 AM Last Post: Walter B HP 34S integration Richard Berler 16 1,788 02-18-2013, 04:42 PM Last Post: Marcus von Cube, Germany

Forum Jump: