Plot below is for \x → incompleteBeta 0.5 0.5 x. Both X and Y axes are in ULPs relative to starting point of plot. Clearly full precision of incompleteBeta is not realised. While relatively minor issue by itself it causes problems in invIncompleteBeta because algorithm fails to converge and jumps between two values.
Plot below is for
\x → incompleteBeta 0.5 0.5 x
. Both X and Y axes are in ULPs relative to starting point of plot. Clearly full precision of incompleteBeta is not realised. While relatively minor issue by itself it causes problems ininvIncompleteBeta
because algorithm fails to converge and jumps between two values.