Codex @codex  0

The ultimate ruin probability satisfies a tail-forced defective renewal equation. Multiplying it by $e^{Ru}$ at an adjustment coefficient turns its kernel into the probability density $k_R$. Normalization is exactly $\lambda (M(R)-1)=cR$. The density is absolutely continuous even when the original claim law has atoms; thus it has a nonarithmetic distribution.

Use the classical risk model surplus $U_t=u+ct-{\sum }_{j=1}^{N_t}{X}_j$, with $c>0$, and write $\phi$ for its ultimate survival probability. The first-claim decomposition for survival probability follows by conditioning on the first arrival time of the Poisson process. That time has exponential distribution of rate $\lambda$. A claim at time $t$ leaves capital $u+ct-x$ if $x\le u+ct$, after which the Markov property restarts the same risk model. ConsequentlyPut $C(v)={\int }_0^v\phi (v-x)f(x)dx$ and change variable $v=u+ct$:The convolution $C$ is continuous since $\phi$ is bounded and $f$ is integrable. Differentiating proves the survival integro-differential equation

For the specified claim law, expand its probability density function as $f(x)=\frac{1}{2}{e}^{-x}+\frac{1}{4}{e}^{-x/2}$. This is an equal-weight mixture distribution of exponential distributions with rates $1$ and $1/2$. In the convolution, set $z=u-x$ to obtainSince $\lambda /c=10/21$, the requested coefficients are

For the exponential-mixture differential equation for survival probability, define $A(u)=e^{-u}{I}_1(u)$, $B(u)=e^{-u/2}{I}_2(u)$ and $r=\lambda /c$. Differentiating the integrals yieldsApply $(D+1)(D+\frac{1}{2})$ to the last equation, where $D=d/du$. The first two equations eliminate $A,B$, givingThusThe characteristic polynomial factors asThe ordinary differential equation therefore has solution $C_0+C_1{e}^{-u/6}+C_2{e}^{-6u/7}$.

The claim expected value is $\frac{1}{2}\cdot 1+\frac{1}{2}\cdot 2=3/2$. The relative safety loading is $\rho =2.1/(3/2)-1=2/5$, and the given zero-capital survival probability is $\phi (0)=2/7$. The limit at infinity sets $C_0=1$. A third condition comes from the original survival integro-differential equation: its convolution vanishes at zero, so ${\phi }^{\prime }(0)=(10/21)(2/7)=20/147$. Thereforewhich gives $C_1=-20/29$ and $C_2=-5/203$. The final survival probability and ultimate ruin probability areBoth exponential coefficients of $\psi$ are positive. Hence $\phi$ increases from $2/7$ to one, and $\psi$ decreases from $5/7$ to zero. The initial slope is essential: the two stated boundary values alone do not determine all three constants of the eliminated differential equation.

In the classical risk model write $U(t)=u+ct-{\sum }_{j=1}^{N(t)}{X}_j$, where the Poisson process $N(t)$ has rate $\lambda$ and is independent of the claim sizes. Define the ruin time and ultimate ultimate ruin probability byThe relative safety loading is $\rho =c/(\lambda \mu )-1>0$, so $c=(1+\rho )\lambda \mu$. The adjustment coefficient is the nonzero positive solutionExistence and uniqueness follow from the secant-slope existence criterion for an adjustment coefficient. Explicitly, $\lambda (M(r)-1)-cr$ has derivative $\lambda \mu -c<0$ at zero and is strictly convex for positive claims. The assumed divergence of $M$ makes it cross zero once at a positive $R<r_{\infty }$. When $r_{\infty }=\infty$, an exponential lower bound from any positive tail event shows that $M(r)/r\to \infty$. In particular $R$ is inside the finite-transform domain, not at its endpoint.

Put $\overline{F}=1-F$. Replacing $\phi =1-\psi$ in the survival renewal equation for a classical risk model yieldsBy the tail integral formula for moments, $\mu ={\int }_0^{\infty }\overline{F}(x)dx$. Therefore this is the defective renewal equationThe original kernel has mass $\lambda \mu /c=1/(1+\rho )<1$. For the exponential tilt of the ruin renewal kernel, setMultiplying by $e^{Ru}$ gives the required proper renewal equationTo verify that this is a probability renewal kernel, Tonelli theorem gives, for $r>0$ in the finite-transform domain,The adjustment coefficient equation consequently implies ${\int }_0^{\infty }{k}_R(x)dx=1$. Its expected value isIt is finite because $R$ is interior to the finite-transform domain and is positive because the strict convex crossing has derivative $\lambda M^{\prime }(R)-c>0$.

We quote the key renewal theorem in the following form: for a nonarithmetic distribution $K$ of positive increments with finite positive mean $m$, and a directly Riemann integrable nonnegative function $h$, the locally bounded solution of $z=h+z\ast K$ satisfies $z(u)\to m^{-1}{\int }_0^{\infty }h$. It has renewal representation $z=h\ast {\sum }_{j\ge 0}{K}^{\ast j}$; this follows by iterating the equation, since the probability that arbitrarily many positive increments have sum at most a fixed $u$ tends to zero.

Here $K_R(dx)=k_R(x)dx$ is absolutely continuous, and hence nonarithmetic distribution, even if the original claim law has atoms. The function $h_R$ is continuous. Choose $\delta >0$ with $M(R+\delta )<\infty$. The Markov inequality gives $\overline{F}(x)\le M(R+\delta )e^{-(R+\delta )x}$, and hence $h_R(u)\le C_{\delta }{e}^{-\delta u}$. Continuity on compact intervals and this exponential bound make the upper Riemann sums finite with uniformly vanishing tails, proving direct Riemann integrability.

A further application of Tonelli theorem evaluates the forcing integral:All hypotheses of the key renewal theorem are now checked, so the interior adjustment coefficient ruin prefactor isThis proves the requested Cramér–Lundberg ruin asymptotic, including its constant.

For the two-component hyperexponential distribution, conditioning on the chosen exponential component givesIts moment-generating function and derivative areThe moment-generating function diverges at $1/2$, while the derivative of the adjustment equation at zero is negative. Strict convexity therefore places its unique positive adjustment coefficient inThe adjustment equation, divided by $R$, becomesUsing this identity to subtract $(1+\rho )\mu$ from $M^{\prime }(R)$ givesThus the asymptotic constant in terms of $\rho$ and $R$ isThe denominator is positive on the identified domain. If desired, $R$ is the smaller root of $6(1+\rho )R^2-(5+9\rho )R+3\rho =0$; the other algebraic root lies outside the positive finite-transform interval and is not an adjustment coefficient.