Solution (source code)

= Solution

Write the surplus in the <classical risk model> as $U_t=u+ct-\sum_{j=1}^{N_t}X_j$, with premium rate $c=(1+\rho)\lambda\mu$. The <adjustment coefficient> is the positive solution of
$$
\lambda\bigl(M_X(R)-1\bigr)=cR.
$$
For the <exponential distribution> of mean $\mu$, $M_X(r)=(1-\mu r)^{-1}$ for $r<1/\mu$. Cancelling the nonzero root $R$ in its equation gives $\lambda\mu/(1-\mu R)=c$, and hence
$$
\boxed{R=\frac{\rho}{(1+\rho)\mu}.}
$$
It is strictly between zero and $1/\mu$, as required by the transform domain.

The <ultimate ruin probability> is $\psi(u)=\mathbb P(\inf_{t\ge0}U_t<0)$. The <Lundberg inequality> states
$$
\boxed{\psi(u)\le e^{-Ru}.}
$$
Ruin can occur only at a claim arrival. Thus the event defining the <finite-claim ruin probability> $\psi_n(u)$ is contained in the ultimate ruin event, giving $\psi_n(u)\le\psi(u)\le e^{-Ru}$ for every $n\ge1$ and $u>0$.

For the explicit calculations, take $\mu=1$, put $d=1+\rho=c/\lambda$ and $h=d+1=2+\rho$, and let $T$ be the first arrival time. It has rate-$\lambda$ <exponential distribution>, independent of the severity. Ruin on the first claim is the event $X_1>u+cT$, so
$$
\boxed{\psi_1(u)=\mathbb E[e^{-(u+cT)}]=e^{-u}\frac{\lambda}{\lambda+c}=\frac{e^{-u}}h.}
$$
For ruin by the second claim, condition on whether the first claim ruins the insurer or leaves capital $u+cT-x$. The independent future arrivals and severities have the same law as a fresh <classical risk model>. Therefore
$$
\begin{aligned}
\psi_2(u)&=\psi_1(u)+\int_0^\infty\lambda e^{-\lambda t}\int_0^{u+ct}e^{-x}\psi_1(u+ct-x)\,dx\,dt\\
&=\frac{e^{-u}}h+\frac{e^{-u}}h\int_0^\infty\lambda e^{-(\lambda+c)t}(u+ct)\,dt\\
&=\boxed{e^{-u}\left(\frac1h+\frac{u}{h^2}+\frac{d}{h^3}\right).}
\end{aligned}
$$
The two terms in the first line refer to disjoint ruin events, so there is no double counting.

For a direct check of the <Lundberg inequality>, now $R=1-1/d$. The ratio $\psi_1(u)/e^{-Ru}=h^{-1}e^{-u/d}$ is less than one. For the second <finite-claim ruin probability>, put $A=h^{-1}+dh^{-3}$ and $B=h^{-2}$. Its ratio is $e^{-u/d}(A+Bu)$. Since
$$
A-dB=\frac{2d+1}{h^3}>0,\qquad 1-A=\frac{d(h^2-1)}{h^3}>0,
$$
the derivative of this ratio is $e^{-u/d}[B-(A+Bu)/d]<0$ for $u\ge0$, and its value at zero is $A<1$. This verifies both required bounds directly from the calculated probabilities.