有需求，就会有代码：
from matplotlib import pyplot as plt from mpl_toolkits.mplot3d import Axes3D import numpy as np fig = plt.figure() ax = plt.axes(projection='3d') xx = np.arange(-0.1, 0.1, 0.01) yy = np.arange(-0.1, 0.1, 0.01) X, Y = np.meshgrid(xx, yy) Z1 = np.arctan((1+Y)/(1+X)) - np.arctan(1) Z2 = np.arctan((1+Y)/(1-X)) - np.arctan(1) Z3 = np.arctan((1-Y)/(1+X)) - np.arctan(1) Z4 = np.arctan((1-Y)/(1-X)) - np.arctan(1) ax.plot_surface(X, Y, Z1, cmap='rainbow') ax.plot_surface(X, Y, Z2, cmap='rainbow') ax.plot_surface(X, Y, Z3, cmap='rainbow') ax.plot_surface(X, Y, Z4, cmap='rainbow') plt.show()
