clc,clear,close all; [X,Y]=meshgrid(-2:0.1:2); Z = -1./sqrt((X-1).^2+Y.^2+0.01)-1./sqrt((X+1).^2+Y.^2+0.01)+1./sqrt(X.^2+(Y-1).^2+0.01); [DX,DY]=gradient(Z); contour(X,Y,Z,13);hold on; quiver(X,Y,DX,DY);

这段代码的作用是绘制一个二维曲面，并在曲面上叠加等高线图和矢量图。首先，代码通过创建一个网格矩阵(X, Y)来定义曲面的坐标范围。然后，使用给定的函数表达式计算曲面的高度(Z)。接下来，使用梯度函数计算曲面的梯度值(DX, DY)。最后，使用contour函数绘制等高线图，并使用quiver函数绘制矢量图。

请注意，在运行这段代码之前，需要确保已经正确安装了MATLAB软件，并在MATLAB命令窗口中运行这段代码。