    ax[0].plot(ts, yerr, k, label='VODE')
    ax[1].plot(ts[1:], dts, k)
    ax[2].plot(ts, ys, k)



    # Load and plot oomph-lib data
    with open("clean-trace") as oomph_file:
        lines = oomph_file.read().splitlines()

    headers = lines[0].split()