<html><head><meta name="color-scheme" content="light dark"></head><body><pre style="word-wrap: break-word; white-space: pre-wrap;">\documentclass[12pt,reqno]{article}

\usepackage[usenames]{color}
\usepackage{amssymb}
\usepackage{graphicx}
\usepackage{amscd}

\let\a\alpha

\usepackage[colorlinks=true,
linkcolor=webgreen,
filecolor=webbrown,
citecolor=webgreen]{hyperref}

\definecolor{webgreen}{rgb}{0,.5,0}
\definecolor{webbrown}{rgb}{.6,0,0}

\usepackage{color}
\usepackage{fullpage}
\usepackage{float}

\usepackage{psfig}
\usepackage{graphics,amsmath,amssymb}
\usepackage{amsthm}
\usepackage{amsfonts}
\usepackage{latexsym}
\usepackage{epsf}

\setlength{\textwidth}{6.5in}
\setlength{\oddsidemargin}{.1in}
\setlength{\evensidemargin}{.1in}
\setlength{\topmargin}{-.1in}
\setlength{\textheight}{8.4in}

\newcommand{\seqnum}[1]{\href{http://oeis.org/#1}{\underline{#1}}}

\begin{document}

\begin{center}
\epsfxsize=4in
\leavevmode\epsffile{logo129.eps}
\end{center}

\theoremstyle{plain}
\newtheorem{theorem}{Theorem}
\newtheorem{corollary}[theorem]{Corollary}
\newtheorem{lemma}[theorem]{Lemma}
\newtheorem{proposition}[theorem]{Proposition}

\theoremstyle{definition}
\newtheorem{definition}[theorem]{Definition}
\newtheorem{example}[theorem]{Example}
\newtheorem{conjecture}[theorem]{Conjecture}

\theoremstyle{remark}
\newtheorem{remark}[theorem]{Remark}

\begin{center}
\vskip 1cm{\LARGE\bf 
Asymptotic Expansions of Central Binomial \\
\vskip .1in
Coefficients and Catalan Numbers
}
\vskip 1cm
\large
Neven Elezovi\'c \\
Faculty of Electrical Engineering and Computing \\
University of Zagreb \\
Unska 3 \\
10000 Zagreb \\
Croatia \\
\href{mailto:neven.elez@fer.hr}{\tt neven.elez@fer.hr} \\
\end{center}

\vskip .2 in

\begin{abstract}
We give a systematic view of the asymptotic expansion of two well-known sequences, the central binomial coefficients and the Catalan numbers.
The main point is explanation of the nature of the best shift in variable $n$, in order to obtain ``nice'' asymptotic expansions. We also give a complete asymptotic expansion of partial sums of these sequences.
\end{abstract}

\vskip .2in

\section{Introduction}

        One of the most beautiful formulas in mathematics is the classical Stirling approximation of the factorial function:
        $$
                n!\approx\sqrt{2\pi n}\biggl(\frac ne\biggr)^n.
        $$
	This is the beginning of the following full asymptotic expansions
	\cite{abram,luke}, Laplace expansion:
        \begin{equation}
	\label{laplace}
                n!\sim \sqrt{2\pi n}\biggl(\frac ne\biggr)^n
                \biggl(
                1+\frac1{12n}+\frac1{288n^2}
                -\frac{139}{51840n^3}
                -\frac{571}{2488320n^4}+\dots
                \biggr)
        \end{equation}
        and Stirling series
        \begin{equation}
	\label{stirling}
                n!\sim \sqrt{2\pi n}\biggl(\frac ne\biggr)^n
                \exp\biggl(
                \frac1{12n}-\frac1{360n^3}+\frac1{1260n^5}+\dots
                \biggr).
        \end{equation}


        The central binomial coefficient has a well-known asymptotic approximation; see e.g., \cite[p.\ 35]{luke}:
        \begin{equation}
        \label{binom}
                \binom{2n}n\sim \frac{2^{2n}}{\sqrt{n\pi}}\biggl[1-\frac1{8n}+\frac1{128n^2}
                +\frac{5}{1024n^3}-\frac{21}{32768n^4}+\dots\biggr].
        \end{equation}
	

	Luschny \cite{luschny} gives the following nice expansions:
	\begin{equation}
	\label{shiftbinom}
		\binom{2n}n\sim\frac{4^n}{\sqrt{N\pi/2}}
		\left(2-\frac2{N^2}+\frac{21}{N^4}-\frac{671}{N^6}
		+\frac{45081}{N^8}\right)
	\end{equation}
	where $N=8n+2$, and for the Catalan numbers
	\begin{equation}
	\label{shiftCatalan}
		\frac1{n+1}\binom{2n}n\sim\frac{4^{n-2}}{M\sqrt{M\pi}}
		\left(128+\frac{160}{M^2}+\frac{84}{M^4}+\frac{715}{M^6}
		-\frac{10180}{M^8}\right)
	\end{equation}
	where $M=4n+3$.
	Here, for the sake of the beauty, the exact value $45080\frac34$ is replaced by 45081, and $10179\frac{13}{16}$ is replaced by 10180.

We would like to thank the anonymous referee who brought to our
attention the existence of the manuscript  \cite{KS} where similar
problems are treated. D.~Kessler and J.~Schiff proved that expansion
mentioned above contains only odd powers of $n+\frac14$ (for the
central binomial coefficient) and $n+\frac34$ (for Catalan numbers). In
this paper we explain why this happens.

	The main subject of this paper is to explain why $N=8n+2$ and $M=4n+3$ are the best choices in such expansions, and also to obtain general form of these expansions, especially in the case of the Laplace expansions. In the last section, the asymptotic expansion of the partial sums of binomial coefficients and Catalan numbers are derived, using a simple and efficient recursive algorithm.



 		
\section{Central binomial coefficients}



	Although the central binomial coefficient is expressed as $\Gamma(2n+1)/\Gamma(n+1)^2$, expansions (\ref{laplace}) or (\ref{stirling}) cannot be used for direct derivation of (\ref{binom}). Instead, one should use the asymptotic expansion of the ratio of two gamma functions. In the standard reference \cite{luke}, the connection with generalized Bernoulli polynomials 
	is used. This approach is improved in a series of recent papers
	\cite{be1}--\cite{BES}.
	Namely, from the duplication formula for the gamma function we have   
	\begin{equation}\label{veza-2}
                \binom{2n}{n}=\frac{\Gamma(2n+1)}{\Gamma(n+1)^2}
                =\frac{4^{n}}{\sqrt\pi}\cdot
                \frac{\Gamma(n+\frac12)}{\Gamma(n+1)}.
        \end{equation}
        
        In \cite{be3}, the following general asymptotic expansion of the quotient of two gamma functions is given:
	\begin{equation}\label{be-m}
		\frac{\Gamma(x+t)}{\Gamma(x+s)}\sim
		x^{t-s}\left(\sum_{m=0}^{\infty}P_m(t,s,r)x^{-m}\right)^{1/r}.
	\end{equation}
	Here, $s$, $t$ and $r\ne0$ are real numbers.
	Coefficients $P_m=P_m(t,s,r)$ are polynomials defined by
	\begin{align}
		&amp;P_0(t,s,r)=1,\\
		&amp;P_m(t,s,r)=\frac{r}{m}\sum_{k=1}^{m}(-1)^{k+1}\frac{B_{k+1}(t)-B_{k+1}(s)}{k+1}P_{m-k}(t,s,r)
	\label{P1}
	\end{align}
	and $B_k(t)$ stands for the Bernoulli polynomials. 
	
	In the sequel, we shall use the following properties of Bernoulli polynomials and Bernoulli numbers:
	\begin{equation}
	\begin{aligned}
		&amp;(-1)^nB_n(-x)=B_n(x)+nx^{n-1},\\
		&amp;B_n(1+x)=B_n(x)+nx^{n-1},\\
		&amp;B_{2n+1}=0,  \qquad (n\ge1),\\
		&amp;B_n(0)=(-1)^nB_n(1)=B_n,\\
		&amp;B_n(\tfrac{1}{2})=-(1-2^{1-n})B_n,\\
		&amp;B_n(-\tfrac{1}{2})=-(1-2^{1-n})B_n+(-1)^n\frac n{2^{n-1}},\\
		&amp;B_n(\tfrac{1}{4})=-2^{-n}(1-2^{1-n})B_n-n4^{-n}E_{n-1},\\
		&amp;B_n(\tfrac{3}{4})=(-1)^{n+1}2^{-n}(1-2^{1-n})B_n+n4^{-n}E_{n-1}.
	\end{aligned}
	\label{Bern}
	\end{equation}	


	Denote $x=n+\alpha$, $t=1/2-\alpha$, $s=1-\alpha$. Applying \eqref{veza-2}, we have
        \begin{equation}
        \label{binom-m}
		\binom{2n}{n}\sim\frac{2^{2n}}{\sqrt{x\pi}}\left(\sum_{k=0}^{\infty}\frac{P_k}{x^k}\right)^{1/r}
	\end{equation}
	where sequence $(P_n)$ is defined by $P_0=1$ and 
	$$
	      P_m=\frac{r}{m} \sum_{k=1}^{m} \frac{B_{k+1}(\frac12+\alpha)-B_{k+1}(\alpha)}{k+1}P_{m-k},\qquad m\geq1.
	 $$ 
	 
	
	In order to obtain a useful formula, the parameter $\alpha$ should be chosen in such a way that the values of Bernoulli polynomials can be         (easily) calculated.
	Some simplifications are also possible if these coefficients are connected in a way which reduces complexity of this expression. Therefore, the following choices are indicated:

	1) $\a=0$: this gives ``natural'' expansion in terms of powers of $n$. Although natural, this choice usually is not the best one.

	2) $\a=\frac12$: this value leads to easily computable coefficients.

	3) $\frac12-\a=1-(1-\a)$: wherefrom it follows $\a=\frac14$. Here, the symmetry property of Bernoulli polynomials is used, and this is the best choice for $\a$.

	4) $\frac12-\a=-(1-\a)$, i.e., $\a=\frac34$: This choice will also reduce computation.

	The value of the Bernoulli polynomials may be calculated explicitly (in terms of Bernoulli and Euler numbers) for some other constants $\a$, for example $\a=\frac16$, but the values will be ``complicated'' compared to the ones chosen above.




	The expansions of the central binomial coefficients are given in
	the following theorem.

