Stirling's approximation carried out to more terms gives: b = n/e + ln(2πn)/2e + (3(ln(2πn))^2 + 2)/(24en) + o(1/n) which becomes a good approximation for

Stirling's approximation carried out to more terms gives: b = n/e + ln(2πn)/2e + (3(ln(2πn))^2 + 2)/(24en) + o(1/n) which becomes a good approximation for large n: https://colab.research.google.com/drive/17cC8YUqKjEBLsNJkJ4IxA_LKDdg9Hflv https://colab.research.google.com/drive/17cC8YUqKjEBLsNJkJ4I...