(* 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[ 10213, 347] NotebookOptionsPosition[ 9159, 306] NotebookOutlinePosition[ 9520, 322] CellTagsIndexPosition[ 9477, 319] WindowFrame->Normal ContainsDynamic->False*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"lhs", "=", RowBox[{"Expand", "[", RowBox[{ RowBox[{"(", RowBox[{"1", "-", RowBox[{"0.57", " ", "B"}]}], ")"}], " ", RowBox[{"(", RowBox[{"1", "-", RowBox[{"1.53", " ", "B"}], "+", RowBox[{"0.63", " ", RowBox[{"B", "^", "2"}]}]}], ")"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.456828184003887*^9, 3.4568282013946238`*^9}}], Cell[BoxData[ RowBox[{"1", "-", RowBox[{"2.1`", " ", "B"}], "+", RowBox[{"1.5021`", " ", SuperscriptBox["B", "2"]}], "-", RowBox[{"0.3591`", " ", SuperscriptBox["B", "3"]}]}]], "Output", CellChangeTimes->{{3.4568281854882717`*^9, 3.4568282022852545`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"rhs", "=", RowBox[{"Expand", "[", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"0.53", "+", RowBox[{"0.37", " ", "B"}], "+", RowBox[{"0.51", " ", RowBox[{"B", "^", "2"}]}]}], ")"}]}], " ", RowBox[{"(", RowBox[{"1", "-", RowBox[{"1.53", " ", "B"}], "+", RowBox[{"0.63", " ", RowBox[{"B", "^", "2"}]}]}], ")"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.4568282034415116`*^9, 3.4568282215822525`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"-", "0.53`"}], "+", RowBox[{"0.44090000000000007`", " ", "B"}], "-", RowBox[{"0.2778`", " ", SuperscriptBox["B", "2"]}], "+", RowBox[{"0.5472`", " ", SuperscriptBox["B", "3"]}], "-", RowBox[{"0.32130000000000003`", " ", SuperscriptBox["B", "4"]}]}]], "Output", CellChangeTimes->{3.45682822347289*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"u", "=", FractionBox[ RowBox[{"-", " ", RowBox[{"(", RowBox[{"0.53", "+", RowBox[{"0.37", " ", "B"}], " ", "+", RowBox[{"0.51", " ", SuperscriptBox["B", "2"]}]}], ")"}]}], RowBox[{ RowBox[{"(", RowBox[{"1", "-", RowBox[{"0.57", " ", "B"}]}], ")"}], RowBox[{"(", RowBox[{"1", "-", RowBox[{"1.97", " ", "B"}], "+", RowBox[{"1.37", " ", SuperscriptBox["B", "2"]}], "-", RowBox[{"0.34", " ", SuperscriptBox["B", "3"]}]}], ")"}]}]]}], ";"}], "\[IndentingNewLine]", RowBox[{"uc", "=", RowBox[{"Series", "[", RowBox[{"u", ",", RowBox[{"{", RowBox[{"B", ",", "0", ",", "5"}], "}"}]}], "]"}]}], "\[IndentingNewLine]", RowBox[{"Map", "[", RowBox[{ RowBox[{ RowBox[{"Coefficient", "[", RowBox[{"uc", ",", "B", ",", "#"}], "]"}], "&"}], ",", RowBox[{"Range", "[", "5", "]"}]}], "]"}]}], "Input", CellChangeTimes->{{3.456851876729271*^9, 3.4568519429644938`*^9}, { 3.4568525291594963`*^9, 3.4568526075198746`*^9}}], Cell[BoxData[ InterpretationBox[ RowBox[{ RowBox[{"-", "0.53`"}], "-", RowBox[{"1.7162`", " ", "B"}], "-", RowBox[{"3.547911`", " ", SuperscriptBox["B", "2"]}], "-", RowBox[{"5.327455959999999`", " ", SuperscriptBox["B", "3"]}], "-", RowBox[{"6.5081253865`", " ", SuperscriptBox["B", "4"]}], "-", RowBox[{"6.894077398925997`", " ", SuperscriptBox["B", "5"]}], "+", InterpretationBox[ SuperscriptBox[ RowBox[{"O", "[", "B", "]"}], "6"], SeriesData[$CellContext`B, 0, {}, 0, 6, 1], Editable->False]}], SeriesData[$CellContext`B, 0, {-0.53, -1.7162, -3.547911, -5.327455959999999, -6.5081253865, \ -6.894077398925997}, 0, 6, 1], Editable->False]], "Output", CellChangeTimes->{{3.4568525354877024`*^9, 3.456852557159855*^9}, { 3.4568525944572077`*^9, 3.45685260792613*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", "1.7162`"}], ",", RowBox[{"-", "3.547911`"}], ",", RowBox[{"-", "5.327455959999999`"}], ",", RowBox[{"-", "6.5081253865`"}], ",", RowBox[{"-", "6.894077398925997`"}]}], "}"}]], "Output", CellChangeTimes->{{3.4568525354877024`*^9, 3.456852557159855*^9}, { 3.4568525944572077`*^9, 3.45685260792613*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"p", "=", FractionBox["1", RowBox[{"(", RowBox[{"1", "-", RowBox[{"1.97", " ", "B"}], "+", RowBox[{"1.37", " ", SuperscriptBox["B", "2"]}], "-", RowBox[{"0.34", " ", SuperscriptBox["B", "3"]}]}], ")"}]]}], ";"}], "\[IndentingNewLine]", RowBox[{"pc", "=", RowBox[{"Series", "[", RowBox[{"p", ",", RowBox[{"{", RowBox[{"B", ",", "0", ",", "5"}], "}"}]}], "]"}]}], "\[IndentingNewLine]", RowBox[{"Map", "[", RowBox[{ RowBox[{ RowBox[{"Coefficient", "[", RowBox[{"pc", ",", "B", ",", "#"}], "]"}], "&"}], ",", RowBox[{"Range", "[", "5", "]"}]}], "]"}]}], "Input", CellChangeTimes->{{3.456851876729271*^9, 3.4568519429644938`*^9}, { 3.4568525291594963`*^9, 3.4568526542079725`*^9}}], Cell[BoxData[ InterpretationBox[ RowBox[{"1", "+", RowBox[{"1.97`", " ", "B"}], "+", RowBox[{"2.5109`", " ", SuperscriptBox["B", "2"]}], "+", RowBox[{"2.587573`", " ", SuperscriptBox["B", "3"]}], "+", RowBox[{"2.327385809999999`", " ", SuperscriptBox["B", "4"]}], "+", RowBox[{"1.8936810356999985`", " ", SuperscriptBox["B", "5"]}], "+", InterpretationBox[ SuperscriptBox[ RowBox[{"O", "[", "B", "]"}], "6"], SeriesData[$CellContext`B, 0, {}, 0, 6, 1], Editable->False]}], SeriesData[$CellContext`B, 0, { 1, 1.97, 2.5109, 2.587573, 2.327385809999999, 1.8936810356999985`}, 0, 6, 1], Editable->False]], "Output", CellChangeTimes->{3.4568526546767282`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ "1.97`", ",", "2.5109`", ",", "2.587573`", ",", "2.327385809999999`", ",", "1.8936810356999985`"}], "}"}]], "Output", CellChangeTimes->{3.4568526546767282`*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"xya", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "2.378"}], ",", "52.4"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "2.499"}], ",", " ", "53.5"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "2.473"}], " ", ",", "55.6"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "2.330"}], ",", "58.0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "2.053"}], " ", ",", "59.5"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1.739"}], ",", " ", "60.0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1.261"}], ",", "60.4"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", ".569"}], ",", " ", "60.5"}], "}"}]}], "}"}]}], ";"}]], "Input", CellChangeTimes->{{3.456845799285078*^9, 3.456845883567407*^9}, { 3.45684622168111*^9, 3.4568462281968184`*^9}, 3.456852236968257*^9}], Cell[BoxData[{ RowBox[{ RowBox[{"xa", "=", RowBox[{"Map", "[", RowBox[{"First", ",", "xya"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"ya", "=", RowBox[{"Map", "[", RowBox[{"Last", ",", "xya"}], "]"}]}], ";"}]}], "Input", CellChangeTimes->{{3.45684623183749*^9, 3.4568462645566587`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ListPlot", "[", RowBox[{"xa", ",", RowBox[{"PlotJoined", "\[Rule]", "True"}]}], "]"}]], "Input", CellChangeTimes->{{3.456851851525823*^9, 3.4568518588384166`*^9}}], Cell[BoxData[ GraphicsBox[{{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[{{1., -2.378}, {2., -2.499}, {3., -2.473}, {4., -2.33}, { 5., -2.053}, {6., -1.739}, {7., -1.261}, {8., -0.569}}]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, -2.5}, PlotRange->{{0, 8.}, {-2.499, -0.569}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output", CellChangeTimes->{3.456851859182171*^9, 3.45685223998392*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ListPlot", "[", RowBox[{"ya", ",", RowBox[{"PlotJoined", "\[Rule]", "True"}]}], "]"}]], "Input", CellChangeTimes->{{3.456851862354087*^9, 3.456851868635417*^9}}], Cell[BoxData[ GraphicsBox[{{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[{{1., 52.4}, {2., 53.5}, {3., 55.6}, {4., 58.}, {5., 59.5}, {6., 60.}, {7., 60.4}, {8., 60.5}}]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 54.}, PlotRange->{{0, 8.}, {52.4, 60.5}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output", CellChangeTimes->{3.4568518689635463`*^9, 3.456852242702705*^9}] }, Open ]] }, WindowSize->{616, 750}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, ShowSelection->True, FrontEndVersion->"6.0 for Microsoft Windows (32-bit) (April 20, 2007)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[590, 23, 406, 12, 31, "Input"], Cell[999, 37, 274, 7, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[1310, 49, 504, 15, 31, "Input"], Cell[1817, 66, 360, 10, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[2214, 81, 1115, 35, 99, "Input"], Cell[3332, 118, 847, 23, 30, "Output"], Cell[4182, 143, 376, 9, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[4595, 157, 813, 24, 93, "Input"], Cell[5411, 183, 730, 21, 30, "Output"], Cell[6144, 206, 207, 5, 30, "Output"] }, Open ]], Cell[6366, 214, 1008, 31, 52, "Input"], Cell[7377, 247, 329, 9, 52, "Input"], Cell[CellGroupData[{ Cell[7731, 260, 195, 4, 31, "Input"], Cell[7929, 266, 500, 13, 239, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8466, 284, 193, 4, 31, "Input"], Cell[8662, 290, 481, 13, 234, "Output"] }, Open ]] } ] *) (* End of internal cache information *)