\begin{theorem}
	\label{th1}
	The following asymptotic expansion is valid:
	\begin{equation}
		\binom{2n}n\sim\frac{4^{n}}{\sqrt{\pi(n+\alpha)}}
		\left( \sum_{m=0}^{\infty}P_m(\alpha)(n+\alpha)^{-m}\right)^{1/r},
	\end{equation}
	where $P_0=1$ and
	\begin{enumerate}
	\item for $\alpha=0$
		\begin{equation}\label{alpha-0}
			P_m=\frac{r}{m}\sum_{k=1}^{\lfloor(m+1)/2\rfloor}\frac{(2^{-2k}-1)B_{2k}}{k} P_{m-2k+1};
		\end{equation}
	\item for $\alpha=\tfrac{1}{4}$
		\begin{equation}\label{alpha-14}
			P_m=\frac{r}{m}\sum_{k=1}^{\lfloor m/2\rfloor}2^{-2k-1}E_kP_{m-2k};
		\end{equation}
	\item for $\alpha=\tfrac{1}{2}$ 
		\begin{equation}\label{alpha-12}
			P_m=\frac{r}{m}\sum_{k=1}^{\lfloor(m+1)/2\rfloor}\frac{(1-2^{-2k})B_{2k}}{k} P_{m-2k+1};
		\end{equation}
	\item for $\alpha=\tfrac{3}{4}$ 
		\begin{equation}\label{alpha-34}
			P_m=\frac{r}{m}\sum_{k=1}^{\lfloor m/2\rfloor}2^{-2k-1}(2-E_k)P_{m-2k};
		\end{equation}
	\end{enumerate}
