decrease abs_eps and rel_eps in gsl_integration_qag

Author: acavelan

model/PkPd/Drug/LSTMDrugThreeComp.cpp

@@ -134,7 +134,7 @@ double LSTMDrugThreeComp::calculateFactor(const Params_fC& p, double duration) c
     F.params = static_cast<void*>(const_cast<Params_fC*>(&p));
 
     // NOTE: tolerances are arbitrary, but seem to be sufficient
-    const double abs_eps = 1e-2, rel_eps = 1e-2;
+    const double abs_eps = 1e-4, rel_eps = 1e-4;
     // NOTE: 1 through 6 are different algorithms of increasing complexity
     const int qag_rule = 1;     // alg 1 seems to be good enough
     double intfC, err_eps;