load monthly_mlo.csv
Keel = monthly_mlo(3:640, [4 9]);
yK = Keel(:,1); TK = Keel(:,2);
hK = plot(yK, TK, 'Color', [0 .7 .7]);
title(['AGREEMENT OF {\color[rgb]{.7 0 0}HOFMANN''S LAW} WITH ' ...
    '{\color[rgb]{0 .7 .7}KEELING CURVE}'], ...
    'FontSize', 16);
text(1981, 333, 'Keeling Curve', 'FontSize', 14, 'Color', [0 .7 .7]);
xlabel('Year', 'FontSize', 14);
grid on
ylabel('Atmospheric CO2 (ppmv)', 'FontSize', 16, 'Color', 'k');
set(gca, 'XLim', [1850 2010], 'YLim', [180 400], 'FontSize', 16, ...
	'YTick', 180:20:400, 'YGrid', 'on');
y = 1850:2010;
base = 287;
% Draw CO2 estimates for base 287
%~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
%Hofmann's raised exponential
line(y, 287 + 1.0245.^(y - 1818), 'Color', [.7 0 0], 'LineWidth', 1.5);
%Goodman's exponential
line(y, 1.00405.^(y - 535), 'Color', [0 .7 0], 'LineWidth', 1.5);
%~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
text(1920, 350, 'Hofmann''s Law (AMU 2009)', 'FontSize', 14, ...
    'Color', [.7 0 0], 'BackgroundColor', 'w', 'Margin', .01);
text(1920, 339, 'CO_2(Y) = 287 + 1.0245^{Y - 1818}', ...
    'FontSize', 14, 'Color', [.7 0 0], 'BackgroundColor', 'w');
text(1920, 250, 'Goodman''s Law (Climate Etc 2012)', 'FontSize', 14, ...
    'Color', [0 .7 0], 'BackgroundColor', 'w', 'Margin', .01);
text(1920, 239, 'CO_2(Y) = 1.00405^{Y - 535}', ...
    'FontSize', 14, 'Color', [0 .7 0], 'BackgroundColor', 'w');
