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

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