(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 9.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[     26570,        820]
NotebookOptionsPosition[     24859,        756]
NotebookOutlinePosition[     25203,        771]
CellTagsIndexPosition[     25160,        768]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[TextData[StyleBox["C2H6 Molecule after Kimball, non overlapping, \
parameterized with MP2/6-31G(d)", "Subsubtitle",
 FontWeight->"Bold"]], "Title",
 CellChangeTimes->{
  3.5792362808724213`*^9, {3.5792373471186934`*^9, 3.5792373542323065`*^9}, {
   3.5792374368032513`*^9, 3.5792374830417323`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"Clear", "[", 
   RowBox[{"k1", ",", "k2", ",", "k3", ",", "s1", ",", "s2", ",", "s3"}], 
   "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{" ", 
  RowBox[{
   RowBox[{"n", "=", "2."}], ";", " ", 
   RowBox[{"cw", "=", 
    RowBox[{"1", "/", "3"}]}], ";", " ", 
   RowBox[{"sw", "=", 
    RowBox[{
     RowBox[{"Sqrt", "[", "2", "]"}], "/", "3"}]}], ";"}]}]}], "Input",
 CellChangeTimes->{3.5792363517745457`*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"c", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"k1", "->", "1.02669238"}], ",", 
     RowBox[{"k2", "->", "1.36467050"}], ",", 
     RowBox[{"k3", "->", "1.20088548"}], ",", 
     RowBox[{"s1", "->", "0.30102681"}], ",", 
     RowBox[{"s2", "->", "0.34863843"}], ",", 
     RowBox[{"s3", "->", "0.34858036"}]}], "}"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.579236366797372*^9, 3.5792363668285723`*^9}}],

Cell["Terms from ChemEdu/Alkan.sys", "Text",
 CellChangeTimes->{3.5792364030830355`*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"T", " ", "=", " ", 
    RowBox[{
     RowBox[{
      RowBox[{"2.25", "*", "n", "*", 
       RowBox[{"(", 
        RowBox[{"k1", "/", 
         RowBox[{"P", "^", "2"}]}], ")"}]}], " ", "+", " ", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"n", " ", "-", " ", "1."}], ")"}], "*", "2.25", "*", 
       RowBox[{"(", 
        RowBox[{"k2", "/", 
         RowBox[{"Q", "^", "2"}]}], ")"}]}], " ", "+", " ", "\n", "     ", 
      RowBox[{"4.5", "*", 
       RowBox[{"(", 
        RowBox[{"n", " ", "+", " ", "1."}], ")"}], "*", 
       RowBox[{"(", 
        RowBox[{"k3", "/", 
         RowBox[{"R", "^", "2"}]}], ")"}]}]}], " ", "/.", " ", "c"}]}], ";"}],
   " "}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"p", " ", "=", " ", "X"}], ";", " ", 
   RowBox[{"p2", " ", "=", " ", 
    RowBox[{"p", "*", "p"}]}], ";", " ", 
   RowBox[{"z", " ", "=", " ", "6."}], ";", " ", 
   RowBox[{"ad", " ", "=", " ", 
    RowBox[{"8.", "/", "3."}]}], ";", " ", 
   RowBox[{"zd", " ", "=", " ", 
    RowBox[{"2.", "/", "3."}]}], ";"}], " "}], "\n", 
 RowBox[{"  ", 
  RowBox[{
   RowBox[{"vd", " ", "=", " ", 
    RowBox[{"4.", "/", "3."}]}], ";", " ", 
   RowBox[{"q", " ", "=", " ", 
    RowBox[{"P", " ", "+", " ", "Q"}]}], ";", " ", 
   RowBox[{"r", " ", "=", " ", 
    RowBox[{"P", " ", "+", " ", "R"}]}], ";", " ", 
   RowBox[{"q2", " ", "=", " ", 
    RowBox[{"q", "*", "q"}]}], ";", " ", 
   RowBox[{"r2", " ", "=", " ", 
    RowBox[{"r", "*", "r"}]}], ";"}], " "}], "\n", 
 RowBox[{"  ", 
  RowBox[{
   RowBox[{"r2p2", " ", "=", " ", 
    RowBox[{"r2", "*", "p2"}]}], ";", " ", 
   RowBox[{"r2p", " ", "=", " ", 
    RowBox[{"r2", "*", "p"}]}], ";", " ", 
   RowBox[{"qr", " ", "=", " ", 
    RowBox[{"q", "*", "r"}]}], ";", " ", 
   RowBox[{"pr", " ", "=", " ", 
    RowBox[{"p", "*", "r"}]}], ";"}], " "}], "\n", 
 RowBox[{"  ", 
  RowBox[{
   RowBox[{"qrp", " ", "=", " ", 
    RowBox[{"qr", "*", "p"}]}], ";"}], " "}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"Vee", " ", "=", " ", 
    RowBox[{
     RowBox[{
      RowBox[{"3.", "*", "n", "*", 
       RowBox[{"(", 
        RowBox[{"s1", "/", "P"}], ")"}]}], " ", "+", " ", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"n", " ", "-", " ", "1"}], ")"}], "*", "3.", "*", 
       RowBox[{"(", 
        RowBox[{"s2", "/", "Q"}], ")"}]}], " ", "+", " ", "\n", "     ", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"n", " ", "+", " ", "1"}], ")"}], "*", "6.", "*", 
       RowBox[{"(", 
        RowBox[{"s3", "/", "R"}], ")"}]}], " ", "+", " ", 
      RowBox[{"8.", "/", "q"}], " ", "+", " ", "\n", "     ", 
      RowBox[{"24.", "/", 
       RowBox[{"Sqrt", "[", 
        RowBox[{
         RowBox[{"zd", "*", "qr"}], " ", "+", " ", "q2", " ", "+", " ", 
         "r2"}], "]"}]}], " ", "+", " ", 
      RowBox[{"24.", "/", "r"}], " ", "+", " ", 
      RowBox[{"4.", "/", 
       RowBox[{"(", 
        RowBox[{"2", "*", "q"}], ")"}]}], " ", "+", " ", "\n", "     ", 
      RowBox[{"24.", "/", 
       RowBox[{"Sqrt", "[", 
        RowBox[{
         RowBox[{"vd", "*", "qr"}], " ", "+", " ", 
         RowBox[{"4.", "*", "q2"}], " ", "+", " ", "r2"}], "]"}]}], " ", "+", 
      " ", 
      RowBox[{"24.", "/", 
       RowBox[{"Sqrt", "[", 
        RowBox[{"ad", "*", "r2"}], "]"}]}], " ", "+", " ", "\n", "     ", 
      RowBox[{"12.", "/", 
       RowBox[{"Sqrt", "[", 
        RowBox[{
         RowBox[{"ad", "*", "qr"}], " ", "+", " ", 
         RowBox[{"4.", "*", "q2"}], " ", "+", " ", 
         RowBox[{"4.", "*", "r2"}]}], "]"}]}], " ", "+", " ", "\n", "     ", 
      RowBox[{"24.", "/", 
       RowBox[{"Sqrt", "[", 
        RowBox[{
         RowBox[{"ad", "*", "qr"}], " ", "+", " ", 
         RowBox[{"4.", "*", "q2"}], " ", "+", " ", 
         RowBox[{"vd", "*", "r2"}]}], "]"}]}]}], " ", "/.", " ", "c"}]}], 
   ";"}], " "}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"Vne", " ", "=", " ", 
    RowBox[{
     RowBox[{
      RowBox[{"-", "2."}], "*", 
      RowBox[{"(", 
       RowBox[{"n", " ", "+", " ", "1."}], ")"}], "*", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"3.", " ", "-", " ", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"p", " ", "-", " ", "1."}], ")"}], "*", 
             RowBox[{"(", 
              RowBox[{"1.", " ", "+", " ", 
               RowBox[{"P", "/", "R"}]}], ")"}]}], ")"}], "^", "2"}]}], ")"}],
         "/", "R"}], ")"}]}], " ", "-", " ", "\n", "    ", 
     RowBox[{"3.", "*", "n", "*", 
      RowBox[{"(", 
       RowBox[{"z", "/", "P"}], ")"}]}], " ", "-", " ", 
     RowBox[{"4.", "*", 
      RowBox[{"(", 
       RowBox[{"z", "/", "q"}], ")"}]}], " ", "-", " ", "\n", "    ", 
     RowBox[{"12.", "/", 
      RowBox[{"Sqrt", "[", 
       RowBox[{
        RowBox[{"zd", "*", "qrp"}], " ", "+", " ", "r2p2", " ", "+", " ", 
        "q2"}], "]"}]}], " ", "-", " ", 
     RowBox[{"12.", "*", 
      RowBox[{"(", 
       RowBox[{"z", "/", 
        RowBox[{"Sqrt", "[", "r2", "]"}]}], ")"}]}], " ", "-", " ", "\n", 
     "    ", 
     RowBox[{"12.", "/", "pr"}], " ", "-", " ", 
     RowBox[{"4.", "*", 
      RowBox[{"(", 
       RowBox[{"z", "/", 
        RowBox[{"(", 
         RowBox[{"2", "*", "q"}], ")"}]}], ")"}]}], " ", "-", " ", "\n", 
     "    ", 
     RowBox[{"12.", "*", 
      RowBox[{"(", 
       RowBox[{"z", "/", 
        RowBox[{"Sqrt", "[", 
         RowBox[{
          RowBox[{"vd", "*", "qr"}], " ", "+", " ", 
          RowBox[{"4.", "*", "q2"}], " ", "+", " ", "r2"}], "]"}]}], ")"}]}], 
     " ", "-", " ", "\n", "    ", 
     RowBox[{"12.", "/", 
      RowBox[{"Sqrt", "[", 
       RowBox[{
        RowBox[{"vd", "*", "qrp"}], " ", "+", " ", "r2p2", " ", "+", " ", 
        RowBox[{"4.", "*", "q2"}]}], "]"}]}], " ", "-", " ", "\n", "    ", 
     RowBox[{"24.", "/", 
      RowBox[{"Sqrt", "[", 
       RowBox[{
        RowBox[{"zd", "*", "r2p"}], " ", "+", " ", "r2p2", " ", "+", " ", 
        "r2"}], "]"}]}], " ", "-", " ", "\n", "    ", 
     RowBox[{"12.", "/", 
      RowBox[{"Sqrt", "[", 
       RowBox[{
        RowBox[{"vd", "*", "qr"}], " ", "+", " ", 
        RowBox[{"vd", "*", "qrp"}], " ", "+", " ", 
        RowBox[{"2.", "*", "r2p"}], " ", "+", " ", "r2p2", " ", "+", " ", 
        RowBox[{"4.", "*", "q2"}], " ", "+", " ", "\n", "       ", "r2"}], 
       "]"}]}], " ", "-", " ", 
     RowBox[{"24.", "/", 
      RowBox[{"Sqrt", "[", 
       RowBox[{
        RowBox[{"vd", "*", "qr"}], " ", "+", " ", 
        RowBox[{"vd", "*", "qrp"}], " ", "-", " ", 
        RowBox[{"zd", "*", "r2p"}], " ", "+", " ", "r2p2", " ", "+", " ", 
        "\n", "       ", 
        RowBox[{"4.", "*", "q2"}], " ", "+", " ", "r2"}], "]"}]}]}]}], ";"}], 
  " "}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"Vnn", " ", "=", " ", 
    RowBox[{
     RowBox[{"6.", "*", 
      RowBox[{"(", 
       RowBox[{"z", "/", 
        RowBox[{"(", 
         RowBox[{"r", "*", "p"}], ")"}]}], ")"}]}], " ", "+", " ", 
     RowBox[{"z", "*", 
      RowBox[{"(", 
       RowBox[{"z", "/", 
        RowBox[{"(", 
         RowBox[{"2", "*", "q"}], ")"}]}], ")"}]}], " ", "+", " ", "\n", 
     "    ", 
     RowBox[{"6.", "*", 
      RowBox[{"(", 
       RowBox[{"z", "/", 
        RowBox[{"Sqrt", "[", 
         RowBox[{
          RowBox[{"vd", "*", "qrp"}], " ", "+", " ", "r2p2", " ", "+", " ", 
          RowBox[{"4.", "*", "q2"}]}], "]"}]}], ")"}]}], " ", "+", " ", "\n", 
     "    ", 
     RowBox[{"6.", "/", 
      RowBox[{"Sqrt", "[", 
       RowBox[{"ad", "*", "r2p2"}], "]"}]}], " ", "+", " ", "\n", "    ", 
     RowBox[{"3.", "/", 
      RowBox[{"Sqrt", "[", 
       RowBox[{
        RowBox[{"ad", "*", "qrp"}], " ", "+", " ", 
        RowBox[{"4.", "*", "r2p2"}], " ", "+", " ", 
        RowBox[{"4.", "*", "q2"}]}], "]"}]}], " ", "+", " ", "\n", "    ", 
     RowBox[{"6.", "/", 
      RowBox[{"Sqrt", "[", 
       RowBox[{
        RowBox[{"ad", "*", "qrp"}], " ", "+", " ", 
        RowBox[{"vd", "*", "r2p2"}], " ", "+", " ", 
        RowBox[{"4.", "*", "q2"}]}], "]"}]}]}]}], ";"}], " "}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"func", " ", "=", " ", 
    RowBox[{
    "T", " ", "+", " ", "Vne", " ", "+", " ", "Vee", " ", "+", " ", "Vnn"}]}],
    ";"}], " "}]}], "Input",
 CellChangeTimes->{{3.579236416436659*^9, 3.579236487775584*^9}, 
   3.57923667291671*^9, 3.5792367206527934`*^9},
 FormatType->"InputForm"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"t", " ", "=", " ", 
  RowBox[{"FindMinimum", "[", 
   RowBox[{"func", ",", " ", 
    RowBox[{"{", 
     RowBox[{"P", ",", " ", "0.25", ",", " ", "0.27"}], "}"}], ",", " ", 
    RowBox[{"{", 
     RowBox[{"Q", ",", " ", "1.1", ",", " ", "1.2"}], "}"}], ",", " ", "\n", 
    "   ", 
    RowBox[{"{", 
     RowBox[{"R", ",", " ", "1.2", ",", " ", "1.3"}], "}"}], ",", " ", 
    RowBox[{"{", 
     RowBox[{"X", ",", " ", "1.3", ",", " ", "1.4"}], "}"}], ",", " ", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"Method", " ", "->", " ", "Automatic"}], ",", " ", "\n", "    ", 
      RowBox[{"MaxIterations", " ", "->", " ", "500"}]}], "}"}]}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.579236416436659*^9, 3.579236487775584*^9}, 
   3.57923667291671*^9, 3.5792367206527934`*^9},
 FormatType->"InputForm"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "79.49412996948107`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"P", "\[Rule]", "0.26224005195523936`"}], ",", 
     RowBox[{"Q", "\[Rule]", "1.2000710223961832`"}], ",", 
     RowBox[{"R", "\[Rule]", "1.261969408118396`"}], ",", 
     RowBox[{"X", "\[Rule]", "1.356989234003253`"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{3.579236756704457*^9, 3.579237297276606*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"u", " ", "=", " ", 
   RowBox[{"t", "[", 
    RowBox[{"[", "2", "]"}], "]"}]}], ";"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   StyleBox[
    RowBox[{"Total", " ", "Energy"}], "Subsubtitle"], 
   StyleBox[" ", "Text"], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"Etot", "=", 
   RowBox[{"t", "[", 
    RowBox[{"[", "1", "]"}], "]"}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{"-", "79.49412996948107`"}]], "Output",
 CellChangeTimes->{3.5792372972922063`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"(*", 
    RowBox[{"Kinetic", " ", "energy"}], "*)"}], "Subsubtitle"], 
  "\[IndentingNewLine]", 
  RowBox[{"T", " ", "/.", "u"}]}]], "Input"],

Cell[BoxData["79.49409996365905`"], "Output",
 CellChangeTimes->{3.579237297307806*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   StyleBox[
    RowBox[{"NUCLEUS", "-", 
     RowBox[{"ELECTRON", " ", "POTENTIAL", " ", "ENERGY"}]}], "Subsubtitle"], 
   StyleBox[
    RowBox[{
     StyleBox[" ", "Subsubtitle"], " "}]], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"Vne", " ", "/.", "u"}]}]], "Input"],

