%一个绘制马鞍面的小程序
x=-4:0.1:4;
y=-3:0.2:3;
[X,Y]=meshgrid(x,y);
Z=X.^2-Y.^2;%%马鞍面方程
mesh(X,Y,Z);
hold on
title('心型线套马鞍面')
i=1
for t=1:0.1:100
a=10;
xx=a*(2*cos(t)-cos(2*t));
yy=a*(2*sin(t)-sin(2*t));
xx_drow(i)=xx;
yy_drow(i)=yy;
i=i+1
end
plot(xx_drow,yy_drow)