diff --git a/src/seismic_hazard_forecasting.py b/src/seismic_hazard_forecasting.py
index 271a1b8..dc7299f 100644
--- a/src/seismic_hazard_forecasting.py
+++ b/src/seismic_hazard_forecasting.py
@@ -185,7 +185,7 @@ def compute_IMT_exceedance(rx_lat, rx_lon, r, fr, p, lambdas, D, percentages_D,
     start = timer()       
         
     try: 
-        method='brentq'
+        method='bisect'
         logger.debug("Now trying Scipy " + method)
         output = root_scalar(exceedance_root_function, bracket=[IMT_min, IMT_max], rtol=rtol, method=method)
         gm_est = output.root