I need to calculate the following integral several thousands times per time step:
where and:
, and
So far I have implemented in Julia as:
