Function beta_continued_fraction_jet
fn beta_continued_fraction_jet(
a: Jet2<f64>,
b: Jet2<f64>,
x: Jet2<f64>,
) -> Jet2<f64>Expand description
beta_continued_fraction in Jet2 on its parameters: the termwise
derivative of the fraction converges to the incomplete beta’s (Boik &
Robison-Cox 1998), and the loop runs until the value and both
derivative components of each factor have settled.