[VA] SRC#001 - Spiky Integral

[Thomas Klemm]

(07-15-2018 09:22 PM)John Keith Wrote:  Or equivalently, using only built-in commands:

<< 0 m NDUPN ->LIST >>

That's what I was originally looking for but it appears to be missing on the HP-48.
Thanks for all your suggestions.

(07-16-2018 12:31 AM)Gerson W. Barbosa Wrote:  We now can use it to compute the integrals much faster than I was able to before.

If you want to create a list of all values you better do that within the FOR-loop:

Code:

«
  « 0 SWAP NDUPN →LIST
  » → n ZEROS
  « { } { 0 1 } 2 n
    FOR k DUP SIZE → a s
      « π a HEAD * 2 k 3 - ^ / +
        k ZEROS EVAL a +
        a 1 k SUB 0 + REVLIST s 1 - ZEROS EVAL +
        ADD
        a k 1 + s SUB k 2 * ZEROS EVAL +
        ADD
      »
    NEXT
    π SWAP HEAD * 2 n 2 - ^ / +
  »
»

Kind regards
Thomas