\end{theorem}

\begin{proof}
	Let us write
	\begin{equation*}
		b_k(\alpha)=[B_{k+1}(\tfrac12+\alpha)-B_{k+1}(\alpha)].
	\end{equation*}

	We have
	$$
		b_{k}(0)=B_{k+1}(\tfrac12)-B_{k+1}.
	$$
	This value is equal to 0 for even $k$, and equal to $(2^{-k}-2)B_{k+1}$ for odd $k$, and hence \eqref{alpha-0} follows.

	For $\a=\frac14$,
	$$
		b_k(\tfrac14)=B_{k+1}(\tfrac14)[(-1)^{k+1}-1]
	$$
	This is equal to 0 for odd $k$, and equal to $(k+1)2^{-2k-1}E_k$ for even $k$. 

	Further,
	$$
		b_k(\tfrac12)=(-1)^{k+1}[B_{k+1}(0)-B_{k+1}(\tfrac12)]
	$$
	and \eqref{alpha-12} follows similarly as in the first case.

	Finally,
	$$
		b_k(\tfrac34)=B_{k+1}(\tfrac14)[1-(-1)^{k+1}]+(k+1)4^{-k}
	$$
	As before, this is equal to 0 for odd $k$, and equal to
	$(k+1)2^{-2k-1}(2-E_k)$ for even $k$. Thus \eqref{alpha-34} follows, completing the proof of the theorem.
\end{proof}

	It is obvious that the choice $\alpha=\frac14$ is superior to others.
	In fact, the equation
	$$
		B_{k+1}(1/2-\alpha)=B_{k+1}(1-\alpha)
	$$
	is an identity for each odd $k$ only if $\alpha=\frac14$. Hence, 
	this value of $\a$ is unique with the property that the asymptotic expansion reduces to even terms.

	We shall give the first few terms of asymptotic expansions of $P_m(\alpha)$ for the values of the 
shift $\alpha$ observed in Theorem \ref{th1}. Using $r=1$ we get:
\begin{multline}\label{A-0}
	\binom{2n}{n}\sim\frac{4^n}{\sqrt{\pi n}}\bigg[1
	-\frac{1}{8n}+
	\frac{1}{128n^2}
	+\frac{5}{1024n^3}
	-\frac{21}{32768n^4}\\
	-\frac{399}{262144n^5}
	+\frac{869}{4194304n^6}+
	\cdots\bigg],
\end{multline}
\begin{multline}\label{Asymptotic-Wn-alpha-14}
	\binom{2n}{n}\sim\frac{4^n}{\sqrt{\pi(n+\frac14)}}\bigg[1
	-\frac{1}{64\left(n+\frac{1}{4}\right)^2}
	+\frac{21}{8192\left(n+\frac{1}{4}\right)^4}
	-\frac{671}{524288\left(n+\frac{1}{4}\right)^6}\\
	+\frac{180323}{134217728\left(n+\frac{1}{4}\right)^8}
	-\frac{20898423}{8589934592\left(n+\frac{1}{4}\right)^{10}}
	+\cdots\bigg],
\end{multline}
\begin{multline}\label{Asymptotic-Wn-alpha-12}
	\binom{2n}{n}\sim\frac{4^n}{\sqrt{\pi(n+\frac{1}{2})}}\bigg[1
	+\frac{1}{8\left(n+\frac{1}{2}\right)}
	+\frac{1}{128\left(n+\frac{1}{2}\right)^2}
	-\frac{5}{1024\left(n+\frac{1}{2}\right)^3}\\
	-\frac{21}{32768\left(n+\frac{1}{2}\right)^4}
	+\frac{399}{262144\left(n+\frac{1}{2}\right)^5}
	+\frac{869}{4194304\left(n+\frac{1}{2}\right)^6}
	+\cdots\bigg],
