Let $X$ and $Y$ be Poisson random variables with means $\lambda$ and $1$, respectively. The difference of $X$ and $Y$ is a Skellam random variable, with probability density function $$\mathbb P(X - Y = k) = \mathrm e^{-\lambda - 1} \lambda^{k/2} I_k(2\sqrt{\lambda}) =: S(\lambda, k),$$ where $I_k$ denotes the modified Bessel function of the first kind. Let $F(\lambda)$ denote the probability that $X$ is larger than $Y$: $$F(\lambda) := \mathbb P(X > Y) = \sum_{k=1}^{\infty} S(\lambda, k) = \mathrm e^{-\lambda - 1} \sum_{k=1}^\infty \lambda^{k/2} I_k(2\sqrt{\lambda}).$$ According to Mathematica, the graph of the function $F$ looks like

_{(source: nyu.edu)}

My question: