Z = np.sinc(np.sqrt(X * X + Y * Y))
V = np.sin(Y)
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
ax.view_init(45, 60)
ax.plot_surface(X, Y, Z, facecolors=cm.Oranges(V))
How to make a 4d plot using Python with matplotlib