  The numerical simulation of domain wall quarks with the two-step multi-boson
(TSMB) algorithm is considered. The inclusion of single quark flavours, as
required for strange quarks, is discussed. The usage of computer memory can be
kept relatively low, independently of the order of polynomial approximations.
Tests are performed with two flavours (N_f=2) of degenerate quarks near the
N_t=4 thermodynamical cross over.
