\documentclass[12pt,reqno]{article}

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

\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}{-.5in}
\setlength{\textheight}{8.9in}

\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}
\newtheorem{result}[theorem]{Result}

\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 
Remodified Bessel Functions via\\ \vskip .1in
Coincidences and Near Coincidences
}
\vskip 1cm
\large
Martin Griffiths\\
School of Education\\
University of Manchester\\
M13 9PL\\
United Kingdom\\
\url{martin.griffiths@manchester.ac.uk}\\
\end{center}

\vskip .2 in

\begin{abstract}
By considering a particular probabilistic scenario associated with
coincidences, we are led to a family of functions akin to the modified
Bessel function of the first kind.  These are in turn solutions to a
certain family of linear differential equations possessing structural
similarities to the modified Bessel differential equation.  The
Stirling number triangle of the second kind arises quite naturally from
these differential equations, as do more complicated, yet related,
truncated number triangles, none of which appear in Sloane's
\textit{On-Line Encyclopedia of Integer Sequences}.
\end{abstract}

\section{Introduction}
Let $X\sim\textup{Po}(\lambda)$ denote a discrete random variable having the Poisson distribution \cite{grimmett} with parameter $\lambda$.  The mass function of $X$ is given by
\[\textup{P}(X=k)=\frac{e^{-\lambda}\lambda^k}{k!}, k\geq 0.\]
Furthermore, suppose that, for some $n\in\mathbb{N}$, the random variables $X_1,X_2,\ldots,X_n$ are independently and identically distributed as $X$.  We shall term the event
\[C_n=\bigcup_{k=0}^{\infty}\{X_1=X_2=\cdots=X_n=k\}\]
a \textit{coincidence}.  In order to place such an event into some of context, we might imagine $n$ call centers each receiving an average of $\lambda$ telephone calls per hour.  In this situation a coincidence is said to occur when, in a given hour, all of the centers receive exactly the same number of calls.

Next, let $r\in\{1,2,\ldots,n\}$ and $S_r$ denote the set $\{1,2,\ldots,n\}\backslash\{r\}$.  We define the event $A_r$ by
\[A_r=\{X_{a_1}=X_{a_2}=\cdots=X_{a_{n-1}}=k\},\]
where $\{a_1,a_2,\ldots,a_{n-1}\}=S_r$.  The event
\[N_n=\bigcup_{k=0}^{\infty}\bigcup_{r=1}^{n}\left\{A_r\cap\{X_r=k+1\}\right\}\]
is known as a \textit{near coincidence}.  Returning to the call-center scenario, a near coincidence is said to occur when, in a given hour, all but one of the centers receive exactly the same number of calls, with the remaining center receiving exactly one more call than all of the others.

In this paper we show how probabilities associated with coincidences, near coincidences and beyond, give rise to functions which may be regarded as extended versions of certain Bessel functions.  Via the linear differential equations these functions satisfy, this leads first to the triangle of Stirling numbers of the second kind and then on to rather more complicated, yet related, truncated number triangles.

\section{A connection with Bessel functions and Stirling numbers of the second kind}
Bessel functions arise as solutions to certain linear differential equations. They come in several varieties, and we will be concerned here with a particular Bessel function that appears, amongst other things, in connection with special relativity \cite{lavenda1, lavenda2} and the Skellam distribution \cite{wiki}.

