  L\"uscher's local bosonic algorithm for Monte Carlo simulations of quantum
field theories with fermions is applied to the simulation of a possibly
supersymmetric Yang-Mills theory with a Majorana fermion in the adjoint
representation. Combined with a correction step in a two-step polynomial
approximation scheme, the obtained algorithm seems to be promising and could be
competitive with more conventional algorithms based on discretized classical
(``molecular dynamics'') equations of motion. The application of the considered
polynomial approximation scheme to optimized hopping parameter expansions is
also discussed.
