clear all; close all; clc;
load a1.csv;


k=0;

for i=1986:2016;
if i==1986
    quarters=[12];
else
     quarters=[3 6 9 12];
end

        for j=quarters
            k=k+1;
            T(k,1)=datenum(i,j,1);
        end
end
	           
plot(T,a1,'b','LineWidth',2);
set(gca,'XTick',[T(1):365*1:T(end)]);
datetick('x','yy','keepticks');
recessionplot;
ax = gca;
ax.YAxis.Exponent=0; 
set(gca,'FontSize',14, 'FontName','Times');

