Numerical solution to a system of differential equations for probability measures

A system of ordinary differential equations describing a stationary distribution of a Markov process with the phase space $\mathbf{R}\times\{1, 2,\dots, M\}$ is considered. A numerical method for finding and calculating its solution being a probability density function is proposed.