图3.1（心型线）

马鞍线代码如下：

clear; clc; a=10; b=10;

x=-1.5:0.1:1.5; y=-1.5:0.1:1.5;

[X,Y]=meshgrid(x,y);%x行y列的矩阵 Z=X.^2/a-Y.^2/b; figure(2);%新建窗口

subplot(221);%分割当前窗口，作图显示在该分割块上 surf (X,Y,Z); title('马鞍面'); xlabel('x'); ylabel('y'); subplot(222); surf (X,Y,Z); title('马鞍面');

xlabel('x'); ylabel('y');

view(20,30);%改变视角到（20，30) subplot(223); surf (X,Y,Z); title('马鞍面'); xlabel('x'); ylabel('y');

view(50,20);%改变视角到（50,20） subplot(224); surf (X,Y,Z); title('马鞍面'); xlabel('x'); ylabel('y');

view(60,10);%改变视角到（60,10）

马鞍线运行结果显示如下图：

图3.2（马鞍面）

7.以不同的视角观察球面x2+y2+z2=r2

和圆柱面 x2+y2=rx所围区域. 解：

代码如下：

clear; clc; r=2;

[x0,y0,z0]=sphere(50); x=r*x0;y=r*y0;z=r*z0; surf(x,y,z);hold on; f=@(x,y)x.^2+y.^2-r*x;

ezsurf(f,[-1.5,2.5],[-2,1.5]); axis equal; view(30,20);

运行结果显示如下图：

图3.3