(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 9.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 157, 7] NotebookDataLength[ 69818, 1833] NotebookOptionsPosition[ 67547, 1752] NotebookOutlinePosition[ 67905, 1768] CellTagsIndexPosition[ 67862, 1765] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["\<\ Effective area for 132 path, notation as in sketch\ \>", "Subtitle", CellChangeTimes->{{3.626152249443516*^9, 3.6261522590926456`*^9}, { 3.626333424940516*^9, 3.6263334515070972`*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"Pn", "[", RowBox[{"n_", ",", "\[Delta]_"}], "]"}], ":=", RowBox[{"{", RowBox[{ RowBox[{"-", RowBox[{"Sin", "[", RowBox[{"\[Delta]", "-", RowBox[{ RowBox[{"(", RowBox[{"n", "-", "2"}], ")"}], " ", RowBox[{"\[Pi]", "/", "3"}]}]}], "]"}]}], ",", RowBox[{"Cos", "[", RowBox[{"\[Delta]", "-", RowBox[{ RowBox[{"(", RowBox[{"n", "-", "2"}], ")"}], " ", RowBox[{"\[Pi]", "/", "3"}]}]}], "]"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"G", "=", RowBox[{"Simplify", "[", RowBox[{"Solve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"y2", "\[Equal]", RowBox[{ RowBox[{"a", " ", "x2"}], "+", "b"}]}], ",", RowBox[{"y1", "\[Equal]", RowBox[{ RowBox[{"a", " ", "x1"}], "+", "b"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"a", ",", "b"}], "}"}]}], "]"}], "]"}]}]}], "Input", CellChangeTimes->{{3.626151445653192*^9, 3.626151502956397*^9}, { 3.626170791564748*^9, 3.626170792118301*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"a", "\[Rule]", FractionBox[ RowBox[{"y1", "-", "y2"}], RowBox[{"x1", "-", "x2"}]]}], ",", RowBox[{"b", "\[Rule]", FractionBox[ RowBox[{ RowBox[{ RowBox[{"-", "x2"}], " ", "y1"}], "+", RowBox[{"x1", " ", "y2"}]}], RowBox[{"x1", "-", "x2"}]]}]}], "}"}], "}"}]], "Output", CellChangeTimes->{{3.626151470133321*^9, 3.626151480897431*^9}, 3.626151568166727*^9, {3.626170769080858*^9, 3.626170792546183*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"Clear", "[", RowBox[{"Gn", ",", "F1", ",", "x1", ",", "p1y", ",", "p1x", ",", "x1s"}], "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Gn", "[", RowBox[{"n_", ",", "\[Delta]_"}], "]"}], ":=", RowBox[{"FullSimplify", "[", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", " ", "x"}], "+", "b"}], ")"}], "/.", RowBox[{"(", RowBox[{ RowBox[{"G", "[", RowBox[{"[", "1", "]"}], "]"}], "/.", RowBox[{"{", RowBox[{ RowBox[{"x1", "\[Rule]", RowBox[{ RowBox[{"Pn", "[", RowBox[{"n", ",", "\[Delta]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], ",", RowBox[{"y1", "\[Rule]", RowBox[{ RowBox[{"Pn", "[", RowBox[{"n", ",", "\[Delta]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}]}], ",", RowBox[{"x2", "\[Rule]", RowBox[{ RowBox[{"Pn", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "\[Delta]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], ",", RowBox[{"y2", "\[Rule]", RowBox[{ RowBox[{"Pn", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "\[Delta]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}]}]}], "}"}]}], ")"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"F1", "[", RowBox[{"\[Delta]_", ",", "x1_"}], "]"}], ":=", RowBox[{"FullSimplify", "[", RowBox[{"(", RowBox[{ RowBox[{ FractionBox[ RowBox[{ RowBox[{ RowBox[{"Pn", "[", RowBox[{"2", ",", "\[Delta]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}], "-", RowBox[{ RowBox[{"Pn", "[", RowBox[{"1", ",", "\[Delta]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}]}], RowBox[{ RowBox[{ RowBox[{"Pn", "[", RowBox[{"2", ",", "\[Delta]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], "-", RowBox[{ RowBox[{"Pn", "[", RowBox[{"1", ",", "\[Delta]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}]], RowBox[{"(", RowBox[{"x", "-", "x1"}], ")"}]}], "+", RowBox[{ RowBox[{"Pn", "[", RowBox[{"1", ",", "\[Delta]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}]}], ")"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"F2", "[", RowBox[{"\[Delta]_", ",", "x1_"}], "]"}], ":=", RowBox[{"FullSimplify", "[", RowBox[{"(", RowBox[{ RowBox[{ FractionBox[ RowBox[{ RowBox[{ RowBox[{"Pn", "[", RowBox[{"6", ",", "\[Delta]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}], "-", RowBox[{ RowBox[{"Pn", "[", RowBox[{"1", ",", "\[Delta]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}]}], RowBox[{ RowBox[{ RowBox[{"Pn", "[", RowBox[{"6", ",", "\[Delta]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], "-", RowBox[{ RowBox[{"Pn", "[", RowBox[{"1", ",", "\[Delta]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}]], RowBox[{"(", RowBox[{"x", "-", "x1"}], ")"}]}], "+", RowBox[{ RowBox[{"Pn", "[", RowBox[{"1", ",", "\[Delta]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}]}], ")"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"F3", "[", RowBox[{"\[Delta]_", ",", "x1_"}], "]"}], ":=", RowBox[{"FullSimplify", "[", RowBox[{"(", RowBox[{ RowBox[{ FractionBox[ RowBox[{ RowBox[{ RowBox[{"Pn", "[", RowBox[{"5", ",", "\[Delta]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}], "-", RowBox[{ RowBox[{"Pn", "[", RowBox[{"6", ",", "\[Delta]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}]}], RowBox[{ RowBox[{ RowBox[{"Pn", "[", RowBox[{"5", ",", "\[Delta]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], "-", RowBox[{ RowBox[{"Pn", "[", RowBox[{"6", ",", "\[Delta]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}]], RowBox[{"(", RowBox[{"x", "-", "x1"}], ")"}]}], "+", RowBox[{ RowBox[{"Pn", "[", RowBox[{"6", ",", "\[Delta]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}]}], ")"}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"x1s", "=", RowBox[{ RowBox[{"FullSimplify", "[", RowBox[{"Solve", "[", RowBox[{ RowBox[{ RowBox[{"F1", "[", RowBox[{"\[Delta]", ",", "x1"}], "]"}], "\[Equal]", RowBox[{"Gn", "[", RowBox[{"2", ",", "\[Delta]"}], "]"}]}], ",", "x"}], "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"{", RowBox[{"\"\\"", ",", RowBox[{"p1x", "=", RowBox[{"FullSimplify", "[", RowBox[{"(", RowBox[{"x", "/.", "x1s"}], ")"}], "]"}]}], ",", RowBox[{"p1y", "=", RowBox[{"FullSimplify", "[", RowBox[{ RowBox[{"Gn", "[", RowBox[{"2", ",", "\[Delta]"}], "]"}], "/.", "x1s"}], "]"}]}]}], "}"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"x2s", "=", RowBox[{ RowBox[{"FullSimplify", "[", RowBox[{"Solve", "[", RowBox[{ RowBox[{ RowBox[{"F1", "[", RowBox[{"\[Delta]", ",", "x1"}], "]"}], "\[Equal]", RowBox[{"Gn", "[", RowBox[{"3", ",", "\[Delta]"}], "]"}]}], ",", "x"}], "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"{", RowBox[{"\"\\"", ",", RowBox[{"p2x", "=", RowBox[{"FullSimplify", "[", RowBox[{"x", "/.", "x2s"}], "]"}]}], ",", RowBox[{"p2y", "=", RowBox[{"FullSimplify", "[", RowBox[{ RowBox[{"Gn", "[", RowBox[{"3", ",", "\[Delta]"}], "]"}], "/.", "x2s"}], "]"}]}]}], "}"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"x3s", "=", RowBox[{ RowBox[{"FullSimplify", "[", RowBox[{"Solve", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"Pn", "[", RowBox[{"1", ",", "\[Delta]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}], "\[Equal]", RowBox[{"Gn", "[", RowBox[{"3", ",", "\[Delta]"}], "]"}]}], ",", "x"}], "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"{", RowBox[{"\"\\"", ",", RowBox[{"p3x", "=", RowBox[{"FullSimplify", "[", RowBox[{"TrigExpand", "[", RowBox[{"x", "/.", "x3s"}], "]"}], "]"}]}], ",", RowBox[{"p3y", "=", RowBox[{"FullSimplify", "[", RowBox[{"TrigExpand", "[", RowBox[{ RowBox[{"Pn", "[", RowBox[{"1", ",", "\[Delta]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}], "]"}], "]"}]}]}], "}"}], "\[IndentingNewLine]", "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"x4s", "=", RowBox[{ RowBox[{"FullSimplify", "[", RowBox[{"Solve", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"Pn", "[", RowBox[{"6", ",", "\[Delta]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}], "\[Equal]", RowBox[{"Gn", "[", RowBox[{"4", ",", "\[Delta]"}], "]"}]}], ",", "x"}], "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"{", RowBox[{"\"\\"", ",", RowBox[{"p4x", "=", RowBox[{"FullSimplify", "[", RowBox[{"TrigExpand", "[", RowBox[{"x", "/.", "x4s"}], "]"}], "]"}]}], ",", RowBox[{"p4y", "=", RowBox[{"FullSimplify", "[", RowBox[{"TrigExpand", "[", RowBox[{ RowBox[{"Pn", "[", RowBox[{"6", ",", "\[Delta]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}], "]"}], "]"}]}]}], "}"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"x5s", "=", RowBox[{ RowBox[{"FullSimplify", "[", RowBox[{"Solve", "[", RowBox[{ RowBox[{ RowBox[{"F2", "[", RowBox[{"\[Delta]", ",", "x1"}], "]"}], "\[Equal]", RowBox[{"Gn", "[", RowBox[{"4", ",", "\[Delta]"}], "]"}]}], ",", "x"}], "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"{", RowBox[{"\"\\"", ",", RowBox[{"p5x", "=", RowBox[{"FullSimplify", "[", RowBox[{"x", "/.", "x5s"}], "]"}]}], ",", RowBox[{"p5y", "=", RowBox[{"FullSimplify", "[", RowBox[{ RowBox[{"Gn", "[", RowBox[{"4", ",", "\[Delta]"}], "]"}], "/.", "x5s"}], "]"}]}]}], "}"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"x7s", "=", RowBox[{ RowBox[{"FullSimplify", "[", RowBox[{"Solve", "[", RowBox[{ RowBox[{ RowBox[{"F3", "[", RowBox[{"\[Delta]", ",", "x1"}], "]"}], "\[Equal]", RowBox[{"Gn", "[", RowBox[{"4", ",", "\[Delta]"}], "]"}]}], ",", "x"}], "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"{", RowBox[{"\"\\"", ",", RowBox[{"p7x", "=", RowBox[{"FullSimplify", "[", RowBox[{"x", "/.", "x7s"}], "]"}]}], ",", RowBox[{"p7y", "=", RowBox[{"FullSimplify", "[", RowBox[{ RowBox[{"Gn", "[", RowBox[{"4", ",", "\[Delta]"}], "]"}], "/.", "x7s"}], "]"}]}]}], "}"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.626151504222093*^9, 3.626151668026886*^9}, { 3.626151786334043*^9, 3.626151826872504*^9}, {3.626152163201243*^9, 3.626152226363241*^9}, {3.626152312438203*^9, 3.6261523495117903`*^9}, { 3.6261523873149033`*^9, 3.626152534095221*^9}, {3.626153149175365*^9, 3.626153190575666*^9}, {3.626153329659944*^9, 3.626153423055274*^9}, { 3.626153454871627*^9, 3.6261534554952927`*^9}, {3.6261535045232687`*^9, 3.6261535146665497`*^9}, {3.626153744224627*^9, 3.6261538594526367`*^9}, 3.626154192707402*^9, {3.6261542438550653`*^9, 3.626154247732071*^9}, { 3.626154288430451*^9, 3.626154331264625*^9}, {3.626154381939871*^9, 3.626154382060169*^9}, {3.626154423218289*^9, 3.6261544714537373`*^9}, 3.626154819863147*^9, {3.626154850481331*^9, 3.626154898707108*^9}, { 3.6261550389986553`*^9, 3.6261550919982357`*^9}, {3.626155174451037*^9, 3.626155177957548*^9}, {3.626155208452577*^9, 3.6261552509894543`*^9}, { 3.6261553528332977`*^9, 3.626155375672233*^9}, {3.6261554218486433`*^9, 3.626155421968625*^9}, {3.626155546966175*^9, 3.6261555696218987`*^9}, { 3.626156675367404*^9, 3.6261566809724083`*^9}, 3.626165974108486*^9, { 3.6261665072352123`*^9, 3.6261665079483232`*^9}, {3.626166542071294*^9, 3.6261665469033947`*^9}, {3.626166697867516*^9, 3.6261667544934*^9}, { 3.626167592813059*^9, 3.626167593045122*^9}, {3.626167636694997*^9, 3.6261678382205067`*^9}, {3.62616787882432*^9, 3.626167908087606*^9}, { 3.626167948393189*^9, 3.626167968056625*^9}, {3.6261679988097477`*^9, 3.626168137422888*^9}, {3.6261685636500387`*^9, 3.6261685773539886`*^9}, { 3.626169011439687*^9, 3.626169031727374*^9}, {3.626169098010456*^9, 3.626169115506217*^9}, {3.626169513629236*^9, 3.626169531461501*^9}, { 3.6261695838880672`*^9, 3.626169650064364*^9}, {3.62616968082687*^9, 3.6261697631081667`*^9}, {3.626169822159769*^9, 3.626169898463022*^9}, { 3.626169964451097*^9, 3.626170005211237*^9}, {3.626170040327117*^9, 3.626170057574769*^9}, {3.626170250685051*^9, 3.626170434356526*^9}, { 3.6261705948191633`*^9, 3.6261706020205812`*^9}, {3.626170862347999*^9, 3.626170878583396*^9}, {3.6261709454207153`*^9, 3.6261709485083942`*^9}, { 3.626171214909965*^9, 3.626171247083808*^9}, {3.62617128386655*^9, 3.626171311119772*^9}, {3.6261713630069857`*^9, 3.626171386501504*^9}, { 3.6261714945332947`*^9, 3.62617149720669*^9}, {3.62617154658932*^9, 3.626171631268448*^9}, {3.626172040530222*^9, 3.626172047468039*^9}, { 3.626172241928832*^9, 3.626172309519906*^9}, {3.6261723639099607`*^9, 3.626172390586672*^9}, {3.62617296042419*^9, 3.626172978610218*^9}, { 3.626173022437262*^9, 3.626173022548375*^9}, {3.626173204723126*^9, 3.62617328292845*^9}, {3.626173383511672*^9, 3.626173448959387*^9}, { 3.626173595817882*^9, 3.626173596009038*^9}, {3.626173648791729*^9, 3.62617366984578*^9}, {3.62617431910483*^9, 3.626174337855749*^9}, { 3.6261743689755383`*^9, 3.626174422507181*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"\<\"p1\"\>", ",", FractionBox[ RowBox[{ RowBox[{"4", " ", RowBox[{"Cos", "[", "\[Delta]", "]"}]}], "-", RowBox[{"Cos", "[", RowBox[{"3", " ", "\[Delta]"}], "]"}], "-", RowBox[{"2", " ", SqrtBox["3"], " ", RowBox[{"Sin", "[", "\[Delta]", "]"}]}], "+", RowBox[{"4", " ", "x1", " ", RowBox[{"Sin", "[", RowBox[{"2", " ", "\[Delta]"}], "]"}]}], "+", RowBox[{ SqrtBox["3"], " ", RowBox[{"(", RowBox[{ RowBox[{"2", " ", "x1"}], "+", RowBox[{"Sin", "[", RowBox[{"3", " ", "\[Delta]"}], "]"}]}], ")"}]}]}], RowBox[{"4", " ", SqrtBox["3"]}]], ",", RowBox[{"-", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "7"}], "-", RowBox[{"3", " ", RowBox[{"Cos", "[", RowBox[{"2", " ", "\[Delta]"}], "]"}]}], "+", RowBox[{"Cos", "[", RowBox[{"4", " ", "\[Delta]"}], "]"}], "+", RowBox[{"2", " ", "x1", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"2", " ", RowBox[{"Cos", "[", RowBox[{"2", " ", "\[Delta]"}], "]"}]}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{ SqrtBox["3"], " ", RowBox[{"Cos", "[", "\[Delta]", "]"}]}], "+", RowBox[{"Sin", "[", "\[Delta]", "]"}]}], ")"}]}], "-", RowBox[{"3", " ", SqrtBox["3"], " ", RowBox[{"Sin", "[", RowBox[{"2", " ", "\[Delta]"}], "]"}]}], "+", RowBox[{ SqrtBox["3"], " ", RowBox[{"Sin", "[", RowBox[{"4", " ", "\[Delta]"}], "]"}]}]}], ")"}], "/", RowBox[{"(", RowBox[{"4", " ", RowBox[{"(", RowBox[{ RowBox[{"3", " ", RowBox[{"Cos", "[", "\[Delta]", "]"}]}], "+", RowBox[{ SqrtBox["3"], " ", RowBox[{"Sin", "[", "\[Delta]", "]"}]}]}], ")"}]}], ")"}]}]}]}], "}"}]], "Output", CellChangeTimes->{ 3.626170949105298*^9, {3.6261712176105833`*^9, 3.626171248242259*^9}, 3.626171314675005*^9, 3.626171633235409*^9, 3.626172272709485*^9, 3.626172310740309*^9, 3.6261723954619007`*^9, 3.62617303564985*^9, 3.626173242954562*^9, 3.626173459065835*^9, 3.6261735988138866`*^9, 3.626173672753909*^9, {3.6261743269537783`*^9, 3.626174340546698*^9}, { 3.6261743744408493`*^9, 3.62617442542183*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"\<\"p2\"\>", ",", RowBox[{ FractionBox["1", "6"], " ", RowBox[{"(", RowBox[{ RowBox[{"3", " ", "x1"}], "+", RowBox[{"4", " ", SqrtBox["3"], " ", RowBox[{"Cos", "[", "\[Delta]", "]"}]}], "-", RowBox[{"3", " ", "x1", " ", RowBox[{"Cos", "[", RowBox[{"2", " ", "\[Delta]"}], "]"}]}], "-", RowBox[{ SqrtBox["3"], " ", RowBox[{"Cos", "[", RowBox[{"3", " ", "\[Delta]"}], "]"}]}], "-", RowBox[{"6", " ", RowBox[{"Sin", "[", "\[Delta]", "]"}]}], "+", RowBox[{ SqrtBox["3"], " ", "x1", " ", RowBox[{"Sin", "[", RowBox[{"2", " ", "\[Delta]"}], "]"}]}]}], ")"}]}], ",", RowBox[{ FractionBox["1", "6"], " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", RowBox[{"Cos", "[", "\[Delta]", "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "3"}], "+", RowBox[{ SqrtBox["3"], " ", "x1", " ", RowBox[{"Cos", "[", "\[Delta]", "]"}]}], "+", RowBox[{"3", " ", "x1", " ", RowBox[{"Sin", "[", "\[Delta]", "]"}]}]}], ")"}]}], "-", RowBox[{ SqrtBox["3"], " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", RowBox[{"Sin", "[", "\[Delta]", "]"}]}], "+", RowBox[{"Sin", "[", RowBox[{"3", " ", "\[Delta]"}], "]"}]}], ")"}]}]}], ")"}]}]}], "}"}]], "Output", CellChangeTimes->{ 3.626170949105298*^9, {3.6261712176105833`*^9, 3.626171248242259*^9}, 3.626171314675005*^9, 3.626171633235409*^9, 3.626172272709485*^9, 3.626172310740309*^9, 3.6261723954619007`*^9, 3.62617303564985*^9, 3.626173242954562*^9, 3.626173459065835*^9, 3.6261735988138866`*^9, 3.626173672753909*^9, {3.6261743269537783`*^9, 3.626174340546698*^9}, { 3.6261743744408493`*^9, 3.62617442643155*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"\<\"p3\"\>", ",", RowBox[{ RowBox[{ RowBox[{"-", FractionBox["1", "2"]}], " ", SqrtBox["3"], " ", RowBox[{"Cos", "[", "\[Delta]", "]"}]}], "+", RowBox[{ SqrtBox["3"], " ", RowBox[{"Sec", "[", "\[Delta]", "]"}]}], "-", FractionBox[ RowBox[{"Sin", "[", "\[Delta]", "]"}], "2"]}], ",", RowBox[{ FractionBox["1", "2"], " ", RowBox[{"(", RowBox[{ RowBox[{"Cos", "[", "\[Delta]", "]"}], "-", RowBox[{ SqrtBox["3"], " ", RowBox[{"Sin", "[", "\[Delta]", "]"}]}]}], ")"}]}]}], "}"}]], "Output",\ CellChangeTimes->{ 3.626170949105298*^9, {3.6261712176105833`*^9, 3.626171248242259*^9}, 3.626171314675005*^9, 3.626171633235409*^9, 3.626172272709485*^9, 3.626172310740309*^9, 3.6261723954619007`*^9, 3.62617303564985*^9, 3.626173242954562*^9, 3.626173459065835*^9, 3.6261735988138866`*^9, 3.626173672753909*^9, {3.6261743269537783`*^9, 3.626174340546698*^9}, { 3.6261743744408493`*^9, 3.6261744266788073`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"\<\"p4\"\>", ",", RowBox[{"-", FractionBox[ RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", SqrtBox["3"]}], "+", RowBox[{ SqrtBox["3"], " ", RowBox[{"Cos", "[", RowBox[{"2", " ", "\[Delta]"}], "]"}]}], "+", RowBox[{"Sin", "[", RowBox[{"2", " ", "\[Delta]"}], "]"}]}], RowBox[{"2", " ", RowBox[{"(", RowBox[{ RowBox[{"Cos", "[", "\[Delta]", "]"}], "+", RowBox[{ SqrtBox["3"], " ", RowBox[{"Sin", "[", "\[Delta]", "]"}]}]}], ")"}]}]]}], ",", RowBox[{ FractionBox["1", "2"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{"Cos", "[", "\[Delta]", "]"}]}], "-", RowBox[{ SqrtBox["3"], " ", RowBox[{"Sin", "[", "\[Delta]", "]"}]}]}], ")"}]}]}], "}"}]], "Output",\ CellChangeTimes->{ 3.626170949105298*^9, {3.6261712176105833`*^9, 3.626171248242259*^9}, 3.626171314675005*^9, 3.626171633235409*^9, 3.626172272709485*^9, 3.626172310740309*^9, 3.6261723954619007`*^9, 3.62617303564985*^9, 3.626173242954562*^9, 3.626173459065835*^9, 3.6261735988138866`*^9, 3.626173672753909*^9, {3.6261743269537783`*^9, 3.626174340546698*^9}, { 3.6261743744408493`*^9, 3.626174427322147*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"\<\"p5\"\>", ",", FractionBox[ RowBox[{ RowBox[{"Cos", "[", RowBox[{"3", " ", "\[Delta]"}], "]"}], "+", RowBox[{"3", " ", SqrtBox["3"], " ", RowBox[{"Sin", "[", "\[Delta]", "]"}]}], "+", RowBox[{ RowBox[{"Cos", "[", "\[Delta]", "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"2", " ", SqrtBox["3"], " ", "x1", " ", RowBox[{"Cos", "[", "\[Delta]", "]"}]}], "-", RowBox[{"2", " ", "x1", " ", RowBox[{"Sin", "[", "\[Delta]", "]"}]}]}], ")"}]}]}], RowBox[{"2", " ", SqrtBox["3"]}]], ",", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "4"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "9"}], "+", RowBox[{"4", " ", SqrtBox["3"], " ", "x1", " ", RowBox[{"Cos", "[", "\[Delta]", "]"}]}]}], ")"}], " ", RowBox[{"Cos", "[", RowBox[{"2", " ", "\[Delta]"}], "]"}]}], "+", RowBox[{"Cos", "[", RowBox[{"4", " ", "\[Delta]"}], "]"}], "+", RowBox[{ RowBox[{"(", RowBox[{ SqrtBox["3"], "+", RowBox[{"4", " ", "x1", " ", RowBox[{"Cos", "[", "\[Delta]", "]"}]}]}], ")"}], " ", RowBox[{"Sin", "[", RowBox[{"2", " ", "\[Delta]"}], "]"}]}], "+", RowBox[{ SqrtBox["3"], " ", RowBox[{"Sin", "[", RowBox[{"4", " ", "\[Delta]"}], "]"}]}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"12", " ", RowBox[{"Cos", "[", "\[Delta]", "]"}]}], "-", RowBox[{"4", " ", SqrtBox["3"], " ", RowBox[{"Sin", "[", "\[Delta]", "]"}]}]}], ")"}]}]}], "}"}]], "Output",\ CellChangeTimes->{ 3.626170949105298*^9, {3.6261712176105833`*^9, 3.626171248242259*^9}, 3.626171314675005*^9, 3.626171633235409*^9, 3.626172272709485*^9, 3.626172310740309*^9, 3.6261723954619007`*^9, 3.62617303564985*^9, 3.626173242954562*^9, 3.626173459065835*^9, 3.6261735988138866`*^9, 3.626173672753909*^9, {3.6261743269537783`*^9, 3.626174340546698*^9}, { 3.6261743744408493`*^9, 3.6261744287506723`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"\<\"p7\"\>", ",", RowBox[{ FractionBox["1", "12"], " ", RowBox[{"(", RowBox[{ RowBox[{"6", " ", "x1"}], "+", RowBox[{"4", " ", SqrtBox["3"], " ", RowBox[{"Cos", "[", "\[Delta]", "]"}]}], "-", RowBox[{ SqrtBox["3"], " ", RowBox[{"Cos", "[", RowBox[{"3", " ", "\[Delta]"}], "]"}]}], "+", RowBox[{"6", " ", RowBox[{"Sin", "[", "\[Delta]", "]"}]}], "-", RowBox[{"4", " ", SqrtBox["3"], " ", "x1", " ", RowBox[{"Sin", "[", RowBox[{"2", " ", "\[Delta]"}], "]"}]}], "-", RowBox[{"3", " ", RowBox[{"Sin", "[", RowBox[{"3", " ", "\[Delta]"}], "]"}]}]}], ")"}]}], ",", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "7"}], " ", SqrtBox["3"]}], "-", RowBox[{"3", " ", SqrtBox["3"], " ", RowBox[{"Cos", "[", RowBox[{"2", " ", "\[Delta]"}], "]"}]}], "+", RowBox[{"6", " ", "x1", " ", RowBox[{"Cos", "[", RowBox[{"3", " ", "\[Delta]"}], "]"}]}], "+", RowBox[{ SqrtBox["3"], " ", RowBox[{"Cos", "[", RowBox[{"4", " ", "\[Delta]"}], "]"}]}], "+", RowBox[{"9", " ", RowBox[{"Sin", "[", RowBox[{"2", " ", "\[Delta]"}], "]"}]}], "-", RowBox[{"2", " ", SqrtBox["3"], " ", "x1", " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", RowBox[{"Sin", "[", "\[Delta]", "]"}]}], "+", RowBox[{"Sin", "[", RowBox[{"3", " ", "\[Delta]"}], "]"}]}], ")"}]}], "-", RowBox[{"3", " ", RowBox[{"Sin", "[", RowBox[{"4", " ", "\[Delta]"}], "]"}]}]}], ")"}], "/", RowBox[{"(", RowBox[{"12", " ", RowBox[{"(", RowBox[{ RowBox[{ SqrtBox["3"], " ", RowBox[{"Cos", "[", "\[Delta]", "]"}]}], "-", RowBox[{"Sin", "[", "\[Delta]", "]"}]}], ")"}]}], ")"}]}]}], "}"}]], "Output", CellChangeTimes->{ 3.626170949105298*^9, {3.6261712176105833`*^9, 3.626171248242259*^9}, 3.626171314675005*^9, 3.626171633235409*^9, 3.626172272709485*^9, 3.626172310740309*^9, 3.6261723954619007`*^9, 3.62617303564985*^9, 3.626173242954562*^9, 3.626173459065835*^9, 3.6261735988138866`*^9, 3.626173672753909*^9, {3.6261743269537783`*^9, 3.626174340546698*^9}, { 3.6261743744408493`*^9, 3.6261744306354513`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell["Test coordinates", "Subsubtitle", CellChangeTimes->{{3.626195928740439*^9, 3.6261959351128473`*^9}, { 3.626200794425498*^9, 3.626200799891798*^9}, {3.626333397897748*^9, 3.626333419505825*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"Manipulate", "[", RowBox[{ RowBox[{"Show", "[", RowBox[{ RowBox[{"Plot", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Gn", "[", RowBox[{"n", ",", "\[Delta]\[Delta]"}], "]"}], ",", RowBox[{"{", RowBox[{"n", ",", "1", ",", "6"}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"x", ",", RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "1.2"}], ",", "1.2"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1.2"}], ",", "1.2"}], "}"}]}], "}"}]}], ",", RowBox[{"AspectRatio", "\[Rule]", "1"}]}], "]"}], ",", RowBox[{"Plot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"F1", "[", RowBox[{"\[Delta]\[Delta]", ",", "x11"}], "]"}], ",", RowBox[{"F2", "[", RowBox[{"\[Delta]\[Delta]", ",", "x11"}], "]"}], ",", RowBox[{"F3", "[", RowBox[{"\[Delta]\[Delta]", ",", "x11"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"x", ",", RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", " ", RowBox[{"{", RowBox[{"Red", ",", "Green", ",", "Orange"}], "}"}]}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"p1x", "/.", RowBox[{"{", RowBox[{ RowBox[{"\[Delta]", "\[Rule]", "\[Delta]\[Delta]"}], ",", RowBox[{"x1", "\[Rule]", "x11"}]}], "}"}]}], ",", RowBox[{"p1y", "/.", RowBox[{"{", RowBox[{ RowBox[{"\[Delta]", "\[Rule]", "\[Delta]\[Delta]"}], ",", RowBox[{"x1", "\[Rule]", "x11"}]}], "}"}]}]}], "}"}], "}"}], ",", RowBox[{"PlotMarkers", "\[Rule]", RowBox[{"{", "\"\\"", "}"}]}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"p2x", "/.", RowBox[{"{", RowBox[{ RowBox[{"\[Delta]", "\[Rule]", "\[Delta]\[Delta]"}], ",", RowBox[{"x1", "\[Rule]", "x11"}]}], "}"}]}], ",", RowBox[{"p2y", "/.", RowBox[{"{", RowBox[{ RowBox[{"\[Delta]", "\[Rule]", "\[Delta]\[Delta]"}], ",", RowBox[{"x1", "\[Rule]", "x11"}]}], "}"}]}]}], "}"}], "}"}], ",", RowBox[{"PlotMarkers", "\[Rule]", RowBox[{"{", "\"\\"", "}"}]}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"p3x", "/.", RowBox[{"{", RowBox[{ RowBox[{"\[Delta]", "\[Rule]", "\[Delta]\[Delta]"}], ",", RowBox[{"x1", "\[Rule]", "x11"}]}], "}"}]}], ",", RowBox[{"p3y", "/.", RowBox[{"{", RowBox[{ RowBox[{"\[Delta]", "\[Rule]", "\[Delta]\[Delta]"}], ",", RowBox[{"x1", "\[Rule]", "x11"}]}], "}"}]}]}], "}"}], "}"}], ",", RowBox[{"PlotMarkers", "\[Rule]", RowBox[{"{", "\"\\"", "}"}]}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"p4x", "/.", RowBox[{"{", RowBox[{ RowBox[{"\[Delta]", "\[Rule]", "\[Delta]\[Delta]"}], ",", RowBox[{"x1", "\[Rule]", "x11"}]}], "}"}]}], ",", RowBox[{"p4y", "/.", RowBox[{"{", RowBox[{ RowBox[{"\[Delta]", "\[Rule]", "\[Delta]\[Delta]"}], ",", RowBox[{"x1", "\[Rule]", "x11"}]}], "}"}]}]}], "}"}], "}"}], ",", RowBox[{"PlotMarkers", "\[Rule]", RowBox[{"{", "\"\\"", "}"}]}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"p5x", "/.", RowBox[{"{", RowBox[{ RowBox[{"\[Delta]", "\[Rule]", "\[Delta]\[Delta]"}], ",", RowBox[{"x1", "\[Rule]", "x11"}]}], "}"}]}], ",", RowBox[{"p5y", "/.", RowBox[{"{", RowBox[{ RowBox[{"\[Delta]", "\[Rule]", "\[Delta]\[Delta]"}], ",", RowBox[{"x1", "\[Rule]", "x11"}]}], "}"}]}]}], "}"}], "}"}], ",", RowBox[{"PlotMarkers", "\[Rule]", RowBox[{"{", "\"\\"", "}"}]}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"p7x", "/.", RowBox[{"{", RowBox[{ RowBox[{"\[Delta]", "\[Rule]", "\[Delta]\[Delta]"}], ",", RowBox[{"x1", "\[Rule]", "x11"}]}], "}"}]}], ",", RowBox[{"p7y", "/.", RowBox[{"{", RowBox[{ RowBox[{"\[Delta]", "\[Rule]", "\[Delta]\[Delta]"}], ",", RowBox[{"x1", "\[Rule]", "x11"}]}], "}"}]}]}], "}"}], "}"}], ",", RowBox[{"PlotMarkers", "\[Rule]", RowBox[{"{", "\"\\"", "}"}]}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"{", RowBox[{"Pn", "[", RowBox[{"1", ",", "\[Delta]\[Delta]"}], "]"}], "}"}], ",", RowBox[{"PlotMarkers", "\[Rule]", RowBox[{"{", RowBox[{"Style", "[", RowBox[{"\"\<1\>\"", ",", "Large", ",", "Bold", ",", "Black"}], "]"}], "}"}]}]}], "]"}], ",", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"{", RowBox[{"Pn", "[", RowBox[{"2", ",", "\[Delta]\[Delta]"}], "]"}], "}"}], ",", RowBox[{"PlotMarkers", "\[Rule]", RowBox[{"{", RowBox[{"Style", "[", RowBox[{"\"\<2\>\"", ",", "Large", ",", "Bold", ",", "Black"}], "]"}], "}"}]}]}], "]"}], ",", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"{", RowBox[{"Pn", "[", RowBox[{"3", ",", "\[Delta]\[Delta]"}], "]"}], "}"}], ",", RowBox[{"PlotMarkers", "\[Rule]", RowBox[{"{", RowBox[{"Style", "[", RowBox[{"\"\<3\>\"", ",", "Large", ",", "Bold", ",", "Black"}], "]"}], "}"}]}]}], "]"}], ",", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"{", RowBox[{"Pn", "[", RowBox[{"4", ",", "\[Delta]\[Delta]"}], "]"}], "}"}], ",", RowBox[{"PlotMarkers", "\[Rule]", RowBox[{"{", RowBox[{"Style", "[", RowBox[{"\"\<4\>\"", ",", "Large", ",", "Bold", ",", "Black"}], "]"}], "}"}]}]}], "]"}], ",", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"{", RowBox[{"Pn", "[", RowBox[{"5", ",", "\[Delta]\[Delta]"}], "]"}], "}"}], ",", RowBox[{"PlotMarkers", "\[Rule]", RowBox[{"{", RowBox[{"Style", "[", RowBox[{"\"\<5\>\"", ",", "Large", ",", "Bold", ",", "Black"}], "]"}], "}"}]}]}], "]"}], ",", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"{", RowBox[{"Pn", "[", RowBox[{"6", ",", "\[Delta]\[Delta]"}], "]"}], "}"}], ",", RowBox[{"PlotMarkers", "\[Rule]", RowBox[{"{", RowBox[{"Style", "[", RowBox[{"\"\<6\>\"", ",", "Large", ",", "Bold", ",", "Black"}], "]"}], "}"}]}]}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"\[Delta]\[Delta]", ",", "0.0001", ",", RowBox[{"\[Pi]", "/", "6"}]}], "}"}], ",", RowBox[{"{", RowBox[{"x11", ",", RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "]"}], "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{{3.626170350295298*^9, 3.626170350309078*^9}, { 3.626170448615747*^9, 3.6261704720560303`*^9}, {3.62617056203384*^9, 3.6261705629061203`*^9}, {3.626170733000451*^9, 3.6261707347568893`*^9}, { 3.6261708167936993`*^9, 3.626170849464717*^9}, {3.626170890083704*^9, 3.626170899611965*^9}, {3.626170956181273*^9, 3.626170971589892*^9}, { 3.626171066866288*^9, 3.626171071973748*^9}, {3.626171108061475*^9, 3.626171169412956*^9}, {3.626171323255369*^9, 3.626171335680759*^9}, { 3.626171646105955*^9, 3.626171651482892*^9}, {3.626171857992269*^9, 3.626171858447365*^9}, {3.626172407513191*^9, 3.626172453991397*^9}, { 3.626173040905837*^9, 3.626173046098222*^9}, {3.6261734654132*^9, 3.62617350014363*^9}}], Cell[BoxData[ TagBox[ StyleBox[ DynamicModuleBox[{$CellContext`x11$$ = -0.040000000000000036`, \ $CellContext`\[Delta]\[Delta]$$ = 0.20426452248333657`, Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{ Hold[$CellContext`\[Delta]\[Delta]$$], 0.0001, Rational[1, 6] Pi}, { Hold[$CellContext`x11$$], -1, 1}}, Typeset`size$$ = { 360., {178., 182.}}, Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ = True, $CellContext`\[Delta]\[Delta]$237639$$ = 0, $CellContext`x11$237640$$ = 0}, DynamicBox[Manipulate`ManipulateBoxes[ 1, StandardForm, "Variables" :> {$CellContext`x11$$ = -1, \ $CellContext`\[Delta]\[Delta]$$ = 0.0001}, "ControllerVariables" :> { Hold[$CellContext`\[Delta]\[Delta]$$, \ $CellContext`\[Delta]\[Delta]$237639$$, 0], Hold[$CellContext`x11$$, $CellContext`x11$237640$$, 0]}, "OtherVariables" :> { Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$, Typeset`skipInitDone$$}, "Body" :> Show[ Plot[ Table[ $CellContext`Gn[$CellContext`n, $CellContext`\[Delta]\[Delta]$$], \ {$CellContext`n, 1, 6}], {$CellContext`x, -1, 1}, PlotRange -> {{-1.2, 1.2}, {-1.2, 1.2}}, AspectRatio -> 1], Plot[{ $CellContext`F1[$CellContext`\[Delta]\[Delta]$$, \ $CellContext`x11$$], $CellContext`F2[$CellContext`\[Delta]\[Delta]$$, \ $CellContext`x11$$], $CellContext`F3[$CellContext`\[Delta]\[Delta]$$, \ $CellContext`x11$$]}, {$CellContext`x, -1, 1}, PlotStyle -> {Red, Green, Orange}], ListPlot[{{ ReplaceAll[$CellContext`p1x, {$CellContext`\[Delta] -> \ $CellContext`\[Delta]\[Delta]$$, $CellContext`x1 -> $CellContext`x11$$}], ReplaceAll[$CellContext`p1y, {$CellContext`\[Delta] -> \ $CellContext`\[Delta]\[Delta]$$, $CellContext`x1 -> $CellContext`x11$$}]}}, PlotMarkers -> {"p1"}], ListPlot[{{ ReplaceAll[$CellContext`p2x, {$CellContext`\[Delta] -> \ $CellContext`\[Delta]\[Delta]$$, $CellContext`x1 -> $CellContext`x11$$}], ReplaceAll[$CellContext`p2y, {$CellContext`\[Delta] -> \ $CellContext`\[Delta]\[Delta]$$, $CellContext`x1 -> $CellContext`x11$$}]}}, PlotMarkers -> {"p2"}], ListPlot[{{ ReplaceAll[$CellContext`p3x, {$CellContext`\[Delta] -> \ $CellContext`\[Delta]\[Delta]$$, $CellContext`x1 -> $CellContext`x11$$}], ReplaceAll[$CellContext`p3y, {$CellContext`\[Delta] -> \ $CellContext`\[Delta]\[Delta]$$, $CellContext`x1 -> $CellContext`x11$$}]}}, PlotMarkers -> {"p3"}], ListPlot[{{ ReplaceAll[$CellContext`p4x, {$CellContext`\[Delta] -> \ $CellContext`\[Delta]\[Delta]$$, $CellContext`x1 -> $CellContext`x11$$}], ReplaceAll[$CellContext`p4y, {$CellContext`\[Delta] -> \ $CellContext`\[Delta]\[Delta]$$, $CellContext`x1 -> $CellContext`x11$$}]}}, PlotMarkers -> {"p4"}], ListPlot[{{ ReplaceAll[$CellContext`p5x, {$CellContext`\[Delta] -> \ $CellContext`\[Delta]\[Delta]$$, $CellContext`x1 -> $CellContext`x11$$}], ReplaceAll[$CellContext`p5y, {$CellContext`\[Delta] -> \ $CellContext`\[Delta]\[Delta]$$, $CellContext`x1 -> $CellContext`x11$$}]}}, PlotMarkers -> {"p5"}], ListPlot[{{ ReplaceAll[$CellContext`p7x, {$CellContext`\[Delta] -> \ $CellContext`\[Delta]\[Delta]$$, $CellContext`x1 -> $CellContext`x11$$}], ReplaceAll[$CellContext`p7y, {$CellContext`\[Delta] -> \ $CellContext`\[Delta]\[Delta]$$, $CellContext`x1 -> $CellContext`x11$$}]}}, PlotMarkers -> {"p7"}], ListPlot[{ $CellContext`Pn[1, $CellContext`\[Delta]\[Delta]$$]}, PlotMarkers -> { Style["1", Large, Bold, Black]}], ListPlot[{ $CellContext`Pn[2, $CellContext`\[Delta]\[Delta]$$]}, PlotMarkers -> { Style["2", Large, Bold, Black]}], ListPlot[{ $CellContext`Pn[3, $CellContext`\[Delta]\[Delta]$$]}, PlotMarkers -> { Style["3", Large, Bold, Black]}], ListPlot[{ $CellContext`Pn[4, $CellContext`\[Delta]\[Delta]$$]}, PlotMarkers -> { Style["4", Large, Bold, Black]}], ListPlot[{ $CellContext`Pn[5, $CellContext`\[Delta]\[Delta]$$]}, PlotMarkers -> { Style["5", Large, Bold, Black]}], ListPlot[{ $CellContext`Pn[6, $CellContext`\[Delta]\[Delta]$$]}, PlotMarkers -> { Style["6", Large, Bold, Black]}]], "Specifications" :> {{$CellContext`\[Delta]\[Delta]$$, 0.0001, Rational[1, 6] Pi}, {$CellContext`x11$$, -1, 1}}, "Options" :> {}, "DefaultOptions" :> {}], ImageSizeCache->{405., {236., 241.}}, SingleEvaluation->True], Deinitialization:>None, DynamicModuleValues:>{}, SynchronousInitialization->True, UnsavedVariables:>{Typeset`initDone$$}, UntrackedVariables:>{Typeset`size$$}], "Manipulate", Deployed->True, StripOnInput->False], Manipulate`InterpretManipulate[1]]], "Output", CellChangeTimes->{ 3.626170478689666*^9, 3.6261705698030853`*^9, 3.626170639915949*^9, 3.626170741559943*^9, {3.626170809247448*^9, 3.62617085008678*^9}, { 3.626170886915401*^9, 3.626170900076673*^9}, 3.6261709718797817`*^9, 3.626171072655197*^9, {3.6261711579587584`*^9, 3.626171170210498*^9}, 3.6261713360969467`*^9, 3.6261716520567703`*^9, 3.626171859143468*^9, { 3.6261724385557547`*^9, 3.626172456457045*^9}, 3.626173046502178*^9, 3.626173500577083*^9, 3.626173606167823*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Area of a convex polygon", "Subsubtitle", CellChangeTimes->{{3.626195928740439*^9, 3.6261959351128473`*^9}, { 3.626200794425498*^9, 3.626200799891798*^9}, {3.626333397897748*^9, 3.626333402217677*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"AreaPolygon", "[", "VertexList_", "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"VL", "=", "VertexList"}], ",", "nmax", ",", "Area"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"nmax", "=", RowBox[{"Length", "[", "VL", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Area", "=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ RowBox[{"VL", "[", RowBox[{"[", "nmax", "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], RowBox[{ RowBox[{"VL", "[", RowBox[{"[", "1", "]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}]}], "-", RowBox[{ RowBox[{ RowBox[{"VL", "[", RowBox[{"[", "nmax", "]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}], RowBox[{ RowBox[{"VL", "[", RowBox[{"[", "1", "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], "+", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"VL", "[", RowBox[{"[", "i", "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], RowBox[{ RowBox[{"VL", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}]}], "-", RowBox[{ RowBox[{ RowBox[{"VL", "[", RowBox[{"[", "i", "]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}], RowBox[{ RowBox[{"VL", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"nmax", "-", "1"}]}], "}"}]}], "]"}]}], ")"}], "/", "2"}]}], ";", "\[IndentingNewLine]", "Area"}]}], "\[IndentingNewLine]", "]"}]}], "\[IndentingNewLine]", RowBox[{"AreaPolygon", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "2"}], "}"}]}], "}"}], "]"}]}], "Input", CellChangeTimes->{{3.626156701676239*^9, 3.626156943264291*^9}, { 3.626156999066063*^9, 3.626157024976871*^9}, {3.626157069204137*^9, 3.626157075290553*^9}, {3.626157127622734*^9, 3.626157130877825*^9}, { 3.626157166808185*^9, 3.626157182829893*^9}, {3.626165996512331*^9, 3.626166166677297*^9}, {3.6261662714815893`*^9, 3.626166347066684*^9}, { 3.626170786540167*^9, 3.626170787566409*^9}, {3.6261885245947638`*^9, 3.626188577742188*^9}, {3.62633338610756*^9, 3.626333386450478*^9}}], Cell[BoxData[ FractionBox["3", "2"]], "Output", CellChangeTimes->{3.626333386932111*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Blue Spot, delfection is \[Pi]-2\[Phi]11\ \>", "Subsubtitle", CellChangeTimes->{{3.626195928740439*^9, 3.6261959351128473`*^9}, { 3.626200794425498*^9, 3.626200799891798*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"\[Pi]", "*", RowBox[{"35", "/", "180.0"}]}], "\[IndentingNewLine]", RowBox[{"\[Pi]", "*", RowBox[{"25", "/", "180.0"}]}], "\[IndentingNewLine]", RowBox[{"\[Pi]", "*", RowBox[{"10", "/", "180.0"}]}]}], "Input", CellChangeTimes->{{3.62620457632329*^9, 3.6262045857259493`*^9}, { 3.6262054117181473`*^9, 3.626205413974619*^9}, {3.626206657164296*^9, 3.6262066600670033`*^9}}], Cell[BoxData["0.6108652381980153`"], "Output", CellChangeTimes->{{3.6262045836654778`*^9, 3.6262045860343103`*^9}, 3.626205414331036*^9, 3.626206660520076*^9}], Cell[BoxData["0.4363323129985824`"], "Output", CellChangeTimes->{{3.6262045836654778`*^9, 3.6262045860343103`*^9}, 3.626205414331036*^9, 3.626206660521532*^9}], Cell[BoxData["0.17453292519943298`"], "Output", CellChangeTimes->{{3.6262045836654778`*^9, 3.6262045860343103`*^9}, 3.626205414331036*^9, 3.626206660522414*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Manipulate", "[", RowBox[{ RowBox[{"Plot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"ArcSin", "[", FractionBox["1", "1.31"], "]"}], ",", RowBox[{"ArcCos", "[", RowBox[{ FractionBox[ RowBox[{"Cos", "[", "e", "]"}], "1.31"], " ", RowBox[{"Cos", "[", RowBox[{ FractionBox["\[Pi]", "3"], "-", "\[Delta]"}], "]"}]}], "]"}], ",", RowBox[{"ArcCos", "[", RowBox[{ FractionBox[ RowBox[{"Cos", "[", "e", "]"}], "1.31"], " ", RowBox[{"Cos", "[", "\[Delta]", "]"}]}], "]"}], ",", RowBox[{"ArcCos", "[", RowBox[{ FractionBox[ RowBox[{"Cos", "[", "e", "]"}], "1.31"], " ", RowBox[{"Cos", "[", RowBox[{ FractionBox["\[Pi]", "3"], "+", "\[Delta]"}], "]"}]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"\[Delta]", ",", "0", ",", RowBox[{"\[Pi]", "/", "6"}]}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{"Gray", ",", "Black", ",", "Blue", ",", "Red"}], "}"}]}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{"0.5", ",", "1.6"}], "}"}]}]}], "]"}], ",", RowBox[{"{", RowBox[{"e", ",", "0.001", ",", RowBox[{"\[Pi]", "/", "2"}]}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.626067425293017*^9, 3.626067484237707*^9}, { 3.62620420470054*^9, 3.6262042468291273`*^9}, {3.626204305983295*^9, 3.626204321943532*^9}, {3.626204503497703*^9, 3.626204534738824*^9}}], Cell[BoxData[ TagBox[ StyleBox[ DynamicModuleBox[{$CellContext`e$$ = 0.6571748646002668, Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{ Hold[$CellContext`e$$], 0.001, Rational[1, 2] Pi}}, Typeset`size$$ = { 360., {115., 120.}}, Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ = True, $CellContext`e$326687$$ = 0}, DynamicBox[Manipulate`ManipulateBoxes[ 1, StandardForm, "Variables" :> {$CellContext`e$$ = 0.001}, "ControllerVariables" :> { Hold[$CellContext`e$$, $CellContext`e$326687$$, 0]}, "OtherVariables" :> { Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$, Typeset`skipInitDone$$}, "Body" :> Plot[{ ArcSin[1/1.31], ArcCos[(Cos[$CellContext`e$$]/1.31) Cos[Pi/3 - $CellContext`\[Delta]]], ArcCos[(Cos[$CellContext`e$$]/1.31) Cos[$CellContext`\[Delta]]], ArcCos[(Cos[$CellContext`e$$]/1.31) Cos[Pi/3 + $CellContext`\[Delta]]]}, {$CellContext`\[Delta], 0, Pi/ 6}, PlotStyle -> {Gray, Black, Blue, Red}, PlotRange -> {0.5, 1.6}], "Specifications" :> {{$CellContext`e$$, 0.001, Rational[1, 2] Pi}}, "Options" :> {}, "DefaultOptions" :> {}], ImageSizeCache->{405., {159., 164.}}, SingleEvaluation->True], Deinitialization:>None, DynamicModuleValues:>{}, SynchronousInitialization->True, UnsavedVariables:>{Typeset`initDone$$}, UntrackedVariables:>{Typeset`size$$}], "Manipulate", Deployed->True, StripOnInput->False], Manipulate`InterpretManipulate[1]]], "Output", CellChangeTimes->{{3.626067464914267*^9, 3.626067484602378*^9}, 3.626204248973016*^9, 3.626204322702071*^9, {3.626204519700961*^9, 3.6262045352398*^9}}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Incident angle for internal reflection", "Subsubtitle", CellChangeTimes->{{3.626333297696978*^9, 3.626333309123588*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"normalA1", "=", RowBox[{"{", RowBox[{ RowBox[{"-", RowBox[{"Sin", "[", "\[Phi]11", "]"}]}], ",", RowBox[{"Cos", "[", "\[Phi]11", "]"}], ",", "0"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"normalA2", "=", RowBox[{"{", RowBox[{ RowBox[{"Sin", "[", "\[Phi]31", "]"}], ",", RowBox[{"Cos", "[", "\[Phi]31", "]"}], ",", "0"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"normalA3", "=", RowBox[{"{", RowBox[{ RowBox[{"Sin", "[", "\[Phi]41", "]"}], ",", RowBox[{"Cos", "[", "\[Phi]41", "]"}], ",", "0"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"inc", "=", RowBox[{"{", RowBox[{"0", ",", RowBox[{"Cos", "[", "ep", "]"}], ",", RowBox[{"-", RowBox[{"Sin", "[", "ep", "]"}]}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Cos", "[", "\[Alpha]A1", "]"}], "==", RowBox[{"normalA1", ".", "inc"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Cos", "[", "\[Alpha]A2", "]"}], "==", RowBox[{"normalA2", ".", "inc"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Cos", "[", "\[Alpha]A3", "]"}], "==", RowBox[{"normalA3", ".", "inc"}]}]}], "Input", CellChangeTimes->{{3.626067768578032*^9, 3.626067768694241*^9}, { 3.6260682729450083`*^9, 3.626068317762583*^9}, {3.626068417887683*^9, 3.626068438017857*^9}, 3.626068558337459*^9, {3.6260686083160973`*^9, 3.626068612418909*^9}, {3.626195729716609*^9, 3.6261957313871737`*^9}, { 3.626196438546402*^9, 3.626196476825663*^9}, {3.626196524277553*^9, 3.626196544566341*^9}, 3.626198937511675*^9, 3.626199127137103*^9, { 3.626204881534009*^9, 3.626204884072157*^9}, {3.626204978668023*^9, 3.626205019835634*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"Cos", "[", "\[Alpha]A1", "]"}], "\[Equal]", RowBox[{ RowBox[{"Cos", "[", "ep", "]"}], " ", RowBox[{"Cos", "[", "\[Phi]11", "]"}]}]}]], "Output", CellChangeTimes->{ 3.626068438388665*^9, 3.6260685586735487`*^9, 3.626068613354025*^9, 3.626195731856443*^9, 3.626196450513282*^9, {3.6261965276463823`*^9, 3.62619654510605*^9}, {3.6261989463615522`*^9, 3.626198948163027*^9}, 3.6261991385215597`*^9, {3.626204979061944*^9, 3.6262050202416277`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"Cos", "[", "\[Alpha]A2", "]"}], "\[Equal]", RowBox[{ RowBox[{"Cos", "[", "ep", "]"}], " ", RowBox[{"Cos", "[", "\[Phi]31", "]"}]}]}]], "Output", CellChangeTimes->{ 3.626068438388665*^9, 3.6260685586735487`*^9, 3.626068613354025*^9, 3.626195731856443*^9, 3.626196450513282*^9, {3.6261965276463823`*^9, 3.62619654510605*^9}, {3.6261989463615522`*^9, 3.626198948163027*^9}, 3.6261991385215597`*^9, {3.626204979061944*^9, 3.62620502024284*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"Cos", "[", "\[Alpha]A3", "]"}], "\[Equal]", RowBox[{ RowBox[{"Cos", "[", "ep", "]"}], " ", RowBox[{"Cos", "[", "\[Phi]41", "]"}]}]}]], "Output", CellChangeTimes->{ 3.626068438388665*^9, 3.6260685586735487`*^9, 3.626068613354025*^9, 3.626195731856443*^9, 3.626196450513282*^9, {3.6261965276463823`*^9, 3.62619654510605*^9}, {3.6261989463615522`*^9, 3.626198948163027*^9}, 3.6261991385215597`*^9, {3.626204979061944*^9, 3.626205020243781*^9}}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["derive Blue spot azimuth", "Subsubtitle", CellChangeTimes->{{3.626333297696978*^9, 3.626333320376446*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"ArcCos", "[", RowBox[{ FractionBox[ RowBox[{"Cos", "[", "e", "]"}], "n"], " ", RowBox[{"Cos", "[", "\[Phi]11", "]"}]}], "]"}], "\[Equal]", RowBox[{"ArcSin", "[", FractionBox["1", "n"], "]"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ FractionBox[ RowBox[{"Cos", "[", "\[Alpha]TIR", "]"}], RowBox[{"Cos", "[", "e", "]"}]], "n"}], "\[Equal]", " ", RowBox[{"Cos", "[", "\[Phi]11", "]"}]}], "\[IndentingNewLine]", RowBox[{"FullSimplify", "[", RowBox[{"\[Pi]", "-", RowBox[{"2", RowBox[{"ArcCos", "[", RowBox[{ FractionBox[ RowBox[{"Cos", "[", "\[Alpha]TIR", "]"}], RowBox[{"Cos", "[", "e", "]"}]], "n"}], "]"}]}]}], "]"}]}], "Input", CellChangeTimes->{{3.626067534365548*^9, 3.626067552641326*^9}, { 3.626067624624234*^9, 3.626067656169186*^9}, {3.626068623732922*^9, 3.62606875051991*^9}, {3.626195749961978*^9, 3.626195750313542*^9}, 3.626195880014627*^9, {3.6261960698863373`*^9, 3.626196110102858*^9}, { 3.626333323210884*^9, 3.6263333238972187`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"ArcCos", "[", FractionBox[ RowBox[{ RowBox[{"Cos", "[", "e", "]"}], " ", RowBox[{"Cos", "[", "\[Phi]11", "]"}]}], "n"], "]"}], "\[Equal]", RowBox[{"ArcSin", "[", FractionBox["1", "n"], "]"}]}]], "Output", CellChangeTimes->{{3.626067543311644*^9, 3.626067553092992*^9}, { 3.6260676378115473`*^9, 3.626067656610134*^9}, {3.6260686356774187`*^9, 3.62606875102989*^9}, 3.6261957508083076`*^9, 3.626195880962776*^9, { 3.6261960765330687`*^9, 3.6261961104684143`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"n", " ", RowBox[{"Cos", "[", "\[Alpha]TIR", "]"}], " ", RowBox[{"Sec", "[", "e", "]"}]}], "\[Equal]", RowBox[{"Cos", "[", "\[Phi]11", "]"}]}]], "Output", CellChangeTimes->{{3.626067543311644*^9, 3.626067553092992*^9}, { 3.6260676378115473`*^9, 3.626067656610134*^9}, {3.6260686356774187`*^9, 3.62606875102989*^9}, 3.6261957508083076`*^9, 3.626195880962776*^9, { 3.6261960765330687`*^9, 3.626196110469441*^9}}], Cell[BoxData[ RowBox[{"2", " ", RowBox[{"ArcSin", "[", RowBox[{"n", " ", RowBox[{"Cos", "[", "\[Alpha]TIR", "]"}], " ", RowBox[{"Sec", "[", "e", "]"}]}], "]"}]}]], "Output", CellChangeTimes->{{3.626067543311644*^9, 3.626067553092992*^9}, { 3.6260676378115473`*^9, 3.626067656610134*^9}, {3.6260686356774187`*^9, 3.62606875102989*^9}, 3.6261957508083076`*^9, 3.626195880962776*^9, { 3.6261960765330687`*^9, 3.626196110741508*^9}}] }, Open ]] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[TextData[{ "Blue Spot Azimuth\n", StyleBox["M. Sillanp \:0308a \:0308a, J. Moilanen, M. Riikonen, and M. \ Pekkola, \[OpenCurlyDoubleQuote]Blue spot on the parhelic circle,\ \[CloseCurlyDoubleQuote] Appl. Opt. 40, 5275\[Dash]5279 (2001).", FontSize->14] }], "Subtitle", CellChangeTimes->{{3.626333230405569*^9, 3.626333241070181*^9}, { 3.62633327624894*^9, 3.6263332833819733`*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"nr", "=", "1.307"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"nb", "=", "1.317"}], ";"}], "\[IndentingNewLine]", RowBox[{"Plot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"2", RowBox[{"ArcSin", "[", RowBox[{"nr", " ", RowBox[{ RowBox[{"Cos", "[", RowBox[{"ArcSin", "[", FractionBox["1", "nr"], "]"}], "]"}], "/", RowBox[{"Cos", "[", RowBox[{"e", "*", RowBox[{"\[Pi]", "/", "180"}]}], "]"}]}]}], "]"}], "*", RowBox[{"180", "/", "\[Pi]"}]}], ",", RowBox[{"2", RowBox[{"ArcSin", "[", RowBox[{"nb", " ", RowBox[{ RowBox[{"Cos", "[", RowBox[{"ArcSin", "[", FractionBox["1", "nb"], "]"}], "]"}], "/", RowBox[{"Cos", "[", RowBox[{"e", "*", RowBox[{"\[Pi]", "/", "180"}]}], "]"}]}]}], "]"}], "*", RowBox[{"180", "/", "\[Pi]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"e", ",", "0", ",", "35"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{"Red", ",", "Blue"}], "}"}]}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{"All", ",", RowBox[{"{", RowBox[{"110", ",", "180"}], "}"}]}], "}"}]}]}], "]"}]}], "Input", CellChangeTimes->{{3.626064689740572*^9, 3.626064762136664*^9}, { 3.6260648789318857`*^9, 3.626064948838785*^9}, {3.626064980439773*^9, 3.626064986687381*^9}, {3.626067180197136*^9, 3.6260671806850567`*^9}}], Cell[BoxData[ GraphicsBox[{{}, {}, {RGBColor[1, 0, 0], LineBox[CompressedData[" 1:eJwV1Hk4llkfB3B7RG9olK2QNUyISuI5X7IvFRNZHnp4nkcZLVrfsZQKpQ1j iUQYiazFkHiRSRQqu0ZGVJZskURS8x5/3Nd9fa5zf8/53dfvnKPkc8SJy8fD w1NHn6V348brg7lz2YwuleJcy2wm7B+m+NgvRpCyNd0Plqw4aXlcfjGBnBC1 eLPk5zVp7iKn04kW76KiFfVl4ZRggcUsUrHwLWTJVo5JKf8G5ZEHC+YzS+Zz aPT6HnKPBAp0h1lTs45afwwyLiGy8sUbbahrEupD57+VEj/Lrpklr6swFz9V WU6+R+5ot6UO6fsrfSaokqiOfm22o+7hM9U7alRNxgK+9ttTp8S9Yc2FPCJK qmYrd1I3Za3bUNZbSxxWduzdRb3wwHP6hPFjUmFYVLGb2q339blP354Q7WOm L3+hjpyUtb3v0UCCkubPOlOX87hLBlQ+JcmaX+33Uq9WffXHRFATEW7o0veg tti6xj+/p5lETZTs9KQ+YeOi72/0griU94Tto/6wZdT7c0gLcQoWNOFQ8zwS qW/XbiWlCSK1vtTSNhs0S3pbSW7S3n1+1NYeBz4dNW4ng7dbRw9T3wkdOv/x Wyd5vAYmwdTVwkLDL/K6yHBxQtwZ6s7fVe0KPbrJ17ebBc5TC9zmrDpU+Yps 9XUxuETNfvo2czToNXnSL9uSRB3kyCfcqNlLegSfzqdQx/YoHbzb00ue3Xy0 PYP6r3GWwQGjPrLzerx4LnVM5JG3n0L6ieieptn/Ucs7hMa/7egn05HdU7XU ORLRlm3aA2STka5Qw1L/kgtz7/cOkKyshtBW6sn8iWMBxu+ITmds4TC1Q6s/ /+S3QWJYVL0gncNE16CXxQ6PIRK2PFRTkdprwfFiUsUQeUy0AtWpA5QNRXcE DpM+PevjW6hjTwqsSvoyQtTNbvY4L+VlUtebTY8RnYGBk8lL+Y0xnMTd40TC NDQ+k3rYLOzOeNE4MflbrC2fet7fTzPx8ASRm/x6t4ZarnrzpvGxSWIXofJ8 iFpXb1jRbmaK7MgWHDC6S/tZ2V87u36asPKfT1pQl1v0+KQ7TRPl3VfkHKlN 3Z/fnr0/TfZMj7Tup3YKK9FID/hEos0Ni68v5TtDdWcnZkjB+1DTBeqHgdJI G/lC2EpFqa9ymYhqsd+uJDVHQmbrrw5R+6if25JpOkck34xc/0wt1vlBOzuZ OsF+UTyPnhfdSpki+3nCq5kmakctPOw5U33vK5EIaTGuo3bdk3mn77dF8iP8 0nhjPhM/53VnsLIWyd6iKpNeaj4+sdS3rYskLWJ57gR1ftGJuCHN70SWMdUu XkD3p6jl2cne7+T+joWDrtTZtSNuP/Av4Rw48+c4dfCaRKFoAV4UWyikbihi 4tHzxB0sWV6479UfNaEWCk86q6fLC1LEdnSijpm6sdDuwQvvELnDwdS3n6VM SpfwYkXNXP8L6ubgP7r+8OZD3tRQftA9et77C+6UVfPD5szwwOB9un+vF767 2MGPznOFx79T59gXKbqN8iPySpaCVDETmx7eS/4mJYA4D7V6c2qr2JIoHBLA 1+oIhyzqAPOHpxrlBCEUu8fSv4SJ2pw6y77fhND8WWlweSmtV796ODJGCNeW CwhpUNtWPYjUzxHClcRWCwvqttbcxshuIcx/3bT+LPX7hZid+luWQXun5s65 pbyDp0vkzDK8W9N4fKyM5qdnuZsOiyCmY3FxtJyJ9m3qERfYYrDMUTK3qqL9 02S0pwSK4cy2r/f8qMdlnZVKosXAclk0vEottHi+6k2lGG5uu5TcSm1Y3Tdr +NMKVPL2ZHpW0/vN7LrvWN0KiGqbtATVMOFrJ2i9W20lnAR8Ep/VMnG68Kjg tc/iULTLlH9Yz8Qc+dUq7NMqZJ7g9xlsY+KJtsGpc29X44KsjlxnPxOpb3Ot Z9NkkGH1pTjjIxMe0sbJCUnyMPhdqzGf1xMqa4Ki87wUcCdtPP7QT55IL68r 54gq4eef9lhnangiI8vtg0jEejwNyP+QDU+oz9WN+Q0q45rqtZJVbvT7yKqf WWxVZBnozzSf9MTjhnW+7j1qsJHljXCN98RBf5eAlYYaOC+Ohz+KPSG8S67C vWQDtvZ5xxV0eWJDYlknS0oLeoFKFR4LnjBReKP7PV4bJ++K7StT9sK3PQHn U1dsRPfzq/vVd3mB8bE7a/6gDhZudB9yDvJCJcdXYPaELiafiPDp5HtBRqB5 XZOaHiTa9BLK//GC/QLXQLlTD/wyx8MYEvsQGFPvZxm7CfK26Udu2+/DiKuk cIOZPl7ffDkSfXEfHGtV6iX5DJBmPe8c3LQP0Wf9Q/urDGDFyn5kLcLCymOp GaKHNsPV6ahHlDYLl/ivv+PV2gKzrFXPdjmzYNkQGOXcuwXh+s5eo6dZeDXw u4HW5a1QVbO+4prDAjslYKHJ0hCpfOs7Mv5mQe+S0YbGL4aQeNnqxyvqjXmp pv/GlW6DAMv5RbG5N/g6LXUq9hvhiR6REAn1xsDolwGf1dtRuu2ovmSVNx6m H/qzq307lqUeO/VIyAcdkS8/ZocbY8ufATv/cvLBB/F5W3NzEwQJK6sLF/mg /eOeDB8eBpRnk0ZaVrMhLe4Uo8jPwI3AzI5j0mwwdXeF9gkysPVLVI2ULBtD ATae7qIMyOsOJnqtZWN+2kTGaTUDEkcOhCwqs7Hus9rvptoMPBeODvhVjw2/ +fmziq50/PLjelkHNv7lSWX1FTFw2F8/2DiSjbrtMiHdxQyknPxwxeUyG5Gn EhJbShnQNuvxPXaVDcnxay/+qmRAsO/H5aIYNlS6zxhnNzBgHBe/SG6wYV3A kg54wwCrbMqoKpeNWDeVFt6VBBu67aJsX7LhEp8+tiBBoNPFX3qnlQ25l/LL Pv9EEOXDtRDsYOO2uRRjSJYgJ0ar6tkrNh7oCOU3qhFIv85+EvCWjV7BkYtx DALHufkoky9sqBbnEdXDBBM8w0/FFDkovtHUMhZA0G5xXyZxPQc4N+ZdfJzA 6n9bQlRUOWDu1gongQREePCWtSYHsZO5T93C6fyzeecfGHDwQzPXMSqZoOxQ btc2Ww5eZeaw5+oJ1EZSC/R+48D3ytPPVc8IXF7zbN0azMHMsZGI8GaCZ8q6 903PcPAfM40ciTYCAa5OOzecA4v+7HGtfwhW3Bg5MRFD65PPPsmaIWgrbtaQ zOOgvrHt8fuVwMcCGa+C9xwYqbuJm0oAodUH/w4Z4aAo7A0zVRLoyOYTchnn IMl4fNZZCqjs25ug9pmDXwsFNRpkARVLq19sBLjoW371hLI8ULcqb/kRYS5+ 2S9ZG7oWCLkQopIhxoWxooK7oSJgbvfKUFuKi//EGl69q0rXD1I8M6bGxfnJ 6leC6kBmim7NOS0u5mwtVH00AK2pzDQNXS76+Z2qZbWA/5ockMnbxkXxSf+p K7rAgbwwmWQHLtTapo1H9IAXBtnGzU5c3Nz42yVzfTpfVglntSsX4cPh679v BoSs9ZTnfbhwcbu15/B2oMbmHr9CEBeNpSoZjcZAfef9o4yzXBDJvAk1BnCk ok4h4gIXGk0PLrwB0G8Z1xIfx0WKOqN9uxnQVOlmwk3mQiK8TiFpB8A8ffyj fwYXC8at5bstgdQN/k6b7nFx+MZewQIr4KL/rYJ/H3DxbvYfRxEbQOfU1FaF R1x6b3BucW0B3vkPvTefctFcODpaawe8bebbFNbKBUSPbl3nABRF6QYO9XBR un8uLGgnoFq9Rq/3PReadadbunYBkt1yaec/cnFLUWCtviNwrPywnsA3LiRP X/aLdgL2pvcXpi3zxcW/xcvGfgFuzIiZPJTyxeLmRD5rZ/q/Dd6P+dV8ERC7 dtdtF+Af9+pGPiNfvJ/MvMnjCgR218w2uvhi7drh2kLqn2X+qFE+5Yv/A6LP l4o= "]]}, {RGBColor[0, 0, 1], LineBox[CompressedData[" 1:eJwV0nk01N8bB3DGEl98FRESaRUqe2Lmc58oWSuFimYMM6PFEpIiRcmShBY7 hbKkovK1/9CGshRNlpA1WbJkncmS3+2Pe+55nXuee557n7eS09nDLBIPD88r vP7utTtiB3I4WUS57Avr6ctUsChJdrJYDEYTW0uc/nr9uPE5+cUYtIvWGPHX DZUP7IQvp6KB8rnPfx0ulHyJfzEDyZqq6cxg77eKT172e4JaSZ75f02yrKUt +T9HnPF3lrPYdE+TCT9yPtKVUOafw66MqQ7gLhSgTV6JLX+tULp3pU9ZMTKQ kq/iYPt3vUmd9itDr0jPPnGx20l7NDz1K9BXQ8u539jJd7vpHP9XKKJrUW8B uy5DYVth52tU0F0Wt4g9X0Sd9Ca/RdJW4RJ/sI93dlydWqhCWfk2J3ivUCFs XM7shX0Niq84rMqHXcxjJ+FR9h6Bm4O0ALb05rb0Mb86lD6VBcLY+3atcXna Xo9II/0BItjeprZaLvofkUX1zg4x7GHdEccZ/0bke5EzIIHN80q4mq3WhL6c PxcjhS1juk0lv7MJKWYt02WwTexPTXmS2WhmycRCATsz4Me1iYVmVLAm/Zwq doWQ4ODHJy3omWk6307s5tubzXPtW9EHn/w8TWz+R0xJt7I2ZKUrfEIfm/G+ 7+GIXwfSzwwWMMf2syIJ1ap0ok+Cq3wOYt9pV3J93N6Jiui5JGvsN6N07VP6 XWg0StGDih0ddrZvyr8HzXK6kzyx5S0D7vV96UHau5dHfbCzV0UZf1brRdxs DTt/7MrE3JwXnb1oqaIxLBR7/OmYlwe5H7WmLeo+wLZscuEbXxhApQvarWzs lgHaPiP7H6ikRy2oA5s2bxUaX/oDZVzfeqgf22OjnoiR7yAa2bbFcvpvv+f5 JePnhtDlVVbs1QG4XjZlg+HkTyRwIX0lDZu2I5oZd2gUpcaM1J/EHjQMyhzN G0VTAfo5nthcl9Mqce5jKHHTdE0w9toKHc3Rn+NI62dTVB62usbgevPpXyh2 T5KVcCCeZ1nP69kNk+j40VsKUtjF+9qdUg9PolGfW6JK2HvsGh7NvphEgxEl 5ruxDwflK6d6TKHcUw22Ln/rmwPUZ8emUY2HWWozdomvDDwYmkMnyGL8b69S IbLRwkBJioMc1lt+b8J22npV9+EeDop+4trbgy3aPKyWlchBtCQLzWVsunqZ bJ4FF1keTNlBuUYFoUHqdMXz32hFwK/bFdjHrB9mdl1cRLe/865oDaLC9iet afSMRcTnrHl9GJtEEk3pa1pEdW5HFBaxn+Z53/2hsoR4zpx+rHQd51PEOHC8 cwkllc3KuWNnvR46/geWkSNJxlU0mAqX1sQJRvHzQkTBbC0thAqvGuKM6HK8 wGe69qI3tuD1+EANdV7wt928Jxw7+lfCPNueFwxChrcUYj/6kDwuk88LJZy9 e8RDqVB/Kb0l3ZEErhM9L6qxFXqeZRZW8AGDJy3U7AbOb2xuf+gXPvDty33J wM62yFt/fIQP/pEr/u2PrVnyPHFBih90HlR25GHvv5MfCW78UPqheEI6HOdn b4lP7VoB0K/gPzaC/Tr7nXHXRUG4I1m6MzMC96tVMRgWLQj3nFXNq7DNyovC tLIFoW0gO/Q79uemnNqwVkFQHK5nbLxFhe/z0Qe0dFdAc3Vbdjq2oCXVNmx6 BYisIJ3MisT1k7MsTXdhOOg5WvYhmgrs3VuDQxiisLGw8l+vGDw/FYKd7CsK Xg3xX2KxR+VslPKjRKE18WpRGbbg4rXy7jJRqCKFfhGIpYJeRdes3moxWD18 1SUJO9kw1vnnOzEIMTH8tz6OCs7mAiaHtoiDaY/nTqNEKlzO9RS4NbMSnGz7 t954QAUOOrM/aEoSUgutrt18SoUqNW2fq33ScEjmRIdnGRVS+nJMZh/IQkKD 0BbjOirYy5ATY+LloUVSzjymiwqb1vhFPaEpwsegjT8eT1EhtfhdMVNECQR8 N1laCdEgLeP4sHDwBvg5f1RGYgMNtnLe/Tw9sBHEq4smjck0SA0r305nbAaD RK+59fY0eFuj4GzXvgVuHfO6y71EA1cXWw9xPWWc+21r55JpIHRwbald/jYQ 8sh+e+YtDbbFFTbTpVQhiggWHB2mAUWxW33pnhq0X6sVuy/tAAvWHtdSxHaA AEveKXOfAxATrRlc152g4Wq10t/bAcqYzvyz3uogziunmpPjALL89Qp1WzQg y+9rXWqPA1jMs7Q3NmvADLdE0kuUDr7R1aeN72iCb8e3eo4aHYaOSQjVGGpB rNBEgrUlHaxeb6qWIGmDysUliyovOkQFugT0lGtDTM6MoWIiHcS9UtJE3HTA 5bxyoN1rOtzgi+3nVdWFczb9F8fG6WBc4xtp06kL5pF6lSnyjtDWe1tbNXwX yMYO976wcARGssd8nbEefOvqeykQ7AgaN/S31c7pwePY9yLCFY7Alaq7cLdg NwS/2zsexOsEpGbjnaUn9UE5Ql6109gJekfmep2kDWD2c8gvqTtOUJLq9l8L 2wCupB14NT7mBF/CPk1kXSeD0Sd3Pa9DDJhRITtmhpKhsSGWgg4zYHVD9udH 4WSQPjIVLmbNAJtV1wrSosnwtEltRe5RBrQkaF9KSiZDaO+Nqj80BnzNSRCM +o8MJY16Hr1uDOiqZ6zz+U4GC5/T7PIIBvxxb4z0HiRDlE5WZXEkAxRXUXjO jZBBY+qRWUE0A+g20v1nf5HhToGBWf49BvR+e599epEM3Xncg2+TGfB9fLs2 dTUFLp4dS9V9xoDhlVyzvXspoF/V+DLkEwOY163eae+nwIeezJSXTQzo4eRQ NptR4LrhjsQeNgPauqjqgoco8E29tdOwjQHVT99KvbengJjcunVKvQx4ZBLV Y36OAjpK//M3nWYALXCLz5F0CjgtsWQuyDKhfTpgwiiDAiZlJPnJtUywOfn1 lHY2Bf5o5Ii4KzDB4kCEvVQuBW4P7/7isZEJBvJT0FpCAYPmW7JJ25kgU1wu Yt9EgYe03Mc3DZnAnrBOc+IhwFQ5qV/7LD5feTh6PR8B4k1Nyb88mXBC/WBA lwAB5sEJx557M+GHhynVToSAjfVXZMl+TOBOUmQPSxOw71dfRkgwExRmttze o0bAP6vW1HGTmcBYvSlweQcBcpLvXZceMCFbW+lsuQYBujputiseMkHj/FrL 3bsImHJMOaP6mAl758SFNQwJ2Mwxzy4qYMJpLjdw/TECwtLrLy1+YkKuzNzZ LjsCxNi6wkw2E6b0pmnJVALymi6Q2C1M8PcdI8swCEjkN5ip+caEyPlerrg7 AQ0upWckR5nw32Ktx3IQAR8tC9NY/7Dgt/x7h/IQAg4sXuPS/mUBQak6cOkG AbETt/sYEiyouVypxokk4OWKDN47cixo/5M/NJ6Azy8fKr+pyoJlnhR6Vx4B wiUzWV6HWPDOQNa/9SUBxV0Oio02LAjziYlrLCDA0PtP8H57FkiM3vr4poyA IxeStMKdWbCp9Qo5qwb/n06Cyt0rLBhatXA0tZaA1YNH7k9fZ8EziwvnEhoI KFLsjgi4yQLdN2dzbrIJGBYfq1BOYIHJM7qMRzd+r97dLMsiFogOdWqd7sP3 nUm2ja1kQdOG4wedBgiw8nd+jd6zIJbafMZ+iABe5jfzY00ssIu3CrH+ScD+ QiPrqXYWKLIb0izHCRCasJnYPMCC72Jm5caTeN4SL15L/mJBtkl1G5ohgHLw /uLwAgvcggxn9DgEeN3v6BgTcgbNigpxzXncr1X1UNQaZ5jj6quqLhFgljNu 9EDZGcq0iow38SBAJ6m1cxRnCHTXclrHh8A7KkB5pYMzFMuiIl5+BCbse1Ij 553h/+77jQI= "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesLabel->{None, None}, AxesOrigin->{0, 110.}, Method->{}, PlotRange->{All, {110, 180}}, PlotRangeClipping->True, PlotRangePadding->{Automatic, Automatic}]], "Output", CellChangeTimes->{{3.6260647358350782`*^9, 3.626064762848359*^9}, { 3.626064880976041*^9, 3.6260649494869223`*^9}, 3.626064987341391*^9, 3.626067181310852*^9}] }, Open ]] }, Open ]] }, WindowSize->{1649, 1005}, WindowMargins->{{0, Automatic}, {Automatic, -1}}, FrontEndVersion->"9.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (November 20, \ 2012)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[579, 22, 197, 4, 51, "Subtitle"], Cell[CellGroupData[{ Cell[801, 30, 1143, 36, 46, "Input"], Cell[1947, 68, 543, 16, 47, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[2527, 89, 13030, 345, 508, "Input"], Cell[15560, 436, 2482, 71, 57, "Output"], Cell[18045, 509, 1968, 55, 48, "Output"], Cell[20016, 566, 1066, 29, 47, "Output"], Cell[21085, 597, 1319, 38, 61, "Output"], Cell[22407, 637, 2222, 63, 59, "Output"], Cell[24632, 702, 2487, 71, 75, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[27156, 778, 205, 3, 29, "Subsubtitle"], Cell[CellGroupData[{ Cell[27386, 785, 9103, 235, 216, "Input"], Cell[36492, 1022, 5995, 124, 494, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[42536, 1152, 213, 3, 29, "Subsubtitle"], Cell[CellGroupData[{ Cell[42774, 1159, 2993, 79, 114, "Input"], Cell[45770, 1240, 90, 2, 47, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[45909, 1248, 188, 4, 29, "Subsubtitle"], Cell[CellGroupData[{ Cell[46122, 1256, 416, 9, 63, "Input"], Cell[46541, 1267, 164, 2, 28, "Output"], Cell[46708, 1271, 164, 2, 28, "Output"], Cell[46875, 1275, 165, 2, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[47077, 1282, 1627, 45, 79, "Input"], Cell[48707, 1329, 2043, 40, 340, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[50799, 1375, 127, 1, 29, "Subsubtitle"], Cell[CellGroupData[{ Cell[50951, 1380, 1796, 47, 131, "Input"], Cell[52750, 1429, 503, 10, 28, "Output"], Cell[53256, 1441, 500, 10, 28, "Output"], Cell[53759, 1453, 501, 10, 28, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[54309, 1469, 113, 1, 29, "Subsubtitle"], Cell[CellGroupData[{ Cell[54447, 1474, 1088, 27, 128, "Input"], Cell[55538, 1503, 534, 12, 47, "Output"], Cell[56075, 1517, 466, 9, 28, "Output"], Cell[56544, 1528, 458, 9, 28, "Output"] }, Open ]] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[57063, 1544, 394, 8, 69, "Subtitle"], Cell[CellGroupData[{ Cell[57482, 1556, 1563, 43, 84, "Input"], Cell[59048, 1601, 8471, 147, 249, "Output"] }, Open ]] }, Open ]] } ] *) (* End of internal cache information *)