#define max_s_n 20

double RFactorNApprox(double tau, int s, int k);