slow loop in calculation of definite integral