text: from matplotlib import pyplot as plt
  rate=5.10e+3
  
  plt.axhline(y=rate, label='Old Data') #plots the old rate as a horizontal line
  plt.plot([ 0.01099996,  0.01199995,  0.01299995,  0.01399995,  0.01499994,
          0.01599994,  0.01799993,  0.01999992,  0.02499991,  0.02999989,
          0.03999985,  0.04999981,  0.05999977,  0.06999974,  0.0799997 ,
          0.08999966,  0.09999962,  0.10999958,  0.11999955,  0.12999951,
          0.13999947,  0.14999943,  0.1599994 ,  0.17999932,  0.19999925,
          0.24999906,  0.29999887,  0.34999868,  0.39999849,  0.4499983 ,
          0.49999811,  0.59999774,  0.69999736,  0.79999698,  0.8999966 ,
          0.99999623,  1.24999586,  1.49999434,  1.74999397,  1.99999245],[5.106e+3, 5.106e+3, 5.106e+3, 5.106e+3, 5.106e+3, 5.106e+3, 5.106e+3, 5.106e+3, 5.106e+3, 5.106e+3, 5.106e+3, 5.106e+3, 5.106e+3, 5.106e+3, 5.106e+3,5.106e+3, 5.106e+3, 5.106e+3, 5.106e+3, 5.106e+3, 5.106e+3, 5.106e+3, 5.106e+3, 5.106e+3, 5.106e+3, 5.106e+3, 5.106e+3, 5.107e+3, 5.108e+3, 5.111e+3, 5.115e+3, 5.13e+3, 5.153e+3, 5.184e+3, 5.223e+3, 5.269e+3, 5.409e+3, 5.572e+3, 5.75e+3, 5.934e+3], 'ro', label='New Data') #Retrieved from ENDF 2011, converted eV to t9 
  
  
  plt.xlabel('T9')
  plt.ylabel('Reaction Rates (cm3/s/mol)')
  plt.title('Li6 (n,g) Li7')
  plt.legend()
  plt.axis([0,2.5,5000,5600])
  plt.show()
  