{
   assert (verbosity < 4);
 
-  _integrate( integrand, IntMethod, Limits, ncomp, 
-             nregions, neval, fail, integral,
-             error, prob, verbosity, userdata, Options );
+  return _integrate( integrand, IntMethod, Limits, ncomp, 
+                    nregions, neval, fail, integral,
+                    error, prob, verbosity, userdata, Options );
 }
 
 int cubaint::integrate (integrand_tc integrand, const intmethod& IntMethod, 
 {
   assert (verbosity < 4);
 
-  _integrate( integrand, IntMethod, Limits, 2*ncomp, 
-             nregions, neval, fail, integral,
-             error, prob, verbosity, userdata, Options );
+  return _integrate( integrand, IntMethod, Limits, 2*ncomp, 
+                    nregions, neval, fail, integral,
+                    error, prob, verbosity, userdata, Options );
 }
 
 template <typename integrandtype, typename compreal>
 
   std::cout << "int_0^1 dx (e-x^2+i0): " << std::endl;
   for (int itype=0; itype<4; itype++) {
     cubaint::integrate (&test_integrand_c, cubaint::intmethod(itype),
-                       {{0,1}}, 2,
+                       {{0,1}}, 1,
                        nregions,
                        neval, fail, &integral_c, 
                        &error_c, &prob_c, 0);