The convergence of an iterative procedure for approximate solution of an equation \(ν A= μ\) is studied where \(A\) is a Markov operator, \(ν\) and \(μ\) are probability measures. The iterative procedure is based on Bayes formula. The convergence is studied with the help of the coupling method.