#Script to illustrate using Sage to analyze data concerning the decomposition of butadiene.

#Times at which data was taken (seconds)
times = [0, 600, 1200, 1800, 2400, 3000, 3600, 4200, 4800, 6000]; #Times in seconds

#Butadiene concentration, moles per liter at each time above
data = [0.31, 0.254, 0.208, 0.172, 0.141, 0.116, 0.0964, 0.0812, 0.0669, 0.0464];

#Plot the raw data versus time. Call the plot "plot1".
pdata = list(zip(times,data))
plt1 = scatter_plot(pdata)
show(plt1)

#Does not look 0th order. Is it first order? Try a logarithmic transformation of the data (as was done for H2O2)
#or butadiene.  See those notebooks.