# This import registers the 3D projection, but is otherwise unused. from mpl_toolkits.mplot3d import Axes3D # noqa: F401 unused import
import matplotlib.pyplot as plt import numpy as np import time
defgenerate(X, Y, phi): ''' Generates Z data for the points in the X, Y meshgrid and parameter phi. ''' R = 1 - np.sqrt(X**2 + Y**2) return np.cos(2 * np.pi * X + phi) * R
# Make the X, Y meshgrid. xs = np.linspace(-1, 1, 50) ys = np.linspace(-1, 1, 50) X, Y = np.meshgrid(xs, ys)
# Set the z axis limits so they aren't recalculated each frame. ax.set_zlim(-1, 1)
# Begin plotting. wframe = None tstart = time.time() for phi in np.linspace(0, 180. / np.pi, 100): # If a line collection is already remove it before drawing. if wframe: ax.collections.remove(wframe)
# Plot the new wireframe and pause briefly before continuing. Z = generate(X, Y, phi) wframe = ax.plot_wireframe(X, Y, Z, rstride=2, cstride=2) plt.pause(.001)