double RFactorKappaApprox(double tau, double k);