plotcubedsphere.m

n = 16;
[Vertices Faces FaceID] = cubedsphere(n, 'equiangular');

color = [1 0.5 0.5;
         0.5 1 0.5;
         0.5 0.5 1;
         1 1 0;
         0 1 1;
         1 0 1];
CData = color(FaceID,:);
fig = figure(1); clf(fig);
pos=get(fig, 'Position');
pos([3 4]) = 400;
set(fig, 'Position', pos);
ax=axes('Parent',fig,'pos',[0 0 1 1]);

patch('Vertices',Vertices,'Faces',Faces, ...
      'FaceVertexCdata',CData,'FaceColor','flat','parent',ax);
axis(ax,'equal')
axis(ax,'off')
view(ax,3)