Slow integral evaluation on HP50g

04222014, 01:56 AM
Post: #1




Slow integral evaluation on HP50g
When trying to evaluate the definite integral
int(0,6,abs(2*sin(e^(x/4))+1),x) on the HP50. The integration never seems to complete (Regardless of whether I am in approximate or exact mode). Can anyone explain the slow evaluation? Appreciatively, Chris 

04222014, 03:50 AM
(This post was last modified: 04222014 03:55 AM by Manolo Sobrino.)
Post: #2




RE: Slow integral evaluation on HP50g
In EMU48 via EQW it works fine and gets this answer (with 10 accurate digits, pay attention to the graph and check if your Number Format is Standard, because then your goal is 12), but only after a little while and using emulation speed. I guess it's Romberg integration, which suffers if the derivatives of the integrand are not continuous, you see that with a small tolerance. Use Approx mode, and change effectively the tolerance via Number Format > FIX 5, for instance.
When you have an abs function, look if there are points in the interval with a value that could permit two slopes (for a global abs any root will be ). There's the trouble: at some point it won't converge for a subinterval there. 

04222014, 03:56 AM
Post: #3




RE: Slow integral evaluation on HP50g
My 50g comes up with an answer of about 12.5733 in STD display after the maximum 64K1 function evaluations. It comes up with 12.5726 in 4 FIX after only 63 evaluations.
I suspect the noncontinuous derivatives at about x=5.2 has a lot to do with it. Removing the ABS, it comes up with a result after 255 evaluations in STD. 

« Next Oldest  Next Newest »

User(s) browsing this thread: 1 Guest(s)