(* Content-type: application/vnd.wolfram.mathematica *)
(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)
(* CreatedBy='Mathematica 8.0' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 157, 7]
NotebookDataLength[ 20343, 585]
NotebookOptionsPosition[ 18813, 527]
NotebookOutlinePosition[ 19225, 545]
CellTagsIndexPosition[ 19182, 542]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"SetDirectory", "[", "\"\<~/KappaLib/\>\"", "]"}],
";"}], "\[IndentingNewLine]",
RowBox[{"<<", "kappaLib-1.2.m"}]}], "Input",
CellChangeTimes->{{3.540187180485169*^9, 3.540187180769973*^9}, {
3.54073721266546*^9, 3.5407372157409163`*^9}}],
Cell[BoxData["\<\"Loading KappaLib v1.2\"\>"], "Print",
CellChangeTimes->{
3.532166858781642*^9, 3.532185546518886*^9, {3.5321916427282057`*^9,
3.532191659074267*^9}, 3.532258220199937*^9, {3.5327883038360443`*^9,
3.532788318439725*^9}, 3.533101929005726*^9, 3.534234782932765*^9,
3.535031367840541*^9, 3.535032110065617*^9, 3.535184820127049*^9,
3.540185710489208*^9, {3.540187181218997*^9, 3.540187190411971*^9},
3.540187586799612*^9, 3.5407372163453836`*^9, 3.540737608593027*^9,
3.5407379323734818`*^9, 3.5407394916248837`*^9}]
}, Open ]],
Cell["\<\
Verify that the Fresnel polynomial for the PDCM medium in
Lindell, Bergamin, Favaro: Decomposable medium conditions in \
four-dimensional representation
factorises into the product of two quadratic forms.\
\>", "Subsection",
CellChangeTimes->{{3.532183469269575*^9, 3.532183474884473*^9}, {
3.5321911799271517`*^9, 3.5321911807574377`*^9}, {3.540187060250759*^9,
3.540187087593748*^9}, {3.540735826543047*^9, 3.5407358354447393`*^9}, {
3.540737750095323*^9, 3.540737801709317*^9}}],
Cell[CellGroupData[{
Cell["Define kappa", "Subsection",
CellChangeTimes->{{3.532183480804378*^9, 3.532183483412216*^9}}],
Cell[BoxData[
RowBox[{
RowBox[{"(*", " ",
RowBox[{
RowBox[{"general", " ",
RowBox[{"(",
RowBox[{"1", ",", "1"}], ")"}]}], "-", "tensor"}], " ", "*)"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"Pmat", "=",
RowBox[{"emMatrix", "[",
RowBox[{"\"\
\"", ",", "4", ",",
RowBox[{"Structure", "\[Rule]", "\"\\""}]}], "]"}]}], ";"}],
"\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{"two", " ", "arbitrary", " ", "bivectors"}], " ", "*)"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"CC", "=",
RowBox[{"emMatrix", "[",
RowBox[{"\"\\"", ",", "4", ",",
RowBox[{"Structure", "\[Rule]", "\"\\""}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"DD", "=",
RowBox[{"emMatrix", "[",
RowBox[{"\"\\"", ",", "4", ",",
RowBox[{"Structure", "\[Rule]", "\"\\""}]}], "]"}]}],
";"}], "\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{"C1", " ", "=", " ", "constant"}], " ", "*)"}],
"\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{"C2", " ", "=", " ", "constant"}], " ", "*)"}],
"\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{"rho", " ", "=", " ",
RowBox[{"scalar", " ", "density", " ", "of", " ", "weight", " ", "1"}]}],
" ", "*)"}]}]}]], "Input",
CellChangeTimes->{{3.540735898943542*^9, 3.5407359333080072`*^9}, {
3.540736914515017*^9, 3.540736993815291*^9}}],
Cell[BoxData[
RowBox[{
RowBox[{"kappa", "=",
RowBox[{
RowBox[{"C1", " ",
RowBox[{"emIdentityKappa", "[", "]"}]}], "+", " ",
RowBox[{"C2", " ",
RowBox[{"emPQToKappa", "[",
RowBox[{"Pmat", ",", "Pmat"}], "]"}]}], "+",
RowBox[{"emBiProduct", "[",
RowBox[{"rho", ",", "CC", ",", "DD"}], "]"}]}]}], ";"}]], "Input",
CellChangeTimes->{{3.540735949960291*^9, 3.540736033505999*^9}, {
3.5407361040329313`*^9, 3.540736109255342*^9}, {3.540736276821391*^9,
3.5407363013873367`*^9}, {3.540736944329253*^9, 3.540736947321188*^9}}]
}, Open ]],
Cell[CellGroupData[{
Cell["Compute Fresnel polynomial", "Subsection",
CellChangeTimes->{{3.532185731773658*^9, 3.5321857350935783`*^9}}],
Cell[BoxData[{
RowBox[{
RowBox[{"vars", "=",
RowBox[{"{",
RowBox[{"xi0", ",", "xi1", ",", "xi2", ",", "xi3"}], "}"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"fresnel", "=",
RowBox[{"emKappaToFresnel", "[",
RowBox[{"kappa", ",", "vars"}], "]"}]}], ";"}]}], "Input",
CellChangeTimes->{{3.532165876043063*^9, 3.532165895943246*^9}, {
3.5321857188312683`*^9, 3.5321857396228228`*^9}, 3.532185775526833*^9, {
3.5407372290173473`*^9, 3.540737232128345*^9}}]
}, Open ]],
Cell[CellGroupData[{
Cell["\<\
Check that the Fresnel polynomial factorises into two quadratic forms\
\>", "Subsection",
CellChangeTimes->{{3.5407376798418617`*^9, 3.540737689897423*^9}}],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"(*", " ",
RowBox[{"define", " ", "adjugate", " ", "of", " ", "P"}], " ", "*)"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"adjP", " ", "=", " ",
RowBox[{"Table", "[", "\n", "\t\t\t",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"(",
RowBox[{"-", "1"}], ")"}], "^",
RowBox[{"(",
RowBox[{"i", "+", "j"}], ")"}]}],
RowBox[{"Det", "[",
RowBox[{"Drop", "[",
RowBox[{"Pmat", ",",
RowBox[{"{",
RowBox[{"i", ",", "i"}], "}"}], ",",
RowBox[{"{",
RowBox[{"j", ",", "j"}], "}"}]}], "]"}], "]"}]}], ",", "\n",
"\t\t\t",
RowBox[{"{",
RowBox[{"j", ",", "1", ",", "4"}], "}"}], ",",
RowBox[{"{",
RowBox[{"i", ",", "1", ",", "4"}], "}"}]}], "\n", "\t\t\t", "]"}]}],
";"}], "\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"(*", " ", "check", " ", "*)"}], "\[IndentingNewLine]",
RowBox[{"Simplify", "[",
RowBox[{
RowBox[{"adjP", "/",
RowBox[{"Det", "[", "Pmat", "]"}]}], "-",
RowBox[{"Inverse", "[", "Pmat", "]"}]}], "]"}]}]}]], "Input",
CellChangeTimes->{{3.5401870171758223`*^9, 3.540187034112459*^9}, {
3.5401873692631397`*^9, 3.540187382253614*^9}, {3.540737021135269*^9,
3.540737021390382*^9}}],
Cell[BoxData[
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"0", ",", "0", ",", "0", ",", "0"}], "}"}], ",",
RowBox[{"{",
RowBox[{"0", ",", "0", ",", "0", ",", "0"}], "}"}], ",",
RowBox[{"{",
RowBox[{"0", ",", "0", ",", "0", ",", "0"}], "}"}], ",",
RowBox[{"{",
RowBox[{"0", ",", "0", ",", "0", ",", "0"}], "}"}]}], "}"}]], "Output",
CellChangeTimes->{
3.5401870346443043`*^9, 3.54018719625863*^9, {3.54018737271117*^9,
3.5401873829715977`*^9}, 3.540187476048547*^9, 3.5401875934857817`*^9,
3.5407360543354397`*^9, 3.5407361339838047`*^9, 3.540736286649625*^9,
3.540736317798861*^9, 3.540737025846758*^9, 3.540737615197083*^9,
3.540737938960113*^9, 3.5407394979462433`*^9}]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"fresnelExp", "=",
RowBox[{
RowBox[{"-", "2"}], " ",
RowBox[{"C2", "^", "2"}], " ", "rho", " ",
RowBox[{"(",
RowBox[{"vars", ".", "Pmat", ".", "CC", ".", "vars"}], ")"}],
RowBox[{"(", " ",
RowBox[{"vars", ".", "adjP", ".", "DD", ".", "vars"}], ")"}]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{"Simplify", "[",
RowBox[{"fresnel", "-", "fresnelExp"}], "]"}]}], "Input",
CellChangeTimes->{{3.540187573104498*^9, 3.540187573258079*^9}, {
3.540736066737051*^9, 3.540736070857925*^9}, {3.540736186133387*^9,
3.540736207724115*^9}, {3.5407370372074633`*^9, 3.540737112050027*^9}, {
3.540737664062924*^9, 3.540737669556992*^9}}],
Cell[BoxData["0"], "Output",
CellChangeTimes->{{3.540187578990944*^9, 3.540187595062821*^9},
3.540736062090716*^9, 3.540736100637714*^9, 3.5407361651524963`*^9, {
3.540736209297638*^9, 3.540736213201111*^9}, 3.540736289493989*^9,
3.540736322609099*^9, {3.5407370551382732`*^9, 3.5407371149994907`*^9}, {
3.540737659559926*^9, 3.540737673459835*^9}, 3.540737940841277*^9,
3.540739499854574*^9}]
}, Open ]]
}, Open ]],
Cell["\<\
Note: the above formula is valid also when P is not invertible\
\>", "Subsection",
CellChangeTimes->{{3.54073788967749*^9, 3.5407379000105886`*^9}}],
Cell[CellGroupData[{
Cell["\<\
Extra: Verify some basic identities for P/\\ Q\
\>", "Section",
CellChangeTimes->{{3.5407371687138567`*^9, 3.540737177056037*^9}}],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"(*", " ",
RowBox[{"Write", " ", "identity", " ", "kappa", " ", "as", " ",
RowBox[{
RowBox[{"id", "/"}], "\\", "id"}]}], " ", "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"id", "=",
RowBox[{"IdentityMatrix", "[", "4", "]"}]}], ";"}],
"\[IndentingNewLine]",
RowBox[{"Union", "[",
RowBox[{"Flatten", "[",
RowBox[{
RowBox[{"emIdentityKappa", "[", "]"}], "-",
RowBox[{"emPQToKappa", "[",
RowBox[{"id", ",", "id"}], "]"}]}], "]"}], "]"}]}]}]], "Input",
CellChangeTimes->{{3.5407337530056057`*^9, 3.540733766097044*^9}, {
3.540735060877563*^9, 3.540735094788509*^9}, {3.540737700422338*^9,
3.54073770849887*^9}}],
Cell[BoxData[
RowBox[{"{", "0", "}"}]], "Output",
CellChangeTimes->{{3.54073507711222*^9, 3.540735096674111*^9},
3.540735298374405*^9, 3.5407379409396057`*^9, 3.540739499954134*^9}]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"(*", " ",
RowBox[{
RowBox[{
RowBox[{"If", " ", "S", " ", "is", " ", "a", " ", "trace"}], "-",
RowBox[{"free", " ",
RowBox[{"(",
RowBox[{"1", ",", "1"}], ")"}]}], "-", "tensor"}], ",", " ",
RowBox[{
RowBox[{
"then", " ", "\[IndentingNewLine]", "\[IndentingNewLine]", " ",
RowBox[{
RowBox[{"S", " ", "/"}], "\\", " ", "Id"}]}], " ", "+", " ",
RowBox[{
RowBox[{
RowBox[{"Id", " ", "/"}], "\\", " ", "S"}], "\[IndentingNewLine]",
"\[IndentingNewLine]", "has", " ", "only", " ", "a", " ", "skewon", " ",
"part"}]}]}], " ", "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"S", "=",
RowBox[{"emMatrix", "[",
RowBox[{"\"\\"", ",", "4", ",",
RowBox[{"Structure", "\[Rule]", "\"\\""}]}], "]"}]}], ";"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"S", "=",
RowBox[{"S", "-",
RowBox[{
RowBox[{"1", "/", "4"}],
RowBox[{"Tr", "[", "S", "]"}],
RowBox[{"IdentityMatrix", "[", "4", "]"}]}]}]}], ";"}],
"\[IndentingNewLine]",
RowBox[{"Tr", "[", "S", "]"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"kappa", "=",
RowBox[{"Simplify", "[",
RowBox[{
RowBox[{"emPQToKappa", "[",
RowBox[{"id", ",", "S"}], "]"}], "+",
RowBox[{"emPQToKappa", "[",
RowBox[{"S", ",", "id"}], "]"}]}], "]"}]}], ";"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"emKappaToMatrix", "[", "kappa", "]"}], "//", "MatrixForm"}],
"\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
"medium", " ", "has", " ", "only", " ", "a", " ", "skewon", " ", "part"}],
" ", "*)"}], "\[IndentingNewLine]",
RowBox[{"Union", "[",
RowBox[{"Flatten", "[",
RowBox[{"Simplify", "[",
RowBox[{"kappa", "+",
RowBox[{"emPoincare", "[", "kappa", "]"}]}], "]"}], "]"}],
"]"}]}]}]], "Input",
CellChangeTimes->{{3.540735100012807*^9, 3.5407352837054358`*^9}, {
3.5407353808891068`*^9, 3.5407354028348227`*^9}, {3.540735518936508*^9,
3.540735525206901*^9}, {3.540735575809033*^9, 3.540735671055735*^9}, {
3.540736334730239*^9, 3.5407363926249533`*^9}, {3.540737915072391*^9,
3.54073791551315*^9}}],
Cell[BoxData["0"], "Output",
CellChangeTimes->{
3.540735187691566*^9, {3.540735236274173*^9, 3.540735284223765*^9},
3.540735347226202*^9, {3.5407353815268927`*^9, 3.5407354270606318`*^9},
3.540735535026779*^9, 3.54073559757342*^9, 3.540735633804201*^9,
3.540735671526009*^9, 3.540737940985022*^9, 3.540739499997499*^9}],
Cell[BoxData[
TagBox[
RowBox[{"(", "\[NoBreak]", GridBox[{
{
RowBox[{
FractionBox["1", "2"], " ",
RowBox[{"(",
RowBox[{"sh11", "+", "sh22", "-", "sh33", "-", "sh44"}], ")"}]}],
"sh32", "sh42", "0", "sh41",
RowBox[{"-", "sh31"}]},
{"sh23",
RowBox[{
FractionBox["1", "2"], " ",
RowBox[{"(",
RowBox[{"sh11", "-", "sh22", "+", "sh33", "-", "sh44"}], ")"}]}],
"sh43",
RowBox[{"-", "sh41"}], "0", "sh21"},
{"sh24", "sh34",
RowBox[{
FractionBox["1", "2"], " ",
RowBox[{"(",
RowBox[{"sh11", "-", "sh22", "-", "sh33", "+", "sh44"}], ")"}]}],
"sh31",
RowBox[{"-", "sh21"}], "0"},
{"0",
RowBox[{"-", "sh14"}], "sh13",
RowBox[{
FractionBox["1", "2"], " ",
RowBox[{"(",
RowBox[{
RowBox[{"-", "sh11"}], "-", "sh22", "+", "sh33", "+", "sh44"}],
")"}]}],
RowBox[{"-", "sh23"}],
RowBox[{"-", "sh24"}]},
{"sh14", "0",
RowBox[{"-", "sh12"}],
RowBox[{"-", "sh32"}],
RowBox[{
FractionBox["1", "2"], " ",
RowBox[{"(",
RowBox[{
RowBox[{"-", "sh11"}], "+", "sh22", "-", "sh33", "+", "sh44"}],
")"}]}],
RowBox[{"-", "sh34"}]},
{
RowBox[{"-", "sh13"}], "sh12", "0",
RowBox[{"-", "sh42"}],
RowBox[{"-", "sh43"}],
RowBox[{
FractionBox["1", "2"], " ",
RowBox[{"(",
RowBox[{
RowBox[{"-", "sh11"}], "+", "sh22", "+", "sh33", "-", "sh44"}],
")"}]}]}
},
GridBoxAlignment->{
"Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}},
"RowsIndexed" -> {}},
GridBoxSpacings->{"Columns" -> {
Offset[0.27999999999999997`], {
Offset[0.7]},
Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
Offset[0.2], {
Offset[0.4]},
Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
Function[BoxForm`e$,
MatrixForm[BoxForm`e$]]]], "Output",
CellChangeTimes->{
3.540735187691566*^9, {3.540735236274173*^9, 3.540735284223765*^9},
3.540735347226202*^9, {3.5407353815268927`*^9, 3.5407354270606318`*^9},
3.540735535026779*^9, 3.54073559757342*^9, 3.540735633804201*^9,
3.540735671526009*^9, 3.540737940985022*^9, 3.5407395000013227`*^9}],
Cell[BoxData[
RowBox[{"{", "0", "}"}]], "Output",
CellChangeTimes->{
3.540735187691566*^9, {3.540735236274173*^9, 3.540735284223765*^9},
3.540735347226202*^9, {3.5407353815268927`*^9, 3.5407354270606318`*^9},
3.540735535026779*^9, 3.54073559757342*^9, 3.540735633804201*^9,
3.540735671526009*^9, 3.540737940985022*^9, 3.540739500268915*^9}]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"(*", " ", "\[IndentingNewLine]", " ",
RowBox[{"Conversely", ",", " ",
RowBox[{
RowBox[{
"the", " ", "skewon", " ", "part", " ", "of", " ", "any", " ", "kappa",
" ", "can", " ", "be", " ", "written", " ", "as", " ",
"\[IndentingNewLine]", "\[IndentingNewLine]", " ",
RowBox[{
RowBox[{"S", "/"}], "\\", "Id"}]}], " ", "+", " ",
RowBox[{
RowBox[{
RowBox[{"Id", " ", "/"}], "\\", " ", "S"}], "\[IndentingNewLine]",
"\[IndentingNewLine]", "where", " ", "S", " ", "is", " ", "the", " ",
"trace"}], "-",
RowBox[{
"free", " ", "part", " ", "of", " ", "the", " ", "first", " ", "trace",
" ", "of", " ", "kappa"}]}]}], " ", "\[IndentingNewLine]", "*)"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"kappa", "=",
RowBox[{"emGeneralKappa", "[", "\"\\"", "]"}]}], ";"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"kappaII", " ", "=", " ",
RowBox[{
RowBox[{"1", "/", "2"}],
RowBox[{"(",
RowBox[{"kappa", " ", "-", " ",
RowBox[{"emPoincare", "[", "kappa", "]"}]}], ")"}]}]}], ";"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"S", "=",
RowBox[{"Table", "[", "\[IndentingNewLine]",
RowBox[{
RowBox[{"Sum", "[", "\[IndentingNewLine]",
RowBox[{
RowBox[{"emReadNormal", "[",
RowBox[{"kappa", ",", "i", ",", "j", ",", "i", ",", "r"}], "]"}],
",",
RowBox[{"{",
RowBox[{"i", ",", "1", ",", "4"}], "}"}]}], "]"}],
"\[IndentingNewLine]", ",", "\[IndentingNewLine]",
RowBox[{"{",
RowBox[{"j", ",", "1", ",", "4"}], "}"}], ",",
RowBox[{"{",
RowBox[{"r", ",", "1", ",", "4"}], "}"}]}], "\[IndentingNewLine]",
"]"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"S", "=",
RowBox[{"(",
RowBox[{"S", "-",
RowBox[{
RowBox[{"1", "/", "2"}], " ",
RowBox[{"emTrace", "[", "kappa", "]"}],
RowBox[{"IdentityMatrix", "[", "4", "]"}]}]}], ")"}]}], ";"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"kappaIIalt", "=",
RowBox[{
RowBox[{"1", "/", "2"}],
RowBox[{"Simplify", "[", "\[IndentingNewLine]",
RowBox[{
RowBox[{"emPQToKappa", "[",
RowBox[{"id", ",", "S"}], "]"}], "+",
RowBox[{"emPQToKappa", "[",
RowBox[{"S", ",", "id"}], "]"}]}], "\[IndentingNewLine]", "]"}]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{"Union", "[",
RowBox[{"Flatten", "[",
RowBox[{"Simplify", "[",
RowBox[{"kappaII", "-", " ", "kappaIIalt"}], "]"}], "]"}],
"]"}]}]}]], "Input",
CellChangeTimes->{{3.540736395879915*^9, 3.540736503422307*^9}, {
3.540736556731851*^9, 3.5407366202859364`*^9}, {3.5407366584551077`*^9,
3.540736668685248*^9}, {3.540736701551873*^9, 3.540736775969623*^9}, {
3.540736879071015*^9, 3.5407368826872187`*^9}, 3.5407371921505003`*^9, {
3.5407377209754343`*^9, 3.54073773431851*^9}}],
Cell[BoxData[
RowBox[{"{", "0", "}"}]], "Output",
CellChangeTimes->{
3.540735187691566*^9, {3.540735236274173*^9, 3.540735284223765*^9},
3.540735347226202*^9, {3.5407353815268927`*^9, 3.5407354270606318`*^9},
3.540735535026779*^9, 3.54073559757342*^9, 3.540735633804201*^9,
3.54073567152789*^9, {3.540736453057632*^9, 3.540736467023405*^9},
3.540736504561767*^9, {3.540736595991777*^9, 3.540736622782304*^9}, {
3.540736659927359*^9, 3.540736669280085*^9}, 3.540736708389797*^9, {
3.5407367587880993`*^9, 3.540736767576867*^9}, 3.5407368839088383`*^9,
3.540737941514192*^9, 3.540739500557535*^9}]
}, Open ]]
}, Open ]]
},
WindowSize->{1211, 940},
WindowMargins->{{15, Automatic}, {Automatic, 17}},
PrintingCopies->1,
PrintingPageRange->{1, Automatic},
FrontEndVersion->"8.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (February 23, \
2011)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)
(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 281, 6, 43, "Input"],
Cell[863, 30, 558, 8, 20, "Print"]
}, Open ]],
Cell[1436, 41, 508, 11, 95, "Subsection"],
Cell[CellGroupData[{
Cell[1969, 56, 100, 1, 26, "Subsection"],
Cell[2072, 59, 1552, 41, 163, "Input"],
Cell[3627, 102, 570, 13, 27, "Input"]
}, Open ]],
Cell[CellGroupData[{
Cell[4234, 120, 116, 1, 34, "Subsection"],
Cell[4353, 123, 496, 12, 43, "Input"]
}, Open ]],
Cell[CellGroupData[{
Cell[4886, 140, 167, 3, 35, "Subsection"],
Cell[CellGroupData[{
Cell[5078, 147, 1370, 37, 133, "Input"],
Cell[6451, 186, 728, 16, 27, "Output"]
}, Open ]],
Cell[CellGroupData[{
Cell[7216, 207, 706, 16, 43, "Input"],
Cell[7925, 225, 412, 6, 27, "Output"]
}, Open ]]
}, Open ]],
Cell[8364, 235, 159, 3, 34, "Subsection"],
Cell[CellGroupData[{
Cell[8548, 242, 141, 3, 37, "Section"],
Cell[CellGroupData[{
Cell[8714, 249, 729, 19, 58, "Input"],
Cell[9446, 270, 187, 3, 27, "Output"]
}, Open ]],
Cell[CellGroupData[{
Cell[9670, 278, 2289, 61, 193, "Input"],
Cell[11962, 341, 334, 5, 27, "Output"],
Cell[12299, 348, 2392, 71, 175, "Output"],
Cell[14694, 421, 356, 6, 27, "Output"]
}, Open ]],
Cell[CellGroupData[{
Cell[15087, 432, 3070, 79, 313, "Input"],
Cell[18160, 513, 625, 10, 27, "Output"]
}, Open ]]
}, Open ]]
}
]
*)
(* End of internal cache information *)