https://keroseneheaterreview.blogspot.com/2024/01/integral-of-intensity-function-in-python.html