calabi yau manifold

calabi yau manifold

matlab code:

t=pi*2/3;  %转角控制参数
n=6;    %指数控制变量
k=2;    %[k] 和 [g] 遍历 [0, n-1] ，因此共有 [n*n] 种不同的组合
g=4;    %每一个组合都确定一个平面片，所有的平面片组成一个奇异的曲面

a=-1:0.05:1;
b=0:0.05:pi/2;
[m,s] = meshgrid(a,b);

%u1=exp(2*k*pi*1i/n)*cosh(m+s*1i).^(2/n);
%u2=exp(2*g*pi*1i/n)*sinh(m+s*1i).^(2/n);

bn=n-1;
for k=0:bn
    for g=0:bn
        x=real(exp(2*k*pi*1i/n)*cosh(m+s*1i).^(2/n));
        y=real(exp(2*g*pi*1i/n)*sinh(m+s*1i).^(2/n));
        z=cos(t)*imag(exp(2*k*pi*1i/n)*cosh(m+s*1i).^(2/n))+sin(t)*imag(exp(2*g*pi*1i/n)*sinh(m+s*1i).^(2/n));
        surface(x,y,z);
    end
end