diff --git a/convoys/regression.py b/convoys/regression.py
index 6a2f2a2..4a23854 100644
--- a/convoys/regression.py
+++ b/convoys/regression.py
@@ -219,7 +219,7 @@ def callback(LL, value_history=[]):
         jac = autograd.grad(lambda x: -generalized_gamma_loss(x, *args))
 
         # Find the maximum a posteriori of the distribution
-        res = scipy.optimize.minimize(f, x0, jac=jac, method='SLSQP',
+        res = scipy.optimize.minimize(f, x0, jac=jac, method='L-BFGS-B',
                                       options={'maxiter': 9999})
         if not res.success:
             raise Exception('Optimization failed with message: %s' %