(* 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[     11454,        331]
NotebookOptionsPosition[      9954,        271]
NotebookOutlinePosition[     10291,        286]
CellTagsIndexPosition[     10248,        283]
WindowFrame->Normal
ContainsDynamic->False*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Earth Sun effective potential", "Section",
 CellChangeTimes->{{3.433601803782052*^9, 3.433601808282585*^9}}],

Cell[CellGroupData[{

Cell["G ", "Subsection",
 CellChangeTimes->{{3.433601929495891*^9, 3.433601929651519*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"G", "=", " ", 
   RowBox[{"6.67", " ", 
    SuperscriptBox["10", 
     RowBox[{"-", "11"}]]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.433601933373866*^9, 3.433601944600246*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Radius of earths orbit", "Subsection",
 CellChangeTimes->{{3.433602187955874*^9, 3.433602192794828*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"RE", "=", 
  RowBox[{"149598000", " ", "1000.0"}]}]], "Input",
 CellChangeTimes->{{3.43360220028607*^9, 3.4336022179156017`*^9}}],

Cell[BoxData["1.49598`*^11"], "Output",
 CellChangeTimes->{{3.433602204523162*^9, 3.433602218307089*^9}, 
   3.433602290578883*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"RE", "=", 
  RowBox[{"1.5", " ", 
   SuperscriptBox["10", "11"]}]}]], "Input",
 CellChangeTimes->{{3.433602223686689*^9, 3.43360222885917*^9}}],

Cell[BoxData["1.5`*^11"], "Output",
 CellChangeTimes->{3.4336022294254503`*^9, 3.433602290645473*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Reduced mass = earth mass", "Subsection",
 CellChangeTimes->{{3.433601847244211*^9, 3.4336018514519863`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"mu", "=", 
   RowBox[{"6", " ", 
    SuperscriptBox["10", "24"]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.433601825486911*^9, 3.433601830948823*^9}, {
  3.43360188226405*^9, 3.433601882955936*^9}, {3.4336022779469633`*^9, 
  3.433602282139187*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Total mass = suns mass ", "Subsection",
 CellChangeTimes->{{3.433601861573824*^9, 3.433601865307724*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"M", " ", "=", " ", 
   RowBox[{"2", " ", 
    SuperscriptBox["10", "30"]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.433601867807991*^9, 3.433601885365899*^9}, {
  3.4336022844347486`*^9, 3.433602287099272*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Angular momentum of earth about sun", "Subsection",
 CellChangeTimes->{{3.433601837094783*^9, 3.433601842971939*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"l", "=", 
  RowBox[{"mu", " ", 
   SuperscriptBox["RE", "2"], " ", 
   FractionBox[
    RowBox[{"2", " ", "Pi"}], 
    RowBox[{"Pi", " ", 
     SuperscriptBox["10", "7"]}]]}]}]], "Input",
 CellChangeTimes->{{3.433601811775181*^9, 3.433601818967599*^9}, 
   3.433601947077173*^9, {3.4336020454722147`*^9, 3.433602140380824*^9}, {
   3.433602233799102*^9, 3.433602236660194*^9}}],

Cell[BoxData["2.6999999999999996`*^40"], "Output",
 CellChangeTimes->{
  3.433602140848755*^9, {3.433602236969019*^9, 3.433602259391671*^9}, 
   3.433602290753587*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"u", "[", "r_", "]"}], "=", 
  RowBox[{
   FractionBox[
    SuperscriptBox["l", "2"], 
    RowBox[{"2", " ", "mu", " ", 
     SuperscriptBox["r", "2"]}]], "-", 
   FractionBox[
    RowBox[{"G", " ", "mu", " ", "M"}], "r"]}]}]], "Input",
 CellChangeTimes->{{3.433601888492647*^9, 3.433601910605774*^9}, {
  3.43360195116537*^9, 3.433601966181568*^9}}],

Cell[BoxData[
 RowBox[{
  FractionBox["6.074999999999999`*^55", 
   SuperscriptBox["r", "2"]], "-", 
  FractionBox["8.004`*^44", "r"]}]], "Output",
 CellChangeTimes->{3.433601967256713*^9, 3.433602108436104*^9, 
  3.4336021437502747`*^9, 3.43360229079132*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"u", "[", 
    RowBox[{"r", " ", "RE"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"r", ",", "0.01", ",", "10"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.433601970777974*^9, 3.433602001733595*^9}, {
  3.4336023151996307`*^9, 3.433602339139729*^9}}],

Cell[BoxData[
 GraphicsBox[{{}, {}, 
   {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwV1nk0Vd3/B3AkZCaEyHBNRaZUSD4n4ilDcVUaJUoyVE8TjzxfhJSQZnk0
iUalopJk7mPMPA/3HvN43XsNmfM7vz/OOuu19tprn8/+vM9eW83jDP0EHw8P
TyP1/P9bcjZoo+fFTkDvxZjw4Epb/rmqmsSZTnCvtZfLTkJb3Yxnwb66DODU
uyX+Ciq2pftd0LFwY8DTw1bRSV6Fts8YCqGdRQzIpzmMfqnIsbUsOm6gepMJ
+XZcn6Dut7ZewZvbOQVMIGQG6uoGU2xjNwpHFUwygeeG9OOPG5Nsx+pdzxit
J4HHKVqlOOWorYXtItbdJoG8Hm4ff/wl1iT3K21JI0H1vTvTT/Udei5Vn3v+
k4QP0bGcpMMfMTorec2FGRLyfy7eeVP4BVvX7bgk69YF7rHl0pMn89Evyqjy
34Au8OY9UjF0vwB5ehXV++O7ICuwha8ktxB1HrF+fSnqglCLxBWc2WK8JHFH
c//abpCXjzfJVi9DYb/gy/lW3eDkk/o7Q7kcH5WeqNU53A01BasneqQrsDjM
9N/ZuG4w7RBPCxitROnJzobEyW54JXfmgJxhDb5wKlm3TLwH5LP99+5OqEGz
dx9CfbV74FWka6n2nxp09wrXszjYA4Ft9/M4hbWY3qIT0ZnXA4bNDzvdjevR
7qKZNbb0wIJT4ODJqHrsl7LjS+dS8xWUnG+01eMae9/QUFovBAaotNIvN2DM
j3fBtKheiGgRPyvwqhF1DuZuEX3aC/HpbRtl2I1Y9LtqbjKrFwj5mfg4kyac
M+AE4jDlEcFTo9lN6J1sfPHUrj5ouXBcwPRLM/KDtQn9ZB+Y3jgU9YDbjE/a
XSbMQ/sgXjCzM0W3BZtkLv4t+qkPPuTGc/r+a8HtUV/902X7wV1DiCj0b0WS
VqqXoN8Pftlau+hPWvFyfstI6F/9MMq88/hndSt+mp09Rf+nH/Re30P/9W2o
6mfhNdXRD6pnb3/y6GjDHCFHTcZUP/Rq1rmlLW9H19QjvSg+APyKQ7uX67dj
HON/HgnEAIxLlCXPXG7HBacCty0pAxDfZGnpLdaBCaxaZY0fAzD3zXH1B4MO
3BDd3SHaNACiggGfdJ070Kdo2SGG4CCQRg7uVbc6sGWjrWuY7yBcuS08oC/c
iRfq9sn5RAwCff7lmStanShx5mQj/dEgVPVHi+tYdaLtq2suGtWDEDyx2jAg
oBMzFSt2lxgPwajRml2hbZ3o9LVd/IP9ECT5fzmWz+3EEZfRXwnHh6D0k+mE
rxAD1WPFHHzuU+MS75YkTRh4a2n3DrG5IfASeHyGGcHA9zd9wqRlh2HYrlno
yEMGlqtEZq8yHIZd7obbbN4xkI/I1lP3GgaBso66Y/UMXFPTcEI7bBhUHSvU
kvoYaO7OfqyXNAw8zOl5h2kGngvVkNpcNwwNt/0jlBWYGCcJdhasYVj3XTzF
ai0T3zw9EL5NaAQKNf2Mh0yZ2JMfN2VvOQL3TXc0F+9j4h7e2Vb3NyNw1WG6
fC6GiWdvrVzp9XME9gSG7Dn7kIkxavoOvuQIHL7p6BaUysSf2zxzL8qNgvmr
R1f35DCxu/bf6SCjUXjPm1SrU8LEP8cSDEMdRsHU1ff3w1ombrryKzn6yiiU
j3LeHOxjIl16sP3mo1G4mpEm/n6MiaeT+WTvZY1CVmTJZOo0E18Wbo56MjYK
y4VECRt+Eovp9PyUFSww9o/as1WIRLLbb/a1Bgt6VfX1RURJlF+W7Jt5kAUd
8u1huitJNLmTk/LtIguMxB0ehsiR6ERr7syNZ4Gh9/iHNAUSr1uLOpUhCwpv
+LPuqpD4ol7relUXC2Z2hMc5q5NY6LmtsH6BBa6fQnSGNUicC79kwjAegyui
MbKf1pIoJ3PLv8dxDMR+v8ZuXRKNU96+GPQegzD1Kb3x9STuMkEmK3wMRIo/
jpEGJPoWk/ITj8fAxK/kfboRidf2zDvPfBuDzPAG0aMbSEzplb2x2DAGa623
pbJNSMy/YFjMxxmDWEO6lMcmEjv47RcFRdgwrGE+8GUziTN3T2wS02LD3U7d
7DFTEmU0Q89Ib2ODnZKvorA5iYafE1+tOswGrwNH7q3YQqKjzecupQA26Iok
fmZR9mmsVlS/zYZNkXN8Xy1IvHpi2EX7HRvmz73m89pKYvIUf6xeKRvm2Eab
5ynnRqqgUQ8b3hFxqgGWJLbLmi9t+sMGH/1NFm2UF+jrlB6LcmDzivXLNYFE
4f++ZqkqcuC1s8SkK2X5nu17n2tzIFRMNeQ8ZS3dOq7mRg7IlO1VDKRscv5o
3CsrDkjtOLnZi7LV99F1uk4cCBtotrGk7LQsqOTdEQ54ZizL4KEcYelxd6cP
B6yOnl+VSa2f9Y/dsb5LHHBa4tu/l/JoprF+WDgHmiS/lfZS36/KUZxXiueA
o9XqHA/Ke3SXlWYlceCZxKbwaqr+a14jd/e85gCxIcVfl3LOs/pjnM8cEGXd
K7tE7Re747t+TCEHrFVbej5S+6sunzKvXc0By53tgp1mJO5ziSktaqfGHbJD
Zqh+5JYd9pif5IDhXxMxvFS/xvltDB7wcuGAtfQH9kaqfmL9grE4F7pynltV
U/2O+7J4z0ebC06G+sWexlT+uH0eAiZcMH2+dbMClZffelUGyQQXAqVn1Aqo
PLk9f1zWeoAL9MfqB4f0SNSPh8UdMVwo3zg490yLRI8K7fLeBC6UanS1cKn8
3heQfBCaygVaY4TwJhqJf4JJw6xcLlS6FlknUfn/dSrsuDaXC9vPizT3rSKR
J/WUUdEiF8rEVdQZslR/SOc/bsLj0HTA8mgl9X8l7VNPuE8bh2Ou507fkCDR
z7qoYvm+cTg9XDuSIECiqPLyDb3Z47C5yfuafj8Tif2spZCScYi3b+EIdDHx
wp3GytUN4/D5hEtUTTsTO1a88HJhjYM57eUxU+q8ePvbNrFQZQKM2z6HeHxn
ol3NNd5nERPg+V7S6WAsEy2XXxKxvz0BI/tPbei7ykRjc0/ZqScTIHC79YJn
KBNXp2xdu/P7BLzJ+Lx+y3kmjgaMO3G5ExCRZT/ptp86H1UOP912dBIUvzBH
01WZWHfaALrNpmBeoaSW9y0DD4g1BdPYv+Gec0DrbHonnrAXUc7YOQfvvxnq
jmS0Yw99sT7OZxHMfDaWvmK1okWfXH/wSR4iaHlFMf/aFsyPjPVvWeAlIuSW
FszeNmGRVmCUefQyIltr8vz5A41Yu8RrddxsOXFHrfvBSskGfLB1WDCzQoAI
j/pv5sN8HXpn/H26IkSI2JfBvy+NVYtZa/1PDukKEwJfhlpyqPtFAX+J7/rv
IkSLbeq0g0kNNq7KlzTzEyNCtQrOnG+pwuSbqwp8eCSI+m9vtPXDf+EPFcEV
xockCdfIpd6dhpU4FJAzZ2QsRagZ93cbLpXj33xq48VzUkRhYlT/uZkyXNJp
/LalQZq4YZpg1cQpxaEw0yibJyuJ5jdsvQeipbiyvfER3UuG4ErNv+vZUoJW
NYUR1payxOHIPO+P3ohV025/lwjKETWvb76M+vATj3630zfoliPSGmnx76eK
0SDkrEli1iqCwwraP7WzGA9ftCESQuWJu1PbB3dnFGGcoI9S0AEFwjF+lH5a
oQiVK7f+GaMpEldbL/01easQR+Lt154cViTe/zDnM5ApxDzuK42igtXENZXI
PbO3C5Ad8VnzU7QSsU6Xb9dX1QLUzJSajjiuTCy9sHE1epuPX2WVq4YN1hBB
5R7l01b5eD3AY8xleg2xbZdDg25RHrpfnGY8rVYhjvyT5TIwm4vldY988h1U
CS3eh3OVtrk4kDXgqJ6tSpjk3JK4H/sDaU+qGiM11YjHYd+uqZA5GLcmlvx0
XY0Qc9HVd7TJwXnf8/XOv9WI4Nzca/Ep3/EB3Whn9iF14p8uV8dyme9IN4vx
1MhTJzZ0Vlcdi81GyTFxLx5dGqHg7Ku9mycbF7acGkm9SSPS5BKCGsO+4Yjr
yP4fCzQi4MJbO9E/Wejm6OdZTllqtr5tYTELa6xY/s2Up0IkfFmUM9ezr3Ap
C97Nu1FF+X/LJtI0F2nE5dplVfGUpdPneOIoK4Ur7ZOlbCog/NrtD434tzL2
wpqFLHwzH53hS3nbPjlBCcrKXJHcQMqszo8Pl+azkK9DrP425XP9ankk5cqP
UotI+URLn/BzykePKDgbLNEIW7XFVG3KdfTEQxaUDf/5bCZPefuO1V47KTdU
hv4SoqyzQfmyJ+U8L6epobksHBdSS31A2fpYvU0a5eN/nqWnUF57h9aWRLl5
Qj37I+X0vKunYynnMDSqKyjLC119cJqyfkNqawvlDJq63lHKz8q0evsoe5nV
5u+mLJP3cmycssjOO3sJylczdWaXKIfTvYcNKf8fGUIYig==
     "]]}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesOrigin->{0, 0},
  PlotRange->{{0.01, 10}, {-2.636376783481485*^33, 1.1555449480732167`*^33}},
  PlotRangeClipping->True,
  PlotRangePadding->{
    Scaled[0.02], 
    Scaled[0.02]}]], "Output",
 CellChangeTimes->{
  3.433602002238459*^9, 3.4336021122051086`*^9, 3.433602146835227*^9, 
   3.433602290856399*^9, {3.433602328862468*^9, 3.4336023396142883`*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Frequency of small oscillations", "Subsection",
 CellChangeTimes->{{3.4336029923109293`*^9, 3.4336029968234167`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Sqrt", "[", 
  RowBox[{"G", " ", 
   FractionBox["M", 
    SuperscriptBox["RE", "3"]]}], "]"}]], "Input",
 CellChangeTimes->{{3.433603030403572*^9, 3.4336030420724163`*^9}}],

Cell[BoxData["1.988112821897337`*^-7"], "Output",
 CellChangeTimes->{3.433603042479271*^9}]
}, Open  ]]
}, Open  ]]
}, Open  ]]
},
WindowSize->{640, 750},
WindowMargins->{{36, Automatic}, {Automatic, 84}},
FrontEndVersion->"6.0 for Mac OS X x86 (32-bit) (March 13, 2008)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[590, 23, 114, 1, 67, "Section"],
Cell[CellGroupData[{
Cell[729, 28, 90, 1, 34, "Subsection"],
Cell[822, 31, 214, 6, 30, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[1073, 42, 110, 1, 34, "Subsection"],
Cell[CellGroupData[{
Cell[1208, 47, 153, 3, 27, "Input"],
Cell[1364, 52, 131, 2, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[1532, 59, 167, 4, 30, "Input"],
Cell[1702, 65, 101, 1, 30, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[1852, 72, 115, 1, 34, "Subsection"],
Cell[1970, 75, 285, 7, 30, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[2292, 87, 111, 1, 34, "Subsection"],
Cell[2406, 90, 246, 6, 30, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[2689, 101, 123, 1, 34, "Subsection"],
Cell[CellGroupData[{
Cell[2837, 106, 401, 10, 49, "Input"],
Cell[3241, 118, 167, 3, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[3445, 126, 384, 11, 50, "Input"],
Cell[3832, 139, 259, 6, 49, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[4128, 150, 315, 8, 27, "Input"],
Cell[4446, 160, 4980, 87, 210, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[9475, 253, 123, 1, 34, "Subsection"],
Cell[CellGroupData[{
Cell[9623, 258, 197, 5, 46, "Input"],
Cell[9823, 265, 91, 1, 30, "Output"]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)