rgb = ls.shade(z, plt.cm.RdYlBu)
surf = ax.plot_surface(x, y, z, rstride=1, cstride=1, linewidth=0,
    antialiased=False, facecolors=rgb)
plt.show()
mlab.show()
Phong shading for shiny Python 3D surface plots