We consider the numerical solution of one type integro-differential equation by a
probability method based on fundamental martingale mixed Gaussian processes. As an
application, we try to simulate estimation ruin with an unknown
parameter driven not classical Levy process, but fractional Brownian
motion.