您可以先为X和Y创建一个meshgrid。然后通过执行Z=np.sin(np.sqrt(X**2 + Y**2))来计算Z。最后，您可以使用matplotlib函数ax.plot_surface(X, Y, Z)绘制曲面。
你可以在下面找到代码：

import numpy as np
from matplotlib import pyplot as plt
from mpl_toolkits.mplot3d import Axes3D

N_points=100
x = np.linspace(-10, 10, N_points)  
y = np.linspace(-10, 10, N_points)  
  
X, Y = np.meshgrid(x, y)
Z=np.sin(np.sqrt(X**2 + Y**2))

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

ax.plot_surface(X, Y, Z)

ax.set_xlabel('X')
ax.set_ylabel('Y')
ax.set_zlabel('Z')
plt.show()

这段代码的输出是：