\end{multline}
\begin{multline}\label{Asymptotic-Wn-alpha-34}
	\binom{2n}{n}\sim\frac{4^n}{\sqrt{\pi(n+\frac{3}{4})}}\bigg[1
	+\frac{1}{4\left(n+\frac{3}{4}\right)}
	+\frac{5}{64\left(n+\frac{3}{4}\right)^2}
	+\frac{5}{256\left(n+\frac{3}{4}\right)^3}\\
	+\frac{21}{8192\left(n+\frac{3}{4}\right)^4}
	+\frac{21}{32768\left(n+\frac{3}{4}\right)^5}
	+\frac{715}{524288\left(n+\frac{3}{4}\right)^6}
	+\cdots\bigg],
\end{multline}


\section{The role of exponent $r$}

	The coefficients $P_m(t,s,r)$ are polynomials in $r$ of degree $m$,
which follows directly from the recursive formula \eqref{P1}.
	
\begin{theorem}
	Let $\alpha=0$. Then 
	\begin{equation}
		P_m(t,s,-r)=(-1)^mP_m(t,s,r).
	\label{r}
	\end{equation}
\end{theorem}

\begin{proof}
	By induction. \eqref{r} holds for $m=0$ and $m=1$. The rest is obvious
	from \eqref{alpha-0}.
\end{proof}

	As a corollary, we get that the coefficients of the expansion for $r=-1$ are identical, up to the sign of odd powers, to the coefficients of the expansion for $r=1$. Therefore,
	from \eqref{A-0} it follows immediately that
\begin{align}\label{A-0r}
	\binom{2n}{n}^{-1}&amp;\sim\frac{\sqrt{\pi n}}{4^n}\bigg[1
	+\frac{1}{8n}+
	\frac{1}{128n^2}
	-\frac{5}{1024n^3}
	-\frac{21}{32768n^4}\nonumber\\
	&amp;\qquad
	+\frac{399}{262144n^5}+
	\frac{869}{4194304n^6}+
	\cdots\bigg].
\end{align}




	
	
	Various choices of $r$ may give useful expansions.
	For example, $r=4$ and $N=4n$ leads to a good approximation with the first two terms:
	\begin{equation*}
	\binom{2n}{n}\sim\frac{2^{2n+1}}{\sqrt{\pi N}}\sqrt[4]{1-\frac{2}{N}+\frac{2}{N^2}-\frac{2}{N^4}-\frac{4}{N^5}-\frac{12}{N^6}+\ldots}
	\end{equation*}
	while $r=2$ and $N=8n+2$ gives a good square root analogue of the formula \eqref{shiftbinom}:
	$$
	\binom{2n}{n}\sim\frac{4^{n+1}}{\sqrt{\pi N}}\sqrt{\frac12-\frac{1}{N^2}+\frac{11}{N^4}-\frac{346}{N^6}
	+\frac{22931}{N^8}+\ldots}.
	$$



\section{Catalan numbers} 

	The standard definition of Catalan numbers is given by a recurrence relation, $C_0=1$ and
	$$
		C_{n+1}=\sum_{k=1}^n C_kC_{n-k}.
	$$
	Catalan numbers occur in various situations. For instance, 
	Stanley \cite[p.\ 219]{Stan} explains 66 such situations.

	The starting point for us will be the following explicit formula:
	\begin{equation}
		C_n=\frac1{n+1}\binom{2n}n=\frac{\Gamma(2n+1)}{\Gamma(n+1)\Gamma(n+2)}
	\end{equation}
	Hence, Catalan numbers can be expressed as a ratio of two gamma functions
	\begin{equation}
	\label{catalan}
                C_n
                =\frac{4^{n}}{\sqrt\pi}\cdot
                \frac{\Gamma(n+\frac12)}{\Gamma(n+2)}.
        \end{equation}
	
	Putting $x=n+\alpha$, $t=\frac12-\a$, $s=2-\a$, from \eqref{P1} we get
	\begin{equation}
	\label{catexp}
		C_n\sim\frac{4^n}{\sqrt\pi}x^{-3/2}
		\left(\sum_{m=0}^\infty P_mx^{-m}\right)^{1/r},
	\end{equation}
	with $P_0=1$ and
	\begin{equation}
		P_m=\frac rm\sum_{k=1}^mc_k(\a)P_{m-k},
	\end{equation}
	where we denote
	\begin{equation}
		c_k(\a)=\frac{B_{n+1}(\a+\frac12)-B_{k+1}(\a-1)}{k+1}.
		\label{cal}
	\end{equation}


	As before, 0 and $\frac12$ are the natural choice for $\a$. Two other good values follow from $\alpha+\frac12=1-(\alpha-1)$ and $\alpha+\frac12=-(\alpha-1)$, wherefrom one gets $\a=\frac34$ and $\a=\frac14$, respectively. 

