Function theta_quadrature 

pub fn theta_quadrature(s: Screening, eps: f64, beta: f64) -> f64

theta_quadrature_n with DEFAULT_NODES nodes.