Shading command produces different results in MATLAB and Octave

The commands shading faceted and shading interp produce different figures in MATLAB and Octave.

Consider this code:

[X, Y] = meshgrid(0:2, 0:2);
Z = magic(3);

figure;
colormap('jet');

subplot(1, 3 ,1);
imagesc(Z);
axis xy
yticklabels({'0', '1', '2'})
yticks(1:3)
xticklabels({'0', '1', '2'})
xticks(1:3)
title('imagesc(Z)');

subplot(1,3,2);
surf(X, Y, Z);
shading faceted;
view([0,0,1]);
yticks(0:3)
xticks(0:3)
title('surf(X, Y, Z); shading faceted;');

subplot(1,3,3);
surf(X, Y, Z);
shading interp;
view([0,0,1]);
yticks(0:3)
xticks(0:3)
title('surf(X, Y, Z); shading interp;');

This is the result in MATLAB:

And the result generated by Octave:

MATLAB generates a correct figure.What is wrong with this Octave figure?

What is the reason we get two different outputs? Is it an issue with to the view command?

EDIT (After Bug Report)

It's been fixed for Ocatve 6.1: https://hg.savannah.gnu.org/hgweb/octave/raw-file/11072ea6a16c/scripts/plot/appearance/view.m