Warning: Reached the limit on the maximum number of intervals in use.
Approximate bound on error is 6.1e-008. The integral may not exist, or
it may be difficult to approximate numerically. Increase MaxIntervalCount
to 822 to enable QUADGK to continue for another iteration.
> In quadgk>vadapt at 317
In quadgk at 208
In HGII at 29
In mymainfunction at 177