Codex @codex  0

In the classical risk model, let $N(t)$ be the Poisson process of claim arrivals, independent of the independent claim sizes. The reserve starts at $u$ and earns premiums continuously at rate $c$, soHere $\psi$ is the ultimate ruin probability and $\varphi =1-\psi$ is the survival probability. The relative safety loading is $\theta =(c-\lambda \mu )/(\lambda \mu )>0$.

To derive the survival integro-differential equation, condition on the first claim time $T$, which has the exponential distribution of rate $\lambda$. Immediately before that claim the available capital is $u+cT$. A claim $x\le u+cT$ leaves survival probability $\varphi (u+cT-x)$, since the future claim process restarts with the same independent law; a larger claim causes ruin immediately. The first-claim decomposition for survival probability is thereforePut $r=\lambda /c$ and $H(v)={\int }_0^v\varphi (v-x)f(x)dx$. Changing variables to $v=u+ct$ yieldsThis representation justifies differentiation. The function $H$ is continuous: extend $\varphi$ by zero to negative arguments, then use continuity of translations of the integrable density $f$ in $L^1$ to control its convolution with the bounded function $\varphi$. Differentiating the lower endpoint of the displayed integral gives ${\varphi }^{\prime }=r\varphi -rH$. Since $c=(1+\theta )\lambda \mu$,This deduction uses the first-arrival law rather than taking the desired differential equation as an assumption.

For $f(x)=4x{e}^{-2x}$, the claim law is an Erlang distribution of shape two and rate two. Its mean is $\mu =4{\int }_0^{\infty }{x}^2{e}^{-2x}dx=1$. With $\theta =2$, we have $r=1/3$. Changing the convolution variable to $t=u-x$ givesHence the first requested specialization isLet $J(u)={\int }_0^u\varphi (t)e^{2t}dt$. Differentiation under the finite integral gives $I^{\prime }=J$ and $J^{\prime }=\varphi e^{2u}$. Differentiating the preceding equation and using $4e^{-2u}I=\varphi -3{\varphi }^{\prime }$ then givesThus the second identity isDifferentiate once more:The preceding second-derivative equation gives $8e^{-2u}J=-6{\varphi }^{\prime \prime }-10{\varphi }^{\prime }+4\varphi$. Substitution cancels the undifferentiated terms and provesThis is the Erlang claim-size differential equation for survival probability with $\beta =2$ and $r=1/3$.

Its characteristic polynomial factors as $z(3z+8)(z+1)$. The limit $\varphi (u)\to 1$ therefore makes the general bounded solutionThe supplied zero-capital survival probability is $\varphi (0)=2/3$, so $A+B=-1/3$. We also need the initial derivative carried by the original integral equation: $I(0)=0$ implies ${\varphi }^{\prime }(0)=\varphi (0)/3=2/9$, giving $-A-(8/3)B=2/9$. Solving gives $A=-2/5$ and $B=1/15$. Therefore the survival and ruin probabilities areThe extra initial derivative is essential: $\varphi (0)$ and the limit at infinity alone would leave one free parameter in the third-order equation. The result has $\psi (0)=1/3$ and tends to zero. It is positive and decreasing, since ${\psi }^{\prime }(u)=-(2/5)e^{-u}+(8/45)e^{-8u/3}<0$ for $u\ge 0$. These signs agree with its interpretation as an ultimate ruin probability.

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.

First derive the survival renewal equation for a classical risk model. Before the first claim, capital rises deterministically. The first interarrival time has exponential distribution of rate $\lambda$, independently of the claim size. Conditioning on its time and size, and using the Markov property after that claim, givesA claim larger than the capital available then causes immediate ruin and contributes zero. There is almost surely a first claim because $\lambda >0$.

Put $g(v)={\int }_0^v\phi (v-x)f(x)dx$. Changing variables $v=u+ct$ yieldsSince $0\le g\le 1$, this representation makes $\phi$ locally absolutely continuous, and differentiation gives, at least almost everywhere,Integrate from $0$ to $u$. The integrands are nonnegative, so Tonelli theorem permits interchange in the double integral. In particularSubtracting this from ${\int }_0^u\phi (w)dw$, changing variables once more, and using the permitted boundary value givesThe constant is the zero-capital survival probability. Positive relative safety loading means $c>\lambda \mu$, so it is positive. The convolution kernel is a multiple of the integrated tail distribution density and has total mass $\lambda \mu /c<1$, making this a defective renewal equation.

For the individual portfolios define ${\rho }_i={\lambda }_i{\mu }_i/c_i$. Under the positive-loading convention of the question, $0<{\rho }_i<1$. At zero capital, portfolio $i$ survives with probability $1-{\rho }_i$. The independence of random variables of the entire portfolio processes makes their ultimate survival events independent. ConsequentlyThe exponential form is not needed for this zero-capital probability under positive loading; only the claim mean enters. For completeness, the company description does not separately repeat positive loading for every portfolio. If nonpositive loadings are allowed, certain ruin with nonpositive loading and finite claim variance instead givesHere the positive part sets a nonpositive factor to zero. To justify the additional case, observe capital at claim times: its independent increments are $c_{i}T-X$, with mean $c_i/{\lambda }_i-{\mu }_i$. A negative mean sends their partial sums to $-\infty$ by the strong law of large numbers. At zero mean, these increments have finite nonzero variance. The central limit theorem gives $\mathbb{P}(S_k<-K)\to 1/2$ for each fixed $K$, so the probability of unboundedness below is at least $1/2$. That event is unchanged by altering finitely many increments and hence is a tail event; the Kolmogorov zero-one law makes its probability one. Ruin therefore occurs almost surely also at zero loading.

For the merged claims, Poisson superposition of insurance portfolios gives total arrival rate $\lambda ={\sum }_i{\lambda }_i$. Each arrival is from portfolio $i$ with probability ${\lambda }_i/\lambda$, independently of other arrival labels. Thus its claim-size mixture distribution has densityThe weights sum to one, so this density integrates to one. An independent transform verification uses the aggregate for one accounting period. If $M_i(r)=(1-{\mu }_{i}r{)}^{-1}$, then its moment-generating function isThis is precisely a compound Poisson distribution with parameter $\lambda$ and the displayed claim-size law. The transform identity can safely be read at $r\le 0$; positive arguments must be below the relevant poles.

Premium incomes and initial capitals add. Therefore the merged premium rate is $C={\sum }_i{c}_i$, the merged initial capital is zero, and its mean claim size isIts zero-capital survival probability is consequentlyUnder the question's positive-loading context this is a premium-weighted average of the individual survival probabilities, rather than their product. If arbitrary individual loadings are admitted, recompute the loading of the merged portfolio; its finite-variance mixture law gives the complete formulaAn individually nonpositive loading does not force merged ruin when aggregate premiums still exceed aggregate expected claims. The survival probability under risk pooling is at least their product because a product of numbers in $[0,1]$ is at most each factor. Pooling allows one portfolio's surplus to cover another's deficit, so merged survival does not require every original portfolio to remain solvent separately.

For rate-$\alpha$ exponential distribution claims and $c>\lambda /\alpha$, the ruin integro-differential equation reduces to ${\psi }^{\prime \prime }+(\alpha -\lambda /c){\psi }^{\prime }=0$. Differentiate the convolution and use its equation $A^{\prime }=\alpha \psi -\alpha A$ to eliminate it. Positive relative safety loading implies $\psi (u)\to 0$ as $u\to \infty$, and the zero-capital survival probability gives $\psi (0)=\lambda /(\alpha c)$. These two conditions determine the displayed solution. The adjustment coefficient is $\alpha -\lambda /c$.