> Do each function separately. The routine that scipy.integrate.quad uses > adapts to the local conditions of the function (when the function is > flat, it uses fewer samples; when steep, more). It's not so easy to do it separately. I want to integrate function that includes linalg.expm - as you see calculation the function is quite expensive. T.Kaz.