Fast_integration.ipynb gives a working version of the integration without alpha. Takes 15 seconds per integral right now