  A new algorithm for simulation of theories with dynamical fermions is
presented. The algorithm is based on obtaining the new configuration U' from
the old one U by solving the equation M(U')\eta= \omega M(U)\eta, where M is
fermionic operator, \eta is random Gaussian vector, and \omega is random real
number close to unity. This algorithm can be used for acceleration of current
simulations in theories with fermions.
