def main(output_directory=DEFAULT_OUTPUT_DIRECTORY):
    """ Run simulation and plot results

    Args:
        output_directory (:obj:`str`, optional): directory to save plots