|
一般的程序代码:
Clear["Global`"];
n = 8; \[Omega] = Exp[I*(2\[Pi])/n]; Subscript[A, 0] = -1 - I; Subscript[A, 1] = 1 + I; Z = x + y I;
Subscript[Z, 0] = Subscript[A, 0] + \[Omega]^k * (Z - Subscript[A, 0]) -(\[Omega]^k - 1)/(1 - \[Omega]^-1) * (Subscript[A, 1] - Subscript[A, 0]);
trans = (ReIm[#]//ComplexExpand//FullSimplify//Factor)&/@Table[Subscript[Z, 0], {k, 0, n-1}];
curve[x_, y_]:= x^2 - x y + y^2 -1;
newcurves = (curve[#[[1]], #[[2]]] // Factor)&/@trans;
Show[ContourPlot[# == 0, {x, -(n+2), (n+2)}, {y, -(n+2), (n+2)}]&/@ newcurves]
复制到Mathematica中运行,修改 n 值即可.
|
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有帐号?注册
x
|