\begin{theorem}
	\label{th2}
	The following asymptotic expansion holds:
	\begin{equation}
		C_n\sim\frac{4^n}{\sqrt\pi}(n+\alpha)^{-3/2}
		\left( \sum_{m=0}^{\infty}P_m(\alpha)(n+\alpha)^{-m}\right)^{1/r},
	\end{equation}
	where $P_0=1$ and
	\begin{enumerate}
	\item for $\alpha=0$
		\begin{equation}\label{calpha-0}
			P_m=\frac{r}{m}\sum_{k=1}^{m}
			\left[\frac{(2^{-k}-2)B_{k+1}}{k+1}+(-1)^k\right] P_{m-k};
		\end{equation}
	\item for $\alpha=\tfrac{1}{2}$ 
		\begin{equation}\label{calpha-12}
			P_m=\frac{r}{m}\sum_{k=1}^{m}
			\left[\frac{(2-2^{-k})B_{k+1}}{k+1}+\frac{(-1)^{k+1}}{2^k}\right] P_{m-k};
		\end{equation}
	\item for $\alpha=\tfrac{3}{4}$
		\begin{equation}\label{calpha-34}
			P_m=\frac{r}{m}\sum_{k=1}^{\lfloor m/2\rfloor}2\cdot4^{-2k-1}(4-E_{2k}) P_{m-2k};
		\end{equation}
	\item for $\alpha=\tfrac{1}{4}$ 
		\begin{equation}\label{calpha-14}
			P_m=\frac{r}{m}\sum_{k=1}^{m}[2^{-2k-1}E_k+(-\tfrac34)^k]P_{m-k}.
		\end{equation}
	\end{enumerate}
\end{theorem}

\begin{proof}

	We need to compute explicit coefficients in formula (\ref{catexp}).

	1) $\a=0$:
	\begin{align*}
		c_k(0)&amp;=\frac{B_{k+1}(\frac12)-B_{k+1}(-1)}{k+1}
	\end{align*}
	Using (\ref{Bern}), 
	we have
	$$
		c_k(0)=\frac{-(2-2^{-k})B_{k+1}}{k+1}-(-1)^{k+1}
	$$
	and \eqref{calpha-0} follows.

	2) $\a=\frac12$:
	\begin{align*}
		c_k(\tfrac12)&amp;=\frac{B_{k+1}(1)-B_{k+1}(-\frac12)}{k+1}
	\end{align*}
	Using (\ref{Bern}), 
	we have
	$$
		c_k(\tfrac12)=\frac{(2-2^{-k})B_{k+1}}{k+1}+\frac{(-1)^k}{2^k}.
	$$

	3) $\a=\frac34$:
	$$
		c_k(\tfrac34)=\frac{B_{k+1}(\frac54)-B_{k+1}(-\frac14)}{k+1}
	$$
	Since
	\begin{align*}
		B_{k+1}(\tfrac54)&amp;=B_{k+1}(\tfrac14)+(k+1)\frac1{4^k},\\
		B_{k+1}(-\tfrac14)&amp;=(-1)^{k+1}\biggl[B_{k+1}(\tfrac14)+\frac{k+1}{4^k}\biggr]
	\end{align*}
	it follows
	$$
		c_k(\tfrac34)=\frac{1+(-1)^k}{k+1}\biggl[B_{k+1}(\tfrac14)
		+\frac{k+1}{4^k}.
	$$
	Therefore, for odd $k$ we have $c_k(\tfrac34)=0$. For even $k$ it follows
	\begin{align*}
		c_k(\tfrac34)&amp;=\frac{2B_{k+1}(\tfrac14)}{k+1}+\frac2{4^k}\\
		&amp;=-2\cdot4^{-k-1}E_{k}+2\cdot4^{-k}.
	\end{align*}
	This proves \eqref{calpha-34}.

	4) $\a=\tfrac14$:
	$$
		c_k(\tfrac14)=\frac{B_{k+1}(\frac34)-B_{k+1}(-\frac34)}{k+1}
	$$
	Since
	$$
		B_{k+1}(-\tfrac34)=(-1)^{k+1}\left[
		B_{k+1}(\tfrac34)+(k+1)(\tfrac34)^k\right]
	$$
	it follows
	$$
		c_k(\tfrac54)=\frac1{k+1}B_{k+1}(\tfrac34)
		[(-1)^{k}+1]+(-1)^{k}(\tfrac34)^k.
	$$
	Therefore, for odd $k$ it holds
	$$
		c_k(\tfrac14)=(-\tfrac34)^k.
	$$
	For even $k$, after reducing in a similar way as before, we get
	$$
		c_k(\tfrac14)=-2\cdot 4^{-k-1}E_k+(-\tfrac34)^k.
	$$
	For odd $k$ this values coincides with previous one, since $E_{2n+1}=0$. Hence, \eqref{calpha-14} is proved.
\end{proof}


	For the convenience of the reader, here is the short list of the observed coefficients, for $r=1$:
\begin{multline}\label{A-0c}
	C_n\sim\frac{4^n}{\sqrt{\pi n^3}}\bigg[1
	-\frac{9}{8n}+
	\frac{145}{128n^2}
	-\frac{1155}{1024n^3}
	+\frac{36939}{32768n^4}\\
	-\frac{295911}{262144n^5}+
	-\frac{4735445}{4194304n^6}+
	\cdots\bigg],
