#Originally developed with Python 2.7.3 using Matplotlib version 1.1.1rc
#Code also works fine with Anaconda3 / matplotlib version 1.4.0!
#Code also works fine with CPython 3.4.2 + Scipy-stack-14.8.27.win32-py3.4 from
#   Christoph Gohlke's unofficial libraries:
#    http://www.lfd.uci.edu/~gohlke/pythonlibs/


from mpl_toolkits.mplot3d import Axes3D
from matplotlib import cm
from matplotlib.ticker import LinearLocator, FormatStrFormatter
import matplotlib.pyplot as plt
import numpy as np

fig = plt.figure()
ax = fig.gca(projection='3d')

X = np.arange(-5, 5, 0.1)
Y = np.arange(-5, 5, 0.1)

X, Y = np.meshgrid(X, Y)

D = np.sqrt(X**2 + Y**2)

Z = np.sin(D)

surf = ax.plot_surface(X, Y, Z, rstride=1, cstride=1, cmap=cm.coolwarm,
        linewidth=0, antialiased=False)  #Try cool vs coolwarm
ax.set_zlim(-1.01, 1.01)

ax.zaxis.set_major_locator(LinearLocator(10))
ax.zaxis.set_major_formatter(FormatStrFormatter('%.02f'))

fig.colorbar(surf, shrink=0.5, aspect=5)

plt.show()