Cell[CellGroupData[{Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"\[Mu]", "[", "q_", "]"}], ":=", RowBox[{"M1", " ", RowBox[{"q", "/", RowBox[{"(", RowBox[{"M1", " ", "+", RowBox[{"q", " ", "M1"}]}], ")"}]}]}]}], ";"}], " ", RowBox[{ RowBox[{"a", "=", "1"}], ";"}], " ", RowBox[{ RowBox[{"M1", "=", "1"}], ";"}], " ", RowBox[{ RowBox[{"G", "=", "1"}], ";"}], " ", RowBox[{ RowBox[{ RowBox[{"Porb", "[", "q_", "]"}], ":=", RowBox[{"Sqrt", "[", RowBox[{"4", " ", SuperscriptBox["\[Pi]", "2"], " ", RowBox[{ SuperscriptBox["a", "3"], "/", RowBox[{"(", RowBox[{"G", " ", RowBox[{"(", RowBox[{"1", "+", "q"}], ")"}]}], ")"}]}]}], "]"}]}], ";"}], " ", RowBox[{ RowBox[{ RowBox[{"\[CapitalOmega]orb", "[", "q_", "]"}], ":=", RowBox[{"2", " ", RowBox[{"\[Pi]", "/", RowBox[{"Porb", "[", "q", "]"}]}]}]}], ";"}], " ", RowBox[{ RowBox[{ RowBox[{"M2", "[", "q_", "]"}], ":=", RowBox[{"q", " ", "M1"}]}], ";"}], " ", RowBox[{ RowBox[{ RowBox[{"\[CapitalPhi]", "[", RowBox[{"q_", ",", "x_", ",", "y_", ",", "z_"}], "]"}], ":=", RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"G", " ", "M1"}], SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"x", "+", RowBox[{"q", "/", RowBox[{"(", RowBox[{"1", "+", "q"}], ")"}]}]}], ")"}], "2"], "+", SuperscriptBox["y", "2"], "+", SuperscriptBox["z", "2"]}]]]}], "-", FractionBox[ RowBox[{"G", " ", RowBox[{"M2", "[", "q", "]"}]}], SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"x", "-", "a", "+", RowBox[{"q", "/", RowBox[{"(", RowBox[{"1", "+", "q"}], ")"}]}]}], ")"}], "2"], "+", SuperscriptBox["y", "2"], "+", SuperscriptBox["z", "2"]}]]], "-", RowBox[{ FractionBox["1", "2"], " ", SuperscriptBox[ RowBox[{"\[CapitalOmega]orb", "[", "q", "]"}], "2"], " ", RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"x", "-", RowBox[{"\[Mu]", "[", "q", "]"}], " ", "+", RowBox[{ RowBox[{"q", "/", RowBox[{"(", RowBox[{"1", "+", "q"}], ")"}]}], "a"}]}], ")"}], "2"], "+", SuperscriptBox["y", "2"]}], ")"}]}]}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"deriv", "[", "q_", "]"}], ":=", RowBox[{"D", "[", RowBox[{ RowBox[{"\[CapitalPhi]", "[", RowBox[{"q", ",", "x", ",", "0", ",", "0"}], "]"}], ",", "x"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"posAtL1", "[", "q_", "]"}], ":=", RowBox[{"FindRoot", "[", RowBox[{ RowBox[{ RowBox[{"Evaluate", "@", RowBox[{"deriv", "[", "q", "]"}]}], "==", "0"}], ",", RowBox[{"{", RowBox[{"x", ",", ".3"}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"\[CapitalPhi]AtL1", "[", "q_", "]"}], ":=", RowBox[{"\[CapitalPhi]", "[", RowBox[{"q", ",", RowBox[{"x", "/.", RowBox[{"posAtL1", "[", "q", "]"}]}], ",", "0", ",", "0"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"posAtL2", "[", "q_", "]"}], ":=", RowBox[{"FindRoot", "[", RowBox[{ RowBox[{ RowBox[{"Evaluate", "@", RowBox[{"deriv", "[", "q", "]"}]}], "==", "0"}], ",", RowBox[{"{", RowBox[{"x", ",", RowBox[{"-", "1"}]}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"\[CapitalPhi]AtL2", "[", "q_", "]"}], ":=", RowBox[{"\[CapitalPhi]", "[", RowBox[{"q", ",", RowBox[{"x", "/.", RowBox[{"posAtL2", "[", "q", "]"}]}], ",", "0", ",", "0"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"posAtL3", "[", "q_", "]"}], ":=", RowBox[{"FindRoot", "[", RowBox[{ RowBox[{ RowBox[{"Evaluate", "@", RowBox[{"deriv", "[", "q", "]"}]}], "==", "0"}], ",", RowBox[{"{", RowBox[{"x", ",", "1.5"}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"\[CapitalPhi]AtL3", "[", "q_", "]"}], ":=", RowBox[{"\[CapitalPhi]", "[", RowBox[{"q", ",", RowBox[{"x", "/.", RowBox[{"posAtL3", "[", "q", "]"}]}], ",", "0", ",", "0"}], "]"}]}], ";"}]}], "Input", CellGroupingRules->{"GroupTogetherGrouping", 10000.}, CellChangeTimes->{{3.384385377953125*^9, 3.38438551671875*^9}, 3.384386055703125*^9, {3.38438609646875*^9, 3.384386104578125*^9}, {3.38438615090625*^9, 3.38438616753125*^9}, 3.384386539828125*^9, 3.384387105734375*^9, 3.760739909476824*^9}, CellID->88106804], Cell[BoxData[ RowBox[{"Manipulate", "[", " ", RowBox[{ RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"p1", ">", RowBox[{"\[CapitalPhi]AtL1", "[", "q", "]"}]}], ",", RowBox[{"p1", "=", RowBox[{"\[CapitalPhi]AtL1", "[", "q", "]"}]}]}], "]"}], ";", " ", RowBox[{"If", "[", RowBox[{ RowBox[{"p2", ">", RowBox[{"\[CapitalPhi]AtL1", "[", "q", "]"}]}], ",", RowBox[{"p2", "=", RowBox[{"\[CapitalPhi]AtL1", "[", "q", "]"}]}]}], "]"}], ";", " ", RowBox[{"If", "[", RowBox[{ RowBox[{"p3", ">", RowBox[{"\[CapitalPhi]AtL3", "[", "q", "]"}]}], ",", RowBox[{"p3", "=", RowBox[{"\[CapitalPhi]AtL3", "[", "q", "]"}]}]}], "]"}], ";", " ", RowBox[{"Show", "[", " ", RowBox[{ RowBox[{"ContourPlot3D", "[", RowBox[{ RowBox[{"\[CapitalPhi]", "[", RowBox[{"q", ",", "x", ",", "y", ",", "z"}], "]"}], ",", RowBox[{"{", RowBox[{"x", ",", RowBox[{"x", "/.", RowBox[{"posAtL1", "[", "q", "]"}]}], ",", RowBox[{"x", "/.", RowBox[{"posAtL3", "[", "q", "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"y", ",", RowBox[{"ControlActive", "[", RowBox[{"0", ",", RowBox[{"-", ".6"}]}], "]"}], ",", ".6"}], "}"}], ",", RowBox[{"{", RowBox[{"z", ",", RowBox[{"-", ".6"}], ",", ".6"}], "}"}], ",", RowBox[{"Boxed", "->", "False"}], ",", RowBox[{"Contours", "->", RowBox[{"{", "p1", "}"}]}], ",", RowBox[{"ContourStyle", "->", RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"RGBColor", "[", RowBox[{".5", ",", ".5", ",", "1"}], "]"}], ",", RowBox[{"Specularity", "[", RowBox[{"White", ",", "20"}], "]"}]}], "}"}], "}"}]}], ",", RowBox[{"ViewPoint", "->", RowBox[{"{", RowBox[{ RowBox[{"-", "0.011"}], ",", RowBox[{"-", "3.199"}], ",", "1.104"}], "}"}]}], ",", RowBox[{"Axes", "->", "False"}], ",", RowBox[{"Mesh", "->", "False"}], ",", RowBox[{"MaxRecursion", "->", "0"}], ",", RowBox[{"PlotPoints", "->", RowBox[{"ControlActive", "[", RowBox[{"7", ",", RowBox[{"If", "[", RowBox[{"quality", ",", "30", ",", "15"}], "]"}]}], "]"}]}], ",", RowBox[{"BoxRatios", "->", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"x", "/.", RowBox[{"posAtL3", "[", "q", "]"}]}], ")"}], "-", RowBox[{"(", RowBox[{"x", "/.", RowBox[{"posAtL2", "[", "q", "]"}]}], ")"}]}], ",", "1.2", ",", "1.2"}], "}"}]}], ",", RowBox[{"Background", "->", "Black"}], ",", RowBox[{"PerformanceGoal", "->", "\"\\""}]}], "]"}], ",", RowBox[{"ContourPlot3D", "[", RowBox[{ RowBox[{"\[CapitalPhi]", "[", RowBox[{"q", ",", "x", ",", "y", ",", "z"}], "]"}], ",", RowBox[{"{", RowBox[{"x", ",", RowBox[{"x", "/.", RowBox[{"posAtL2", "[", "q", "]"}]}], ",", RowBox[{"x", "/.", RowBox[{"posAtL1", "[", "q", "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"y", ",", RowBox[{"ControlActive", "[", RowBox[{"0", ",", RowBox[{"-", ".6"}]}], "]"}], ",", ".6"}], "}"}], ",", RowBox[{"{", RowBox[{"z", ",", RowBox[{"-", ".6"}], ",", ".6"}], "}"}], ",", RowBox[{"Boxed", "->", "False"}], ",", RowBox[{"Contours", "->", RowBox[{"{", "p2", "}"}]}], ",", RowBox[{"ContourStyle", "->", RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"RGBColor", "[", RowBox[{"1", ",", ".47", ",", "0"}], "]"}], ",", RowBox[{"Specularity", "[", RowBox[{"White", ",", "20"}], "]"}]}], "}"}], "}"}]}], ",", RowBox[{"ViewPoint", "->", RowBox[{"{", RowBox[{ RowBox[{"-", "0.011"}], ",", RowBox[{"-", "3.199"}], ",", "1.104"}], "}"}]}], ",", RowBox[{"Axes", "->", "False"}], ",", RowBox[{"Mesh", "->", "False"}], ",", RowBox[{"MaxRecursion", "->", "0"}], ",", RowBox[{"PlotPoints", "->", RowBox[{"ControlActive", "[", RowBox[{"7", ",", RowBox[{"If", "[", RowBox[{"quality", ",", "30", ",", "15"}], "]"}]}], "]"}]}], ",", RowBox[{"BoxRatios", "->", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"x", "/.", RowBox[{"posAtL3", "[", "q", "]"}]}], ")"}], "-", RowBox[{"(", RowBox[{"x", "/.", RowBox[{"posAtL2", "[", "q", "]"}]}], ")"}]}], ",", "1.2", ",", "1.2"}], "}"}]}], ",", RowBox[{"Background", "->", "Black"}], ",", RowBox[{"PerformanceGoal", "->", "\"\\""}]}], "]"}], ",", RowBox[{"ContourPlot3D", "[", RowBox[{ RowBox[{"\[CapitalPhi]", "[", RowBox[{"q", ",", "x", ",", "y", ",", "z"}], "]"}], ",", RowBox[{"{", RowBox[{"x", ",", RowBox[{"x", "/.", RowBox[{"posAtL3", "[", "q", "]"}]}], ",", RowBox[{"x", "/.", RowBox[{"posAtL2", "[", "q", "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"y", ",", RowBox[{"ControlActive", "[", RowBox[{"0", ",", RowBox[{"-", ".6"}]}], "]"}], ",", ".6"}], "}"}], ",", RowBox[{"{", RowBox[{"z", ",", RowBox[{"-", ".6"}], ",", ".6"}], "}"}], ",", RowBox[{"Boxed", "->", "False"}], ",", RowBox[{"Contours", "->", RowBox[{"{", "p3", "}"}]}], ",", RowBox[{"ContourStyle", "->", RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"RGBColor", "[", RowBox[{"1", ",", "0", ",", "0", ",", RowBox[{"ControlActive", "[", RowBox[{"1", ",", ".3"}], "]"}]}], "]"}], ",", RowBox[{"Specularity", "[", RowBox[{"White", ",", "20"}], "]"}]}], "}"}], "}"}]}], ",", RowBox[{"ViewPoint", "->", RowBox[{"{", RowBox[{ RowBox[{"-", "0.011"}], ",", RowBox[{"-", "3.199"}], ",", "1.104"}], "}"}]}], ",", RowBox[{"Axes", "->", "False"}], ",", RowBox[{"Mesh", "->", "False"}], ",", RowBox[{"MaxRecursion", "->", "0"}], ",", RowBox[{"PlotPoints", "->", RowBox[{"ControlActive", "[", RowBox[{"7", ",", RowBox[{"If", "[", RowBox[{"quality", ",", "30", ",", "15"}], "]"}]}], "]"}]}], ",", RowBox[{"BoxRatios", "->", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"x", "/.", RowBox[{"posAtL3", "[", "q", "]"}]}], ")"}], "-", RowBox[{"(", RowBox[{"x", "/.", RowBox[{"posAtL2", "[", "q", "]"}]}], ")"}]}], ",", "1.2", ",", "1.2"}], "}"}]}], ",", RowBox[{"Background", "->", "Black"}], ",", RowBox[{"PerformanceGoal", "->", "\"\\""}]}], "]"}], ",", RowBox[{"Lighting", "->", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"\"\\"", ",", RowBox[{"GrayLevel", "[", ".75", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"\"\\"", ",", "White", ",", RowBox[{"ImageScaled", "[", RowBox[{"{", RowBox[{"0", ",", "0", ",", "1"}], "}"}], "]"}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{"\"\\"", ",", "White", ",", RowBox[{"{", RowBox[{ RowBox[{"-", RowBox[{"\[Mu]", "[", "q", "]"}]}], ",", "0", ",", "0"}], "}"}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{"\"\\"", ",", "White", ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"\[Mu]", "[", "q", "]"}], "/", "q"}], ",", "0", ",", "0"}], "}"}]}], "}"}]}], "}"}]}], ",", RowBox[{"ImageSize", "->", RowBox[{"{", RowBox[{"450", ",", "300"}], "}"}]}], ",", RowBox[{"PlotRange", "->", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "1.2"}], ",", "1.2"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", ".6"}], ",", ".6"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", ".6"}], ",", ".6"}], "}"}]}], "}"}]}], ",", RowBox[{"SphericalRegion", "->", "True"}], ",", RowBox[{"ViewAngle", "->", RowBox[{"\[Pi]", "/", "15"}]}]}], "]"}]}], ",", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"q", ",", ".4", ",", "\"\\""}], "}"}], ",", ".2", ",", ".9", ",", RowBox[{"Appearance", "->", "\"\\""}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"p1", ",", RowBox[{"-", "4"}], ",", "\"\\""}], "}"}], ",", RowBox[{"-", "5"}], ",", RowBox[{"\[CapitalPhi]AtL1", "[", "q", "]"}], ",", RowBox[{"ControlType", "->", "Slider"}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"p2", ",", RowBox[{"\[CapitalPhi]AtL1", "[", ".4", "]"}], ",", "\"\\""}], "}"}], ",", RowBox[{"-", "5"}], ",", RowBox[{"\[CapitalPhi]AtL1", "[", "q", "]"}], ",", RowBox[{"ControlType", "->", "Slider"}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"p3", ",", RowBox[{"\[CapitalPhi]AtL3", "[", ".4", "]"}], ",", "\"\\""}], "}"}], ",", RowBox[{"-", "6"}], ",", RowBox[{"\[CapitalPhi]AtL3", "[", "q", "]"}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"quality", ",", "False", ",", "\"\\""}], "}"}], ",", RowBox[{"{", RowBox[{"True", ",", "False"}], "}"}]}], "}"}], ",", " ", RowBox[{"TrackedSymbols", ":>", RowBox[{"{", RowBox[{"q", ",", "p1", ",", "p2", ",", "p3", ",", "quality"}], "}"}]}], ",", RowBox[{"AutorunSequencing", "->", RowBox[{"{", RowBox[{"2", ",", "3", ",", "4"}], "}"}]}], ",", RowBox[{"SaveDefinitions", "->", "True"}]}], "]"}]], "Input", CellGroupingRules->{"GroupTogetherGrouping", 10000.}, CellChangeTimes->{ 3.35696210375764*^9, {3.3843875190625*^9, 3.3843876118125*^9}, { 3.384388997390625*^9, 3.384389006375*^9}, 3.38438990796875*^9, { 3.384536027061058*^9, 3.384536037108576*^9}, { 3.3845364849653616`*^9, 3.38453649157516*^9}, { 3.3846939198163166`*^9, 3.3846939320663953`*^9}, 3.760739909464436*^9}], Cell[BoxData[ TagBox[ StyleBox[ DynamicModuleBox[{$CellContext`p1$$ = -4, $CellContext`p2$$ = \ -2.7352387492694854`, $CellContext`p3$$ = -2.4912540035827324`, \ $CellContext`q$$ = 0.4, $CellContext`quality$$ = False, Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{{ Hold[$CellContext`q$$], 0.4, "mass ratio"}, 0.2, 0.9}, {{ Hold[$CellContext`p1$$], -4, "evolve star 1"}, -5, Dynamic[ $CellContext`\[CapitalPhi]AtL1[$CellContext`q$$]]}, {{ Hold[$CellContext`p2$$], -2.7352387492694854`, "evolve star 2"}, -5, Dynamic[ $CellContext`\[CapitalPhi]AtL1[$CellContext`q$$]]}, {{ Hold[$CellContext`p3$$], -2.4912540035827324`, "gravitational potential"}, -6, Dynamic[ $CellContext`\[CapitalPhi]AtL3[$CellContext`q$$]]}, {{ Hold[$CellContext`quality$$], False, "high quality"}, { True, False}}}, Typeset`size$$ = {450., {148., 152.}}, Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ = False, $CellContext`q$46409$$ = 0, $CellContext`p1$46410$$ = 0, $CellContext`p2$46411$$ = 0, $CellContext`p3$46412$$ = 0, $CellContext`quality$46413$$ = False}, DynamicBox[Manipulate`ManipulateBoxes[ 1, StandardForm, "Variables" :> {$CellContext`p1$$ = -4, $CellContext`p2$$ = \ -2.7352387492694854`, $CellContext`p3$$ = -2.4912540035827324`, \ $CellContext`q$$ = 0.4, $CellContext`quality$$ = False}, "ControllerVariables" :> { Hold[$CellContext`q$$, $CellContext`q$46409$$, 0], Hold[$CellContext`p1$$, $CellContext`p1$46410$$, 0], Hold[$CellContext`p2$$, $CellContext`p2$46411$$, 0], Hold[$CellContext`p3$$, $CellContext`p3$46412$$, 0], Hold[$CellContext`quality$$, $CellContext`quality$46413$$, False]}, "OtherVariables" :> { Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$, Typeset`skipInitDone$$}, "Body" :> ( If[$CellContext`p1$$ > \ $CellContext`\[CapitalPhi]AtL1[$CellContext`q$$], $CellContext`p1$$ = \ $CellContext`\[CapitalPhi]AtL1[$CellContext`q$$]]; If[$CellContext`p2$$ > \ $CellContext`\[CapitalPhi]AtL1[$CellContext`q$$], $CellContext`p2$$ = \ $CellContext`\[CapitalPhi]AtL1[$CellContext`q$$]]; If[$CellContext`p3$$ > \ $CellContext`\[CapitalPhi]AtL3[$CellContext`q$$], $CellContext`p3$$ = \ $CellContext`\[CapitalPhi]AtL3[$CellContext`q$$]]; Show[ ContourPlot3D[ $CellContext`\[CapitalPhi][$CellContext`q$$, \ $CellContext`x, $CellContext`y, $CellContext`z], {$CellContext`x, ReplaceAll[$CellContext`x, $CellContext`posAtL1[$CellContext`q$$]], ReplaceAll[$CellContext`x, $CellContext`posAtL3[$CellContext`q$$]]}, {$CellContext`y, ControlActive[0, -0.6], 0.6}, {$CellContext`z, -0.6, 0.6}, Boxed -> False, Contours -> {$CellContext`p1$$}, ContourStyle -> {{ RGBColor[0.5, 0.5, 1], Specularity[White, 20]}}, ViewPoint -> {-0.011, -3.199, 1.104}, Axes -> False, Mesh -> False, MaxRecursion -> 0, PlotPoints -> ControlActive[7, If[$CellContext`quality$$, 30, 15]], BoxRatios -> {ReplaceAll[$CellContext`x, $CellContext`posAtL3[$CellContext`q$$]] - ReplaceAll[$CellContext`x, $CellContext`posAtL2[$CellContext`q$$]], 1.2, 1.2}, Background -> Black, PerformanceGoal -> "Speed"], ContourPlot3D[ $CellContext`\[CapitalPhi][$CellContext`q$$, \ $CellContext`x, $CellContext`y, $CellContext`z], {$CellContext`x, ReplaceAll[$CellContext`x, $CellContext`posAtL2[$CellContext`q$$]], ReplaceAll[$CellContext`x, $CellContext`posAtL1[$CellContext`q$$]]}, {$CellContext`y, ControlActive[0, -0.6], 0.6}, {$CellContext`z, -0.6, 0.6}, Boxed -> False, Contours -> {$CellContext`p2$$}, ContourStyle -> {{ RGBColor[1, 0.47, 0], Specularity[White, 20]}}, ViewPoint -> {-0.011, -3.199, 1.104}, Axes -> False, Mesh -> False, MaxRecursion -> 0, PlotPoints -> ControlActive[7, If[$CellContext`quality$$, 30, 15]], BoxRatios -> {ReplaceAll[$CellContext`x, $CellContext`posAtL3[$CellContext`q$$]] - ReplaceAll[$CellContext`x, $CellContext`posAtL2[$CellContext`q$$]], 1.2, 1.2}, Background -> Black, PerformanceGoal -> "Speed"], ContourPlot3D[ $CellContext`\[CapitalPhi][$CellContext`q$$, \ $CellContext`x, $CellContext`y, $CellContext`z], {$CellContext`x, ReplaceAll[$CellContext`x, $CellContext`posAtL3[$CellContext`q$$]], ReplaceAll[$CellContext`x, $CellContext`posAtL2[$CellContext`q$$]]}, {$CellContext`y, ControlActive[0, -0.6], 0.6}, {$CellContext`z, -0.6, 0.6}, Boxed -> False, Contours -> {$CellContext`p3$$}, ContourStyle -> {{ RGBColor[1, 0, 0, ControlActive[1, 0.3]], Specularity[White, 20]}}, ViewPoint -> {-0.011, -3.199, 1.104}, Axes -> False, Mesh -> False, MaxRecursion -> 0, PlotPoints -> ControlActive[7, If[$CellContext`quality$$, 30, 15]], BoxRatios -> {ReplaceAll[$CellContext`x, $CellContext`posAtL3[$CellContext`q$$]] - ReplaceAll[$CellContext`x, $CellContext`posAtL2[$CellContext`q$$]], 1.2, 1.2}, Background -> Black, PerformanceGoal -> "Speed"], Lighting -> {{"Ambient", GrayLevel[0.75]}, {"Directional", White, ImageScaled[{0, 0, 1}]}, { "Point", White, {-$CellContext`\[Mu][$CellContext`q$$], 0, 0}}, { "Point", White, {$CellContext`\[Mu][$CellContext`q$$]/$CellContext`\ q$$, 0, 0}}}, ImageSize -> {450, 300}, PlotRange -> {{-1.2, 1.2}, {-0.6, 0.6}, {-0.6, 0.6}}, SphericalRegion -> True, ViewAngle -> Pi/15]), "Specifications" :> {{{$CellContext`q$$, 0.4, "mass ratio"}, 0.2, 0.9, Appearance -> "Labeled"}, {{$CellContext`p1$$, -4, "evolve star 1"}, -5, Dynamic[ $CellContext`\[CapitalPhi]AtL1[$CellContext`q$$]], ControlType -> Slider}, {{$CellContext`p2$$, -2.7352387492694854`, "evolve star 2"}, -5, Dynamic[ $CellContext`\[CapitalPhi]AtL1[$CellContext`q$$]], ControlType -> Slider}, {{$CellContext`p3$$, -2.4912540035827324`, "gravitational potential"}, -6, Dynamic[ $CellContext`\[CapitalPhi]AtL3[$CellContext`q$$]]}, \ {{$CellContext`quality$$, False, "high quality"}, {True, False}}}, "Options" :> { TrackedSymbols :> {$CellContext`q$$, $CellContext`p1$$, \ $CellContext`p2$$, $CellContext`p3$$, $CellContext`quality$$}, AutorunSequencing -> {2, 3, 4}}, "DefaultOptions" :> {ControllerLinking -> True}], ImageSizeCache->{495., {245., 251.}}, SingleEvaluation->True], Deinitialization:>None, DynamicModuleValues:>{}, Initialization:>({$CellContext`\[CapitalPhi]AtL1[ Pattern[$CellContext`q, Blank[]]] := $CellContext`\[CapitalPhi][$CellContext`q, ReplaceAll[$CellContext`x, $CellContext`posAtL1[$CellContext`q]], 0, 0], $CellContext`\[CapitalPhi][ Pattern[$CellContext`q, Blank[]], Pattern[$CellContext`x, Blank[]], Pattern[$CellContext`y, Blank[]], Pattern[$CellContext`z, Blank[]]] := -(($CellContext`G $CellContext`M1)/ Sqrt[($CellContext`x + $CellContext`q/( 1 + $CellContext`q))^2 + $CellContext`y^2 + \ $CellContext`z^2]) - ($CellContext`G $CellContext`M2[$CellContext`q])/ Sqrt[($CellContext`x - $CellContext`a + $CellContext`q/( 1 + $CellContext`q))^2 + $CellContext`y^2 + \ $CellContext`z^2] - ((1/ 2) $CellContext`\[CapitalOmega]orb[$CellContext`q]^2) \ (($CellContext`x - $CellContext`\[Mu][$CellContext`q] + \ ($CellContext`q/( 1 + $CellContext`q)) $CellContext`a)^2 + \ $CellContext`y^2), $CellContext`G = 1, $CellContext`M1 = 1, $CellContext`M2[ Pattern[$CellContext`q, Blank[]]] := $CellContext`q $CellContext`M1, $CellContext`a = 1, $CellContext`\[CapitalOmega]orb[ Pattern[$CellContext`q, Blank[]]] := 2 (Pi/$CellContext`Porb[$CellContext`q]), $CellContext`Porb[ Pattern[$CellContext`q, Blank[]]] := Sqrt[(4 Pi^2) ($CellContext`a^3/($CellContext`G ( 1 + $CellContext`q)))], $CellContext`\[Mu][ Pattern[$CellContext`q, Blank[]]] := $CellContext`M1 \ ($CellContext`q/($CellContext`M1 + $CellContext`q $CellContext`M1)), \ $CellContext`posAtL1[ Pattern[$CellContext`q, Blank[]]] := FindRoot[Evaluate[ $CellContext`deriv[$CellContext`q]] == 0, {$CellContext`x, 0.3}], $CellContext`deriv[ Pattern[$CellContext`q, Blank[]]] := D[ $CellContext`\[CapitalPhi][$CellContext`q, $CellContext`x, 0, 0], $CellContext`x], $CellContext`\[CapitalPhi]AtL3[ Pattern[$CellContext`q, Blank[]]] := $CellContext`\[CapitalPhi][$CellContext`q, ReplaceAll[$CellContext`x, $CellContext`posAtL3[$CellContext`q]], 0, 0], $CellContext`posAtL3[ Pattern[$CellContext`q, Blank[]]] := FindRoot[Evaluate[ $CellContext`deriv[$CellContext`q]] == 0, {$CellContext`x, 1.5}], $CellContext`posAtL2[ Pattern[$CellContext`q, Blank[]]] := FindRoot[Evaluate[ $CellContext`deriv[$CellContext`q]] == 0, {$CellContext`x, -1}]}; Typeset`initDone$$ = True), SynchronousInitialization->True, UnsavedVariables:>{Typeset`initDone$$}, UntrackedVariables:>{Typeset`size$$}], "Manipulate", Deployed->True, StripOnInput->False], Manipulate`InterpretManipulate[1]]], "Output", CellGroupingRules->{"GroupTogetherGrouping", 10000.}, CellChangeTimes->{3.760739909464703*^9}, CellID->395119005] }, {3}]]