Introduction - If you have any usage issues, please Google them yourself
This submission contains a function that computes face normal data for the VertexNormals property of the surface that fixes the lighting (facenorm.m) as well as a workaround function that you can just use in place of surf (surfo.m).
Only use this if you get erroneous surface plots and you can t go without OpenGL and flat lighting.