(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 7.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     11411,        334]
NotebookOptionsPosition[     10765,        311]
NotebookOutlinePosition[     11110,        326]
CellTagsIndexPosition[     11067,        323]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "Mathematica", " ", "code", " ", "to", " ", "draw", " ", "images", " ", 
    "of", " ", "unit", " ", "disk"}], "*)"}], 
  RowBox[{
   RowBox[{"ctor", "[", "z_", "]"}], ":=", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"output", "=", 
      RowBox[{"{", "}"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"Do", "[", 
      RowBox[{
       RowBox[{"AppendTo", "[", 
        RowBox[{"output", ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"Re", "[", 
            RowBox[{"z", "[", 
             RowBox[{"[", "k", "]"}], "]"}], "]"}], ",", 
           RowBox[{"Im", "[", 
            RowBox[{"z", "[", 
             RowBox[{"[", "k", "]"}], "]"}], "]"}]}], "}"}]}], "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"k", ",", "1", ",", 
         RowBox[{"Length", "[", "z", "]"}]}], "}"}]}], "]"}], ";", 
     "\[IndentingNewLine]", "output"}], ")"}]}]}]], "Input",
 CellChangeTimes->{3.4622732997077665`*^9, 3.4622733362951927`*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"plotpolylistCA", "[", "vert_", "]"}], ":=", 
  RowBox[{"(", 
   RowBox[{
    RowBox[{"vert3", "=", 
     RowBox[{"{", "}"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"Do", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"vert1", "=", 
        RowBox[{"ctor", "[", 
         RowBox[{"vert", "[", 
          RowBox[{"[", "k", "]"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"AppendTo", "[", 
        RowBox[{"vert3", ",", "vert1"}], "]"}], ";"}], ",", 
      RowBox[{"{", 
       RowBox[{"k", ",", "1", ",", 
        RowBox[{"Length", "[", "vert", "]"}]}], "}"}]}], "]"}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"plotpolylist", "[", "vert3", "]"}]}], ")"}]}]], "Input",
 CellChangeTimes->{{3.4622732997077665`*^9, 3.462273305972309*^9}}],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"plotpolylist", "[", "vertlist_", "]"}], ":=", 
   RowBox[{"Show", "[", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"(*", 
          RowBox[{
           RowBox[{"GrayLevel", "[", "1", "]"}], ",", 
           RowBox[{"Thickness", "[", ".01", "]"}], ",", 
           RowBox[{"{", 
            RowBox[{"Polygon", "[", 
             RowBox[{"vertlist", "[", 
              RowBox[{"[", "k", "]"}], "]"}], "]"}], "}"}], ","}], "*)"}], 
         RowBox[{
          RowBox[{"GrayLevel", "[", 
           RowBox[{"If", "[", 
            RowBox[{
             RowBox[{"k", "\[Equal]", 
              RowBox[{"Length", "[", "vertlist", "]"}]}], ",", "0", ",", 
             ".1"}], "]"}], "]"}], ",", 
          RowBox[{"Thickness", "[", 
           RowBox[{"If", "[", 
            RowBox[{
             RowBox[{"k", "\[Equal]", 
              RowBox[{"Length", "[", "vertlist", "]"}]}], ",", ".005", ",", 
             ".005"}], "]"}], "]"}], ",", 
          RowBox[{"{", 
           RowBox[{"Line", "[", 
            RowBox[{"vertlist", "[", 
             RowBox[{"[", "k", "]"}], "]"}], "]"}], "}"}]}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
          RowBox[{"Axes", "\[Rule]", "None"}], ",", 
          RowBox[{"AspectRatio", "\[Rule]", "Automatic"}]}], "}"}]}], "]"}], 
      ",", 
      RowBox[{"{", 
       RowBox[{"k", ",", "1", ",", 
        RowBox[{"Length", "[", "vertlist", "]"}]}], "}"}]}], "]"}], 
    "]"}]}]}]], "Input",
 CellChangeTimes->{{3.4622732997077665`*^9, 3.462273309893506*^9}, 
   3.4622733399508114`*^9}],

Cell[BoxData[
 RowBox[{"\n", 
  RowBox[{"makeradial", ":=", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"num0", "=", "16"}], ";", "\[IndentingNewLine]", 
     RowBox[{"angles", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"2", " ", "Pi", " ", 
         RowBox[{"k", "/", "num0"}]}], ",", 
        RowBox[{"{", 
         RowBox[{"k", ",", "1", ",", "num0"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"out", "=", 
      RowBox[{"{", "}"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"lambda1", "=", ".9"}], ";", "\[IndentingNewLine]", 
     RowBox[{"inner", "=", ".1"}], ";", "\[IndentingNewLine]", 
     RowBox[{"num", "=", 
      RowBox[{"Ceiling", "[", 
       RowBox[{
        RowBox[{"Log", "[", "inner", "]"}], "/", 
        RowBox[{"Log", "[", "lambda1", "]"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"radpts", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"lambda1", "^", "k"}], ",", 
        RowBox[{"{", 
         RowBox[{"k", ",", "1", ",", "num"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"AppendTo", "[", 
      RowBox[{"radpts", ",", "0"}], "]"}], ";", "\[IndentingNewLine]", 
     RowBox[{"PrependTo", "[", 
      RowBox[{"radpts", ",", "1"}], "]"}], ";", "\[IndentingNewLine]", 
     RowBox[{"Do", "[", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"AppendTo", "[", 
          RowBox[{"out", ",", 
           RowBox[{
            RowBox[{"Exp", "[", 
             RowBox[{"I", " ", 
              RowBox[{"angles", "[", 
               RowBox[{"[", "k", "]"}], "]"}]}], "]"}], " ", "radpts"}]}], 
          "]"}], ";"}], ")"}], ",", 
       RowBox[{"{", 
        RowBox[{"k", ",", "1", ",", 
         RowBox[{"Length", "[", "angles", "]"}]}], "}"}]}], "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"lambda2", "=", ".8"}], ";", "\[IndentingNewLine]", 
     RowBox[{"inner2", "=", ".25"}], ";", "\[IndentingNewLine]", 
     RowBox[{"num", "=", 
      RowBox[{"Ceiling", "[", 
       RowBox[{
        RowBox[{"Log", "[", "inner2", "]"}], "/", 
        RowBox[{"Log", "[", "lambda2", "]"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"radpts3", "=", 
      RowBox[{"{", 
       RowBox[{
       ".25", ",", ".5", ",", ".75", ",", ".9", ",", ".95", ",", ".99"}], 
       "}"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"Do", "[", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"num", "=", 
          RowBox[{"5", " ", 
           RowBox[{"Ceiling", "[", 
            RowBox[{"2", " ", 
             RowBox[{"Pi", "/", 
              RowBox[{"(", 
               RowBox[{"1", "-", 
                RowBox[{"radpts3", "[", 
                 RowBox[{"[", "k", "]"}], "]"}]}], ")"}]}]}], "]"}]}]}], ";", 
         "\[IndentingNewLine]", 
         RowBox[{"AppendTo", "[", 
          RowBox[{"out", ",", 
           RowBox[{
            RowBox[{"radpts3", "[", 
             RowBox[{"[", "k", "]"}], "]"}], " ", 
            RowBox[{"Table", "[", 
             RowBox[{
              RowBox[{"Exp", "[", 
               RowBox[{"2", " ", "Pi", " ", "I", " ", 
                RowBox[{"k", "/", "num"}]}], "]"}], ",", 
              RowBox[{"{", 
               RowBox[{"k", ",", "0", ",", "num"}], "}"}]}], "]"}]}]}], "]"}],
          ";"}], ")"}], ",", 
       RowBox[{"{", 
        RowBox[{"k", ",", "1", ",", 
         RowBox[{"Length", "[", "radpts3", "]"}]}], "}"}]}], "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"AppendTo", "[", 
      RowBox[{"out", ",", 
       RowBox[{"Table", "[", 
        RowBox[{
         RowBox[{"Exp", "[", 
          RowBox[{"2", " ", "Pi", " ", "I", " ", 
           RowBox[{"k", "/", "100"}]}], "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"k", ",", "0", ",", "100"}], "}"}]}], "]"}]}], "]"}]}], 
    ")"}]}]}]], "Input",
 CellChangeTimes->{{3.4622732997077665`*^9, 3.462273344403167*^9}, {
  3.462273930017418*^9, 3.4622739413137865`*^9}}],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"makeradial", ";"}]}]], "Input",
 CellChangeTimes->{{3.4622735991495557`*^9, 3.4622736043676376`*^9}, 
   3.4622739068465805`*^9}],

