% Axes settings
set(gca,'FontSize',14);
axis vis3d equal;
view([0,0]);
grid off; axis off;
xlim([-5,10]);
ylim([-2,8]);
zlim([-55,20]);

% FLAG: choose f = 1 if want to show 5 pt stencil, 0 otherwise
f = 1;

% put together the file name
filename='an_nat_ord_7x7'; % base file name
if( f==1 )
    filename = [ filename '_5pt']; % with 5 pt stencil
else
    filename = [ filename '_no5pt']; % without 5 pt stencil
end
filename=[ filename '.avi']; % .avi extension

[~,N] = b700(f,1);

cla
[~,M0] = b700(f,1);

cla
[~,M1] = b701(f,1);

cla
[~,M2] = b702(f,1);

cla
[~,M3] = b703(f,1);

cla
[~,M4] = b704(f,1);

cla
[~,M5] = b705(f,1);

cla
[~,M6] = b706(f,1);

cla
[~,M7] = b707(f,1);

cla
[~,M8] = b708(f,1);

cla
[~,M9] = b709(f,1);

cla
[~,M10] = b710(f,1);

cla
[~,M11] = b711(f,1);

cla
[~,M12] = b712(f,1);

cla
[~,M13] = b713(f,1);

cla
[~,M14] = b714(f,1);

M = [M0, M1, M2,M3,M4,M5,M6,M7,M8,M9, M10,M11,M12,M13,M14];

movie2avi(M, filename,'quality',10);
