Plan for the experiments

• Plot $u^n$ against $\uex = Ie^{-at}$ for various choices of the parameters $I$, $a$, $\Delta t$, and $\theta$
• How does the discrete solution compare with the exact solution when $\Delta t$ is varied and $\theta=0,0.5,1$?
• Use the decay_mod.py module (little modification of the plotting, see experiments/decay_mod.py)
• Make separate program for running (automating) the experiments (script)

1. python decay_mod.py --I 1 --a 2 --makeplot --T 5 --dt 0.5 0.25 0.1 0.05
2. Combine generated figures FE_*.png, BE_*.png, and CN_*.png to new figures with multiple plots
3. Run script as python decay_exper0.py 0.5 0.25 0.1 0.05 ($\Delta t$ values on the command line)