Simpson's Rule on HP17b/19b/200LX anyone? Peter A. Gebhardt Unregistered Posts: 174 Threads: 20 Joined: Sep 2006 07-09-2007, 11:07 AM Does anyone have available a routine which allows to integrate functions with Simpson's Rule (or even Romberg's algorithm ...) on one of the Solver-equipped calculator series (and it's newer siblings) above? Hints for implemention or code snippets would be welcome! Reason: To enter the function in question only once - as demonstrated in Ed Keefe's articles on Palmptop Papers, where he used an implementation of the Trapezoid Rule. TIA. Best regards Peter A. Gebhardt couldn't be reached - temporarily shutdown or gone forever ??? Any mirrors / backups available? Edited: 9 July 2007, 11:19 a.m. bill platt Unregistered Posts: 2,448 Threads: 90 Joined: Jul 2005 07-09-2007, 12:19 PM 1424241 Edited: 9 July 2007, 12:19 p.m. Bill (Smithville, NJ) Unregistered Posts: 614 Threads: 66 Joined: Jul 2006 07-09-2007, 01:28 PM Peter, If you have the HPMuseum DVD or CD's, then check out pages 36 thru 41 of the "HP-27S and HP-19B Technical Applications Step-by-Step Solutions" book. It's a Solver example of Numerical Integration using Simpson's Rule. This should also be usable for the HP-200LX and the HP-17BII. Bill Peter A. Gebhardt Unregistered Posts: 174 Threads: 20 Joined: Sep 2006 07-09-2007, 05:51 PM Thx. a lot Bill! I've already ordered the DVD from Dave today. What bothers me conc. the implementation on the Solver is, how to get rid of the nuissance to enter the function several times (as demonstrated by y(0),y(i) and y(n) in the following code snippet: Quote:``` h = (xn - x0) / n Integral = y(0) For i = 1 To n - 1 If i Mod 2 = 0 Then Integral = Integral + 2 * y(i) Else Integral = Integral + 4 * y(i) End If Next Integral = Integral + y(n) ``` An additional requirement: Either achieving minimal error for a given number of function evaluations and?/or computing an approximation having a given error bound, at minimal cost. Best regards Peter A. Gebhardt PS: @Bill Platt, thank you too, for reminding me to express the problem in a way better to understand ;-) Edited: 9 July 2007, 5:58 p.m. Bill (Smithville, NJ) Unregistered Posts: 614 Threads: 66 Joined: Jul 2006 07-09-2007, 07:44 PM Hi Peter, Quote: I've already ordered the DVD from Dave today. Great - be sure to set aside some time (a lot of time) for perusing them. They are a great resource. Quote: how to get rid of the nuissance to enter the function several times The secret is to use the SUM function to loop. When you get the DVD and look at the Technical Applications Book, it'll all become clear (or, in my case, a little less murky). Using the Sum let's you enter the function only one time. The Technical Applications book is a Must Read for anyone wanting to make good use of the Solver. It has a great section on the Let and Get functions that I haven't found in any other reference book. Bill Peter A. Gebhardt Unregistered Posts: 174 Threads: 20 Joined: Sep 2006 07-10-2007, 07:27 PM Current status - I have a working solution, BUT: The implementation of an algorithm (as an example this one here) does work reliably only with N == powers of 2. Using other multiples of 2 shows an oscillating behaviour of the results dampening out with increasing N. (clearly visible thanks to the PLOT function of the HPCALC Solver!) Excuse me for asking for help still - I suppose that this is an effect of the used algorithm (the integration) as such - are there any enhancements, I should/can add? It's about 40 years since I was used to calculus, so pls. could you point me to some information where I can find useful info targetted to the 200LX Solver capabilities. For example, solutions using arrays for intermediate storage are a no-no, because of the read-only characteristics of the .STA lists. TIA Peter A. Gebhardt Edited: 10 July 2007, 7:35 p.m. Bill (Smithville, NJ) Unregistered Posts: 614 Threads: 66 Joined: Jul 2006 07-10-2007, 08:39 PM Peter, You might want to try MERCURY from the HP200LX SUPER Site. It's supposed to be a replacement for Borlands Eureka's program and can be run on the HP-200LX. It does Integration as well as many other equations. You can find it here: Bill Peter A. Gebhardt Unregistered Posts: 174 Threads: 20 Joined: Sep 2006 07-10-2007, 09:00 PM Bill, Thx. again for your help offered. Because of the advisory process I use, I can't leave the Solver (data has to be used with other equations later on), so your advice is very welcome, but not what I'm looking for. As an example, part of my work requires solving an Exponential Reciprocal Gamma Function, where part of the function is the Integral over x^(alpha-1)*exp(-x/beta) dx (for what I'm wanting to implement the algorithm). So I'm still looking for a solution which I can "integrate" (no pun intended!) into the Solver environment. Best regards Peter A. Gebhardt Bill (Smithville, NJ) Unregistered Posts: 614 Threads: 66 Joined: Jul 2006 07-11-2007, 01:02 PM Peter, I went through my set of The HP Palmtop Paper On Disk and found the following Solver Equation: ```Integral !CALC Num'cal Integration from TECH APP booklet E.G. Set A=-1,B=1,X=1 N=20 and find the area of a circle(I)! I=L( H:(A-B)/(-2*N))/3* SIGMA(R, 0, 2*N, 1, 0*L(X,A+R*G(H))+(SGN(R)+SGN(2*N-R)+2*MOD(R,2))* ( ! function goes on next line(s) ! ! e.g. area under a line (right triangle)! ! X ! (2*SQRT(1-X^2)) !area in top and bottom of circle! ) ) ``` This is based on the one in the Tech App Guide. You can just replace (2*SQRT(1-X^2)) with your F(X). Note that N=20 means 40 iterations. Have Fun. Bill Edited: 11 July 2007, 1:03 p.m. Peter A. Gebhardt Unregistered Posts: 174 Threads: 20 Joined: Sep 2006 07-11-2007, 01:06 PM Bill, Thx. a lot! Peter A. Gebhardt « Next Oldest | Next Newest »

 Possibly Related Threads… Thread Author Replies Views Last Post OT: Jeppesen E6B Wind-Easy Computer (Slide Rule) Eddie W. Shore 18 6,576 10-12-2013, 03:26 PM Last Post: George Litauszky slide rule magnifiers Al 15 4,524 09-17-2013, 04:49 PM Last Post: aurelio trig scales on the Post Versalog slide rule Al 12 3,625 09-15-2013, 06:01 AM Last Post: John I. First experience with a slide rule aurelio 27 7,004 08-17-2013, 01:16 PM Last Post: Garth Wilson 200LX problem: Insert batteries and get a continuous beep... ideas ? Gene Wright 8 2,686 06-19-2013, 04:52 PM Last Post: gene wright A slide rule to conquer the moon revisited Palmer O. Hanson, Jr. 12 3,525 06-19-2013, 03:05 PM Last Post: Ken Shaw A slide rule to conquer the Moon Gilles Carpentier 11 3,178 06-03-2013, 02:36 PM Last Post: Ingo a new variant of the 200LX? robert rozee 1 1,086 06-03-2013, 10:33 AM Last Post: BShoring OT: My Slide Rule Collection Eddie W. Shore 18 5,517 04-25-2013, 11:47 PM Last Post: Garth Wilson More Slide Rule Stuff... . .... Thomas Chrapkiewicz 1 1,095 09-25-2012, 03:41 PM Last Post: Kiyoshi Akima

Forum Jump: