  We describe an implementation of the Rational Hybrid Monte Carlo (RHMC)
algorithm for dynamical computations with two flavours of staggered quarks. We
discuss several variants of the method, the performance and possible sources of
error for each of them, and we compare the performance and results to the
inexact R algorithm.
