function []=getfunky

format long

fprintf('elapsed time of model2d');
[TMH,TM,C13T,C13U,C12,t2,x] = model2d031210;
% figure;
% plot(t2,TMH);
% figure;
% plot(t2,log(TMH));
for i = t2(1):t2(end);
    ratio = C13T(:,1)+C13U(:,1)./C12(:,1);
end

fprintf('total methane in headspace at time = {0,30,90,150,210}\n');
TMHatTime = [TMH(1);TMH(31:60:211)]'

fprintf('total methane at time = {0,30,90,150,210}\n');
TMatTime = [TM(1),TM(31:60:211)]

fprintf('ratio of methane at time = {0,30,90,150,210}\n');
RatioatTime = [ratio(1);ratio(31:60:211)]'


%plot labeld methane vs time to find K-value
subplot(221)
plot(t2,C13T(:,31));title('C13T vs time');
subplot(222)
plot(t2, log(C13T(:,31)));title('ln(C13T) vs time');
%use polyfit to calculate slope = -(k) => consumption
[PF, S]=polyfit(t2,log(C13T(:,31)),1);
Slope=PF(1)
