[I,abserr]=cubature_logan('f', 10^-10, 1, 1);
fprintf('\n \n[ I ]= %1.15e',I);
fprintf('\n \n');

