t=readtable("US_PREG_03.xlsx");
rawdata = table2array(t(:,2:end));
rawdata = strrep(rawdata, ',', '');
rawdata = str2double(rawdata);
timeVector = datetime(1920, 1, 1):calmonths(1):datetime(2023, 12, 1);
productionVector = reshape(rawdata', [], 1);
plot(timeVector, productionVector);
title('Monthly Oil Production (1920-2023)');
ylabel('Oil Production');
bar(years, rawdata, 'stacked');
title('Yearly Oil Production (1920-2023)');
subtitle('with Monthly Breakdown');
ylabel('Oil Production');
legend('Jan', 'Feb', 'Mar', 'Apr', 'May', 'Jun', 'Jul', 'Aug', 'Sep', 'Oct', 'Nov', 'Dec', ...
'Location', 'BestOutside');
months = {'Jan', 'Feb', 'Mar', 'Apr', 'May', 'Jun', 'Jul', 'Aug', 'Sep', 'Oct', 'Nov', 'Dec'};
years = {'2019', '2020', '2021', '2022'};
data = rand(4, 12) * 100;
dataTable = array2table(data, 'RowNames', years, 'VariableNames', months);
writetable(dataTable, 'data.xlsx', 'WriteRowNames', true);
dataTable = readtable('data.xlsx', 'ReadRowNames', true);
data = table2array(dataTable);
years = dataTable.Properties.RowNames;
months = dataTable.Properties.VariableNames;
heatmap(months, years, data);
title('Heatmap of Monthly Data');
plot(data', 'LineWidth', 2);
legend(years, 'Location', 'northeastoutside');
xticks(1:length(months));
title('Monthly Data Over Years');