We should explain why using `fig, ax = plt.subplots()` and then `ax.XXX` is preferred by us over `plt.XXX`