(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 6.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 7979, 261] NotebookOptionsPosition[ 7423, 238] NotebookOutlinePosition[ 7802, 254] CellTagsIndexPosition[ 7759, 251] WindowFrame->Normal ContainsDynamic->False*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[{ StyleBox[ RowBox[{ RowBox[{"h", " ", "=", " ", "0.1"}], ";"}], FontSize->18], "\[IndentingNewLine]", StyleBox[ RowBox[{ RowBox[{"Ap", " ", "=", " ", RowBox[{"{", RowBox[{"1.8", ",", " ", "0.7", ",", " ", "0"}], "}"}]}], ";"}], FontSize->18], "\[IndentingNewLine]", StyleBox[ RowBox[{ RowBox[{"Bp", " ", "=", " ", RowBox[{"{", RowBox[{"0.6", ",", " ", RowBox[{"h", " ", "+", " ", "0.4"}], ",", " ", "0"}], "}"}]}], ";"}], FontSize->18], "\[IndentingNewLine]", StyleBox[ RowBox[{ RowBox[{"Cp", " ", "=", " ", RowBox[{"{", RowBox[{"0.6", ",", " ", "0", ",", " ", "0"}], "}"}]}], ";"}], FontSize->18], "\[IndentingNewLine]", StyleBox[ RowBox[{ RowBox[{"Dp", " ", "=", " ", RowBox[{"{", RowBox[{"0", ",", " ", "0.4", ",", " ", "0"}], "}"}]}], ";"}], FontSize->18], "\[IndentingNewLine]", StyleBox[ RowBox[{ RowBox[{"eBA", " ", "=", " ", RowBox[{ RowBox[{"(", RowBox[{"Ap", " ", "-", " ", "Bp"}], ")"}], " ", "/", " ", SqrtBox[ RowBox[{ RowBox[{"(", RowBox[{"Ap", " ", "-", " ", "Bp"}], ")"}], ".", RowBox[{"(", RowBox[{"Ap", " ", "-", " ", "Bp"}], ")"}]}]]}]}], ";"}], FontSize->18], "\[IndentingNewLine]", StyleBox[ RowBox[{ RowBox[{"eBD", " ", "=", " ", RowBox[{ RowBox[{"(", RowBox[{"Dp", " ", "-", " ", "Bp"}], ")"}], " ", "/", " ", SqrtBox[ RowBox[{ RowBox[{"(", RowBox[{"Dp", " ", "-", " ", "Bp"}], ")"}], ".", RowBox[{"(", RowBox[{"Dp", " ", "-", " ", "Bp"}], ")"}]}]]}]}], ";"}], FontSize->18], "\[IndentingNewLine]", StyleBox[ RowBox[{ RowBox[{"eBC", " ", "=", " ", RowBox[{ RowBox[{"(", RowBox[{"Cp", " ", "-", " ", "Bp"}], ")"}], " ", "/", " ", SqrtBox[ RowBox[{ RowBox[{"(", RowBox[{"Cp", " ", "-", " ", "Bp"}], ")"}], ".", RowBox[{"(", RowBox[{"Cp", " ", "-", " ", "Bp"}], ")"}]}]]}]}], ";"}], FontSize->18], "\[IndentingNewLine]", StyleBox[ RowBox[{ RowBox[{"eAC", " ", "=", " ", RowBox[{ RowBox[{"(", RowBox[{"Cp", " ", "-", " ", "Ap"}], ")"}], " ", "/", " ", SqrtBox[ RowBox[{ RowBox[{"(", RowBox[{"Cp", " ", "-", " ", "Ap"}], ")"}], ".", RowBox[{"(", RowBox[{"Cp", " ", "-", " ", "Ap"}], ")"}]}]]}]}], ";"}], FontSize->18], "\[IndentingNewLine]", StyleBox[ RowBox[{ RowBox[{"eCD", " ", "=", " ", RowBox[{ RowBox[{"(", RowBox[{"Dp", " ", "-", " ", "Cp"}], ")"}], " ", "/", " ", SqrtBox[ RowBox[{ RowBox[{"(", RowBox[{"Dp", " ", "-", " ", "Cp"}], ")"}], ".", RowBox[{"(", RowBox[{"Dp", " ", "-", " ", "Cp"}], ")"}]}]]}]}], ";"}], FontSize->18], "\[IndentingNewLine]", StyleBox[ RowBox[{ RowBox[{"Solve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"TBA", " ", "eBA"}], " ", "+", " ", RowBox[{"TBC", " ", "eBC"}], " ", "+", " ", RowBox[{"TBD", " ", "eBD"}]}], " ", "\[Equal]", " ", "0"}], ",", " ", RowBox[{ RowBox[{ RowBox[{ RowBox[{"-", "TBA"}], " ", "eBA"}], " ", "+", " ", RowBox[{"TAC", " ", "eAC"}], " ", "+", " ", RowBox[{"{", RowBox[{"0", ",", " ", RowBox[{"-", "1"}], ",", " ", "0"}], "}"}]}], " ", "\[Equal]", " ", "0"}], ",", " ", RowBox[{ RowBox[{ RowBox[{"TCD", " ", "eCD"}], " ", "-", " ", RowBox[{"TBC", " ", "eBC"}], " ", "+", " ", RowBox[{"{", RowBox[{"0", ",", " ", "3", ",", " ", "0"}], "}"}], " ", "-", " ", RowBox[{"TAC", " ", "eAC"}]}], " ", "\[Equal]", " ", "0"}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{ "TBA", ",", " ", "TBC", ",", " ", "TBD", ",", " ", "TAC", ",", " ", "TCD"}], "}"}]}], "]"}], " "}], FontSize->18]}], "Input", CellChangeTimes->{{3.434815125013697*^9, 3.434815540719946*^9}, { 3.4348161920587463`*^9, 3.434816196294259*^9}, {3.4348166096041074`*^9, 3.434816622106347*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"TBA", "\[Rule]", "2.433105012119288`"}], ",", RowBox[{"TBC", "\[Rule]", "0.`"}], ",", RowBox[{"TBD", "\[Rule]", "2.433105012119288`"}], ",", RowBox[{"TAC", "\[Rule]", RowBox[{"-", "2.7784887978899606`"}]}], ",", RowBox[{"TCD", "\[Rule]", RowBox[{"-", "2.8844410203711917`"}]}]}], "}"}], "}"}]], "Output", CellChangeTimes->{{3.4348155328611746`*^9, 3.434815541782364*^9}, 3.434816210860674*^9, 3.4348166253881855`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"TBA", " ", "eBA"}], " ", "+", " ", RowBox[{"TBC", " ", "eBC"}], " ", "+", " ", RowBox[{"TBD", " ", "eBD"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"-", "TBA"}], " ", "eBA"}], " ", "+", " ", RowBox[{"TAC", " ", "eAC"}], " ", "+", " ", RowBox[{"{", RowBox[{"0", ",", " ", RowBox[{"-", "1"}], ",", " ", "0"}], "}"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"TCD", " ", "eCD"}], " ", "-", " ", RowBox[{"TBC", " ", "eBC"}], " ", "+", " ", RowBox[{"{", RowBox[{"0", ",", " ", "3", ",", " ", "0"}], "}"}], " ", "-", " ", RowBox[{"TAC", " ", "eAC"}]}]}], "Input", CellChangeTimes->{{3.4348162874906416`*^9, 3.434816313872731*^9}, 3.4348163584628386`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"0.9863939238321437`", " ", "TBA"}], "+", RowBox[{"0.`", " ", "TBC"}], "-", RowBox[{"0.7682212795973759`", " ", "TBD"}]}], ",", RowBox[{ RowBox[{ RowBox[{"-", "0.16439898730535732`"}], " ", "TBA"}], "-", RowBox[{"1.`", " ", "TBC"}], "-", RowBox[{"0.6401843996644799`", " ", "TBD"}]}], ",", "0"}], "}"}]], "Output", CellChangeTimes->{3.4348163628233852`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"-", "0.8637789008984336`"}], " ", "TAC"}], "-", RowBox[{"0.9863939238321437`", " ", "TBA"}]}], ",", RowBox[{ RowBox[{"-", "1"}], "-", RowBox[{"0.5038710255240861`", " ", "TAC"}], "+", RowBox[{"0.16439898730535732`", " ", "TBA"}]}], ",", "0"}], "}"}]], "Output", CellChangeTimes->{3.4348163628390145`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"0.8637789008984336`", " ", "TAC"}], "+", RowBox[{"0.`", " ", "TBC"}], "-", RowBox[{"0.8320502943378436`", " ", "TCD"}]}], ",", RowBox[{"3", "+", RowBox[{"0.5038710255240861`", " ", "TAC"}], "+", RowBox[{"1.`", " ", "TBC"}], "+", RowBox[{"0.5547001962252291`", " ", "TCD"}]}], ",", "0"}], "}"}]], "Output", CellChangeTimes->{3.4348163628390145`*^9}] }, Open ]] }, ScreenStyleEnvironment->"Presentation", WindowSize->{1016, 668}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, FrontEndVersion->"6.0 for Microsoft Windows (32-bit) (May 21, 2008)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[590, 23, 4181, 131, 402, "Input"], Cell[4774, 156, 524, 12, 36, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[5335, 173, 749, 19, 83, "Input"], Cell[6087, 194, 458, 13, 36, "Output"], Cell[6548, 209, 410, 12, 36, "Output"], Cell[6961, 223, 446, 12, 36, "Output"] }, Open ]] } ] *) (* End of internal cache information *)