[X, Y, Z] = ndgrid(linspace(-8,8,40), linspace(-8,8,40), linspace(-8,8,40));
D = X.*Y.*Z;
D = D * (100 / max(D(:)));

figure
for ii = 10:10:100
    patch(isosurface(D, ii, D), 'FaceColor', 'interp', 'LineStyle', 'none')
end
t = text(30, 20, -10, 't=5');  % step 1
view([3 1 1])                  % step 2
set(t, 'Units', 'normalized')  % step 3

figure
for ii = 10:10:100
    patch(isosurface(D, ii, D), 'FaceColor', 'interp', 'LineStyle', 'none')
end
t = text(30, 20, -10, 't=5');  % step 1
view([1 2 1])                  % step 2
set(t, 'Units', 'normalized')  % step 3