pub const NEWTON_TOLERANCE: f64 = 1e-15;
Convergence tolerance for Newton-Raphson iteration
Iteration stops when |delta| < tolerance * (1 + |x|).
Value: 1e-15 (approximately machine epsilon)