Code for interpreting the grabbed output

• Run through the output string, line by line
• If the current line prints \( \theta \), \( \Delta t \), and \( E \), split the line into these three pieces and store the data
• Store data in a dictionary errors with keys dt and the three \( \theta \) values

errors = {'dt': dt_values, 1: [], 0: [], 0.5: []}
for line in output.splitlines():
    words = line.split()
    if words[0] in ('0.0', '0.5', '1.0'):  # line with E?
        # typical line: 0.0   1.25:    7.463E+00
        theta = float(words[0])
        E = float(words[2])
        errors[theta].append(E)

Next: plot \( E \) versus \( \Delta t \) for \( \theta=0,0.5,1 \)

Complete program: experiments/decay_exper1.py. Fine recipe for

• how to run other programs
• how to extract and interpret output from other programs
• how to automate many manual steps in creating simulations and figures