mat

iceN = reshapetimeseries(t, extent_N, 'bin', 52);
iceS = reshapetimeseries(t, extent_S, 'bin', 52);

ice = cat(3, iceN, iceS);

%% 
% By default, the function calculates the mean and upper/lower bounds of
% the data at each x-location in the data.  We can plot either
% <boundedline_documentation.html bounded lines>:

[A, h] = ensemble2bnd(1:52,ice, 'dims', 'xey', 'plot', 'boundedline')
set(gca,'xlim', [1 52]);
xlabel('Week of year');
ylabel('Ice extent');

%%
% line plots with error bars:

cla
[A, h] =