Simulation of the Beta-Stacy process