\end{multline}
\begin{multline}\label{Asymptotic-Wn-alpha-12c}
	C_n\sim\frac{4^n}{\sqrt{\pi(n+\frac{1}{2})^3}}\bigg[1
	-\frac{3}{8\left(n+\frac{1}{2}\right)}
	+\frac{25}{128\left(n+\frac{1}{2}\right)^2}
	-\frac{105}{1024\left(n+\frac{1}{2}\right)^3}\\
	+\frac{1659}{32768\left(n+\frac{1}{2}\right)^4}
	-\frac{6237}{262144\left(n+\frac{1}{2}\right)^{5}}
	+\frac{50765}{4194304\left(n+\frac{1}{2}\right)^{6}}
	+\cdots\bigg].
\end{multline}
\begin{multline}\label{Asymptotic-Wn-alpha-34c}
	C_n\sim\frac{4^n}{\sqrt{\pi(n+\frac{3}{4})^3}}\bigg[1
	+\frac{5}{64\left(n+\frac{3}{4}\right)^2}
	+\frac{21}{8192\left(n+\frac{3}{4}\right)^4}
	+\frac{715}{524288\left(n+\frac{3}{4}\right)^6}\\
	-\frac{162877}{134217728\left(n+\frac{3}{4}\right)^8}
	+\frac{19840275}{8589934592\left(n+\frac{3}{4}\right)^{10}}
	+\cdots\bigg].
\end{multline}
\begin{multline}\label{Asymptotic-Wn-alpha-14c}
	C_n\sim\frac{4^n}{\sqrt{\pi(n+\frac{1}{4})^3}}\bigg[1
	-\frac{3}{4\left(n+\frac{1}{4}\right)}
	+\frac{35}{64\left(n+\frac{1}{4}\right)^2}
	-\frac{105}{256\left(n+\frac{1}{4}\right)^3}\\
	+\frac{2541}{8192\left(n+\frac{1}{4}\right)^4}
	-\frac{7623}{32768\left(n+\frac{1}{4}\right)^{5}}
	+\frac{90805}{524288\left(n+\frac{1}{4}\right)^{6}}
	+\cdots\bigg].
\end{multline}
	
	As one can see, the expansion in terms of $n+\frac34$ has additional property that it contains only odd terms.		

\begin{corollary}
	The value $\a=\frac34$ is the unique value for which asymptotic expansion of Catalan numbers contains only odd terms.
\end{corollary}

\begin{proof}
	If this is the case, then coefficient $c_1(\a)$ from (\ref{cal}) must vanish:
	$$
		B_2(\alpha+\tfrac12)-B_2(\alpha-1)=0.
	$$
	From this one obtain $\a=\tfrac34$.
\end{proof}
	
	


