set(h, 'FaceVertexCData', rand(3,3))
This is the line which gives color to your surfaces. As you have three surfaces there, it should be a 3x3 matrix. Note that RGB color code has 1X3 vector for each color. So, that line selected three random colors. Say you want to fill Red color for each surface. YOu need to follow like this:
color = repmat([1 0 0],3,1)
set(h, 'FaceVertexCData', color)
