function C4H6_kinetics() time = [0; 1000;1800;2800;3600;4400;5200;6200]; %time in seconds data = [0.01;0.00625; 0.00476; 0.00370; 0.00313; 0.00270; 0.00241; 0.00208]; %C4H6 conc in mol/L end