有趣的surface函数

整理文档时偶然发现了在暑假时玩的一个有趣的surface函数，分享给大家。

在心形图形上贴上图片：
clear
clc
[c,f]=meshgrid(0:0.1:pi,0:0.1:2*pi);

x=sin(c).*cos(f);
y=sin(c).*sin(f);
z=cos(c);
z=z;

y=2.*y;
x=5.*x;
z=1.3.*z;

[x1,y1,z1]=cart2pol(x,y,z); 
   
x1=x1+pi/4;
x4=x1-pi/2;

[x2,y2,z2]=pol2cart(x1,y1,z1);
[x3,y3,z3]=pol2cart(x4,y1,z1);

q=find(y3<0);w=size(q);
for i=1:w
    v(i)=mod(q(i),63);
    if v(i)==0
        v(i)=63;
    end 
    u(i)=(q(i)-v(i))/63;
end

for i=1:w    
    if v(i)==0