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),
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),