This code has been used test a simulation of the first bechmark case.