  In this paper, we are concerned with the numerical solution of one type
integro-differential equation by a probability method based on the fundamental
martingale of mixed Gaussian processes. As an application, we will try to
simulate the estimation of ruin probability with an unknown parameter driven
not by the classical L\'evy process but by the mixed fractional Brownian
motion.
