Method of successive approximations for solving integral equations of actuarial mathematics
Abstract
The dissertation considers various generalizations of the classical risk process describing the stochastic evolution of the capital of an insurance company (Cramer Lundberg model), including processes with variable deterministic or random premiums, non Poisson flows of premiums and claims, and a stochastic Markovian environment. Integral equations for the probability of nonruin as a function of initial capital are derived for these generalizations. For processes in a stochastic Markovian environment, systems of integral equations for nonruin probabilities corresponding to different initial states are obtained. General necessary and sufficient, as well as specific sufficient, conditions for the existence and uniqueness of solutions are established. A successive approximation method for numerical and analytical solution of the integral equations is theoretically and practically validated; its uniform convergence and rate of convergence are established. A technique for estimating the accuracy of approximate solutions is developed by constructing upper and lower approximations to the exact solution. The method is tested on numerical examples and compared with the Monte Carlo method and known solution approximations. The developed method increases the accuracy of actuarial calculations: it allows the probability of ruin to be calculated with any prescribed accuracy, empirical approximations to be verified and improved iteratively, and the accuracy and parameters of Monte Carlo simulations to be estimated and corrected.
Related
- Existence and uniqueness for non-linear forward equations: Addendum to "Mean-field approximations in insurance" (Hornung, 2026)
- Weak solutions to granular medium equations via stochastic particle approximations
- Convergence and error analysis for numerical approximations of stochastic Leray-Lions equations
- Multilevel Picard approximations for McKean-Vlasov stochastic differential equations with jumps
- Structure-Preserving Dynamical Low-Rank Approximations for Stochastic Vlasov--Poisson Equations
BibTeX
@misc{norkin2026method,
title = {{Method of successive approximations for solving integral equations of actuarial mathematics}},
author = {Bogdan Norkin},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.03279},
eprint = {2609.03279},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.03279},
note = {CC BY-NC-ND 4.0},
}
RIS
TY - UNPB
TI - Method of successive approximations for solving integral equations of actuarial mathematics
AU - Bogdan Norkin
PY - 2026
DA - 2026/09/03
PB - arXiv
DO - 10.48550/arxiv.2609.03279
UR - https://arxiv.org/abs/2609.03279
AB - The dissertation considers various generalizations of the classical risk process describing the stochastic evolution of the capital of an insurance company (Cramer Lundberg model), including processes with variable deterministic or random premiums, non Poisson flows of premiums and claims, and a stochastic Markovian environment. Integral equations for the probability of nonruin as a function of initial capital are derived for these generalizations. For processes in a stochastic Markovian environment, systems of integral equations for nonruin probabilities corresponding to different initial states are obtained. General necessary and sufficient, as well as specific sufficient, conditions for the existence and uniqueness of solutions are established. A successive approximation method for numerical and analytical solution of the integral equations is theoretically and practically validated; its uniform convergence and rate of convergence are established. A technique for estimating the accuracy of approximate solutions is developed by constructing upper and lower approximations to the exact solution. The method is tested on numerical examples and compared with the Monte Carlo method and known solution approximations. The developed method increases the accuracy of actuarial calculations: it allows the probability of ruin to be calculated with any prescribed accuracy, empirical approximations to be verified and improved iteratively, and the accuracy and parameters of Monte Carlo simulations to be estimated and corrected.
ER -