I'm trying to model the number of possible chess games after n plies.(perft results A048987): 1, 20, 400, 8902, 197281, 4865609, 119060324, 3195901860, 84998978956, 2439530234167, 69352859712417, 2097651003696806, 62854969236701747, 1981066775000396239, 61885021521585529237, 2015099950053364471960.
I expect this sequence to tend towards an upper asymptote. It starts off with an almost factorial growth which can be seen by the consecutive ratios, which (rounded off to 2 decimal places) are: 20, 20, 22.25, 22.16, 24.66, 24.47, 26.84, 26.6, 28.7, 28.43, 30.25, 29.96, 31.52, 31.24, 32.56. This is approximately $$x+20+\frac{\cos(\pi x)-1}{2}$$ and thus $$\prod_{x=0}^{n}(x+20+\frac{\cos(\pi x)-1}{2})$$ provides a good approximation for the first few terms of the sequence, although it does not tend towards an upper asymptote, and it already starts to deviate from the actual value in the last 3 terms. I'm looking to fit this data to a function, like the generalised logistic function, but one that starts off with an almost factorial growth, rather than an exponential one. What's the right model to use here?