Call 'simulate' instead of 'conc_final' from the main program, to check that it works.