load('data2.mat');
figure;
contourf(K,D,r);
title('Total Methane - exponential rate');xlabel('K');ylabel('D');colorbar; 
figure;
contourf(K,D,m);
title('Total Methane - linear rate');xlabel('K');ylabel('D');colorbar;