Cell[BoxData[
 RowBox[{"-", "270.4620043570094`"}]], "Output",
 CellChangeTimes->{3.579237297307806*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   StyleBox[
    RowBox[{"ELECTRON", "-", 
     RowBox[{"ELECTRON", " ", "POTENTIAL", " ", "ENERGY"}]}], "Subsubtitle"], 
   StyleBox[" ", "Subsubtitle"], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Vee", " ", "/.", "u"}], " ", "/.", "c"}]}]], "Input"],

Cell[BoxData["69.42993493962692`"], "Output",
 CellChangeTimes->{3.579237297323406*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   StyleBox[
    RowBox[{"NUCLEUS", "-", 
     RowBox[{"NUCLEUS", " ", "POTENTIAL", " ", "ENERGY"}]}], "Subsubtitle"], 
   StyleBox[" ", "Text"], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"Vnn", " ", "/.", "u"}]}]], "Input"],

Cell[BoxData["42.04383948424231`"], "Output",
 CellChangeTimes->{3.5792372973390064`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  StyleBox[" ", "Text"], 
  StyleBox[
   RowBox[{"(*", 
    RowBox[{"C", "-", 
     RowBox[{"H", " ", "distance", "  ", "\[CapitalARing]"}]}], "*)"}], 
   "Subsubtitle"], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"ch", " ", "=", " ", 
     RowBox[{
      RowBox[{"X", "*", 
       RowBox[{"(", 
        RowBox[{"P", " ", "+", " ", "R"}], ")"}]}], " ", "/.", "u"}]}], ";"}],
    "\[IndentingNewLine]", 
   RowBox[{"dch", " ", "=", " ", 
    RowBox[{
    "ch", "*", "0.529177", 
     "                                "}]}]}]}]], "Input"],

Cell[BoxData["1.0945157482873067`"], "Output",
 CellChangeTimes->{3.5792372973390064`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  StyleBox[" ", "Text"], 
  RowBox[{"(*", 
   StyleBox[
    RowBox[{"C", "-", 
     RowBox[{"C", " ", "distance", "  ", "\[CapitalARing]"}]}], "Subsubtitle"], 
   StyleBox[" ", "Subsubtitle"], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"cc", " ", "=", " ", 
     RowBox[{
      RowBox[{"2.", "*", 
       RowBox[{"(", 
        RowBox[{"P", " ", "+", " ", "Q"}], ")"}]}], " ", "/.", "u"}]}], ";"}],
    "\[IndentingNewLine]", 
   RowBox[{"dcc", "=", 
    RowBox[{
    "cc", "*", "0.529177", 
     "                                              "}]}]}]}]], "Input"],

Cell[BoxData["1.5476427747841255`"], "Output",
 CellChangeTimes->{3.5792372973546066`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   StyleBox[
    RowBox[{"Virial", " ", "Theorem"}], "Subsubtitle"], 
   StyleBox[" ", "Subsubtitle"], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"-", 
      RowBox[{"(", 
       RowBox[{"Vne", " ", "+", " ", "Vee", " ", "+", " ", "Vnn"}], ")"}]}], 
     "/", "T"}], " ", "/.", "u"}], " ", "/.", "c", 
   "             "}]}]], "Input"],

Cell[BoxData["2.0000003774597364`"], "Output",
 CellChangeTimes->{3.5792372973702064`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   StyleBox[
    RowBox[{
    "Atomization", " ", "Energy", " ", "at", " ", "0", " ", "K", " ", "in", 
     " ", 
     RowBox[{"kcal", "/", "mol"}]}], "Subsubtitle"], 
   StyleBox[" ", "Subsubtitle"], "*)"}], 
  StyleBox["\[IndentingNewLine]", "Text"], 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", 
      RowBox[{"t", "[", 
       RowBox[{"[", "1", "]"}], "]"}]}], " ", "-", " ", 
     RowBox[{"n", "*", "37.78899"}], " ", "-", 
     RowBox[{
      RowBox[{"(", " ", 
       RowBox[{
        RowBox[{"2", "*", "n"}], " ", "+", " ", "2"}], ")"}], "*", 
      "0.50177"}]}], ")"}], "*", "627.5096", " "}]}]], "Input"],

Cell[BoxData["568.2287489370774`"], "Output",
 CellChangeTimes->{3.579237297385806*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"plot1", " ", "=", " ", 
    RowBox[{
     RowBox[{"Graphics", "[", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"-", "q"}], ",", "0"}], "}"}], ",", "P"}], "]"}], ",", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"-", "q"}], ",", "0"}], "}"}], ",", "0.1"}], "]"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"q", ",", "0"}], "}"}], ",", "P"}], "]"}], ",", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"q", ",", "0"}], "}"}], ",", "0.1"}], "]"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"0", ",", "0"}], "}"}], ",", "Q"}], "]"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{
             RowBox[{"-", "q"}], "-", 
             RowBox[{"r", "*", "cw"}]}], ",", 
            RowBox[{
             RowBox[{"-", "r"}], "*", "2", "*", "sw"}]}], "}"}], ",", "R"}], 
         "]"}], ",", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{
             RowBox[{"-", "q"}], "-", " ", 
             RowBox[{"pr", "*", "cw"}]}], ",", 
            RowBox[{
             RowBox[{"-", "pr"}], "*", "2", "*", "sw"}]}], "}"}], ",", 
          "0.1"}], "]"}], ",", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"q", " ", "+", 
             RowBox[{"r", "*", "cw"}]}], ",", 
            RowBox[{"r", "*", "2", "*", "sw"}]}], "}"}], ",", "R"}], "]"}], 
        ",", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"q", "+", " ", 
             RowBox[{"pr", "*", "cw"}]}], ",", 
            RowBox[{"pr", "*", "2", "*", "sw"}]}], "}"}], ",", "0.1"}], "]"}],
         ",", "\[IndentingNewLine]", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Dashing", "[", 
           RowBox[{"{", 
            RowBox[{"0.02", ",", "0.02"}], "}"}], "]"}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"Circle", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"q", "+", 
               RowBox[{"r", "*", "cw"}]}], ",", 
              RowBox[{
               RowBox[{"-", "r"}], "*", "sw"}]}], "}"}], ",", "R"}], "]"}], 
          ",", 
          RowBox[{"Disk", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"q", " ", "+", " ", 
               RowBox[{"pr", "*", "cw"}]}], ",", 
              RowBox[{
               RowBox[{"-", "pr"}], "*", "sw"}]}], "}"}], ",", "0.1"}], "]"}],
           ",", 
          RowBox[{"Circle", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{
               RowBox[{"-", "q"}], "-", 
               RowBox[{"r", "*", "cw"}]}], ",", 
              RowBox[{"r", "*", "sw"}]}], "}"}], ",", "R"}], "]"}], ",", 
          RowBox[{"Disk", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{
               RowBox[{"-", "q"}], "-", 
               RowBox[{"pr", "*", "cw"}]}], ",", 
              RowBox[{"pr", "*", "sw"}]}], "}"}], ",", "0.1"}], "]"}]}], 
         "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Thickness", "[", "0.01", "]"}], ",", 
          RowBox[{"Line", "[", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{"-", "q"}], ",", "0"}], "}"}], ",", 
             RowBox[{"{", 
              RowBox[{"q", ",", "0"}], "}"}]}], "}"}], "]"}], ",", 
          RowBox[{"Line", "[", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{"q", ",", "0"}], "}"}], ",", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"q", " ", "+", " ", 
                RowBox[{"pr", "*", "cw"}]}], ",", 
               RowBox[{"pr", "*", "2", "*", "sw"}]}], "}"}]}], "}"}], "]"}], 
          ",", 
          RowBox[{"Line", "[", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{"-", "q"}], ",", "0"}], "}"}], ",", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{
                RowBox[{"-", "q"}], "-", " ", 
                RowBox[{"pr", "/", "3."}]}], ",", 
               RowBox[{
                RowBox[{"-", "pr"}], "*", "2", "*", "sw"}]}], "}"}]}], "}"}], 
           "]"}], ",", "\[IndentingNewLine]", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"Dashing", "[", 
             RowBox[{"{", 
              RowBox[{"0.01", ",", "0.01"}], "}"}], "]"}], ",", 
            "\[IndentingNewLine]", 
            RowBox[{"Line", "[", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"{", 
                RowBox[{"q", ",", "0"}], "}"}], ",", 
               RowBox[{"{", 
                RowBox[{
                 RowBox[{"q", " ", "+", " ", 
                  RowBox[{"pr", "*", "cw"}]}], ",", 
                 RowBox[{
                  RowBox[{"-", "pr"}], "*", "sw"}]}], "}"}]}], "}"}], "]"}], 
            ",", "\[IndentingNewLine]", 
            RowBox[{"Line", "[", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"{", 
                RowBox[{
                 RowBox[{"-", "q"}], ",", "0"}], "}"}], ",", 
               RowBox[{"{", 
                RowBox[{
                 RowBox[{
                  RowBox[{"-", "q"}], " ", "-", " ", 
                  RowBox[{"pr", "*", "cw"}]}], ",", 
                 RowBox[{"pr", "*", "sw"}]}], "}"}]}], "}"}], "]"}]}], 
           "}"}]}], "}"}]}], "}"}], "]"}], " ", "/.", "u"}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{"Show", "[", 
  RowBox[{"plot1", ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"AspectRatio", " ", "->", " ", "Automatic"}], ",", 
     RowBox[{"PlotRange", " ", "->", " ", "All"}], ",", 
     RowBox[{"Axes", " ", "->", " ", "True"}], ",", 
     RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", 
     RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{
  StyleBox[
   RowBox[{"(*", 
    RowBox[{"Projection", " ", "on", " ", "xy", " ", "plane"}], " ", "*)"}], 
   "Subsubtitle"]}]}], "Input",
 CellChangeTimes->{3.579237288540591*^9}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{-1.4623110743514225`, 0}, 0.26224005195523936`], 
   DiskBox[{-1.4623110743514225`, 0}, 0.1], 
   CircleBox[{1.4623110743514225`, 0}, 0.26224005195523936`], 
   DiskBox[{1.4623110743514225`, 0}, 0.1], 
   CircleBox[{0, 0}, 1.2000710223961832`], 
   CircleBox[{-1.9703808943759675`, -1.4370384602223387`}, 1.261969408118396],
    DiskBox[{-2.1517563502467008`, -1.9500457193703256`}, 0.1], 
   CircleBox[{1.9703808943759675`, 1.4370384602223387`}, 1.261969408118396], 
   DiskBox[{2.1517563502467008`, 1.9500457193703256`}, 0.1], 
   {Dashing[{0.02, 0.02}], 
    CircleBox[{1.9703808943759675`, -0.7185192301111694}, 1.261969408118396], 
    DiskBox[{2.1517563502467008`, -0.9750228596851628}, 0.1], 
    CircleBox[{-1.9703808943759675`, 0.7185192301111694}, 1.261969408118396], 
    DiskBox[{-2.1517563502467008`, 0.9750228596851628}, 0.1]}, 
   {Thickness[0.01], 
    LineBox[{{-1.4623110743514225`, 0}, {1.4623110743514225`, 0}}], 
    LineBox[{{1.4623110743514225`, 0}, {2.1517563502467008`, 
      1.9500457193703256`}}], 
    LineBox[{{-1.4623110743514225`, 
       0}, {-2.1517563502467008`, -1.9500457193703256`}}], 
    {Dashing[{0.01, 0.01}], 
     LineBox[{{1.4623110743514225`, 0}, {
       2.1517563502467008`, -0.9750228596851628}}], 
     LineBox[{{-1.4623110743514225`, 0}, {-2.1517563502467008`, 
       0.9750228596851628}}]}}},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->All]], "Output",
 CellChangeTimes->{3.5792372974014063`*^9}]
}, Open  ]]
}, Open  ]]
},
WindowSize->{741, 817},
WindowMargins->{{24, Automatic}, {Automatic, 30}},
FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (January 25, 2013)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 303, 5, 63, "Title"],
Cell[885, 29, 448, 13, 52, "Input"],
Cell[1336, 44, 446, 11, 52, InheritFromParent],
Cell[1785, 57, 88, 1, 30, "Text"],
Cell[1876, 60, 8378, 234, 491, "Input"],
Cell[CellGroupData[{
Cell[10279, 298, 827, 20, 66, "Input"],
Cell[11109, 320, 444, 11, 31, "Output"]
}, Open  ]],
Cell[11568, 334, 131, 4, 31, "Input"],
Cell[CellGroupData[{
Cell[11724, 342, 262, 8, 52, "Input"],
Cell[11989, 352, 106, 2, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[12132, 359, 189, 6, 52, "Input"],
Cell[12324, 367, 87, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[12448, 373, 309, 9, 52, "Input"],
Cell[12760, 384, 104, 2, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[12901, 391, 304, 8, 52, "Input"],
Cell[13208, 401, 87, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[13332, 407, 265, 7, 52, "Input"],
Cell[13600, 416, 89, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[13726, 422, 576, 19, 76, "Input"],
Cell[14305, 443, 90, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[14432, 449, 611, 19, 76, "Input"],
Cell[15046, 470, 90, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[15173, 476, 414, 13, 52, "Input"],
Cell[15590, 491, 90, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[15717, 497, 671, 21, 52, "Input"],
Cell[16391, 520, 87, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[16515, 526, 6776, 193, 372, "Input"],
Cell[23294, 721, 1537, 31, 320, "Output"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)