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