(* 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[ 23039, 768] NotebookOptionsPosition[ 22792, 755] NotebookOutlinePosition[ 23137, 770] CellTagsIndexPosition[ 23094, 767] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[{ RowBox[{ RowBox[{"sq3", "=", RowBox[{"Sqrt", "[", "3", "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"nc", "=", "9"}], ";"}], "\n", RowBox[{ RowBox[{"a", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", "nc", "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"b", "=", "a"}], ";", " ", RowBox[{"c", "=", "a"}], ";", " ", RowBox[{"g", "=", "a"}], ";", " ", RowBox[{"sig", "=", "a"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"k1", "=", "1.0"}], ";", RowBox[{"k2", "=", "1.0"}], ";", RowBox[{"k3", "=", "1.0"}], ";", RowBox[{"k4", "=", "1.0"}], ";", RowBox[{"k5", "=", "1.0"}], ";"}], "\n", RowBox[{ RowBox[{"k6", "=", "1.0"}], ";", RowBox[{"k7", "=", "1.0"}], ";"}], "\n", RowBox[{ RowBox[{"sig1", "=", "0.31403448"}], ";", RowBox[{"sig2", "=", "0.30821367"}], ";", RowBox[{"sig3", "=", "0.30820426"}], ";", RowBox[{"sig4", "=", "0.3"}], ";", RowBox[{"sig5", "=", "0.3"}], ";"}], "\n", RowBox[{ RowBox[{"sig6", "=", "0.3"}], ";", RowBox[{"sig7", "=", "0.3"}], ";"}], "\n", RowBox[{ RowBox[{"For", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", "10"}], ",", RowBox[{"i", "++"}], ",", " ", RowBox[{ RowBox[{"g", "[", RowBox[{"[", "i", "]"}], "]"}], "=", "1.0"}]}], "]"}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"g", "[", RowBox[{"[", "1", "]"}], "]"}], "=", "2.35018989"}], ";", RowBox[{ RowBox[{"g", "[", RowBox[{"[", "2", "]"}], "]"}], "=", "0.95875665"}], ";", RowBox[{ RowBox[{"g", "[", RowBox[{"[", "6", "]"}], "]"}], "=", RowBox[{"g", "[", RowBox[{"[", "2", "]"}], "]"}]}], ";", RowBox[{ RowBox[{"g", "[", RowBox[{"[", "3", "]"}], "]"}], "=", "1.21283096"}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"g", "[", RowBox[{"[", "4", "]"}], "]"}], "=", RowBox[{"g", "[", RowBox[{"[", "3", "]"}], "]"}]}], ";", RowBox[{ RowBox[{"g", "[", RowBox[{"[", "5", "]"}], "]"}], "=", RowBox[{"g", "[", RowBox[{"[", "3", "]"}], "]"}]}], ";", RowBox[{ RowBox[{"g", "[", RowBox[{"[", "8", "]"}], "]"}], "=", RowBox[{"g", "[", RowBox[{"[", "3", "]"}], "]"}]}], ";", RowBox[{ RowBox[{"g", "[", RowBox[{"[", "9", "]"}], "]"}], "=", RowBox[{"g", "[", RowBox[{"[", "3", "]"}], "]"}]}], ";", RowBox[{ RowBox[{"g", "[", RowBox[{"[", "7", "]"}], "]"}], "=", RowBox[{"g", "[", RowBox[{"[", "3", "]"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"sig", "[", RowBox[{"[", "1", "]"}], "]"}], "=", "0.31403448"}], ";", RowBox[{ RowBox[{"sig", "[", RowBox[{"[", "2", "]"}], "]"}], "=", "0.30821367"}], ";", RowBox[{ RowBox[{"sig", "[", RowBox[{"[", "6", "]"}], "]"}], "=", RowBox[{"sig", "[", RowBox[{"[", "2", "]"}], "]"}]}], ";", RowBox[{ RowBox[{"sig", "[", RowBox[{"[", "3", "]"}], "]"}], "=", "0.30820426"}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{ RowBox[{"sig", "[", RowBox[{"[", "4", "]"}], "]"}], "=", RowBox[{"sig", "[", RowBox[{"[", "3", "]"}], "]"}]}], ";", RowBox[{ RowBox[{"sig", "[", RowBox[{"[", "5", "]"}], "]"}], "=", RowBox[{"sig", "[", RowBox[{"[", "3", "]"}], "]"}]}], ";", RowBox[{ RowBox[{"sig", "[", RowBox[{"[", "8", "]"}], "]"}], "=", RowBox[{"sig", "[", RowBox[{"[", "3", "]"}], "]"}]}], ";", RowBox[{ RowBox[{"sig", "[", RowBox[{"[", "9", "]"}], "]"}], "=", RowBox[{"sig", "[", RowBox[{"[", "3", "]"}], "]"}]}], ";", RowBox[{ RowBox[{"sig", "[", RowBox[{"[", "7", "]"}], "]"}], "=", RowBox[{"sig", "[", RowBox[{"[", "3", "]"}], "]"}]}], ";"}], "\n"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"ss", "=", RowBox[{"Sin", "[", "x4", "]"}]}], ";", RowBox[{"cs", "=", RowBox[{"Cos", "[", "x4", "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"1", ",", "1"}], "]"}], "]"}], "=", RowBox[{"-", "2"}]}], ";", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"1", ",", "2"}], "]"}], "]"}], "=", "0"}], ";", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"1", ",", "3"}], "]"}], "]"}], "=", "x1"}], ";"}], "\n", RowBox[{"(*", RowBox[{"F", " ", "right"}], "*)"}]}], "\n", RowBox[{ RowBox[{ RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}], "=", RowBox[{"-", "2"}]}], ";", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"2", ",", "2"}], "]"}], "]"}], "=", "9"}], ";", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"2", ",", "3"}], "]"}], "]"}], "=", "x2"}], ";"}], "\[IndentingNewLine]"}], "\n", RowBox[{ RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"3", ",", "1"}], "]"}], "]"}], "=", RowBox[{"-", "2"}]}], ";", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"3", ",", "2"}], "]"}], "]"}], "=", "0"}], ";", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"3", ",", "3"}], "]"}], "]"}], "=", "x3"}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"4", ",", "1"}], "]"}], "]"}], "=", RowBox[{"-", "2"}]}], ";", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"4", ",", "2"}], "]"}], "]"}], "=", "0"}], ";", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"4", ",", "3"}], "]"}], "]"}], "=", "x3"}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"5", ",", "1"}], "]"}], "]"}], "=", RowBox[{"-", "2"}]}], ";", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"5", ",", "2"}], "]"}], "]"}], "=", "0"}], ";", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"5", ",", "3"}], "]"}], "]"}], "=", "x3"}], ";"}], "\n", RowBox[{"(*", RowBox[{"F", " ", "left"}], "*)"}]}], "\n", RowBox[{ RowBox[{ RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"6", ",", "1"}], "]"}], "]"}], "=", RowBox[{"-", "2"}]}], ";", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"6", ",", "2"}], "]"}], "]"}], "=", "9"}], ";", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"6", ",", "3"}], "]"}], "]"}], "=", "x2"}], ";"}], "\[IndentingNewLine]"}], "\n", RowBox[{ RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"7", ",", "1"}], "]"}], "]"}], "=", RowBox[{"-", "2"}]}], ";", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"7", ",", "2"}], "]"}], "]"}], "=", "0"}], ";", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"7", ",", "3"}], "]"}], "]"}], "=", "x3"}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"8", ",", "1"}], "]"}], "]"}], "=", RowBox[{"-", "2"}]}], ";", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"8", ",", "2"}], "]"}], "]"}], "=", "0"}], ";", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"8", ",", "3"}], "]"}], "]"}], "=", "x3"}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"9", ",", "1"}], "]"}], "]"}], "=", RowBox[{"-", "2"}]}], ";", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"9", ",", "2"}], "]"}], "]"}], "=", "0"}], ";", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"9", ",", "3"}], "]"}], "]"}], "=", "x3"}], ";"}], "\[IndentingNewLine]", "\n", RowBox[{"(*", RowBox[{ "coord", " ", "with", " ", "zero", " ", "value", " ", "not", " ", "explicitely", " ", RowBox[{"written", "!"}]}], "*)"}], "\n", RowBox[{"(*", "nuclei", "*)"}]}], "\n", RowBox[{ RowBox[{ RowBox[{"b", "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}], "=", RowBox[{"x1", "+", "x2"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{ RowBox[{"b", "[", RowBox[{"[", RowBox[{"6", ",", "1"}], "]"}], "]"}], "=", RowBox[{"-", RowBox[{"b", "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}]}]}], ";"}], "\n", RowBox[{"(*", "clouds", "*)"}]}], "\n", RowBox[{ RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}], "=", RowBox[{"b", "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\n", RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"3", ",", "1"}], "]"}], "]"}], "=", RowBox[{ RowBox[{"b", "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}], "+", RowBox[{ RowBox[{"(", RowBox[{"x2", "+", "x3"}], ")"}], "*", "cs"}]}]}], ";", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"3", ",", "2"}], "]"}], "]"}], "=", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"x2", "+", "x3"}], ")"}]}], "*", "ss"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"4", ",", "1"}], "]"}], "]"}], "=", RowBox[{"a", "[", RowBox[{"[", RowBox[{"3", ",", "1"}], "]"}], "]"}]}], ";", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"4", ",", "2"}], "]"}], "]"}], "=", RowBox[{ RowBox[{"-", "0.5"}], "*", RowBox[{"a", "[", RowBox[{"[", RowBox[{"3", ",", "2"}], "]"}], "]"}]}]}], ";", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"4", ",", "3"}], "]"}], "]"}], "=", RowBox[{"0.5", "*", "sq3", "*", RowBox[{"a", "[", RowBox[{"[", RowBox[{"3", ",", "2"}], "]"}], "]"}]}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"5", ",", "1"}], "]"}], "]"}], "=", RowBox[{"a", "[", RowBox[{"[", RowBox[{"3", ",", "1"}], "]"}], "]"}]}], ";", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"5", ",", "2"}], "]"}], "]"}], "=", RowBox[{"a", "[", RowBox[{"[", RowBox[{"4", ",", "2"}], "]"}], "]"}]}], ";", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"5", ",", "3"}], "]"}], "]"}], "=", RowBox[{"-", RowBox[{"a", "[", RowBox[{"[", RowBox[{"4", ",", "3"}], "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]"}], "\n", RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"6", ",", "1"}], "]"}], "]"}], "=", RowBox[{"-", RowBox[{"b", "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}]}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"7", ",", "1"}], "]"}], "]"}], "=", RowBox[{"-", RowBox[{"a", "[", RowBox[{"[", RowBox[{"3", ",", "1"}], "]"}], "]"}]}]}], ";", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"7", ",", "2"}], "]"}], "]"}], "=", RowBox[{"-", RowBox[{"a", "[", RowBox[{"[", RowBox[{"3", ",", "2"}], "]"}], "]"}]}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"8", ",", "1"}], "]"}], "]"}], "=", RowBox[{"-", RowBox[{"a", "[", RowBox[{"[", RowBox[{"3", ",", "1"}], "]"}], "]"}]}]}], ";", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"8", ",", "2"}], "]"}], "]"}], "=", RowBox[{"-", RowBox[{"a", "[", RowBox[{"[", RowBox[{"4", ",", "2"}], "]"}], "]"}]}]}], ";", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"8", ",", "3"}], "]"}], "]"}], "=", RowBox[{"a", "[", RowBox[{"[", RowBox[{"4", ",", "3"}], "]"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"9", ",", "1"}], "]"}], "]"}], "=", RowBox[{"-", RowBox[{"a", "[", RowBox[{"[", RowBox[{"3", ",", "1"}], "]"}], "]"}]}]}], ";", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"9", ",", "2"}], "]"}], "]"}], "=", RowBox[{"-", RowBox[{"a", "[", RowBox[{"[", RowBox[{"4", ",", "2"}], "]"}], "]"}]}]}], ";", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"9", ",", "3"}], "]"}], "]"}], "=", RowBox[{"a", "[", RowBox[{"[", RowBox[{"5", ",", "3"}], "]"}], "]"}]}], ";"}], "\n", "\n", RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{"CALC", ".", "OF"}], " ", "ENERGY"}], "-", "TERMS"}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"KI", "/", "KJ"}], "*)"}]}], "\n", RowBox[{ RowBox[{ RowBox[{"vnn", "=", RowBox[{"9", "*", RowBox[{"9", "/", RowBox[{"(", RowBox[{"2", "*", RowBox[{"(", RowBox[{"x1", "+", "x2"}], ")"}]}], ")"}]}]}]}], ";", " ", RowBox[{"vee", "=", "0.0"}], ";", RowBox[{"vne", "=", "0.0"}], ";", RowBox[{"ekin", "=", "0.0"}], ";"}], "\n", RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{"WI", "/", "WJ"}], " ", RowBox[{"EXCL", ".", "I"}]}], "=", "J"}], "*)"}]}], "\n", RowBox[{ RowBox[{"For", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", " ", "nc"}], " ", ",", RowBox[{"i", "++"}], ",", "\n", RowBox[{"For", "[", " ", RowBox[{ RowBox[{"j", "=", RowBox[{"i", "+", "1"}]}], ",", RowBox[{"j", "<", RowBox[{"nc", "+", "1"}]}], ",", RowBox[{"j", "++"}], ",", "\n", RowBox[{ RowBox[{"d", "=", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "-", RowBox[{"a", "[", RowBox[{"[", RowBox[{"j", ",", "1"}], "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "-", RowBox[{"a", "[", RowBox[{"[", RowBox[{"j", ",", "2"}], "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "-", RowBox[{"a", "[", RowBox[{"[", RowBox[{"j", ",", "3"}], "]"}], "]"}]}], ")"}], "^", "2"}]}], "]"}]}], ";", "\n", RowBox[{"If", "[", " ", RowBox[{ RowBox[{"d", ">", "0"}], " ", ",", RowBox[{"vee", "=", RowBox[{"vee", "+", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "*", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"j", ",", "1"}], "]"}], "]"}], "/", "d"}]}]}]}]}], "]"}]}]}], "]"}]}], "]"}], "\n", RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{"KJ", "/", "WI"}], " ", "I"}], "=", "J"}], "*)"}]}], "\n", RowBox[{ RowBox[{"vne", "=", RowBox[{"vne", "+", RowBox[{"1.5", "*", RowBox[{"c", "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}], "*", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"2", ",", "2"}], "]"}], "]"}], "/", RowBox[{"c", "[", RowBox[{"[", RowBox[{"2", ",", "3"}], "]"}], "]"}]}]}]}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"vne", "=", RowBox[{"vne", "+", RowBox[{"1.5", "*", RowBox[{"c", "[", RowBox[{"[", RowBox[{"6", ",", "1"}], "]"}], "]"}], "*", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"6", ",", "2"}], "]"}], "]"}], "/", RowBox[{"c", "[", RowBox[{"[", RowBox[{"6", ",", "3"}], "]"}], "]"}]}]}]}]}], ";"}], "\[IndentingNewLine]", "\n", RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{"KJ", "/", "WI"}], " ", RowBox[{"I", "<>", "J"}], " ", "inc"}], ";"}], "*)"}]}], "\n", RowBox[{ RowBox[{"For", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", RowBox[{"nc", "+", "1"}]}], ",", RowBox[{"i", "++"}], ",", "\n", RowBox[{ RowBox[{"If", "[", " ", RowBox[{ RowBox[{"i", "\[NotEqual]", "2"}], " ", ",", "\n", RowBox[{ RowBox[{"d", "=", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "-", RowBox[{"b", "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], ")"}], "^", "2"}]}], "]"}]}], ";", "\n", RowBox[{"If", "[", " ", RowBox[{ RowBox[{"d", ">", "0"}], ",", RowBox[{"vne", "=", RowBox[{"vne", "+", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "*", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"2", ",", "2"}], "]"}], "]"}], "/", "d"}]}]}]}]}], "]"}]}]}], "]"}], ";", "\n", RowBox[{"If", "[", " ", RowBox[{ RowBox[{"i", "\[NotEqual]", "6"}], ",", "\n", RowBox[{ RowBox[{"d", "=", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "-", RowBox[{"b", "[", RowBox[{"[", RowBox[{"6", ",", "1"}], "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], ")"}], "^", "2"}]}], "]"}]}], ";", "\n", RowBox[{"If", "[", " ", RowBox[{ RowBox[{"d", ">", "0"}], " ", ",", RowBox[{"vne", "=", RowBox[{"vne", "+", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "*", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"6", ",", "2"}], "]"}], "]"}], "/", "d"}]}]}]}]}], "]"}]}]}], "]"}]}]}], "]"}], "\n"}], "\n", RowBox[{ RowBox[{"For", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", " ", RowBox[{"nc", "+", "1"}]}], ",", RowBox[{"i", "++"}], ",", "\n", RowBox[{"(*", "T", "*)"}], "\n", RowBox[{ RowBox[{"ekin", "=", RowBox[{"ekin", "-", RowBox[{"1.125", "*", RowBox[{"g", "[", RowBox[{"[", "i", "]"}], "]"}], "*", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "/", RowBox[{ RowBox[{"c", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "^", "2"}]}]}]}]}], ";", "\n", RowBox[{"(*", RowBox[{"WI", "/", "WI"}], "*)"}], "\n", RowBox[{"vee", "=", RowBox[{"vee", "+", RowBox[{"3.0", "*", RowBox[{ RowBox[{"sig", "[", RowBox[{"[", "i", "]"}], "]"}], "/", RowBox[{"c", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}]}]}]}]}]}]}], "]"}], "\n", "\n"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Epot", "=", RowBox[{"vne", "+", "vee", "+", "vnn"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"func", "=", RowBox[{"Ekin", "+", "Epot"}]}], ";"}], "\n"}], "\n", RowBox[{"t", "=", RowBox[{"FindMinimum", "[", RowBox[{"func", ",", RowBox[{"{", RowBox[{"x1", ",", "1.17843"}], "}"}], ",", RowBox[{"{", RowBox[{"x2", ",", "0.1595"}], "}"}], ",", RowBox[{"{", RowBox[{"x3", ",", "0.8044"}], "}"}], ",", RowBox[{"{", RowBox[{"x4", ",", "1.34"}], "}"}], ",", RowBox[{"Method", "\[Rule]", "\"\\""}], ",", RowBox[{"MaxIterations", "\[Rule]", "500"}]}], "]"}]}]}], "Input", CellChangeTimes->{{3.577942230718297*^9, 3.577942230718297*^9}, { 3.577942264180356*^9, 3.577942265131957*^9}, {3.577942319014452*^9, 3.5779424054386044`*^9}, {3.577942440101865*^9, 3.5779424567158937`*^9}, { 3.5779425225012093`*^9, 3.5779425465408516`*^9}, {3.5779425971161404`*^9, 3.577942597365741*^9}, {3.577942660701852*^9, 3.577942690669505*^9}, { 3.5779427231487617`*^9, 3.577942868431817*^9}, {3.5779429100370903`*^9, 3.577943391641136*^9}, 3.577943422747591*^9, {3.577943618964735*^9, 3.5779436491819887`*^9}, {3.5779437045308857`*^9, 3.577943727384926*^9}, { 3.577943802077857*^9, 3.5779438458203335`*^9}, {3.577943880187194*^9, 3.577943888221208*^9}, {3.577943985190978*^9, 3.577944388373687*^9}, { 3.5779445035798893`*^9, 3.5779445527355757`*^9}, {3.5779445893488398`*^9, 3.577944718829067*^9}, {3.577944766611951*^9, 3.5779448631605206`*^9}, { 3.5779675928735037`*^9, 3.577967600143117*^9}, {3.577969013512399*^9, 3.5779690139960003`*^9}}] }, WindowSize->{1014, 849}, WindowMargins->{{19, Automatic}, {19, Automatic}}, 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[557, 20, 22231, 733, 1552, "Input"] } ] *) (* End of internal cache information *)