simulation code, results, plot