Cell[BoxData[
 RowBox[{"plotpolylistCA", "[", "out", "]"}]], "Input",
 CellChangeTimes->{{3.4622744009612103`*^9, 3.4622744074299192`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"f", "[", "z_", "]"}], ":=", 
   RowBox[{"z", "+", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"z", "^", "2"}], ")"}], "/", "2"}]}]}], ";"}], "\n", 
 RowBox[{"plotpolylistCA", "[", 
  RowBox[{"f", "[", "out", "]"}], "]"}]}], "Input",
 CellChangeTimes->{{3.4622732997077665`*^9, 3.4622733242660217`*^9}, 
   3.4622735323143616`*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"f", "[", "z_", "]"}], ":=", 
   RowBox[{"z", "+", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"z", "^", "8"}], ")"}], "/", "8"}]}]}], ";"}], "\n", 
 RowBox[{"plotpolylistCA", "[", 
  RowBox[{"f", "[", "out", "]"}], "]"}]}], "Input",
 CellChangeTimes->{{3.4622732997077665`*^9, 3.4622733295776043`*^9}, 
   3.4622735332829876`*^9}],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"plotpolylistCA", "[", 
   RowBox[{"Exp", "[", "out", "]"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.4622732997077665`*^9, 3.4622733295776043`*^9}, {
  3.4622735339391537`*^9, 3.462273535110879*^9}}],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"Manipulate", "[", 
   RowBox[{
    RowBox[{"plotpolylistCA", "[", 
     RowBox[{"out", " ", "+", " ", 
      RowBox[{"a", " ", 
       RowBox[{"out", "^", "2"}]}], " ", "+", " ", 
      RowBox[{"b", "  ", 
       RowBox[{"out", "^", "3"}]}]}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"a", ",", 
      RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"b", ",", 
      RowBox[{"-", "2"}], ",", "2"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.4622740600132694`*^9, 3.4622741210277233`*^9}}],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"koebefcn", "[", "z_", "]"}], ":=", 
     RowBox[{"z", " ", "+", " ", 
      RowBox[{"2", " ", 
       RowBox[{"z", "^", "2"}]}], " ", "+", " ", 
      RowBox[{"3", 
       RowBox[{"z", "^", "3"}]}], " ", "+", " ", 
      RowBox[{"4", " ", 
       RowBox[{"z", "^", "4"}]}], " ", "+", 
      RowBox[{"5", " ", 
       RowBox[{"z", "^", "5"}]}], " ", "+", 
      RowBox[{"6", " ", 
       RowBox[{"z", "^", "6"}]}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"plotpolylistCA", "[", 
    RowBox[{"koebefcn", "[", "out", "]"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.462274230994362*^9, 3.4622743070555615`*^9}}]
},
WindowSize->{614, 750},
WindowMargins->{{89, Automatic}, {Automatic, 36}},
FrontEndVersion->"7.0 for Microsoft Windows (32-bit) (November 10, 2008)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[545, 20, 1006, 28, 92, "Input"],
Cell[1554, 50, 800, 21, 92, "Input"],
Cell[2357, 73, 1748, 46, 212, "Input"],
Cell[4108, 121, 3990, 106, 412, "Input"],
Cell[8101, 229, 187, 4, 52, "Input"],
Cell[8291, 235, 139, 2, 31, "Input"],
Cell[8433, 239, 381, 11, 52, "Input"],
Cell[8817, 252, 381, 11, 52, "Input"],
Cell[9201, 265, 257, 5, 52, "Input"],
Cell[9461, 272, 585, 16, 72, "Input"],
Cell[10049, 290, 712, 19, 72, "Input"]
}
]
*)

(* End of internal cache information *)