fn finish<const DERIVATIVE: bool>( value: Scaled, slope: Scaled, y: f64, ) -> Result<(f64, f64)>
Converts (P, (1 - x²) P') to doubles, rejecting results beyond the double range.
(P, (1 - x²) P')