\section{Expansions of Stirling's type}

	Stirling expansion of the factorial function (\ref{stirling}) includes the exponential function. Using the known asymptotic expansion 
	$$
		\ln\Gamma(x+t)=(x+t-\tfrac12)\ln x-x
		+\tfrac12\ln(2\pi)+\sum_{k=1}^\infty
		\frac{(-1)^{k+1}B_{k+1}(t)}{k(k+1}\,x^{-k}
	$$
	we immediately get
        \begin{equation}
        \label{exp1}
                \frac{\Gamma(x+t)}{\Gamma(x+s)} 
		\sim x^{t-s}\cdot \exp\left(\sum_{k=1}^{\infty}Q_{k}(t,s)x^{-k}\right)
        \end{equation}
        where $Q_k(t,s)$ are polynomials of order $k$ defined by
        \begin{equation}
	\label{exp2}
                Q_k(t,s)=\frac{(-1)^{k+1}B_{k+1}(t)-B_{k+1}(s)}{k(k+1)}.
        \end{equation}


	Hence, we obtain:


\begin{theorem}
	The binomial coefficient has the following asymptotic expansions of Stirling's type:
	\begin{equation}
	\label{39}
		\binom{2n}n\sim\frac{4^n}{\sqrt{\pi n}}
		\exp\biggl(\sum_{k=1}^\infty \frac{(2^{-2k}-1)B_{2k}}{k(2k-1)}n^{-2k+1}\biggr) 
	\end{equation}
	\begin{equation}
		\binom{2n}n\sim\frac{4^n}{\sqrt{\pi (n+\frac14)}}
		\exp\biggl(\sum_{k=1}^\infty \frac{2^{-4k-2}E_{2k}} {k}n^{-2k}\biggr). 		
	\end{equation}
\end{theorem}

	The proof is already carried out in the Theorem~\ref{th1}.

\bigskip

	A similar result can be stated for Catalan numbers:

\begin{theorem}
	
	\begin{equation}
		C_n\sim\frac{4^n}{n\sqrt{\pi n}}
		\exp\biggl(\sum_{k=1}^\infty (-1)^{k+1}\frac{(2^{-k}-2)B_{k+1}-k-1}{k(k+1)}n^{-k}\biggr) 
	\end{equation}
	\begin{equation}
	\label{42}
		C_n\sim\frac{4^n}{\sqrt{\pi(n+\frac34)^3}}
		\exp\biggl(\sum_{k=1}^\infty \frac{2^{-4k-2}(4-E_{2k})} {k}(n+\tfrac34)^{-2k}\biggr). 		
	\end{equation}
\end{theorem}
	
	Formulae (\ref{39})--(\ref{42}) are derived in the manuscript \cite{KS}.


\section{The sum of binomial coefficients and Catalan numbers}

	In a recent paper, Mattarei \cite{Mat} proves the following asymptotic expansions:
	\begin{gather}
	\label{bins}
		\sum_{k=0}^n\binom{2k}k=\frac{4^{n+1}}{3\sqrt{\pi n}}
		\bigl(1+\frac1{24n}+\frac{59}{384n^2}+\frac{2425}{9216n^3}+O(n^{-4})\bigr)\\
		\sum_{k=0}^nC_n=\frac{4^{n+1}}{3n\sqrt{\pi n}}
		\bigl(1-\frac5{8n}+\frac{475}{384n^2}+\frac{1225}{9216n^3}+O(n^{-4})\bigr)
	\label{cats}
	\end{gather}
	The calculation was tedious. For example, it relies on a computer
	algebra system,
	since it is based on the following formulas:
	\begin{gather*}
		\sum_{k=0}^n\binom{2k}k=\frac43\binom{2n}n\sum_{j=0}^m
		\frac{3^{-j}}{(2n-1)(2n/3-1)\cdot(2n/(2j-1)-1)}+O(4^nn^{-m-3/2}),\\
		\sum_{k=0}^nC_n=\frac23\binom{2n}n\sum_{j=0}^m
		\frac{(-3)^j+3^{-j}}{(2n-1)(2n/3-1)\cdot(2n/(2j+1)-1)}+O(4^nn^{-m-5/2})
	\end{gather*}
	The final calculation of the coefficients \eqref{bins} and \eqref{cats} was carried out using {\sc Maple}, with $m=4$.

	We shall derive an efficient algorithm for recursive calculations of asymptotic expansions of this and similar sums, which enables an easy calculation of the arbitrary coefficient in these expansions.

	The theorem will be formulated in such a way that it may be easily applied to both binomial and Catalan sums. It is evident that a similar statement is valid for the asymptotic expansion of the sum of more general functions.

\begin{theorem}
	Suppose that $a(n)$ has the following expansion, $P_0(\a)=1$ and
	\begin{equation}
	\label{coef-a}
		a(n)\sim\frac{4^n}{\sqrt{\pi}}\sum_{k=0}^\infty P_k(\a)(n+\a)^{-k-r},
	\end{equation}
	where $r&gt;0$ is a real number. Then
	\begin{equation}
		\sum_{k=0}^n a(k)\sim\frac 43\cdot\frac{4^{n+1}}{\sqrt\pi}\sum_{k=0}^\infty
		S_k(\a)(n+\a)^{-k-r}
	\end{equation}
	where the coefficients of this expansion satisfy $S_0(\a)=1$ and
	\begin{equation}
	\label{coefs}
		S_k(\a)=P_k(\a)+\frac13\sum_{j=0}^{k-1}(-1)^{k-j}\binom{-j-r}{k-j}S_j(\a)
	\end{equation}
\end{theorem}



\begin{proof}
	Denote $\Sigma(n)=\sum_{k=0}^n a(k)$. Suppose that
	$$
		\Sigma(n)\sim C\cdot\frac{4^n}{\sqrt\pi}n^{-r}+O(n^{-r-1}).
	$$
	Then
	\begin{align*}
		\Sigma(n)&amp;\sim a(n)+\Sigma(n-1)\\
		&amp;\sim\frac{4^n}{\sqrt\pi}\,n^{-r}+C\cdot\frac{4^{n-1}}{\sqrt\pi}\,(n-1)^{-r}
		+O(n^{-r-1})\\
		&amp;\sim\frac{4^n}{\sqrt\pi}\,n^{-r}+C\cdot\frac{4^{n-1}}{\sqrt\pi}\,n^{-r}
		+O(n^{-r-1})\\
		&amp;\sim C\cdot \frac{4^n}{\sqrt\pi}n^{-r}+O(n^{-r-1})
	\end{align*}
	and from here it follows that $C=\frac43$. The fact that $\Sigma(n)$ indeed has the asymptotic behavior of this type may be proved in the same way as it is done for the case $r=1/2$ in \cite{Mat}.

	Hence, we obtain that $\Sigma(n)$ has the asymptotic expansion of the following form:
	\begin{equation}	
		\Sigma(n)=\frac{4^{n+1}}{3\sqrt\pi}\sum_{k=0}^\infty S_k(\a)(n+\a)^{-k-r}.
	\end{equation}
	Then, using the asymptotic expansion \eqref{coef-a}, we get
	\begin{align*}
		\frac{4^{n+1}}{3\sqrt\pi}&amp;\sum_{k=0}^\infty S_k(\a)(n+\a)^{-k-r}\\
		&amp;=\frac{4^n}{\sqrt\pi}\sum_{k=0}^\infty P_k(\a)(n+\a)^{-k-r}+		\frac{4^{n}}{3\sqrt\pi}\sum_{k=0}^\infty S_k(\a)(n+\a-1)^{-k-r}\\
		&amp;=\frac{4^n}{\sqrt\pi}\sum_{k=0}^\infty P_k(\a)(n+\a)^{-k-r}\\
		&amp;\qquad+		\frac{4^{n}}{3\sqrt\pi}\sum_{k=0}^\infty S_k(\a)
		(n+\a)^{-k-r}\sum_{j=0}^\infty (-1)^j\binom{-k-r}{j} (n+\a)^{-j}.
	\end{align*}
	Hence
	\begin{equation*}
		4S_k(\a)=3P_k(\a)+\sum_{j=0}^k (-1)^{k-j}\binom{-j-r}{k-j}S_j(\a).
	\end{equation*}
	Extracting from the right side the member $S_k$,  we get \eqref{coefs}.
\end{proof} 

	Taking $\a=0$ and $r=1/2$ or $r=3/2$, it is easy to obtain the following asymptotics:
	\begin{align*}
		\sum_{k=0}^n\binom{2k}k&amp;\sim\frac{4^{n+1}}{3\sqrt{\pi n}}
		\biggl(1+\frac1{24n}+\frac{59}{384n^2}+\frac{2425}{9216n^3}
		+\frac{576793}{884736n^4}\\
		&amp;\qquad\qquad+\frac{5000317}{2359296n^5}+\frac{953111599}{113246208n^6}
		+\dots\biggr)\\
		\sum_{k=0}^nC_n&amp;\sim\frac{4^{n+1}}{3n\sqrt{\pi n}}
		\biggl(1-\frac5{8n}+\frac{475}{384n^2}+\frac{1225}{9216n^3}
		+\frac{395857}{98304n^4}\\
		&amp;\qquad\qquad
		+\frac{27786605}{2359296n^5}
		+\frac{6798801295}{113246208n^6}
		\biggr).
	\end{align*}
	
	Here, there is no good value of $\a$ which would lead to the expansion similar to \eqref {Asymptotic-Wn-alpha-14} or \eqref{Asymptotic-Wn-alpha-34c}, as it is evident from the formula \eqref{coefs}.


\begin{thebibliography}{99}

\bibitem{abram}
        M. Abramowitz and I. A. Stegun, eds.,
        {\it Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables}, 
        National Bureau of Standards, Applied Mathematics Series \textbf{55}, 9th printing, Washington, 1970.

\bibitem{erd}
        A. Erd\'elyi,
        {\it Asymptotic Expansions}, Dover Publications, New York, 1956.

\bibitem{luke}
        Y. L. Luke, 
        {\it The Special Functions and Their Approximations}, Vol.\ I, Academic Press, New York, 1969.

\bibitem{be1}
	T. Buri\'c and N. Elezovi\'c,  Bernoulli polynomials and
	asymptotic expansions of the quotient of gamma functions, \textit{J.
	Comput. Appl. Math.} {\bf 235} (2011),  3315--3331.

\bibitem{be2}
        T. Buri\'c and N. Elezovi\'c,
        New asymptotic expansions of the gamma function and improvements of Stirling's type formulas, \textit{J. Comput. Anal. Appl.} \textbf{13} (2011), 785--795.

\bibitem{be3}
	T. Buri\'c and N. Elezovi\'c,  New asymptotic expansions of
	the quotient of gamma functions, \textit{Integral Transform. Spec. Funct.}
	{\bf 23} (2012), 355--368.

\bibitem{BES}
    T. Buri\'c, N. Elezovi\'c, and R. \v Simi\'c,
    Asymptotic expansions of
    the multiple quotients of gamma functions and applications, 
    \textit{Math. Inequal. Appl.}, to appear.


\bibitem{EGP}
N. Elezovi\'c, C. Giordano, and J. Pe\v{c}ari\'c, 
The best bounds in Gautschi's inequality, \textit{Math. Inequal. Appl.} \textbf{3} (2000), 239--252.

\bibitem{luschny}
P. Luschny, Approximation formulas for the factorial function $n!$,\newline \url{http://www.luschny.de/math/factorial/approx/SimpleCases.html}.

\bibitem{KS}
D. Kessler and J. Schiff,
The asymptotics of factorials, binomial 
coefficients and Catalan numbers, manuscript,	
\url{http://u.math.biu.ac.il/~schiff/Papers/prepap3.pdf}.

\bibitem{Mat}
S. Mattarei, Asymptotics of partial sums of central binomial
coefficients and Catalan numbers, preprint,
\url{http://arXiv.math.CO/0906.4290}.

\bibitem{Stan}
	Richard P. Stanley,
	{\it Enumerative Combinatorics}, Vol.\ 2, Cambridge University Press, 1999.

\bibitem{tri-erd}
        F. G. Tricomi and A. Erd\'elyi,
        The asymptotic expansion of a ratio of Gamma functions,
        {\it Pacific J. Math.}, {\bf 1} (1951), 133--142.




\end{thebibliography}

\bigskip
\hrule
\bigskip

\noindent 2010 {\it Mathematics Subject Classification}:
Primary 11B65; Secondary 41A60, 33B15.

\noindent \emph{Keywords: }
binomial coefficient; Catalan number; asymptotic expansion; gamma function.

\bigskip
\hrule
\bigskip

\noindent (Concerned with sequences
\seqnum{A001163},
\seqnum{A001164},
\seqnum{A046968}, and
\seqnum{A046969}.)

\bigskip
\hrule
\bigskip

\vspace*{+.1in}
\noindent
Received March 21 2013;
revised versions received  June 6 2013; October 14 2013;
November 26 2013.
Published in {\it Journal of Integer Sequences},
January 3 2014.

\bigskip
\hrule
\bigskip

\noindent
Return to
\htmladdnormallink{Journal of Integer Sequences home page}{http://www.cs.uwaterloo.ca/journals/JIS/}.
\vskip .1in


\end{document}

                                                                                
</pre></body></html>