\begin{definition}
The \textit{modified Bessel function of the first kind} \cite{weisstein2},
\[I_m(x)=\sum_{k=0}^{\infty}\frac{x^{2k+m}}{2^{2k+m}k!\Gamma(k+m+1)},\]
is one of the solutions to the \textit{modified Bessel differential equation} \cite{weisstein1} given by
\[x^2y''+xy'-\left(m^2+x^2\right)y=0,\]
where $\Gamma(x)$ is the \textit{gamma function} \cite{griffiths, knuth}.
\end{definition}

The probability of a coincidence occurring in any given hour is given by
\begin{align*}
\textup{P}(C_n)
&=\sum_{k=0}^{\infty}\left(\frac{e^{-\lambda}\lambda^k}{k!}\right)^n\\
&=e^{-n\lambda}\sum_{k=0}^{\infty}\frac{\lambda^{nk}}{(k!)^n}.
\end{align*}
Then, noting that
\begin{align*}
e^{-2\lambda}I_0(2\lambda)
&=e^{-2\lambda}\sum_{k=0}^{\infty}\frac{(2\lambda)^{2k}}{2^{2k}k!\Gamma(k+1)}\\
&=e^{-2\lambda}\sum_{k=0}^{\infty}\frac{\lambda^{2k}}{(k!)^2}\\
&=\textup{P}(C_2),
\end{align*}
we are led first to extend the definition of $I_0(x)$ as follows:
\begin{definition}
Let $n\in\mathbb{N}$.  Then
\[I_0(n,x)=\sum_{k=0}^{\infty}\frac{x^{nk}}{n^{nk}(k!)^n}.\]
\end{definition}
The function $I_0(n,x)$ is related to the probability $\textup{P}(C_n)$ by way of
\begin{align*}
\textup{P}(C_n)
&=e^{-n\lambda}\sum_{k=0}^{\infty}\frac{\lambda^{nk}}{(k!)^n}\\
&=e^{-n\lambda}I_0(n,n\lambda).
\end{align*}
 It is clear that $I_0(1,x)=e^x$ and $I_0(2,x)=I_0(x)$.  Note also that $y=I_0(1,x)$ and $y=I_0(2,x)$ satisfy $y'-y=0$ and $xy''+y'-xy=0$, respectively.  We now find a third-order linear differential equation having $I_0(3,x)$ as a solution.

\begin{result}\label{simplediffeq}
The function $y=I_0(3,x)$ is a solution to
\[x^2y'''+3xy''+y'-x^2y=0.\]
\end{result}

\begin{proof}
Let $y=I_0(3,x)$.  We have
\[y'=\sum_{k=1}^{\infty}\frac{x^{3k-1}}{3^{3k-1}(k!)^2(k-1)!}\]
and
\[y''=\sum_{k=1}^{\infty}\frac{(3k-1)x^{3k-2}}{3^{3k-1}(k!)^2(k-1)!},\]
so that
\[x(y'+xy'')=\sum_{k=1}^{\infty}\frac{x^{3k}}{3^{3k-2}k!((k-1)!)^2},\]
and hence
\begin{align*}
\left(x(y'+xy'')\right)'
&=\sum_{k=1}^{\infty}\frac{x^{3k-1}}{3^{3(k-1)}((k-1)!)^3}\\
&=\sum_{k=0}^{\infty}\frac{x^{3k+2}}{3^{3k}(k!)^3}\\
&=x^2y.
\end{align*}
It follows from this that $y=I_0(3,x)$ satisfies
\[x^2y'''+3xy''+y'-x^2y=0,\]
as required.
\end{proof}

Adopting a method similar to that used in Result \ref{simplediffeq}, we may show that $y=I_0(4,x)$ is a solution to
\[\left(x\left(x^2y'''+3xy''+y'\right)\right)'-x^3y=0,\]
and so on.  It is in fact the case that $I_0(n,x)$ satisfies the $n$th-order linear differential equation
\begin{equation}\label{stirlingdiffeq}
\sum_{k=1}^{n}S(n,k)x^{k-1}y^{(k)}-x^{n-1}y=0,
\end{equation}
where $S(n,k)$ is a Stirling number of the second kind, enumerating the partitions of $n$ distinct objects into exactly $k$ non-empty parts, and $y^{(k)}$ denotes the $k$th derivative of $y$ with respect to $x$.  The number triangle associated with $S(n,k)$ appears in Sloane's \textit{On-Line Encyclopedia of Integer Sequences} \cite{sloane} as sequence \seqnum{A008277}.  It is straightforward to prove (\ref{stirlingdiffeq}) by using induction in conjunction with the well-known result
\begin{equation}\label{stirlingrec}
S(n,k)=kS(n-1,k)+S(n-1,k-1),
\end{equation}
which may be found in \cite{cameron} and \cite{knuth}.

\section{Near coincidences}
Let us now consider the probability of the occurrence of a near coincidence, assuming that $n\geq 2$.  We have
\begin{align*}
\textup{P}(N_n)
&={n \choose 1} \sum_{k=0}^{\infty}\frac{e^{-\lambda}\lambda^{k+1}}{(k+1)!}\left(\frac{e^{-\lambda}\lambda^k}{k!}\right)^{n-1}\\
&=ne^{-n\lambda}\sum_{k=0}^{\infty}\frac{\lambda^{nk+1}}{(k!)^n(k+1)}.
\end{align*}
This leads us to generalize $I_0(n,x)$ as follows:
\begin{definition}\label{remodified}
For $n\geq m+1$,
\[I_m(n,x)=\sum_{k=0}^{\infty}\frac{x^{nk+m}}{n^{nk+m}(k!)^n(k+1)^m}.\]
\end{definition}
\noindent Note then that $\textup{P}(N_n)=ne^{-n\lambda}I_1(n,n\lambda)$.  In this section we will indeed consider the special case $m=1$, which is the one associated with near coincidences.

\begin{result}
The function $y=I_1(2,x)$ satisfies the differential equation
\[x^3y'''+2x^2y''-xy'+y-x^2\left(xy'+y\right)=0.\]
\end{result}

\begin{proof}
First,
\begin{align*}
\left(xI_1(2,x)\right)'
&=\sum_{k=0}^{\infty}\frac{2(k+1)x^{2k+1}}{2^{2k+1}(k!)^2(k+1)}\\
&=\sum_{k=0}^{\infty}\frac{x^{2k+1}}{2^{2k}(k!)^2}\\
&=xI_0(2,x).
\end{align*}
Therefore
\begin{align*}
I_0(2,x)
&=\frac{1}{x}\left(xI'_1(2,x)+I_1(2,x)\right)\\
&=I'_1(2,x)+\frac{1}{x}I_1(2,x),
\end{align*}
and so, since $y=I_0(2,x)$ is a solution to $xy''+y'-xy=0$, we may obtain
\[x\left(I'_1(2,x)+\frac{1}{x}I_1(2,x)\right)''+\left(I'_1(2,x)+\frac{1}{x}I_1(2,x)\right)'
-x\left(I'_1(2,x)+\frac{1}{x}I_1(2,x)\right)=0.\]
From this it follows that $y=I_1(2,x)$ does in fact satisfy
\[x^3y'''+2x^2y''-xy'+y-x^2\left(xy'+y\right)=0.\]
\end{proof}

Taking things further,
\begin{align*}
\left(x^2I_1(3,x)\right)'
&=\sum_{k=0}^{\infty}\frac{3(k+1)x^{3k+2}}{3^{3k+1}(k!)^3(k+1)}\\
&=\sum_{k=0}^{\infty}\frac{x^{3k+2}}{3^{3k}(k!)^3}\\
&=x^2I_0(3,x),
\end{align*}
from which we may obtain, using Result \ref{simplediffeq}, that $y=I_1(3,x)$ is a solution to
\[x^4y''''+5x^3y'''+x^2y''+2xy'-2y-x^3\left(xy'+2y\right)=0.\]
More generally,
\begin{align*}
\left(x^{n-1}I_1(n,x)\right)'
&=\sum_{k=0}^{\infty}\frac{n(k+1)x^{nk+n-1}}{n^{nk+1}(k!)^n(k+1)}\\
&=\sum_{k=0}^{\infty}\frac{x^{nk+n-1}}{n^{nk}(k!)^n}\\
&=x^{n-1}I_0(n,x),
\end{align*}
giving
\begin{equation}\label{indstart1}
I_0(n,x)=I'_1(n,x)+\frac{n-1}{x}I_1(n,x).
\end{equation}

\begin{result}\label{nearresult}
The function $y=I_1(n,x)$ is a solution to
\begin{equation}\label{fdiff}
\sum_{k=1}^n f(n,k) x^{k+1} y^{(k+1)}+(n-1)(-1)^{n+1}\left(xy'-y\right)-x^n\left(xy'+(n-1)y\right)=0,
\end{equation}
where
\begin{equation}\label{ftriang}
f(n,k)=S(n,k)+\frac{n-1}{(k+1)!}\sum_{j=k+1}^n (-1)^{j-k+1}j!S(n,j),
\end{equation}
noting that the sum on the far right is defined to be zero when $k\geq n$.
\end{result}

\begin{proof}
Starting with (\ref{indstart1}) and proceeding by induction gives
\[I_0^{(k)}(n,x)=I_1^{(k+1)}(n,x)+k!(n-1)\sum_{j=0}^k\frac{(-1)^{k-j}I_1^{(j)}(n,x)}{x^{k-j+1}j!}.\]
From this we obtain, using (\ref{stirlingdiffeq}) and induction once more, the general result that $I_1(n,x)$ satisfies the differential equation
\begin{align}\label{gendiffeq1}
\sum_{k=1}^{n} f(n,k) x^{k+1} y^{(k+1)}+(n-1)\sum_{k=1}^{n}(-1)^{k+1} k! & S(n,k)\left(xy'-y\right)\nonumber\\
&-x^n\left(xy'+(n-1)y\right)=0,
\end{align}
where $f(n,k)$ is as given in the statement of the result.

Since, by definition,
\[\sum_{k=1}^n S(n,k)(x)_k=x^n,\]
where $(x)_k=x(x-1)(x-2)\cdots(x-k+1)$ denotes the \textit{falling factorial}, on setting $x=-1$ it follows that
\begin{align*}
(-1)^n
&=\sum_{k=1}^n S(n,k)(-1)_k\\
&=\sum_{k=1}^n S(n,k)k!(-1)^k.
\end{align*}
Therefore (\ref{gendiffeq1}) may be simplified somewhat to give the desired result.
\end{proof}

Since in Definition \ref{remodified} we require $n\geq m+1$, the number triangle for $f(n,k)$ is, unlike that for $S(n,k)$, truncated.  Its first few rows may be seen in Table \ref{Table1} of Section \ref{sectiontables}.

\section{Further coincidences}
Next, let $r,s\in\{1,2,\ldots,n\}$ such that $r\neq s$, and let $Q_{r,s}$ denote the set $\{1,2,\ldots,n\}\backslash\{r,s\}$.  We define $B_{r,s}$ by
\[B_{r,s}=\{X_{a_1}=X_{a_2}=\cdots=X_{a_{n-2}}=k\},\]
where $\{a_1,a_2,\ldots,a_{n-2}\}=Q_{r,s}$, and consider the event
\[M_n=\bigcup_{k=0}^{\infty}\bigcup_{r,s}\left\{B_{r,s}\cap\{X_r=X_s=k+1\}\right\},\]
where the inner union is over all possible pairs $(r,s)$ such that $r,s\in\{1,2,\ldots,n\}$ and $r\neq s$.  It is clear that
\begin{align*}
\textup{P}(M_n)
&={n \choose 2} \sum_{k=0}^{\infty}\left(\frac{e^{-\lambda}\lambda^{k+1}}{(k+1)!}\right)^2
\left(\frac{e^{-\lambda}\lambda^k}{k!}\right)^{n-2}\\
&=e^{-n\lambda}{n \choose 2}\sum_{k=0}^{\infty}\frac{\lambda^{nk+2}}{(k!)^n(k+1)^2}\\
&=e^{-n\lambda}{n \choose 2}I_2(n,n\lambda).
\end{align*}
For this new scenario we now obtain, in correspondence to Result \ref{nearresult}, a family of linear differential equations and the associated truncated number triangle.

First,
\begin{align*}
\left(x\left(x^{n-2}I_2(n,x)\right)'\right)'
&=\left(x\sum_{k=0}^{\infty}\frac{n(k+1)x^{nk+n-1}}{n^{nk+2}(k!)^n(k+1)^2}\right)'\\
&=\left(\sum_{k=0}^{\infty}\frac{x^{nk+n}}{n^{nk+1}(k!)^n(k+1)}\right)'\\
&=\sum_{k=0}^{\infty}\frac{n(k+1)x^{nk+n-1}}{n^{nk+1}(k!)^n(k+1)}\\
&=\sum_{k=0}^{\infty}\frac{x^{nk+n-1}}{n^{nk}(k!)^n}\\
&=x^{n-1}I_0(n,x),
\end{align*}
from which we have
\begin{equation}\label{indstart2}
I_0(n,x)=I''_2(n,x)+\frac{2n-3}{x}I'_2(n,x)+\frac{(n-2)^2}{x^2}I_2(n,x).
\end{equation}
It follows from this that
\begin{align*}
I_0^{(k)}(n,x)=&I_2^{(k+2)}(n,x)+k!(2n-3)\sum_{j=0}^k\frac{(-1)^{k-j}I_2^{(j+1)}(n,x)}{x^{k-j+1}j!}\\
&\qquad\qquad\qquad+k!(n-2)^2\sum_{j=0}^k\frac{(-1)^{k-j}(k-j+1)I_2^{(j)}(n,x)}{x^{k-j+2}j!}.
\end{align*}

Then, to obtain a linear differential equation satisfied by $I_2(n,x)$, we may use the fact that $I_0(n,x)$ is a solution to (\ref{stirlingdiffeq}) to give
\begin{align}\label{gdiff}
\sum_{k=1}^n g(n,k)& x^{k+2} y^{(k+2)}+(n-2)^2\sum_{k=1}^n k!S(n,k)\sum_{j=0}^2\frac{(-1)^{k+j}}{j!}x^j y^{(j)}(k-j+1)\nonumber\\
&+(2n-3)(-1)^{n+1}\left(x^2y''-xy'\right)-x^n\left(x^2y''+(2n-3)xy'+(n-2)^2y\right)=0,
\end{align}
where
\begin{align}\label{gtriang}
g(n,k)=&S(n,k)+\frac{2n-3}{(k+1)!}\sum_{j=k+1}^n(-1)^{j-k-1}j!S(n,j)\nonumber\\
&\qquad\qquad+\frac{(n-2)^2}{(k+2)!}\sum_{j=k+2}^n(-1)^{j-k}j!S(n,j)(j-k-1),
\end{align}
and the first and second sums on the right are defined to be zero when $k\geq n$ and $k\geq n-1$, respectively.  The truncated number triangle for $g(n,k)$ is shown in Table \ref{Table2} of Section \ref{sectiontables}.

The sequence of numbers along the $m$th diagonal of the triangle of Stirling numbers of the second kind is given by $\{S(n+m,n): n=1,2,3,\ldots\}$, where the convention is that the upper-most diagonal, consisting solely of 1s, is the zeroth diagonal. The sequences for $m$ = 1, 2, 3 and 4 appear in \cite{sloane} as 
\seqnum{A000127}, \seqnum{A001296}, \seqnum{A001297},
and \seqnum{A001298}, respectively.  Using the recurrence relation
(\ref{stirlingrec}), it is possible to show, by induction, that for
each $m\in\mathbb{N}$ there exists a polynomial $p_m(x)$ such that
$S(n+m,n)=p_m(n)$, $n=1,2,3,\ldots$.  It is in fact reasonably
straightforward to show that $p_m(x)$ has degree $2m$ and leading
coefficient $\frac{1}{2^m m!}$.  It follows from this, in conjunction
with (\ref{ftriang}), that the $m$th diagonal of the table for $f(n,k)$
is a polynomial sequence of degree $2m$.  On using (\ref{gtriang}), a
similar result applies to the table for $g(n,k)$.

It is possible to generalize the results (\ref{indstart1}) and (\ref{indstart2}).  We have
\[x^{n-1}I_0(n,x)=\left(x\cdots\left(x\left(x^{n-m}I_m(n,x)\right)'\right)'\cdots\right)',\]
where the nesting is to a depth of $m$.  It follows from this that
\begin{equation}\label{generalresult1}
x^{n-1}I_0(n,x)=\sum_{k=1}^m S(m,k)x^{k-1}\left(x^{n-m}I_m(n,x)\right)^{(k)}.
\end{equation}
It is also straightforward to show that
\begin{equation}\label{generalresult2}
\left(x^q h(x)\right)^{(k)}=\sum_{j=0}^k x^{q-j}(q)_j {k \choose j}h^{(k-j)}(x),
\end{equation}
where $(q)_0=1$ by definition.  From (\ref{generalresult1}) and (\ref{generalresult2}) we may obtain
\begin{align*}
I_0(n,x)
&=\frac{1}{x^{n-1}}\sum_{k=1}^m S(m,k)x^{k-1}\sum_{j=0}^k x^{n-m-j}(n-m)_j {k \choose j}I_m^{(k-j)}(n,x)\\
&=\sum_{k=1}^m S(m,k)\sum_{j=0}^k x^{k-j-m}(n-m)_j {k \choose j}I_m^{(k-j)}(n,x).
\end{align*}
This result, in conjunction with (\ref{stirlingdiffeq}), allows us to find a linear differential equation satisfied by  $y=I_m(n,x)$ for any $m\in\mathbb{N}$.

\section{Tables}\label{sectiontables}
\renewcommand\arraystretch{1.4}
\begin{table}[H]
\begin{center}
  \begin{tabular}{c c c c c c c c c}
  $n$ & $f(n,1)$ & $f(n,2)$ & $f(n,3)$ & $f(n,4)$ & $f(n,5)$ & $f(n,6)$ & $f(n,7)$ & $f(n,8)$\\
    \hline
    2 & 2 & 1 & & & & & & \\
    3 & 1 & 5 & 1 & & & & & \\
    4 & 4 & 13 & 9 & 1 &  &  & & \\
    5 & 1 & 35 & 45 & 14 & 1 & & & \\
    6 & 6 & 81 & 190 & 110 & 20 & 1 & & \\
    7 & 1 & 189 & 721 & 686 & 224 & 27 & 1 & \\
    8 & 8 & 421 & 2583 & 3759 & 1932 & 406 & 35 & 1
  \end{tabular}\caption{The coefficients $f(n,k)$ of the differential equation (\ref{fdiff}).}\label{Table1}
\end{center}
\end{table}

\renewcommand\arraystretch{1.4}
\begin{table}[H]
\begin{center}
  \begin{tabular}{c c c c c c c c c}
  $n$ & $g(n,1)$ & $g(n,2)$ & $g(n,3)$ & $g(n,4)$ & $g(n,5)$ & $g(n,6)$ & $g(n,7)$ & $g(n,8)$\\
    \hline
    3 & 2 & 6 & 1 & & & & & \\
    4 & $-2$ & 21 & 11 & 1 &  &  & & \\
    5 & 46 & 50 & 69 & 17 & 1 & & & \\
    6 & $-150$ & 201 & 318 & 162 & 24 & 1 & & \\
    7 & 526 & 294 & 1421 & 1141 & 319 & 32 & 1 & \\
    8 & $-1498$ & 1429 & 5481 & 7035 & 3120 & 562 & 41 & 1
  \end{tabular}\caption{The coefficients $g(n,k)$ of the differential equation (\ref{gdiff}).}\label{Table2}
\end{center}
\end{table}

\section{Acknowledgement}

The author is grateful to the anonymous referee for providing a number of constructive comments and suggestions as to how the original version of this article might be improved.



\begin{thebibliography}{99}

\bibitem{cameron} P. J. Cameron, \textit{Combinatorics: Topics, Techniques, Algorithms}, Cambridge University Press, 1994.

\bibitem{griffiths}  M. Griffiths, \textit{The Backbone of Pascal's Triangle}, United Kingdom Mathematics Trust, 2008.

\bibitem{grimmett}  G. Grimmett and D. Stirzaker, \textit{Probability and Random Processes}, 3rd ed., Oxford University Press, 2001.

\bibitem{lavenda1}  B. H. Lavenda, Is relativistic quantum mechanics compatible with special relativity? \textit{Zeitschrift f\"{u}r Naturforschung A} \textbf{56a} (2001), 347--365.

\bibitem{lavenda2}  B. H. Lavenda, Special relativity via modified Bessel functions, \textit{Zeitschrift f\"{u}r Naturforschung A} \textbf{55a} (2000), 745--753.

\bibitem{knuth} D. E. Knuth, 
\textit{The Art of Computer Programming Volume 1}, Addison-Wesley, 1968.

\bibitem{sloane} N. J. A. Sloane (Ed.), The On-Line Encyclopedia of Integer Sequences, 2011.  Available at \url{http://oeis.org/} .

\bibitem{weisstein1}  E. W. Weisstein, Modified Bessel differential equation,
MathWorld, 2011.  
Available at \\ \url{http://mathworld.wolfram.com/ModifiedBesselDifferentialEquation.html} .

\bibitem{weisstein2}  E. W. Weisstein, Modified Bessel function of the first 
kind, MathWorld, 2011.  
Available at \\ \url{http://mathworld.wolfram.com/ModifiedBesselFunctionoftheFirstKind.html} .

\bibitem{wiki}  {\it Skellam distribution}, Wikipedia, The Free Encyclopedia, 2011.  Available at \url{http://en.wikipedia.org/wiki/Skellam_distribution} .

\end{thebibliography}

\bigskip
\hrule
\bigskip

\noindent 2010 {\it Mathematics Subject Classification}:
Primary 33C10; Secondary 11B73, 33C90, 34A05, 60E05.

\noindent \emph{Keywords: } Bessel functions, linear differential
equations, Poisson distributions, Stirling numbers of the second kind.

\bigskip
\hrule
\bigskip

\noindent (Concerned with sequences
\seqnum{A000127},
\seqnum{A001296},
\seqnum{A001297},
\seqnum{A001298}, and
\seqnum{A008277}.)

\bigskip
\hrule
\bigskip

\vspace*{+.1in}
\noindent
Received January 28 2011;
revised version received May 29 2011. 
Published in {\it Journal of Integer Sequences}, June 11 2011.

\bigskip
\hrule
\bigskip

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

\end{document}

.
