[WP34s] Regularized incomplete Beta function

[Paul Dale]

(05-11-2014 01:50 PM)Dieter Wrote:  1. Catching the case p=1/2 at the beginning is not neccessary. The estimate will be exactly zero, and this is what will be returned directly after the first iteration, i.e. almost immediately.

I'm getting "no root found" without that check. That's why I added it.

Quote:2. This is part of the iteration loop:

Code:

...

The stack shuffling command obviously was added because the incomplete Beta function now accepts the three parameters in a different order than before. After this, right before the call to the Beta function, the T-register holds 1/2. After that call, the comment says that both Z and T hold the current estimate x – which was completely lost before the call (was neither in X nor Y nor Z nor T). ?!?

The eight level stack is in play in xrom code. I just didn't document the stack diagrams fully.


- Pauli