old_contents: from pysb.examples.tyson_oscillator import model
  from pysb.integrate import odesolve
  import numpy as np
  import matplotlib.pyplot as plt
  
  t = np.linspace(0, 100, 10001)
  x = odesolve(model, t)
  
  plt.plot(t, x['CT'],  lw=2, label='YT')  # Good validation of mass balance for cdc2, should be constant at 1
  plt.plot(t, x['YT'], lw=2, label='YT')
  plt.plot(t, x['M'], lw=2, label='M')
  
  plt.legend(loc=0)
  plt.xlabel('time')
  plt.ylabel('population')
  
  plt.show()
  
new_contents: from pysb.examples.tyson_oscillator import model
  from pysb.integrate import odesolve
  import numpy as np
  import matplotlib.pyplot as plt
  
  t = np.linspace(0, 100, 10001)
  x = odesolve(model, t)
  
  plt.plot(t, x['CT'],  lw=2, label='CT')  # Good validation of mass balance for cdc2, should be constant at 1
  plt.plot(t, x['YT'], lw=2, label='YT')
  plt.plot(t, x['M'], lw=2, label='M')
  
  plt.legend(loc=0)
  plt.xlabel('time')
  plt.ylabel('population')
  
  plt.show()
  
subject: Update run_tyson_oscillator.py