Subintervals for integration

[erossiter]

Subintervals are currently set at 10, but Alex noted this was a pretty arbitrary number when he set it. Do you want to change it?

This explains the role of subintervals:
https://www.gnu.org/software/gsl/manual/html_node/QAG-adaptive-integration.html