(*********************************************************************** Mathematica-Compatible Notebook This notebook can be used on any computer system with Mathematica 3.0, MathReader 3.0, or any compatible application. The data for the notebook starts with the line of stars above. To get the notebook into a Mathematica-compatible application, do one of the following: * Save the data starting with the line of stars above into a file with a name ending in .nb, then open the file inside the application; * Copy the data starting with the line of stars above to the clipboard, then use the Paste menu command inside the application. Data for notebooks contains only printable 7-bit ASCII and can be sent directly in email or through ftp in text mode. Newlines can be CR, LF or CRLF (Unix, Macintosh or MS-DOS style). NOTE: If you modify the data for this notebook not in a Mathematica- compatible application, you must delete the line below containing the word CacheID, otherwise Mathematica-compatible applications may try to use invalid cache data. For more information on notebooks and Mathematica-compatible applications, contact Wolfram Research: web: http://www.wolfram.com email: info@wolfram.com phone: +1-217-398-0700 (U.S.) Notebook reader applications are available free of charge from Wolfram Research. ***********************************************************************) (*CacheID: 232*) (*NotebookFileLineBreakTest NotebookFileLineBreakTest*) (*NotebookOptionsPosition[ 29274, 1070]*) (*NotebookOutlinePosition[ 55309, 2018]*) (* CellTagsIndexPosition[ 55265, 2014]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell["Multilocus genotype frequencies", "Title", Evaluatable->False, AspectRatioFixed->True], Cell["\<\ Published results using this notebook should cite: Hospital, F., Dillmann, C. and Melchinger, A.E. (1996) A general algorithm to \ compute multilocus genotype frequencies under various mating systems. Comput. \ Appl. Biosci, 12(6):455-462. For more information please contact fred@moulon.inra.fr\ \>", "Text"], Cell[CellGroupData[{ Cell[TextData["Self-Fertilization"], "Subtitle", Evaluatable->False, AspectRatioFixed->True], Cell[CellGroupData[{ Cell[TextData["Initialization"], "Section", Evaluatable->False, AspectRatioFixed->True], Cell["ClearAll[\"Global`*\"];", "Input", InitializationCell->True, AspectRatioFixed->True], Cell[CellGroupData[{ Cell[TextData["Number of loci"], "Subsection", Evaluatable->False, AspectRatioFixed->True], Cell[CellGroupData[{ Cell["n=3", "Input", InitializationCell->True, AspectRatioFixed->True], Cell[BoxData[ \(3\)], "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[TextData["General definitions"], "Subsection", Evaluatable->False, InitializationCell->True, AspectRatioFixed->True], Cell["<<\"definitions.m\"", "Input", InitializationCell->True, AspectRatioFixed->True], Cell["\<\ S1= Table[{i,i},{i,1,ngam/2}]; S2=Flatten[ Table[{i,j}, {i,1,ngam/2-1},{j,i+1,ngam-i}],1]; S3=Table[{ngam/2+1-i,ngam/2+i},{i,1,ngam/2}]; S4=Flatten[ Table[{i,j}, {j,ngam/2+2,ngam},{i,ngam+2-j,j-1}],1]; S5=Table[{i,i},{i,ngam/2+1,ngam}]; L1=Length[S1];L2=Length[S2];L3=Length[S3];\ \>", "Input", InitializationCell->True, AspectRatioFixed->True], Cell[CellGroupData[{ Cell["genlist=Join[S1,S2,S3,S4,S5]", "Input", InitializationCell->True, AspectRatioFixed->True], Cell[BoxData[ \({{1, 1}, {2, 2}, {3, 3}, {4, 4}, {1, 2}, {1, 3}, {1, 4}, {1, 5}, {1, 6}, {1, 7}, {2, 3}, {2, 4}, {2, 5}, {2, 6}, {3, 4}, {3, 5}, {4, 5}, { 3, 6}, {2, 7}, {1, 8}, {4, 6}, {5, 6}, {3, 7}, {4, 7}, {5, 7}, {6, 7}, {2, 8}, {3, 8}, {4, 8}, {5, 8}, {6, 8}, {7, 8}, {5, 5}, {6, 6}, { 7, 7}, {8, 8}}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell["ngen=Length[genlist]", "Input", InitializationCell->True, AspectRatioFixed->True], Cell[BoxData[ \(36\)], "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[TextData["Initial genotype frequencies"], "Subsection", Evaluatable->False, AspectRatioFixed->True], Cell[CellGroupData[{ Cell[TextData["No symmetry"], "Subsubsection", Evaluatable->False, AspectRatioFixed->True], Cell[CellGroupData[{ Cell["ulist=Table[{genlist[[i]]},{i,ngen}]", "Input", AspectRatioFixed->True], Cell[BoxData[ \({{{1, 1}}, {{2, 2}}, {{3, 3}}, {{4, 4}}, {{1, 2}}, {{1, 3}}, {{1, 4}}, {{1, 5}}, {{1, 6}}, {{1, 7}}, {{2, 3}}, {{2, 4}}, {{2, 5}}, {{ 2, 6}}, {{3, 4}}, {{3, 5}}, {{4, 5}}, {{3, 6}}, {{2, 7}}, {{1, 8}}, {{4, 6}}, {{5, 6}}, {{3, 7}}, {{4, 7}}, {{5, 7}}, {{6, 7}}, {{ 2, 8}}, {{3, 8}}, {{4, 8}}, {{5, 8}}, {{6, 8}}, {{7, 8}}, {{5, 5}}, {{6, 6}}, {{7, 7}}, {{8, 8}}}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell["nsym=Length[ulist]", "Input", AspectRatioFixed->True], Cell[BoxData[ \(36\)], "Output"] }, Open ]], Cell["\<\ Clear[index]; index=Table[\"+\",{ngam},{ngam}]; Do[index[[ulist[[i,1,1]],ulist[[i,1,2]]]]=i,{i,1,nsym}];\ \>", "Input", AspectRatioFixed->True], Cell[TextData[ "Index for elements of freq[t], vector of all genotype frequencies at t"], "Text", Evaluatable->False, AspectRatioFixed->True], Cell[CellGroupData[{ Cell["\<\ TableForm[index, TableHeadings->{Table[gamname[i],{i,ngam}], Table[gamname[i],{i,ngam}]}, TableAlignments->Center]\ \>", "Input", AspectRatioFixed->True], Cell[BoxData[ TagBox[GridBox[{ {"\<\"\"\>", "\<\"000\"\>", "\<\"001\"\>", "\<\"010\"\>", "\<\"011\"\>", "\<\"100\"\>", "\<\"101\"\>", "\<\"110\"\>", "\<\"111\"\>"}, {"\<\"000\"\>", "1", "5", "6", "7", "8", "9", "10", "20"}, {"\<\"001\"\>", "\<\"+\"\>", "2", "11", "12", "13", "14", "19", "27"}, {"\<\"010\"\>", "\<\"+\"\>", "\<\"+\"\>", "3", "15", "16", "18", "23", "28"}, {"\<\"011\"\>", "\<\"+\"\>", "\<\"+\"\>", "\<\"+\"\>", "4", "17", "21", "24", "29"}, {"\<\"100\"\>", "\<\"+\"\>", "\<\"+\"\>", "\<\"+\"\>", "\<\"+\"\>", "33", "22", "25", "30"}, {"\<\"101\"\>", "\<\"+\"\>", "\<\"+\"\>", "\<\"+\"\>", "\<\"+\"\>", "\<\"+\"\>", "34", "26", "31"}, {"\<\"110\"\>", "\<\"+\"\>", "\<\"+\"\>", "\<\"+\"\>", "\<\"+\"\>", "\<\"+\"\>", "\<\"+\"\>", "35", "32"}, {"\<\"111\"\>", "\<\"+\"\>", "\<\"+\"\>", "\<\"+\"\>", "\<\"+\"\>", "\<\"+\"\>", "\<\"+\"\>", "\<\"+\"\>", "36"} }, RowSpacings->1, ColumnSpacings->3, RowAlignments->Baseline, ColumnAlignments->{Center}], (TableForm[ #, TableHeadings -> {{"000", "001", "010", "011", "100", "101", "110", "111"}, {"000", "001", "010", "011", "100", "101", "110", "111"}}, TableAlignments -> Center]&)]], "Output"] }, Open ]], Cell["Define all initial frequencies below:", "Text", Evaluatable->False, AspectRatioFixed->True], Cell["\<\ Do[f[1,x]=0,{x,nsym}]; f[1,index[[1,ngam]]]=1;\ \>", "Input", AspectRatioFixed->True] }, Open ]], Cell[CellGroupData[{ Cell[TextData["Symmetry"], "Subsubsection", Evaluatable->False, InitializationCell->True, AspectRatioFixed->True], Cell["rs[k_]:=Sum[10^(n-l-1)*(1+rho[k,l]),{l,n-1}]", "Input", InitializationCell->True, AspectRatioFixed->True], Cell["genrs[i_,j_]:=(-1)^(gam[i][[1]]+gam[j][[1]])*(rs[i]+rs[j])", "Input", InitializationCell->True, AspectRatioFixed->True], Cell["rslist[ens_]:=Apply[genrs,ens,2]", "Input", InitializationCell->True, AspectRatioFixed->True], Cell["\<\ usUnion[ens_]:=Map[ens[[##]]&, Sort[Map[Position[ens,#][[1,1]]&,Union[ens]]]]\ \>", "Input", InitializationCell->True, AspectRatioFixed->True], Cell[CellGroupData[{ Cell["\<\ SS2=Map[Apply[S2[[##]]&,Position[rslist[S2],##],2]&, usUnion[rslist[S2]]]\ \>", "Input", InitializationCell->True, AspectRatioFixed->True], Cell[BoxData[ \({{{1, 2}}, {{1, 3}}}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell["LL2=Length[SS2]", "Input", InitializationCell->True, AspectRatioFixed->True], Cell[BoxData[ \(2\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell["S4r=Reverse[S4]", "Input", InitializationCell->True, AspectRatioFixed->True], Cell[BoxData[ \({{3, 4}, {2, 4}}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell["\<\ SS4=Map[Apply[S4r[[##]]&,Position[rslist[S4r],##],2]&, usUnion[rslist[S4r]]]\ \>", "Input", InitializationCell->True, AspectRatioFixed->True], Cell[BoxData[ \({{{3, 4}}, {{2, 4}}}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell["\<\ ulist=Join[ Table[{S1[[i]],S5[[L1-i+1]]},{i,L1}], Table[Join[SS2[[i]],SS4[[i]]],{i,LL2}], Table[{S3[[i]]},{i,L3}]]; ulist//ColumnForm\ \>", "Input", InitializationCell->True, AspectRatioFixed->True], Cell[BoxData[ InterpretationBox[GridBox[{ {\({{1, 1}, {4, 4}}\)}, {\({{2, 2}, {3, 3}}\)}, {\({{1, 2}, {3, 4}}\)}, {\({{1, 3}, {2, 4}}\)}, {\({{2, 3}}\)}, {\({{1, 4}}\)} }, GridBaseline->{Baseline, {1, 1}}, ColumnAlignments->{Left}], ColumnForm[ {{{1, 1}, {4, 4}}, {{2, 2}, {3, 3}}, {{1, 2}, {3, 4}}, {{1, 3}, {2, 4}}, {{2, 3}}, {{1, 4}}}], Editable->False]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell["nsym=Length[ulist]", "Input", InitializationCell->True, AspectRatioFixed->True], Cell[BoxData[ \(6\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Clear[index]; index=Table[\"+\",{ngam},{ngam}]; Do[Apply[(index[[##]]=i)&,ulist[[i]],2],{i,nsym}]; TableForm[index, TableHeadings->{Table[gamname[i],{i,ngam}], Table[gamname[i],{i,ngam}]}, TableAlignments->Center]\ \>", "Input", InitializationCell->True, AspectRatioFixed->True], Cell[BoxData[ TagBox[GridBox[{ {"\<\"\"\>", "\<\"00\"\>", "\<\"01\"\>", "\<\"10\"\>", "\<\"11\"\>"}, {"\<\"00\"\>", "1", "3", "4", "6"}, {"\<\"01\"\>", "\<\"+\"\>", "2", "5", "4"}, {"\<\"10\"\>", "\<\"+\"\>", "\<\"+\"\>", "2", "3"}, {"\<\"11\"\>", "\<\"+\"\>", "\<\"+\"\>", "\<\"+\"\>", "1"} }, RowSpacings->1, ColumnSpacings->3, RowAlignments->Baseline, ColumnAlignments->{Center}], (TableForm[ #, TableHeadings -> {{"00", "01", "10", "11"}, {"00", "01", "10", "11"}}, TableAlignments -> Center]&)]], "Output"] }, Open ]], Cell[TextData["Define all initial frequencies below:"], "Text", Evaluatable->False, InitializationCell->True, AspectRatioFixed->True], Cell[CellGroupData[{ Cell["\<\ Do[f[1,x]=0,{x,1,nsym-1}]; f[1,nsym]=1\ \>", "Input", InitializationCell->True, AspectRatioFixed->True], Cell[BoxData[ \(1\)], "Output"] }, Open ]] }, Closed]] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[TextData["Tables"], "Section", Evaluatable->False, AspectRatioFixed->True], Cell["\<\ Format[zf[t_,x_]]:=Subscripted[\"f*\"[x]][t]; Format[zf[t_,x_],TeXForm]:=Subscripted[\"f^{\\ast}\"[x]][t]; Format[zf[t,x_],TeXForm]:=Subscripted[\"f^{\\ast}\"[x]]\ \>", "Input", InitializationCell->True, AspectRatioFixed->True], Cell[CellGroupData[{ Cell["tabf=Table[zf[t,x],{x,nsym}]", "Input", InitializationCell->True, AspectRatioFixed->True], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ InterpretationBox[\("f*"\_1\), Subscripted[ "f*"[ 1]], Editable->False], "[", "t", "]"}], ",", RowBox[{ InterpretationBox[\("f*"\_2\), Subscripted[ "f*"[ 2]], Editable->False], "[", "t", "]"}], ",", RowBox[{ InterpretationBox[\("f*"\_3\), Subscripted[ "f*"[ 3]], Editable->False], "[", "t", "]"}], ",", RowBox[{ InterpretationBox[\("f*"\_4\), Subscripted[ "f*"[ 4]], Editable->False], "[", "t", "]"}], ",", RowBox[{ InterpretationBox[\("f*"\_5\), Subscripted[ "f*"[ 5]], Editable->False], "[", "t", "]"}], ",", RowBox[{ InterpretationBox[\("f*"\_6\), Subscripted[ "f*"[ 6]], Editable->False], "[", "t", "]"}], ",", RowBox[{ InterpretationBox[\("f*"\_7\), Subscripted[ "f*"[ 7]], Editable->False], "[", "t", "]"}], ",", RowBox[{ InterpretationBox[\("f*"\_8\), Subscripted[ "f*"[ 8]], Editable->False], "[", "t", "]"}], ",", RowBox[{ InterpretationBox[\("f*"\_9\), Subscripted[ "f*"[ 9]], Editable->False], "[", "t", "]"}], ",", RowBox[{ InterpretationBox[\("f*"\_10\), Subscripted[ "f*"[ 10]], Editable->False], "[", "t", "]"}]}], "}"}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[TextData["Index for genotypes"], "Subsection", Evaluatable->False, AspectRatioFixed->True], Cell[CellGroupData[{ Cell[TextData["Visual"], "Subsubsection", Evaluatable->False, AspectRatioFixed->True], Cell[CellGroupData[{ Cell["\<\ Table[Join[ {zf[t,i]}, Apply[StringJoin[gamname[#1],\"/\",gamname[#2]]&,ulist[[i]],2] ],{i,nsym}]//TableForm\ \>", "Input", AspectRatioFixed->True], Cell[BoxData[ TagBox[GridBox[{ {\(zf[t, 1]\), "\<\"000/000\"\>"}, {\(zf[t, 2]\), "\<\"001/001\"\>"}, {\(zf[t, 3]\), "\<\"010/010\"\>"}, {\(zf[t, 4]\), "\<\"011/011\"\>"}, {\(zf[t, 5]\), "\<\"000/001\"\>"}, {\(zf[t, 6]\), "\<\"000/010\"\>"}, {\(zf[t, 7]\), "\<\"000/011\"\>"}, {\(zf[t, 8]\), "\<\"000/100\"\>"}, {\(zf[t, 9]\), "\<\"000/101\"\>"}, {\(zf[t, 10]\), "\<\"000/110\"\>"}, {\(zf[t, 11]\), "\<\"001/010\"\>"}, {\(zf[t, 12]\), "\<\"001/011\"\>"}, {\(zf[t, 13]\), "\<\"001/100\"\>"}, {\(zf[t, 14]\), "\<\"001/101\"\>"}, {\(zf[t, 15]\), "\<\"010/011\"\>"}, {\(zf[t, 16]\), "\<\"010/100\"\>"}, {\(zf[t, 17]\), "\<\"011/100\"\>"}, {\(zf[t, 18]\), "\<\"010/101\"\>"}, {\(zf[t, 19]\), "\<\"001/110\"\>"}, {\(zf[t, 20]\), "\<\"000/111\"\>"}, {\(zf[t, 21]\), "\<\"011/101\"\>"}, {\(zf[t, 22]\), "\<\"100/101\"\>"}, {\(zf[t, 23]\), "\<\"010/110\"\>"}, {\(zf[t, 24]\), "\<\"011/110\"\>"}, {\(zf[t, 25]\), "\<\"100/110\"\>"}, {\(zf[t, 26]\), "\<\"101/110\"\>"}, {\(zf[t, 27]\), "\<\"001/111\"\>"}, {\(zf[t, 28]\), "\<\"010/111\"\>"}, {\(zf[t, 29]\), "\<\"011/111\"\>"}, {\(zf[t, 30]\), "\<\"100/111\"\>"}, {\(zf[t, 31]\), "\<\"101/111\"\>"}, {\(zf[t, 32]\), "\<\"110/111\"\>"}, {\(zf[t, 33]\), "\<\"100/100\"\>"}, {\(zf[t, 34]\), "\<\"101/101\"\>"}, {\(zf[t, 35]\), "\<\"110/110\"\>"}, {\(zf[t, 36]\), "\<\"111/111\"\>"} }, RowSpacings->1, ColumnSpacings->3, RowAlignments->Baseline, ColumnAlignments->{Left}], (TableForm[ #]&)]], "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["TeX", "Subsubsection", Evaluatable->False, AspectRatioFixed->True], Cell["\<\ Table[Join[ {zf[t,i]}, Apply[StringJoin[\"\\\\frac{\",gamname[#1],\"}{\",gamname[#2],\"}\"]&,ulist[[\ i]],2] ],{i,nsym}]//ColumnForm//TeXForm\ \>", "Input", AspectRatioFixed->True] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[TextData["Recursion"], "Subsection", Evaluatable->False, AspectRatioFixed->True], Cell[CellGroupData[{ Cell[TextData["Definition"], "Subsubsection", Evaluatable->False, InitializationCell->True, AspectRatioFixed->True], Cell["\<\ genrec[x_]:=Sum[ 2^delta[ ulist[[x,1,1]] , ulist[[x,1,2]] ] * Simplify[(Apply[Plus, Apply[P[##][[ ulist[[x,1,1]] ]] * P[##][[ ulist[[x,1,2]] ]]&,ulist[[k]],2]] ) ]* zf[t,k] ,{k,1,nsym}]\ \>", "Input", InitializationCell->True, AspectRatioFixed->True] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[TextData["Frequencies"], "Subsection", Evaluatable->False, AspectRatioFixed->True], Cell[CellGroupData[{ Cell[TextData["Definition"], "Subsubsection", Evaluatable->False, InitializationCell->True, AspectRatioFixed->True], Cell["\<\ f[t_,x_]:=f[t,x]=Sum[ 2^delta[ ulist[[x,1,1]] , ulist[[x,1,2]] ] * Simplify[(Apply[Plus, Apply[P[##][[ ulist[[x,1,1]] ]] * P[##][[ ulist[[x,1,2]] ]]&,ulist[[k]],2]] ) ]* f[t-1,k] ,{k,1,nsym}]\ \>", "Input", InitializationCell->True, AspectRatioFixed->True], Cell["freq[t_]:=Table[f[t,x],{x,nsym}]", "Input", InitializationCell->True, AspectRatioFixed->True], Cell["f[t_,i_,j_]:=freq[t][[index[[i,j]]]]", "Input", InitializationCell->True, AspectRatioFixed->True], Cell[CellGroupData[{ Cell[BoxData[ \(freq[1]\)], "Input"], Cell[BoxData[ \({0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(freq[2]\)], "Input"], Cell[BoxData[ \({1\/4\ \((\(-1\) + r)\)\^2, r\^2\/4, \(-\(1\/2\)\)\ \((\(-1\) + r)\)\ r, \(-\(1\/2\)\)\ \((\(-1\) + r)\)\ r, r\^2\/2, 1\/2\ \((\(-1\) + r)\)\^2, \(-\(1\/2\)\)\ \((\(-1\) + r)\)\ r, \(-\(1\/2\)\)\ \((\(-1\) + r)\)\ r, r\^2\/4, 1\/4\ \((\(-1\) + r)\)\^2} \)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(freq[2] // ColumnForm\)], "Input"], Cell[BoxData[ InterpretationBox[GridBox[{ {\(1\/4\ \((\(-1\) + r\_1)\)\^2\ \((\(-1\) + r\_2)\)\^2\)}, {\(1\/4\ \((\(-1\) + r\_1)\)\^2\ r\_2\%2\)}, {\(1\/4\ r\_1\%2\ r\_2\%2\)}, {\(1\/4\ r\_1\%2\ \((\(-1\) + r\_2)\)\^2\)}, { \(\(-\(1\/2\)\)\ \((\(-1\) + r\_1)\)\^2\ \((\(-1\) + r\_2)\)\ r\_2\)}, {\(1\/2\ \((\(-1\) + r\_1)\)\ r\_1\ \((\(-1\) + r\_2)\)\ r\_2\)}, { \(\(-\(1\/2\)\)\ \((\(-1\) + r\_1)\)\ r\_1\ \((\(-1\) + r\_2)\)\^2\)}, { \(\(-\(1\/2\)\)\ \((\(-1\) + r\_1)\)\ r\_1\ \((\(-1\) + r\_2)\)\^2\)}, {\(1\/2\ \((\(-1\) + r\_1)\)\ r\_1\ \((\(-1\) + r\_2)\)\ r\_2\)}, { \(\(-\(1\/2\)\)\ \((\(-1\) + r\_1)\)\^2\ \((\(-1\) + r\_2)\)\ r\_2\)}, {\(\(-\(1\/2\)\)\ \((\(-1\) + r\_1)\)\ r\_1\ r\_2\%2\)}, {\(1\/2\ \((\(-1\) + r\_1)\)\ r\_1\ \((\(-1\) + r\_2)\)\ r\_2\)}, {\(1\/2\ \((\(-1\) + r\_1)\)\ r\_1\ \((\(-1\) + r\_2)\)\ r\_2\)}, {\(\(-\(1\/2\)\)\ \((\(-1\) + r\_1)\)\ r\_1\ r\_2\%2\)}, {\(\(-\(1\/2\)\)\ r\_1\%2\ \((\(-1\) + r\_2)\)\ r\_2\)}, {\(\(-\(1\/2\)\)\ r\_1\%2\ \((\(-1\) + r\_2)\)\ r\_2\)}, {\(1\/2\ r\_1\%2\ \((\(-1\) + r\_2)\)\^2\)}, {\(1\/2\ r\_1\%2\ r\_2\%2\)}, {\(1\/2\ \((\(-1\) + r\_1)\)\^2\ r\_2\%2\)}, {\(1\/2\ \((\(-1\) + r\_1)\)\^2\ \((\(-1\) + r\_2)\)\^2\)}, {\(\(-\(1\/2\)\)\ r\_1\%2\ \((\(-1\) + r\_2)\)\ r\_2\)}, {\(\(-\(1\/2\)\)\ r\_1\%2\ \((\(-1\) + r\_2)\)\ r\_2\)}, {\(\(-\(1\/2\)\)\ \((\(-1\) + r\_1)\)\ r\_1\ r\_2\%2\)}, {\(1\/2\ \((\(-1\) + r\_1)\)\ r\_1\ \((\(-1\) + r\_2)\)\ r\_2\)}, {\(1\/2\ \((\(-1\) + r\_1)\)\ r\_1\ \((\(-1\) + r\_2)\)\ r\_2\)}, {\(\(-\(1\/2\)\)\ \((\(-1\) + r\_1)\)\ r\_1\ r\_2\%2\)}, { \(\(-\(1\/2\)\)\ \((\(-1\) + r\_1)\)\^2\ \((\(-1\) + r\_2)\)\ r\_2\)}, {\(1\/2\ \((\(-1\) + r\_1)\)\ r\_1\ \((\(-1\) + r\_2)\)\ r\_2\)}, { \(\(-\(1\/2\)\)\ \((\(-1\) + r\_1)\)\ r\_1\ \((\(-1\) + r\_2)\)\^2\)}, { \(\(-\(1\/2\)\)\ \((\(-1\) + r\_1)\)\ r\_1\ \((\(-1\) + r\_2)\)\^2\)}, {\(1\/2\ \((\(-1\) + r\_1)\)\ r\_1\ \((\(-1\) + r\_2)\)\ r\_2\)}, { \(\(-\(1\/2\)\)\ \((\(-1\) + r\_1)\)\^2\ \((\(-1\) + r\_2)\)\ r\_2\)}, {\(1\/4\ r\_1\%2\ \((\(-1\) + r\_2)\)\^2\)}, {\(1\/4\ r\_1\%2\ r\_2\%2\)}, {\(1\/4\ \((\(-1\) + r\_1)\)\^2\ r\_2\%2\)}, {\(1\/4\ \((\(-1\) + r\_1)\)\^2\ \((\(-1\) + r\_2)\)\^2\)} }, GridBaseline->{Baseline, {1, 1}}, ColumnAlignments->{Left}], ColumnForm[ { Times[ Rational[ 1, 4], Power[ Plus[ -1, r[ 1]], 2], Power[ Plus[ -1, r[ 2]], 2]], Times[ Rational[ 1, 4], Power[ Plus[ -1, r[ 1]], 2], Power[ r[ 2], 2]], Times[ Rational[ 1, 4], Power[ r[ 1], 2], Power[ r[ 2], 2]], Times[ Rational[ 1, 4], Power[ r[ 1], 2], Power[ Plus[ -1, r[ 2]], 2]], Times[ Rational[ -1, 2], Power[ Plus[ -1, r[ 1]], 2], Plus[ -1, r[ 2]], r[ 2]], Times[ Rational[ 1, 2], Plus[ -1, r[ 1]], r[ 1], Plus[ -1, r[ 2]], r[ 2]], Times[ Rational[ -1, 2], Plus[ -1, r[ 1]], r[ 1], Power[ Plus[ -1, r[ 2]], 2]], Times[ Rational[ -1, 2], Plus[ -1, r[ 1]], r[ 1], Power[ Plus[ -1, r[ 2]], 2]], Times[ Rational[ 1, 2], Plus[ -1, r[ 1]], r[ 1], Plus[ -1, r[ 2]], r[ 2]], Times[ Rational[ -1, 2], Power[ Plus[ -1, r[ 1]], 2], Plus[ -1, r[ 2]], r[ 2]], Times[ Rational[ -1, 2], Plus[ -1, r[ 1]], r[ 1], Power[ r[ 2], 2]], Times[ Rational[ 1, 2], Plus[ -1, r[ 1]], r[ 1], Plus[ -1, r[ 2]], r[ 2]], Times[ Rational[ 1, 2], Plus[ -1, r[ 1]], r[ 1], Plus[ -1, r[ 2]], r[ 2]], Times[ Rational[ -1, 2], Plus[ -1, r[ 1]], r[ 1], Power[ r[ 2], 2]], Times[ Rational[ -1, 2], Power[ r[ 1], 2], Plus[ -1, r[ 2]], r[ 2]], Times[ Rational[ -1, 2], Power[ r[ 1], 2], Plus[ -1, r[ 2]], r[ 2]], Times[ Rational[ 1, 2], Power[ r[ 1], 2], Power[ Plus[ -1, r[ 2]], 2]], Times[ Rational[ 1, 2], Power[ r[ 1], 2], Power[ r[ 2], 2]], Times[ Rational[ 1, 2], Power[ Plus[ -1, r[ 1]], 2], Power[ r[ 2], 2]], Times[ Rational[ 1, 2], Power[ Plus[ -1, r[ 1]], 2], Power[ Plus[ -1, r[ 2]], 2]], Times[ Rational[ -1, 2], Power[ r[ 1], 2], Plus[ -1, r[ 2]], r[ 2]], Times[ Rational[ -1, 2], Power[ r[ 1], 2], Plus[ -1, r[ 2]], r[ 2]], Times[ Rational[ -1, 2], Plus[ -1, r[ 1]], r[ 1], Power[ r[ 2], 2]], Times[ Rational[ 1, 2], Plus[ -1, r[ 1]], r[ 1], Plus[ -1, r[ 2]], r[ 2]], Times[ Rational[ 1, 2], Plus[ -1, r[ 1]], r[ 1], Plus[ -1, r[ 2]], r[ 2]], Times[ Rational[ -1, 2], Plus[ -1, r[ 1]], r[ 1], Power[ r[ 2], 2]], Times[ Rational[ -1, 2], Power[ Plus[ -1, r[ 1]], 2], Plus[ -1, r[ 2]], r[ 2]], Times[ Rational[ 1, 2], Plus[ -1, r[ 1]], r[ 1], Plus[ -1, r[ 2]], r[ 2]], Times[ Rational[ -1, 2], Plus[ -1, r[ 1]], r[ 1], Power[ Plus[ -1, r[ 2]], 2]], Times[ Rational[ -1, 2], Plus[ -1, r[ 1]], r[ 1], Power[ Plus[ -1, r[ 2]], 2]], Times[ Rational[ 1, 2], Plus[ -1, r[ 1]], r[ 1], Plus[ -1, r[ 2]], r[ 2]], Times[ Rational[ -1, 2], Power[ Plus[ -1, r[ 1]], 2], Plus[ -1, r[ 2]], r[ 2]], Times[ Rational[ 1, 4], Power[ r[ 1], 2], Power[ Plus[ -1, r[ 2]], 2]], Times[ Rational[ 1, 4], Power[ r[ 1], 2], Power[ r[ 2], 2]], Times[ Rational[ 1, 4], Power[ Plus[ -1, r[ 1]], 2], Power[ r[ 2], 2]], Times[ Rational[ 1, 4], Power[ Plus[ -1, r[ 1]], 2], Power[ Plus[ -1, r[ 2]], 2]]}], Editable->False]], "Output"] }, Closed]], Cell[BoxData[ \(\(freq[3]; \)\)], "Input"], Cell[BoxData[ \(\(freq[4]; \)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(freq[4] /. {r[1] -> .1, r[2] -> .2} // ColumnForm\)], "Input"], Cell[BoxData[ InterpretationBox[GridBox[{ {"0.251349329888000028`"}, {"0.0750957958880000298`"}, {"0.0116712838880000036`"}, {"0.0322128178880000026`"}, {"0.034285124224000012`"}, {"0.00961149196800000282`"}, {"0.00743652403200000655`"}, {"0.0180138522240000078`"}, {"0.00185836377600000091`"}, {"0.0227815960320000066`"}, {"0.0013799920320000012`"}, {"0.00961149196800000282`"}, {"0.00185836377600000091`"}, {"0.00630892022400000129`"}, {"0.004846148224000002`"}, {"0.000843420032000000396`"}, {"0.00065256396800000056`"}, {"0.000121095968000000131`"}, {"0.00327091996800000295`"}, {"0.0176263879680000146`"}, {"0.000843420032000000396`"}, {"0.004846148224000002`"}, {"0.00630892022400000129`"}, {"0.00185836377600000091`"}, {"0.00961149196800000282`"}, {"0.0013799920320000012`"}, {"0.0227815960320000066`"}, {"0.00185836377600000091`"}, {"0.0180138522240000078`"}, {"0.00743652403200000655`"}, {"0.00961149196800000282`"}, {"0.034285124224000012`"}, {"0.0322128178880000026`"}, {"0.0116712838880000036`"}, {"0.0750957958880000298`"}, {"0.251349329888000028`"} }, GridBaseline->{Baseline, {1, 1}}, ColumnAlignments->{Left}], ColumnForm[ {0.25134932988800002, 0.075095795888000028, 0.011671283888000004, 0.032212817888000005, 0.034285124224000013, 0.0096114919680000029, 0.0074365240320000067, 0.01801385222400001, 0.0018583637760000008, 0.022781596032000008, 0.0013799920320000012, 0.0096114919680000029, 0.0018583637760000008, 0.0063089202240000017, 0.0048461482240000023, 0.00084342003200000041, 0.00065256396800000057, 0.00012109596800000014, 0.0032709199680000028, 0.017626387968000014, 0.00084342003200000041, 0.0048461482240000023, 0.0063089202240000017, 0.0018583637760000008, 0.0096114919680000029, 0.0013799920320000012, 0.022781596032000008, 0.0018583637760000008, 0.01801385222400001, 0.0074365240320000067, 0.0096114919680000029, 0.034285124224000013, 0.032212817888000005, 0.011671283888000004, 0.075095795888000028, 0.25134932988800002}], Editable->False]], "Output"] }, Closed]], Cell[BoxData[ \(\(freq[5]; \)\)], "Input"], Cell[BoxData[ \(\(freq[6]; \)\)], "Input"], Cell[BoxData[ \(Save["\", freq]\)], "Input"] }, Open ]] }, Open ]] }, Open ]] }, Open ]] }, Open ]] }, FrontEndVersion->"X 3.0", ScreenRectangle->{{0, 1024}, {0, 768}}, AutoGeneratedPackage->Automatic, WindowToolbars->{}, CellGrouping->Automatic, WindowSize->{520, 600}, WindowMargins->{{Automatic, 201}, {Automatic, 39}}, PrivateNotebookOptions->{"ColorPalette"->{RGBColor, 128}}, ShowCellLabel->True, ShowCellTags->False, RenderingOptions->{"ObjectDithering"->True, "RasterDithering"->False}, CharacterEncoding->"XAutomaticEncoding", StyleDefinitions -> Notebook[{ Cell[CellGroupData[{ Cell["Style Definitions", "Subtitle"], Cell["\<\ Modify the definitions below to change the default appearance of \ all cells in a given style. Make modifications to any definition using \ commands in the Format menu.\ \>", "Text"], Cell[CellGroupData[{ Cell["Style Environment Names", "Section"], Cell[StyleData[All, "Working"], PageWidth->WindowWidth, ScriptMinSize->9], Cell[StyleData[All, "Presentation"], PageWidth->WindowWidth, ScriptMinSize->12, FontSize->16], Cell[StyleData[All, "Condensed"], PageWidth->WindowWidth, CellBracketOptions->{"Margins"->{1, 1}, "Widths"->{0, 5}}, ScriptMinSize->8, FontSize->11], Cell[StyleData[All, "Printout"], PageWidth->PaperWidth, ScriptMinSize->5, FontSize->10, PrivateFontOptions->{"FontType"->"Outline"}] }, Open ]], Cell[CellGroupData[{ Cell["Notebook Options", "Section"], Cell["\<\ The options defined for the style below will be used at the \ Notebook level.\ \>", "Text"], Cell[StyleData["Notebook"], PageHeaders->{{Cell[ TextData[ { CounterBox[ "Page"]}], "PageNumber"], None, Cell[ TextData[ { ValueBox[ "FileName"]}], "Header"]}, {Cell[ TextData[ { ValueBox[ "FileName"]}], "Header"], None, Cell[ TextData[ { CounterBox[ "Page"]}], "PageNumber"]}}, CellFrameLabelMargins->6, StyleMenuListing->None] }, Open ]], Cell[CellGroupData[{ Cell["Styles for Headings", "Section"], Cell[CellGroupData[{ Cell[StyleData["Title"], CellMargins->{{7, Inherited}, {8, 40}}, Evaluatable->False, CellGroupingRules->{"TitleGrouping", 0}, CellHorizontalScrolling->False, PageBreakBelow->False, GroupPageBreakWithin->False, TextAlignment->Center, CounterIncrements->"Title", CounterAssignments->{{"Section", 0}, {"Equation", 0}, {"Figure", 0}, { "Subtitle", 0}, {"Subsubtitle", 0}}, AspectRatioFixed->True, FontFamily->"times", FontSize->24, FontWeight->"Bold"], Cell[StyleData["Title", "Presentation"], CellMargins->{{24, 10}, {20, 40}}, LineSpacing->{1, 0}, FontSize->44], Cell[StyleData["Title", "Condensed"], CellMargins->{{8, 10}, {4, 8}}, FontSize->20], Cell[StyleData["Title", "Printout"], CellMargins->{{2, 10}, {15, 30}}, FontSize->24] }, Open ]], Cell[CellGroupData[{ Cell[StyleData["Subtitle"], CellMargins->{{7, Inherited}, {6, 15}}, Evaluatable->False, CellGroupingRules->{"TitleGrouping", 10}, CellHorizontalScrolling->False, PageBreakBelow->False, GroupPageBreakWithin->False, TextAlignment->Center, CounterIncrements->"Subtitle", CounterAssignments->{{"Section", 0}, {"Equation", 0}, {"Figure", 0}, { "Subsubtitle", 0}}, AspectRatioFixed->True, FontFamily->"times", FontSize->18, FontWeight->"Bold"], Cell[StyleData["Subtitle", "Presentation"], CellMargins->{{24, 10}, {15, 20}}, LineSpacing->{1, 0}, FontSize->36], Cell[StyleData["Subtitle", "Condensed"], CellMargins->{{8, 10}, {4, 4}}, FontSize->14], Cell[StyleData["Subtitle", "Printout"], CellMargins->{{2, 10}, {10, 15}}, FontSize->18] }, Open ]], Cell[CellGroupData[{ Cell[StyleData["Subsubtitle"], CellMargins->{{7, Inherited}, {6, 20}}, Evaluatable->False, CellGroupingRules->{"TitleGrouping", 20}, CellHorizontalScrolling->False, PageBreakBelow->False, GroupPageBreakWithin->False, TextAlignment->Center, CounterIncrements->"Subsubtitle", CounterAssignments->{{"Section", 0}, {"Equation", 0}, {"Figure", 0}}, AspectRatioFixed->True, FontFamily->"times", FontSize->14, FontSlant->"Italic"], Cell[StyleData["Subsubtitle", "Presentation"], CellMargins->{{24, 10}, {10, 20}}, LineSpacing->{1, 0}, FontSize->24], Cell[StyleData["Subsubtitle", "Condensed"], CellMargins->{{8, 10}, {8, 12}}, FontSize->12], Cell[StyleData["Subsubtitle", "Printout"], CellMargins->{{2, 10}, {8, 10}}, FontSize->14] }, Open ]], Cell[CellGroupData[{ Cell[StyleData["Section"], CellDingbat->"\[GraySquare]", CellMargins->{{22, Inherited}, {8, 20}}, Evaluatable->False, CellGroupingRules->{"SectionGrouping", 30}, CellHorizontalScrolling->False, PageBreakBelow->False, CounterIncrements->"Section", CounterAssignments->{{"Subsection", 0}, {"Subsubsection", 0}}, AspectRatioFixed->True, FontFamily->"times", FontSize->18, FontWeight->"Bold"], Cell[StyleData["Section", "Presentation"], CellMargins->{{40, 10}, {11, 32}}, LineSpacing->{1, 0}, FontSize->24], Cell[StyleData["Section", "Condensed"], CellMargins->{{18, Inherited}, {6, 12}}, FontSize->12], Cell[StyleData["Section", "Printout"], CellMargins->{{13, 0}, {7, 22}}, FontSize->14] }, Open ]], Cell[CellGroupData[{ Cell[StyleData["Subsection"], CellDingbat->"\[FilledSquare]", CellMargins->{{19, Inherited}, {8, 15}}, Evaluatable->False, CellGroupingRules->{"SectionGrouping", 40}, CellHorizontalScrolling->False, PageBreakBelow->False, CounterIncrements->"Subsection", CounterAssignments->{{"Subsubsection", 0}}, AspectRatioFixed->True, FontFamily->"times", FontSize->14, FontWeight->"Bold"], Cell[StyleData["Subsection", "Presentation"], CellMargins->{{36, 10}, {11, 32}}, LineSpacing->{1, 0}, FontSize->22], Cell[StyleData["Subsection", "Condensed"], CellMargins->{{16, Inherited}, {6, 12}}, FontSize->12], Cell[StyleData["Subsection", "Printout"], CellMargins->{{9, 0}, {7, 22}}, FontSize->12] }, Open ]], Cell[CellGroupData[{ Cell[StyleData["Subsubsection"], CellDingbat->"\[EmptySquare]", CellMargins->{{18, Inherited}, {8, 12}}, Evaluatable->False, CellGroupingRules->{"SectionGrouping", 50}, CellHorizontalScrolling->False, PageBreakBelow->False, CounterIncrements->"Subsubsection", AspectRatioFixed->True, FontFamily->"times", FontSize->12, FontWeight->"Bold"], Cell[StyleData["Subsubsection", "Presentation"], CellMargins->{{34, 10}, {11, 26}}, LineSpacing->{1, 0}, FontSize->18], Cell[StyleData["Subsubsection", "Condensed"], CellMargins->{{17, Inherited}, {6, 12}}, FontSize->10], Cell[StyleData["Subsubsection", "Printout"], CellMargins->{{9, 0}, {7, 14}}, FontSize->11] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Styles for Body Text", "Section"], Cell[CellGroupData[{ Cell[StyleData["Text"], CellMargins->{{7, 10}, {7, 7}}, Evaluatable->False, CellHorizontalScrolling->False, PageBreakWithin->Automatic, LineSpacing->{1, 3}, CounterIncrements->"Text", AspectRatioFixed->True, FontFamily->"times", FontSize->12], Cell[StyleData["Text", "Presentation"], CellMargins->{{24, 10}, {10, 10}}, LineSpacing->{1, 5}], Cell[StyleData["Text", "Condensed"], CellMargins->{{8, 10}, {6, 6}}, LineSpacing->{1, 1}], Cell[StyleData["Text", "Printout"], CellMargins->{{2, 2}, {6, 6}}] }, Open ]], Cell[CellGroupData[{ Cell[StyleData["SmallText"], CellMargins->{{7, 10}, {6, 6}}, Evaluatable->False, CellHorizontalScrolling->False, PageBreakWithin->Automatic, LineSpacing->{1, 3}, CounterIncrements->"SmallText", AspectRatioFixed->True, FontFamily->"times", FontSize->10], Cell[StyleData["SmallText", "Presentation"], CellMargins->{{24, 10}, {8, 8}}, LineSpacing->{1, 5}, FontSize->12], Cell[StyleData["SmallText", "Condensed"], CellMargins->{{8, 10}, {5, 5}}, LineSpacing->{1, 2}, FontSize->9], Cell[StyleData["SmallText", "Printout"], CellMargins->{{2, 2}, {5, 5}}, FontSize->7] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Styles for Input/Output", "Section"], Cell["\<\ The cells in this section define styles used for input and output \ to the kernel. Be careful when modifying, renaming, or removing these \ styles, because the front end associates special meanings with these style \ names.\ \>", "Text"], Cell[CellGroupData[{ Cell[StyleData["Input"], PageWidth->Infinity, CellMargins->{{42, 10}, {5, 7}}, Evaluatable->True, CellGroupingRules->"InputGrouping", CellHorizontalScrolling->True, PageBreakWithin->False, GroupPageBreakWithin->False, CellLabelMargins->{{23, Inherited}, {Inherited, Inherited}}, DefaultFormatType->DefaultInputFormatType, AutoItalicWords->{}, FormatType->InputForm, ShowStringCharacters->True, NumberMarks->True, CounterIncrements->"Input", AspectRatioFixed->True, FontFamily->"courier", FontSize->12, FontWeight->"Bold"], Cell[StyleData["Input", "Presentation"], CellMargins->{{72, Inherited}, {8, 10}}, LineSpacing->{1, 0}], Cell[StyleData["Input", "Condensed"], CellMargins->{{40, 10}, {2, 3}}], Cell[StyleData["Input", "Printout"], CellMargins->{{39, 0}, {4, 6}}, FontSize->9] }, Open ]], Cell[StyleData["InputOnly"], Evaluatable->True, CellGroupingRules->"InputGrouping", CellHorizontalScrolling->True, DefaultFormatType->DefaultInputFormatType, AutoItalicWords->{}, FormatType->InputForm, ShowStringCharacters->True, NumberMarks->True, CounterIncrements->"Input", StyleMenuListing->None, FontWeight->"Bold"], Cell[CellGroupData[{ Cell[StyleData["Output"], PageWidth->Infinity, CellMargins->{{42, 10}, {7, 5}}, CellEditDuplicate->True, Evaluatable->False, CellGroupingRules->"OutputGrouping", CellHorizontalScrolling->True, PageBreakWithin->False, GroupPageBreakWithin->False, GeneratedCell->True, CellAutoOverwrite->True, CellLabelMargins->{{23, Inherited}, {Inherited, Inherited}}, DefaultFormatType->DefaultOutputFormatType, LineSpacing->{1, 0}, AutoItalicWords->{}, FormatType->InputForm, CounterIncrements->"Output", AspectRatioFixed->True, FontFamily->"courier", FontSize->12], Cell[StyleData["Output", "Presentation"], CellMargins->{{72, Inherited}, {10, 8}}, LineSpacing->{1, 0}], Cell[StyleData["Output", "Condensed"], CellMargins->{{41, Inherited}, {3, 2}}], Cell[StyleData["Output", "Printout"], CellMargins->{{39, 0}, {6, 4}}, FontSize->9] }, Open ]], Cell[CellGroupData[{ Cell[StyleData["Message"], PageWidth->Infinity, CellMargins->{{42, Inherited}, {Inherited, Inherited}}, Evaluatable->False, CellGroupingRules->"OutputGrouping", PageBreakWithin->False, GroupPageBreakWithin->False, GeneratedCell->True, CellAutoOverwrite->True, ShowCellLabel->False, CellLabelMargins->{{23, Inherited}, {Inherited, Inherited}}, DefaultFormatType->DefaultOutputFormatType, AutoItalicWords->{}, FormatType->InputForm, CounterIncrements->"Message", AspectRatioFixed->True, StyleMenuListing->None, FontFamily->"courier", FontSize->12, FontColor->RGBColor[0, 0, 1]], Cell[StyleData["Message", "Presentation"], CellMargins->{{72, Inherited}, {Inherited, Inherited}}, LineSpacing->{1, 0}], Cell[StyleData["Message", "Condensed"], CellMargins->{{41, Inherited}, {Inherited, Inherited}}], Cell[StyleData["Message", "Printout"], CellMargins->{{39, Inherited}, {Inherited, Inherited}}, FontSize->8, FontColor->GrayLevel[0]] }, Open ]], Cell[CellGroupData[{ Cell[StyleData["Print"], PageWidth->Infinity, CellMargins->{{42, Inherited}, {Inherited, Inherited}}, Evaluatable->False, CellGroupingRules->"OutputGrouping", CellHorizontalScrolling->True, PageBreakWithin->False, GroupPageBreakWithin->False, GeneratedCell->True, CellAutoOverwrite->True, ShowCellLabel->False, CellLabelMargins->{{23, Inherited}, {Inherited, Inherited}}, DefaultFormatType->DefaultOutputFormatType, AutoItalicWords->{}, FormatType->InputForm, CounterIncrements->"Print", AspectRatioFixed->True, StyleMenuListing->None, FontFamily->"courier", FontSize->12], Cell[StyleData["Print", "Presentation"], CellMargins->{{72, Inherited}, {Inherited, Inherited}}, LineSpacing->{1, 0}], Cell[StyleData["Print", "Condensed"], CellMargins->{{41, Inherited}, {Inherited, Inherited}}], Cell[StyleData["Print", "Printout"], CellMargins->{{39, Inherited}, {Inherited, Inherited}}, FontSize->8] }, Open ]], Cell[CellGroupData[{ Cell[StyleData["Graphics"], PageWidth->Infinity, CellMargins->{{7, Inherited}, {Inherited, Inherited}}, Evaluatable->False, CellGroupingRules->"GraphicsGrouping", CellHorizontalScrolling->True, PageBreakWithin->False, GroupPageBreakWithin->False, GeneratedCell->True, CellAutoOverwrite->True, ShowCellLabel->False, DefaultFormatType->DefaultOutputFormatType, FormatType->InputForm, CounterIncrements->"Graphics", AspectRatioFixed->True, ImageSize->{282, 287}, ImageMargins->{{34, Inherited}, {Inherited, 0}}, StyleMenuListing->None, FontFamily->"courier", FontSize->12], Cell[StyleData["Graphics", "Presentation"], ImageMargins->{{62, Inherited}, {Inherited, 0}}], Cell[StyleData["Graphics", "Condensed"], ImageSize->{175, 175}, ImageMargins->{{38, Inherited}, {Inherited, 0}}], Cell[StyleData["Graphics", "Printout"], ImageSize->{250, 250}, ImageMargins->{{30, Inherited}, {Inherited, 0}}, FontSize->9] }, Open ]], Cell[CellGroupData[{ Cell[StyleData["CellLabel"], StyleMenuListing->None, FontFamily->"Helvetica", FontSize->9, FontColor->RGBColor[0, 0, 1]], Cell[StyleData["CellLabel", "Presentation"], FontSize->12], Cell[StyleData["CellLabel", "Condensed"], FontSize->9], Cell[StyleData["CellLabel", "Printout"], FontFamily->"Courier", FontSize->8, FontSlant->"Italic", FontColor->GrayLevel[0]] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Formulas and Programming", "Section"], Cell[CellGroupData[{ Cell[StyleData["InlineFormula"], CellMargins->{{10, 4}, {0, 8}}, CellHorizontalScrolling->True, ScriptLevel->1, SingleLetterItalics->True], Cell[StyleData["InlineFormula", "Presentation"], CellMargins->{{24, 10}, {10, 10}}, LineSpacing->{1, 5}], Cell[StyleData["InlineFormula", "Condensed"], CellMargins->{{8, 10}, {6, 6}}, LineSpacing->{1, 1}], Cell[StyleData["InlineFormula", "Printout"], CellMargins->{{2, 0}, {6, 6}}] }, Open ]], Cell[CellGroupData[{ Cell[StyleData["DisplayFormula"], CellMargins->{{42, Inherited}, {Inherited, Inherited}}, CellHorizontalScrolling->True, ScriptLevel->0, SingleLetterItalics->True, StyleMenuListing->None, UnderoverscriptBoxOptions->{LimitsPositioning->True}], Cell[StyleData["DisplayFormula", "Presentation"], LineSpacing->{1, 5}], Cell[StyleData["DisplayFormula", "Condensed"], LineSpacing->{1, 1}], Cell[StyleData["DisplayFormula", "Printout"]] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Styles for Headers and Footers", "Section"], Cell[StyleData["Header"], CellMargins->{{7, 0}, {4, 1}}, Evaluatable->False, CellHorizontalScrolling->False, PageBreakWithin->Automatic, AspectRatioFixed->True, StyleMenuListing->None, FontFamily->"times", FontSize->12, FontSlant->"Italic"], Cell[StyleData["Footer"], CellMargins->{{7, 0}, {0, 4}}, Evaluatable->False, CellHorizontalScrolling->False, PageBreakWithin->Automatic, TextAlignment->Center, AspectRatioFixed->True, StyleMenuListing->None, FontFamily->"times", FontSize->12, FontSlant->"Italic"], Cell[StyleData["PageNumber"], CellMargins->{{0, 0}, {4, 1}}, StyleMenuListing->None, FontFamily->"Times", FontSize->10] }, Open ]], Cell[CellGroupData[{ Cell["Palette Styles", "Section"], Cell["\<\ The cells below define styles that define standard \ ButtonFunctions, for use in palette buttons.\ \>", "Text"], Cell[StyleData["Paste"], StyleMenuListing->None, ButtonStyleMenuListing->Automatic, ButtonBoxOptions->{ButtonFunction:>(FrontEndExecute[ { FrontEnd`NotebookApply[ FrontEnd`InputNotebook[ ], #, After]}]&)}], Cell[StyleData["Evaluate"], StyleMenuListing->None, ButtonStyleMenuListing->Automatic, ButtonBoxOptions->{ButtonFunction:>(FrontEndExecute[ { FrontEnd`NotebookApply[ FrontEnd`InputNotebook[ ], #, All], SelectionEvaluate[ FrontEnd`InputNotebook[ ], All]}]&)}], Cell[StyleData["EvaluateCell"], StyleMenuListing->None, ButtonStyleMenuListing->Automatic, ButtonBoxOptions->{ButtonFunction:>(FrontEndExecute[ { FrontEnd`NotebookApply[ FrontEnd`InputNotebook[ ], #, All], FrontEnd`SelectionMove[ FrontEnd`InputNotebook[ ], All, Cell, 1], FrontEnd`SelectionEvaluateCreateCell[ FrontEnd`InputNotebook[ ], All]}]&)}], Cell[StyleData["CopyEvaluate"], StyleMenuListing->None, ButtonStyleMenuListing->Automatic, ButtonBoxOptions->{ButtonFunction:>(FrontEndExecute[ { FrontEnd`SelectionCreateCell[ FrontEnd`InputNotebook[ ], All], FrontEnd`NotebookApply[ FrontEnd`InputNotebook[ ], #, All], FrontEnd`SelectionEvaluate[ FrontEnd`InputNotebook[ ], All]}]&)}], Cell[StyleData["CopyEvaluateCell"], StyleMenuListing->None, ButtonStyleMenuListing->Automatic, ButtonBoxOptions->{ButtonFunction:>(FrontEndExecute[ { FrontEnd`SelectionCreateCell[ FrontEnd`InputNotebook[ ], All], FrontEnd`NotebookApply[ FrontEnd`InputNotebook[ ], #, All], FrontEnd`SelectionEvaluateCreateCell[ FrontEnd`InputNotebook[ ], All]}]&)}] }, Open ]], Cell[CellGroupData[{ Cell["Hyperlink Styles", "Section"], Cell["\<\ The cells below define styles useful for making hypertext \ ButtonBoxes. The \"Hyperlink\" style is for links within the same Notebook, \ or between Notebooks.\ \>", "Text"], Cell[CellGroupData[{ Cell[StyleData["Hyperlink"], StyleMenuListing->None, ButtonStyleMenuListing->Automatic, FontColor->RGBColor[0, 0, 1], FontVariations->{"Underline"->True}, ButtonBoxOptions->{ButtonFunction:>(FrontEndExecute[ { FrontEnd`NotebookLocate[ #2]}]&), Active->True, ButtonNote->ButtonData}], Cell[StyleData["Hyperlink", "Presentation"]], Cell[StyleData["Hyperlink", "Condensed"]], Cell[StyleData["Hyperlink", "Printout"], FontColor->GrayLevel[0], FontVariations->{"Underline"->False}] }, Open ]], Cell["\<\ The following styles are for linking automatically to the on-line \ help system.\ \>", "Text"], Cell[CellGroupData[{ Cell[StyleData["MainBookLink"], StyleMenuListing->None, ButtonStyleMenuListing->Automatic, FontColor->RGBColor[0, 0, 1], FontVariations->{"Underline"->True}, ButtonBoxOptions->{ButtonFunction:>(FrontEndExecute[ { FrontEnd`HelpBrowserLookup[ "MainBook", #]}]&), Active->True, ButtonFrame->"None"}], Cell[StyleData["MainBookLink", "Presentation"]], Cell[StyleData["MainBookLink", "Condensed"]], Cell[StyleData["MainBookLink", "Printout"], FontColor->GrayLevel[0], FontVariations->{"Underline"->False}] }, Open ]], Cell[CellGroupData[{ Cell[StyleData["AddOnsLink"], StyleMenuListing->None, ButtonStyleMenuListing->Automatic, FontFamily->"Courier", FontColor->RGBColor[0, 0, 1], FontVariations->{"Underline"->True}, ButtonBoxOptions->{ButtonFunction:>(FrontEndExecute[ { FrontEnd`HelpBrowserLookup[ "AddOns", #]}]&), Active->True, ButtonFrame->"None"}], Cell[StyleData["AddOnsLink", "Presentation"]], Cell[StyleData["AddOnsLink", "Condensed"]], Cell[StyleData["AddOnLink", "Printout"], FontColor->GrayLevel[0], FontVariations->{"Underline"->False}] }, Open ]], Cell[CellGroupData[{ Cell[StyleData["RefGuideLink"], StyleMenuListing->None, ButtonStyleMenuListing->Automatic, FontFamily->"Courier", FontColor->RGBColor[0, 0, 1], FontVariations->{"Underline"->True}, ButtonBoxOptions->{ButtonFunction:>(FrontEndExecute[ { FrontEnd`HelpBrowserLookup[ "RefGuideLink", #]}]&), Active->True, ButtonFrame->"None"}], Cell[StyleData["RefGuideLink", "Presentation"]], Cell[StyleData["RefGuideLink", "Condensed"]], Cell[StyleData["RefGuideLink", "Printout"], FontColor->GrayLevel[0], FontVariations->{"Underline"->False}] }, Open ]], Cell[CellGroupData[{ Cell[StyleData["GettingStartedLink"], StyleMenuListing->None, ButtonStyleMenuListing->Automatic, FontColor->RGBColor[0, 0, 1], FontVariations->{"Underline"->True}, ButtonBoxOptions->{ButtonFunction:>(FrontEndExecute[ { FrontEnd`HelpBrowserLookup[ "GettingStarted", #]}]&), Active->True, ButtonFrame->"None"}], Cell[StyleData["GettingStartedLink", "Presentation"]], Cell[StyleData["GettingStartedLink", "Condensed"]], Cell[StyleData["GettingStartedLink", "Printout"], FontColor->GrayLevel[0], FontVariations->{"Underline"->False}] }, Open ]], Cell[CellGroupData[{ Cell[StyleData["OtherInformationLink"], StyleMenuListing->None, ButtonStyleMenuListing->Automatic, FontColor->RGBColor[0, 0, 1], FontVariations->{"Underline"->True}, ButtonBoxOptions->{ButtonFunction:>(FrontEndExecute[ { FrontEnd`HelpBrowserLookup[ "OtherInformation", #]}]&), Active->True, ButtonFrame->"None"}], Cell[StyleData["OtherInformationLink", "Presentation"]], Cell[StyleData["OtherInformationLink", "Condensed"]], Cell[StyleData["OtherInformationLink", "Printout"], FontColor->GrayLevel[0], FontVariations->{"Underline"->False}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Placeholder Styles", "Section"], Cell["\<\ The cells below define styles useful for making placeholder \ objects in palette templates.\ \>", "Text"], Cell[CellGroupData[{ Cell[StyleData["Placeholder"], Editable->False, Selectable->False, StyleBoxAutoDelete->True, Placeholder->True, StyleMenuListing->None], Cell[StyleData["Placeholder", "Presentation"]], Cell[StyleData["Placeholder", "Condensed"]], Cell[StyleData["Placeholder", "Printout"]] }, Open ]], Cell[CellGroupData[{ Cell[StyleData["SelectionPlaceholder"], Editable->False, Selectable->False, StyleBoxAutoDelete->True, Placeholder->Primary, StyleMenuListing->None, DrawHighlighted->True], Cell[StyleData["SelectionPlaceholder", "Presentation"]], Cell[StyleData["SelectionPlaceholder", "Condensed"]], Cell[StyleData["SelectionPlaceholder", "Printout"]] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["FormatType Styles", "Section"], Cell["\<\ The cells below define styles that are mixed in with the styles \ of most cells. If a cell's FormatType matches the name of one of the styles \ defined below, then that style is applied between the cell's style and its \ own options.\ \>", "Text"], Cell[StyleData["CellExpression"], PageWidth->Infinity, CellMargins->{{6, Inherited}, {Inherited, Inherited}}, ShowCellLabel->False, ShowSpecialCharacters->False, AllowInlineCells->False, AutoItalicWords->{}, StyleMenuListing->None, FontFamily->"Courier", Background->GrayLevel[1]], Cell[StyleData["InputForm"], AllowInlineCells->False, StyleMenuListing->None, FontFamily->"Courier"], Cell[StyleData["OutputForm"], PageWidth->Infinity, TextAlignment->Left, LineSpacing->{1, -5}, StyleMenuListing->None, FontFamily->"Courier"], Cell[StyleData["StandardForm"], LineSpacing->{1.25, 0}, StyleMenuListing->None, FontFamily->"Courier"], Cell[StyleData["TraditionalForm"], LineSpacing->{1.25, 0}, SingleLetterItalics->True, TraditionalFunctionNotation->True, DelimiterMatching->None, StyleMenuListing->None], Cell["\<\ The style defined below is mixed in to any cell that is in an \ inline cell within another.\ \>", "Text"], Cell[StyleData["InlineCell"], TextAlignment->Left, ScriptLevel->1, StyleMenuListing->None], Cell[StyleData["InlineCellEditing"], StyleMenuListing->None, Background->RGBColor[1, 0.749996, 0.8]] }, Open ]] }, Open ]] }] ] (*********************************************************************** Cached data follows. If you edit this Notebook file directly, not using Mathematica, you must remove the line containing CacheID at the top of the file. The cache data will then be recreated when you save this file from within Mathematica. ***********************************************************************) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[1731, 51, 96, 2, 76, "Title", Evaluatable->False], Cell[1830, 55, 319, 6, 104, "Text"], Cell[CellGroupData[{ Cell[2174, 65, 96, 2, 42, "Subtitle", Evaluatable->False], Cell[CellGroupData[{ Cell[2295, 71, 91, 2, 49, "Section", Evaluatable->False], Cell[2389, 75, 94, 2, 27, "Input", InitializationCell->True], Cell[CellGroupData[{ Cell[2508, 81, 94, 2, 40, "Subsection", Evaluatable->False], Cell[CellGroupData[{ Cell[2627, 87, 74, 2, 27, "Input", InitializationCell->True], Cell[2704, 91, 35, 1, 27, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[2788, 98, 127, 3, 40, "Subsection", Evaluatable->False, InitializationCell->True], Cell[2918, 103, 90, 2, 27, "Input", InitializationCell->True], Cell[3011, 107, 358, 14, 177, "Input", InitializationCell->True], Cell[CellGroupData[{ Cell[3394, 125, 99, 2, 27, "Input", InitializationCell->True], Cell[3496, 129, 358, 5, 27, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3891, 139, 91, 2, 27, "Input", InitializationCell->True], Cell[3985, 143, 36, 1, 27, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[4070, 150, 108, 2, 40, "Subsection", Evaluatable->False], Cell[CellGroupData[{ Cell[4203, 156, 94, 2, 36, "Subsubsection", Evaluatable->False], Cell[CellGroupData[{ Cell[4322, 162, 79, 1, 27, "Input"], Cell[4404, 165, 449, 6, 27, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[4890, 176, 61, 1, 27, "Input"], Cell[4954, 179, 36, 1, 27, "Output"] }, Open ]], Cell[5005, 183, 155, 5, 57, "Input"], Cell[5163, 190, 148, 4, 32, "Text", Evaluatable->False], Cell[CellGroupData[{ Cell[5336, 198, 165, 6, 72, "Input"], Cell[5504, 206, 1431, 27, 183, "Output"] }, Open ]], Cell[6950, 236, 101, 2, 32, "Text", Evaluatable->False], Cell[7054, 240, 97, 4, 42, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[7188, 249, 119, 3, 36, "Subsubsection", Evaluatable->False, InitializationCell->True], Cell[7310, 254, 115, 2, 27, "Input", InitializationCell->True], Cell[7428, 258, 129, 2, 27, "Input", InitializationCell->True], Cell[7560, 262, 103, 2, 27, "Input", InitializationCell->True], Cell[7666, 266, 156, 5, 42, "Input", InitializationCell->True], Cell[CellGroupData[{ Cell[7847, 275, 152, 5, 42, "Input", InitializationCell->True], Cell[8002, 282, 54, 1, 27, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8093, 288, 86, 2, 27, "Input", InitializationCell->True], Cell[8182, 292, 35, 1, 27, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8254, 298, 86, 2, 27, "Input", InitializationCell->True], Cell[8343, 302, 50, 1, 27, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8430, 308, 155, 5, 42, "Input", InitializationCell->True], Cell[8588, 315, 54, 1, 27, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8679, 321, 212, 8, 87, "Input", InitializationCell->True], Cell[8894, 331, 479, 13, 117, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[9410, 349, 89, 2, 27, "Input", InitializationCell->True], Cell[9502, 353, 35, 1, 27, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[9574, 359, 293, 11, 132, "Input", InitializationCell->True], Cell[9870, 372, 635, 14, 98, "Output"] }, Open ]], Cell[10520, 389, 139, 3, 32, "Text", Evaluatable->False, InitializationCell->True], Cell[CellGroupData[{ Cell[10684, 396, 117, 5, 42, "Input", InitializationCell->True], Cell[10804, 403, 35, 1, 27, "Output"] }, Open ]] }, Closed]] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[10912, 412, 83, 2, 33, "Section", Evaluatable->False], Cell[10998, 416, 241, 6, 57, "Input", InitializationCell->True], Cell[CellGroupData[{ Cell[11264, 426, 99, 2, 27, "Input", InitializationCell->True], Cell[11366, 430, 1682, 52, 27, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[13085, 487, 99, 2, 40, "Subsection", Evaluatable->False], Cell[CellGroupData[{ Cell[13209, 493, 89, 2, 36, "Subsubsection", Evaluatable->False], Cell[CellGroupData[{ Cell[13323, 499, 159, 6, 72, "Input"], Cell[13485, 507, 1848, 43, 669, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[15382, 556, 76, 2, 36, "Subsubsection", Evaluatable->False], Cell[15461, 560, 192, 7, 72, "Input"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[15702, 573, 89, 2, 40, "Subsection", Evaluatable->False], Cell[CellGroupData[{ Cell[15816, 579, 121, 3, 36, "Subsubsection", Evaluatable->False, InitializationCell->True], Cell[15940, 584, 266, 10, 117, "Input", InitializationCell->True] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[16255, 600, 91, 2, 40, "Subsection", Evaluatable->False], Cell[CellGroupData[{ Cell[16371, 606, 121, 3, 36, "Subsubsection", Evaluatable->False, InitializationCell->True], Cell[16495, 611, 272, 10, 117, "Input", InitializationCell->True], Cell[16770, 623, 103, 2, 27, "Input", InitializationCell->True], Cell[16876, 627, 107, 2, 27, "Input", InitializationCell->True], Cell[CellGroupData[{ Cell[17008, 633, 40, 1, 27, "Input"], Cell[17051, 636, 149, 2, 27, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[17237, 643, 40, 1, 27, "Input"], Cell[17280, 646, 314, 5, 45, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[17631, 656, 54, 1, 27, "Input"], Cell[17688, 659, 8630, 326, 908, "Output"] }, Closed]], Cell[26333, 988, 46, 1, 24, "Input"], Cell[26382, 991, 46, 1, 27, "Input"], Cell[CellGroupData[{ Cell[26453, 996, 84, 1, 27, "Input"], Cell[26540, 999, 2493, 54, 657, "Output"] }, Closed]], Cell[29048, 1056, 46, 1, 24, "Input"], Cell[29097, 1059, 46, 1, 27, "Input"], Cell[29146, 1062, 64, 1, 27, "Input"] }, Open ]] }, Open ]] }, Open ]] }, Open ]] }, Open ]] } ] *) (*********************************************************************** End of Mathematica Notebook file. ***********************************************************************)