(* ::Package:: *)

(************************************************************************)
(* This file was generated automatically by the Mathematica front end.  *)
(* It contains Initialization cells from a Notebook file, which         *)
(* typically will have the same name as this file except ending in      *)
(* ".nb" instead of ".m".                                               *)
(*                                                                      *)
(* This file is intended to be loaded into the Mathematica kernel using *)
(* the package loading commands Get or Needs.  Doing so is equivalent   *)
(* to using the Evaluate Initialization Cells menu command in the front *)
(* end.                                                                 *)
(*                                                                      *)
(* DO NOT EDIT THIS FILE.  This entire file is regenerated              *)
(* automatically each time the parent Notebook file is saved in the     *)
(* Mathematica front end.  Any changes you make to this file will be    *)
(* overwritten.                                                         *)
(************************************************************************)



cCos[theta_,xi_]:=.5 (\[ExponentialE]^(xi+\[ImaginaryI] theta)+\[ExponentialE]^(-xi-\[ImaginaryI] theta))
cSin[theta_,xi_]:=(-.5 \[ImaginaryI]) (\[ExponentialE]^(xi+\[ImaginaryI] theta)-\[ExponentialE]^(-xi-\[ImaginaryI] theta))
z1[theta_,xi_,n_,k_]:=\[ExponentialE]^((k 2 \[Pi] \[ImaginaryI])/n) cCos[theta,xi]^(2.0/n)
z2[theta_,xi_,n_,k_]:=\[ExponentialE]^((k 2 \[Pi] \[ImaginaryI])/n) cSin[theta,xi]^(2.0/n)
z1p[theta_,xi_,n_,k_]:=\[ExponentialE]^((xi+\[ImaginaryI] theta)/n) \[ExponentialE]^((k 2 \[Pi] \[ImaginaryI])/n)
z2p[theta_,xi_,n_,k_]:=\[ExponentialE]^(-((xi+\[ImaginaryI] theta)/n)) \[ExponentialE]^(-((k 2 \[Pi] \[ImaginaryI])/n))


MakePolygons[vl_List]:=Block[{l=vl,l1=RotateLeft/@vl,mesh},mesh={l,l1,RotateLeft[l1],RotateLeft[l]};mesh=Map[Drop[#1,-1]&,mesh,{1}];mesh=Map[Drop[#1,-1]&,mesh,{2}];Polygon/@Transpose[(Flatten[#1,1]&)/@mesh]]