(* 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[     27894,       1206]
NotebookOptionsPosition[     26280,       1148]
NotebookOutlinePosition[     26793,       1167]
CellTagsIndexPosition[     26750,       1164]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"CH4", " ", "molecule", " ", "with", " ", 
     RowBox[{"C", "[", 
      RowBox[{"1", "s2", " ", "4", "sp"}], "]"}]}], ",", " ", 
    RowBox[{"no", " ", "overlap"}], ",", " ", 
    RowBox[{"parametrized", " ", "ES", " ", "06.11", ".2012"}]}], " ", "*)"}],
   "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"Variables", ":", " ", "P"}], ",", "Q", ",", "R"}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"Clear", "[", 
     RowBox[{
     "k1", ",", "k2", ",", "s1", ",", "s2", ",", " ", "z", ",", " ", "P", ",",
       "Q", ",", "R", ",", "Vee", ",", "Vne", ",", "Vnn", ",", "T", ",", "d1",
       ",", "d2", ",", "d3", ",", "par", ",", "w3", ",", "ad"}], "]"}], ";"}],
    "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Clear", "[", 
     RowBox[{
     "Vnec", ",", "Veec", ",", "Vneb", ",", "Veeb", ",", "Veecb", ",", 
      "Vnecb", ",", "Vnncb", ",", "Vnebb", ",", "Veebb", ",", "Vnnbb"}], 
     "]"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Clear", "[", 
     RowBox[{"Tc", ",", "Tb", ",", "wd"}], "]"}], ";"}]}]}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"z", "=", "6.0"}], ";", " ", 
    RowBox[{"wd", "=", "0.94281"}], ";", " ", 
    RowBox[{"ad", "=", 
     RowBox[{"Sqrt", "[", 
      RowBox[{"8", "/", "3"}], "]"}]}], ";", " ", 
    RowBox[{"w3", "=", 
     RowBox[{"Sqrt", "[", "3", "]"}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]"}], 
  RowBox[{"(*", " ", 
   RowBox[{"C", " ", "core"}], " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Tc", "=", 
   RowBox[{"2.25", "*", 
    RowBox[{"k1", "/", 
     RowBox[{"P", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Vnec", "=", 
   RowBox[{
    RowBox[{"-", "3.0"}], "*", 
    RowBox[{"z", "/", "P"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"Veec", "=", 
     RowBox[{"3.0", "*", 
      RowBox[{"s1", "/", "P"}]}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]"}], 
  RowBox[{"(*", " ", 
   RowBox[{"CH", " ", "bond"}], " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Tb", "=", 
   RowBox[{"2.25", "*", 
    RowBox[{"k2", "/", 
     RowBox[{"Q", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Vneb", "=", 
   RowBox[{
    RowBox[{"-", "2.0"}], "*", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"1.5", "-", 
       RowBox[{"0.5", "*", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"R", "/", "Q"}], ")"}], "^", "2"}]}]}], ")"}], "/", 
     "Q"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"Veeb", "=", 
     RowBox[{"3.0", "*", 
      RowBox[{"s2", "/", "Q"}]}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]"}], 
  RowBox[{"(*", " ", 
   RowBox[{"core", "-", "bond"}], " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Vnecb", "=", 
   RowBox[{
    RowBox[{
     RowBox[{"-", "2.0"}], "*", 
     RowBox[{"z", "/", 
      RowBox[{"(", 
       RowBox[{"P", "+", "Q"}], ")"}]}]}], "-", 
    RowBox[{"2", "/", 
     RowBox[{"(", 
      RowBox[{"P", "+", "Q", "+", "R"}], ")"}]}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Veecb", "=", 
   RowBox[{"4.0", "/", 
    RowBox[{"(", 
     RowBox[{"P", "+", "Q"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"Vnncb", "=", 
     RowBox[{"z", "/", 
      RowBox[{"(", 
       RowBox[{"P", "+", "Q", "+", "R"}], ")"}]}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]"}], 
  RowBox[{"(*", " ", 
   RowBox[{"bond", "-", "bond"}], " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"d1", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"P", "+", "Q", "+", "R"}], ")"}], "*", "ad"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"d2", "=", 
   RowBox[{"Sqrt", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{"P", "+", "Q"}], ")"}], "^", "2"}], "+", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"P", "+", "Q", "+", "R"}], ")"}], "^", "2"}], "+", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"2", "/", "3"}], ")"}], "*", 
      RowBox[{"(", 
       RowBox[{"P", "+", "Q"}], ")"}], "*", 
      RowBox[{"(", 
       RowBox[{"P", "+", "Q", "+", "R"}], ")"}]}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"d3", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"P", "+", "Q"}], ")"}], "*", "ad"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Vnebb", "=", 
   RowBox[{
    RowBox[{"-", "4.0"}], "/", "d2"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Veebb", "=", 
   RowBox[{"4.0", "/", "d3"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"Vnnbb", "=", 
     RowBox[{"1.0", "/", "d1"}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]"}], 
  RowBox[{"(*", " ", 
   RowBox[{"sum", " ", "terms"}], " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"T", " ", "=", " ", 
   RowBox[{"Tc", "+", 
    RowBox[{"4", "*", "Tb"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Vne", " ", "=", " ", 
   RowBox[{"Vnec", "+", 
    RowBox[{"4", "*", "Vneb"}], "+", 
    RowBox[{"4", "*", "Vnecb"}], "+", 
    RowBox[{"6", "*", "Vnebb"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Vee", " ", "=", " ", 
  RowBox[{"Veec", "+", 
   RowBox[{"4", "*", "Veeb"}], "+", 
   RowBox[{"4", "*", "Veecb"}], "+", 
   RowBox[{"6", "*", "Veebb"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Vnn", " ", "=", " ", 
   RowBox[{
    RowBox[{"4", "*", "Vnncb"}], "+", 
    RowBox[{"6", "*", "Vnnbb"}]}]}], ";"}]}], "Input"],

Cell[BoxData[
 RowBox[{
  FractionBox["30.696938456699066`", 
   RowBox[{"P", "+", "Q"}]], "+", 
  FractionBox[
   RowBox[{"3.`", " ", "s1"}], "P"], "+", 
  FractionBox[
   RowBox[{"12.`", " ", "s2"}], "Q"]}]], "Output"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"par", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"k1", " ", "->", " ", "1.0"}], ",", " ", 
       RowBox[{"k2", "->", " ", "1.0"}], ",", " ", 
       RowBox[{"s1", " ", "->", " ", "0.4"}], ",", " ", 
       RowBox[{"s2", " ", "->", " ", "0.4"}]}], "}"}]}], ";"}], " ", "*)"}], 
  " ", "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"par", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"k1", " ", "\[Rule]", " ", "1.01613950"}], ",", " ", 
       RowBox[{"k2", " ", "->", " ", "0.65858068"}], ",", 
       RowBox[{"s1", "->", " ", 
        RowBox[{"0.25742249", "\[IndentingNewLine]", 
         RowBox[{"s2", "->", " ", "0.22958899"}]}]}]}], "}"}]}], ";"}], 
   "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"par", "=", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"k1", "->", "1.05518276"}], ",", 
      RowBox[{"k2", "->", " ", "1.21738169"}], ",", 
      RowBox[{"s1", "->", " ", "0.24406300"}], ",", 
      RowBox[{"s2", "->", " ", "0.30451569"}]}], "}"}]}], ";"}], 
  " "}]], "Input"],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"func", "=", 
    RowBox[{
     RowBox[{"T", "+", "Vne", "+", "Vee", "+", "Vnn"}], " ", "/.", " ", 
     "par"}]}], ";"}]}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"t", "=", 
   RowBox[{"FindMinimum", "[", 
    RowBox[{"func", ",", 
     RowBox[{"{", 
      RowBox[{"P", ",", "0.26"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"Q", ",", "1.11"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"R", ",", "0.476"}], "}"}]}], "]"}], " "}]}]], "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "40.45715398294954`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"P", "\[Rule]", "0.2667945964258979`"}], ",", 
     RowBox[{"Q", "\[Rule]", "1.2442368403134039`"}], ",", 
     RowBox[{"R", "\[Rule]", "0.5370596483068801`"}]}], "}"}]}], 
  "}"}]], "Output"]
}, Open  ]],

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

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Vne", " ", "/.", "u"}]], "Input"],

Cell[BoxData[
 RowBox[{"-", "120.39772695342145`"}]], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Vee", " ", "/.", "u"}], "  ", "/.", " ", "par"}]], "Input"],

Cell[BoxData["25.99650499051755`"], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Vnn", " ", "/.", " ", "u"}]], "Input"],

Cell[BoxData["13.512208913086864`"], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"-", 
     RowBox[{"(", 
      RowBox[{"Vne", "+", "Vee", "+", "Vnn"}], ")"}]}], "/", "T"}], " ", "/.",
    " ", "u"}], " ", "/.", "par"}]], "Input"],

Cell[BoxData["2.000625618427296`"], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"(", 
    RowBox[{"P", "+", "Q", "+", "R"}], ")"}], "*", "0.529177"}], " ", "/.", 
  " ", "u"}]], "Input"],

Cell[BoxData["1.0838026961114835`"], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"plot1", "=", 
  RowBox[{
   RowBox[{
    RowBox[{"Graphics", "[", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{"0", ",", "0"}], "}"}], ",", "P"}], "]"}], ",", 
       RowBox[{"Disk", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{"0", ",", "0"}], "}"}], ",", "0.1"}], "]"}], ",", 
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{"0", ",", "0"}], "}"}], ",", "Q"}], "]"}], ",", "\n", "\t", 
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{
            RowBox[{"-", "wd"}], "*", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"P", "+", "Q"}], ")"}], "/", "2"}]}], ",", 
           RowBox[{"wd", "*", "w3", "*", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"P", "+", "Q"}], ")"}], "/", "2"}]}]}], "}"}], ",", 
         "Q"}], "]"}], ",", 
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{
            RowBox[{"-", "wd"}], "*", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"P", "+", "Q"}], ")"}], "/", "2"}]}], ",", 
           RowBox[{
            RowBox[{"-", "wd"}], "*", "w3", "*", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"P", "+", "Q"}], ")"}], "/", "2"}]}]}], "}"}], ",", 
         "Q"}], "]"}], ",", 
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"wd", "*", 
            RowBox[{"(", 
             RowBox[{"P", "+", "Q"}], ")"}]}], ",", "0"}], "}"}], ",", "Q"}], 
        "]"}], ",", 
       RowBox[{"Disk", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{
            RowBox[{"-", "wd"}], "*", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"P", "+", "Q", "+", "R"}], ")"}], "/", "2"}]}], ",", 
           RowBox[{"wd", "*", "w3", "*", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"P", "+", "Q", "+", "R"}], ")"}], "/", "2"}]}]}], "}"}],
          ",", "0.1"}], "]"}], ",", 
       RowBox[{"Disk", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{
            RowBox[{"-", "wd"}], "*", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"P", "+", "Q", "+", "R"}], ")"}], "/", "2"}]}], ",", 
           RowBox[{
            RowBox[{"-", "wd"}], "*", "w3", "*", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"P", "+", "Q", "+", "R"}], ")"}], "/", "2"}]}]}], "}"}],
          ",", "0.1"}], "]"}], ",", 
       RowBox[{"Disk", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"wd", "*", 
            RowBox[{"(", 
             RowBox[{"P", "+", "Q", "+", "R"}], ")"}]}], ",", "0"}], "}"}], 
         ",", "0.1"}], "]"}], ",", "\n", "\t", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Thickness", "[", "0.012", "]"}], ",", 
         RowBox[{"Line", "[", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"0", ",", "0"}], "}"}], ",", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"wd", "*", 
               RowBox[{"(", 
                RowBox[{"P", "+", "Q", "+", "R"}], ")"}]}], ",", "0"}], 
             "}"}]}], "}"}], "]"}], ",", 
         RowBox[{"Line", "[", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"0", ",", "0"}], "}"}], ",", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{
               RowBox[{"-", "wd"}], "*", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{"P", "+", "Q", "+", "R"}], ")"}], "/", "2"}]}], ",", 
              RowBox[{"wd", "*", "w3", "*", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{"P", "+", "Q", "+", "R"}], ")"}], "/", "2"}]}]}], 
             "}"}]}], "}"}], "]"}], ",", 
         RowBox[{"Line", "[", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"0", ",", "0"}], "}"}], ",", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{
               RowBox[{"-", "wd"}], "*", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{"P", "+", "Q", "+", "R"}], ")"}], "/", "2"}]}], ",", 
              RowBox[{
               RowBox[{"-", "wd"}], "*", "w3", "*", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{"P", "+", "Q", "+", "R"}], ")"}], "/", "2"}]}]}], 
             "}"}]}], "}"}], "]"}]}], "}"}]}], "}"}], "]"}], "/.", " ", "u"}],
    "\n", "       "}]}]], "Input"],

Cell[BoxData[
 TagBox[
  RowBox[{"\[SkeletonIndicator]", "Graphics", "\[SkeletonIndicator]"}],
  False,
  Editable->False]], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 StyleBox[
  RowBox[{"Show", "[", 
   RowBox[{"plot1", ",", " ", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"AspectRatio", "->", "Automatic"}], ",", " ", 
      RowBox[{"PlotRange", "->", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"-", "3"}], ",", "3"}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "}"}]}], ",", 
      RowBox[{"Axes", "->", "True"}], ",", 
      RowBox[{"GridLines", "->", "Automatic"}], ",", 
      RowBox[{"Frame", "->", "True"}]}], "}"}]}], "]"}],
  FontSize->12]], "Input"],

Cell[GraphicsData["PostScript", "\<\
%!
%%Creator: Mathematica
%%AspectRatio: 1 
MathPictureStart
/Mabs {
Mgmatrix idtransform
Mtmatrix dtransform
} bind def
/Mabsadd { Mabs
3 -1 roll add
3 1 roll add
exch } bind def
%% Graphics
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10  scalefont  setfont
% Scaling calculations
0.5 0.166667 0.5 0.166667 [
[.16667 -0.0125 -6 -9 ]
[.16667 -0.0125 6 0 ]
[.33333 -0.0125 -6 -9 ]
[.33333 -0.0125 6 0 ]
[.5 -0.0125 -3 -9 ]
[.5 -0.0125 3 0 ]
[.66667 -0.0125 -3 -9 ]
[.66667 -0.0125 3 0 ]
[.83333 -0.0125 -3 -9 ]
[.83333 -0.0125 3 0 ]
[1 -0.0125 -3 -9 ]
[1 -0.0125 3 0 ]
[ 0 0 -0.125 0 ]
[-0.0125 .16667 -12 -4.5 ]
[-0.0125 .16667 0 4.5 ]
[-0.0125 .33333 -12 -4.5 ]
[-0.0125 .33333 0 4.5 ]
[-0.0125 .5 -6 -4.5 ]
[-0.0125 .5 0 4.5 ]
[-0.0125 .66667 -6 -4.5 ]
[-0.0125 .66667 0 4.5 ]
[-0.0125 .83333 -6 -4.5 ]
[-0.0125 .83333 0 4.5 ]
[-0.0125 1 -6 -4.5 ]
[-0.0125 1 0 4.5 ]
[ 0 0 -0.125 0 ]
[ 0 1 .125 0 ]
[ 1 0 .125 0 ]
[ 0 0 0 0 ]
[ 1 1 0 0 ]
] MathScale
% Start of Graphics
1 setlinecap
1 setlinejoin
newpath
0 0 .5 r
.25 Mabswid
[ ] 0 setdash
.16667 0 m
.16667 1 L
s
.33333 0 m
.33333 1 L
s
.5 0 m
.5 1 L
s
.66667 0 m
.66667 1 L
s
.83333 0 m
.83333 1 L
s
0 .16667 m
1 .16667 L
s
0 .33333 m
1 .33333 L
s
0 .5 m
1 .5 L
s
0 .66667 m
1 .66667 L
s
0 .83333 m
1 .83333 L
s
0 g
.16667 0 m
.16667 .00625 L
s
[(-2)] .16667 -0.0125 0 1 Mshowa
.33333 0 m
.33333 .00625 L
s
[(-1)] .33333 -0.0125 0 1 Mshowa
.5 0 m
.5 .00625 L
s
[(0)] .5 -0.0125 0 1 Mshowa
.66667 0 m
.66667 .00625 L
s
[(1)] .66667 -0.0125 0 1 Mshowa
.83333 0 m
.83333 .00625 L
s
[(2)] .83333 -0.0125 0 1 Mshowa
1 0 m
1 .00625 L
s
[(3)] 1 -0.0125 0 1 Mshowa
.125 Mabswid
.03333 0 m
.03333 .00375 L
s
.06667 0 m
.06667 .00375 L
s
.1 0 m
.1 .00375 L
s
.13333 0 m
.13333 .00375 L
s
.2 0 m
.2 .00375 L
s
.23333 0 m
.23333 .00375 L
s
.26667 0 m
.26667 .00375 L
s
.3 0 m
.3 .00375 L
s
.36667 0 m
.36667 .00375 L
s
.4 0 m
.4 .00375 L
s
.43333 0 m
.43333 .00375 L
s
.46667 0 m
.46667 .00375 L
s
.53333 0 m
.53333 .00375 L
s
.56667 0 m
.56667 .00375 L
s
.6 0 m
.6 .00375 L
s
.63333 0 m
.63333 .00375 L
s
.7 0 m
.7 .00375 L
s
.73333 0 m
.73333 .00375 L
s
.76667 0 m
.76667 .00375 L
s
.8 0 m
.8 .00375 L
s
.86667 0 m
.86667 .00375 L
s
.9 0 m
.9 .00375 L
s
.93333 0 m
.93333 .00375 L
s
.96667 0 m
.96667 .00375 L
s
.25 Mabswid
0 0 m
1 0 L
s
0 .16667 m
.00625 .16667 L
s
[(-2)] -0.0125 .16667 1 0 Mshowa
0 .33333 m
.00625 .33333 L
s
[(-1)] -0.0125 .33333 1 0 Mshowa
0 .5 m
.00625 .5 L
s
[(0)] -0.0125 .5 1 0 Mshowa
0 .66667 m
.00625 .66667 L
s
[(1)] -0.0125 .66667 1 0 Mshowa
0 .83333 m
.00625 .83333 L
s
[(2)] -0.0125 .83333 1 0 Mshowa
0 1 m
.00625 1 L
s
[(3)] -0.0125 1 1 0 Mshowa
.125 Mabswid
0 .03333 m
.00375 .03333 L
s
0 .06667 m
.00375 .06667 L
s
0 .1 m
.00375 .1 L
s
0 .13333 m
.00375 .13333 L
s
0 .2 m
.00375 .2 L
s
0 .23333 m
.00375 .23333 L
s
0 .26667 m
.00375 .26667 L
s
0 .3 m
.00375 .3 L
s
0 .36667 m
.00375 .36667 L
s
0 .4 m
.00375 .4 L
s
0 .43333 m
.00375 .43333 L
s
0 .46667 m
.00375 .46667 L
s
0 .53333 m
.00375 .53333 L
s
0 .56667 m
.00375 .56667 L
s
0 .6 m
.00375 .6 L
s
0 .63333 m
.00375 .63333 L
s
0 .7 m
.00375 .7 L
s
0 .73333 m
.00375 .73333 L
s
0 .76667 m
.00375 .76667 L
s
0 .8 m
.00375 .8 L
s
0 .86667 m
.00375 .86667 L
s
0 .9 m
.00375 .9 L
s
0 .93333 m
.00375 .93333 L
s
0 .96667 m
.00375 .96667 L
s
.25 Mabswid
0 0 m
0 1 L
s
0 .99375 m
0 1 L
s
.16667 .99375 m
.16667 1 L
s
.33333 .99375 m
.33333 1 L
s
.5 .99375 m
.5 1 L
s
.66667 .99375 m
.66667 1 L
s
.83333 .99375 m
.83333 1 L
s
.125 Mabswid
.03333 .99625 m
.03333 1 L
s
.06667 .99625 m
.06667 1 L
s
.1 .99625 m
.1 1 L
s
.13333 .99625 m
.13333 1 L
s
.2 .99625 m
.2 1 L
s
.23333 .99625 m
.23333 1 L
s
.26667 .99625 m
.26667 1 L
s
.3 .99625 m
.3 1 L
s
.36667 .99625 m
.36667 1 L
s
.4 .99625 m
.4 1 L
s
.43333 .99625 m
.43333 1 L
s
.46667 .99625 m
.46667 1 L
s
.53333 .99625 m
.53333 1 L
s
.56667 .99625 m
.56667 1 L
s
.6 .99625 m
.6 1 L
s
.63333 .99625 m
.63333 1 L
s
.7 .99625 m
.7 1 L
s
.73333 .99625 m
.73333 1 L
s
.76667 .99625 m
.76667 1 L
s
.8 .99625 m
.8 1 L
s
.86667 .99625 m
.86667 1 L
s
.9 .99625 m
.9 1 L
s
.93333 .99625 m
.93333 1 L
s
.96667 .99625 m
.96667 1 L
s
.25 Mabswid
0 1 m
1 1 L
s
.99375 0 m
1 0 L
s
.99375 .16667 m
1 .16667 L
s
.99375 .33333 m
1 .33333 L
s
.99375 .5 m
1 .5 L
s
.99375 .66667 m
1 .66667 L
s
.99375 .83333 m
1 .83333 L
s
.125 Mabswid
.99625 .03333 m
1 .03333 L
s
.99625 .06667 m
1 .06667 L
s
.99625 .1 m
1 .1 L
s
.99625 .13333 m
1 .13333 L
s
.99625 .2 m
1 .2 L
s
.99625 .23333 m
1 .23333 L
s
.99625 .26667 m
1 .26667 L
s
.99625 .3 m
1 .3 L
s
.99625 .36667 m
1 .36667 L
s
.99625 .4 m
1 .4 L
s
.99625 .43333 m
1 .43333 L
s
.99625 .46667 m
1 .46667 L
s
.99625 .53333 m
1 .53333 L
s
.99625 .56667 m
1 .56667 L
s
.99625 .6 m
1 .6 L
s
.99625 .63333 m
1 .63333 L
s
.99625 .7 m
1 .7 L
s
.99625 .73333 m
1 .73333 L
s
.99625 .76667 m
1 .76667 L
s
.99625 .8 m
1 .8 L
s
.99625 .86667 m
1 .86667 L
s
.99625 .9 m
1 .9 L
s
.99625 .93333 m
1 .93333 L
s
.99625 .96667 m
1 .96667 L
s
.25 Mabswid
1 0 m
1 1 L
s
0 .5 m
1 .5 L
s
.5 0 m
.5 1 L
s
0 0 m
1 0 L
1 1 L
0 1 L
closepath
clip
newpath
.5 Mabswid
newpath
.5 .5 .04447 0 365.73 arc
s
.5 .5 m
.5 .5 .01667 0 365.73 arc
F
newpath
.5 .5 .20737 0 365.73 arc
s
newpath
.38128 .70563 .20737 0 365.73 arc
s
newpath
.38128 .29437 .20737 0 365.73 arc
s
newpath
.73744 .5 .20737 0 365.73 arc
s
.33909 .77871 m
.33909 .77871 .01667 0 365.73 arc
F
.33909 .22129 m
.33909 .22129 .01667 0 365.73 arc
F
.82183 .5 m
.82183 .5 .01667 0 365.73 arc
F
.012 w
.5 .5 m
.82183 .5 L
s
.5 .5 m
.33909 .77871 L
s
.5 .5 m
.33909 .22129 L
s
% End of Graphics
MathPictureEnd
\
\>"], "Graphics",
 ImageSize->{288, 288},
 ImageMargins->{{43, 0}, {0, 0}},
 ImageRegion->{{0, 1}, {0, 
  1}},ImageCache->GraphicsData["CompressedBitmap", "\<\
eJztXclyI7kRLYmLqKUX9bRa6pbUpPal7YsdYftg+2jLB3+Do2PCETOHCTvG
8lz1I/oXfkp/CkfIqsqs5QHM2qvIYow4bBSW9xJZQAJIAP/8+vjDv3/6+vjj
919nf//5639/+PH7/83+9p+fX4IGG563MTN/M8/8Xnhe+LV4+W9xZr5MiPz+
M/1v7C3mj1788+jNF4vr8PFi7sUjPHjP3xZX9HjofXt+kAfP327dqW7cZd6C
TL3H+aYP+C/0v4HnP6QE3n0QZjJ8KZ7+5SMfec+Pz/TL/D0/+PHvQPw7kO9d
CPTbs/fwEOZjOBGPDR/RXyO5+VBxCVdRqTwIqvkj5YaS3IOsw7BQeCYskM0f
EtHNZ5O+xxQtLDHkMZ8HaoCSjOJJAlTX2Ul8CcP8EuPEItlsJDQyJXSqlBuu
p5d/BqXchHX9QCCi0b1bEN9er5EiQvVNZbkRVUSL7otABd/QL48VOiasMOo9
zvUO57CZfqGL/H5j/kftxRvzNYiEvvzDfMLqHrcobNQiLH1YsbANu8rdc9gT
h92BsHtlWndYD6mH1ENqHhLZIl9qAEL9NiV1dHttkEgPqYfUQ1pzSH2z2IFK
6iFVCInU8a1JedgCcGMT6635NTK/3jULieYb/Fmc+IMj803ptVmV2jhdcvGb
4OkxP90rAC53g3XGCUfg6WcmttNs3U4ZyAQ8nTGJrWZhnjCQHfD0gkmgyZsa
YX5iIK/A0yt+OmgW5kcQDwGOvfhUQpAAvY0VYT0C8V6DeHeM9Tf1gfvAhe7b
CrC2XhVBOgDxrOAi8cyn8hb7PUj4zlrAMBbv5eO2JUpotN8x+gNrUaNYLqHo
bpXQSqjhfRDPDjcebxOEXdUH/S2I90FZFGozL+uDjuIdKosagrCL+qC/BvGO
lEWhDvW8PuivQLyPyqLQgsBZcejjMHO83Kzt2pFR+kmJBBljs+LMwmyja9Za
g6oInQkImxanE2YbXe2mN5Fi/4myGQRV90R/DyBsG4T9Yr6j/YuXaZEcNQeI
GxoanChlgNKelqj8eKlbq/zbmdFVzsxX/tRKubaukA5/ViLYzVzPpbaxRaCj
1/64PuioIXS3HCsBHfWK7sa2VOioY3V3P62Bjrozt02ARp1aW0QP3Tpiagpw
7iEUajDdBqPEQ2W4DdVSlaPD0FHHS9BpJvkWPEVDMfd4pnLApNcTE3IDnqJh
r3vsuN6A0VCfAG/bAKNlBvfEQuWAz50Srh3wKJYFRX75XDUGMtJtYGjkX7tl
g/YdCHuvhJahg8DQbhjadcnQclXtMJaFFzaawyT2QecQD1lHW4oYNAODDmIm
OV868SEeTWEesJxbgVmyMB/rEj0J+aIbgIcM+LxmwHA4YYM5YkUoG2auQYQT
JknzDOSClr3sK0u1VP+Y2XcJMEl41g3AWyzhLgEmCU+7AXjCgD93A/A2qwTN
dOyaX8gCRzB3OYXboigfMEmYWlVaXUaTPtSC0mgeTVwiOmPOb5d/vclBLFNn
ssPyf8WA0Yz9NkNCE6GYjhcsy9OvbRCl1P5mhytmh6WAJk1ec6Fo/ijKJPyc
cX5boNyK9GyXYUqriV6MfQaHxuNvOe0VoF0PiSEXmvDUSBgoA66xBB3J79r5
tCISe0zCE68q80GepwfMZI+xIoHf8tN6SLxiEsiZzK+TlP/ONoMcgyzv+emX
JjAjQxe9Au79ofZOz3w0K6t4k7tmZdWpcKMsvAdcFfILNbh4Gs/DM85gWT7c
061ZZbWxQ6unVMnIcxCRlfHPrppiOhc0WbGMMcX+LWU9CCrbdzD4nfkO3nUv
cDU4AvHeg7BRKsycKRDmd8H5uVwX8GZ9uy+IasJhK/72Ux2hJq9shRxxCjRv
kPJ9kX39BfmiZVIZAHwAT9HkFnrx8zZBjCCRFo2eEnKJHipQUC4ouZjtWrmg
dxBxQ973Wi+IqV0uBSffUIMnBqfVpLE4n+ImUMcbueKgsV7OtZEbUCRS6T1W
ADv5lD2BnAkGoEDEGw1KTkFYibyRyopN4p5o8zhnhMRtvbr7yhMQVog0frbP
FWwlGnGkR3r+BD373W48uexHW11dshbG3u93XIWonRKXcqzET4iU20smwzA9
I5XvmIp7zsSlobiakC/KIatERh8TO6O4uSiDJ+tINxEWbXzCnJBM3a41ubUO
tYxRR/Ew/IDrCZFwOx64Ebj9WCoiloA5sBFz+6y4EWjTluozlCBGNYaALHPG
CcMbJhZdxwjjHjExlPMytygXMa0al0AsaQCb+J9YFVHOGt8680EI3B4huYkh
u3pqI0Y1hkhQ10N9yXEgrkQfgcbg7lX53IzQCOqUUcV6Udnti0gljYQa6Fm7
Y0TqJKjAaFhYXcdcXcgqkFlByeMpkQcegiACYqRl7JHRyOvYQUr281ZPCtVe
Cc7ayB5Ddri8bshU3WWqaA+DdrueeyEuN200JkZWjnX7PdFGgw7tloeGBYAG
x6hvRYNOUXFE1j38fAI5203dCgWAYKI+2CoAWUDLKQDt3kmEClkQGQWA6g7J
GgnglDXAugJKb/4UpG2YNmqqUDuKZplOud7R6zNh2jOQVrvdFL2EyBCukfZn
rm1EO9mDNCEAq0mC+mCkRNaTTKi2UW8hc/5o00oRsrLsntFUKUJ2ymSRRZC0
7IrSRuNuxLKE7eBIdZAAZix1+1x7apq2hTLQDunQqsqMlQCZ8kMsA8StYRmI
a4SWuZxphJhvcn7IaEGNnXaHPrI5EeaMApCpWdTuWwVAdTsKajCaIqhqWK8t
pI/GSqjrQoI4Z02guqZtHO4ZI0mrPb0ADbBQ413CATmo67LSJvez5MCEw6yu
Uy0UBTIe3R2bpL3gV0FW5JOKbz7IaK2SttXIOeSn1oV5VNujeM2KQVM2MXHQ
y2jQoOYBtcCIIvJM9xu1+IFfiFgLFVrOX8smgHh3HdYDeovHXALqqRumf8Lg
EtVnPdsvxW551ZsPqnrtKS4VcZcTFZFf0DRRzcFH9kB1k/SU8SPdnrJIqGGm
yvXbPmy3IE1FBGVGk6X2JPqXiIxMqzwsE/Gk50lEpnlHcSyPWCbFSH9kEaIC
9aQLnjODZRHR4lTPYG0OSSkmZcnniHNDXjRLBYJvVtKeJ4QcliYslQmrhLiZ
IRudwv6Vg0LMdzd9ixGxIGlG78HRJCpMXQ6otF51k4cw1xm4IUkLHGnyJtdU
IjJp1xZXZmlvdUGtTd2glPc1Np+ECY26Nip1LPFKIL8FpJ2xSSvnrGHpy4bM
chxIOZqzp++JZQnfujDXjCGS1PoQYYdNkRy2Z8rqto48rAZ3w6SzjzeQ1S1V
37kRxxGDQx2qdsQpDXviQd1DTPfkAeo9FYyT7mB+RVt0Bm1WKId77vmFVZwl
UtAuPE84ElnHB54Lx9vewkkV8c1Gi+Hu2WGiuIVqvaVkZfsrykyzJpDq1RJ7
Nt20Gl4TQOuAmQRB0RAJ/0XABivq2BoWRL9A2C8U25R/dbwFrPZOYUcYt0dI
gjZqCBvxDanO/6fVbk+oX0erw9m9vdy7C1stgDXw8lt778619+9dew/vNfXs
7/J+jfXaWdPlPVAKequ3aW1l9xeu1Y7QldjDu7K7rld2n3w3Tzboj55AjFpy
LggktYLnt6zlOTz6k96RNUhBXTtvaV3P11qfs9VyH4WpP5ZPvJCsh/FZzwOV
RWyrL5N1IJfw40hRz3Y4porvJI4507GLLhe7HVCW1bqm8tE9E7mPyPw9CDOl
lHRYZiLeOec8yFNRceOtm6eYZjshczljNB9BaqI9pDafci6bxrHy1p2AuZx3
zlOI5Q2mh/aTw4cxBDEQuYYlTh4jm4paToDeiVNBBKo8AVqO3kbKj8wCGWgh
4wfZe1FxpwdapdLRHof+np8qjkP3k8Y3dlSJX9Q5MhqK57efxuXJA3Tq+SU/
qPF0ejEwxeTKeWMABSGHGj2JllzbQJEPOVnl1zYk7tOYiK6nHLPomfY2DWqv
6OpE+fUaxKtIt0SjTriCUF+BoL9m6bsv+CwVcOdujunc3Tydu/2oc/dLde4G
L/cdaf2Nc9mk2dz9fRCwfq6vPZckdvEmypW98bMbN6t27/7aFt8RvHxRJPwM
m8Dc5furm7v1e0l3Zv46cJf6csBkw7T1tnq0CnHmlLDb3aARwOfO16utgGnw
i2Zsl7ngVAzY6kdEbSuahl3mDNUUYEVRVUK3+pNoNxGiAjQedTbAuYYVTQIu
QTmQ4+lMWZTbT6Vy6Ghucrr60N2uNJVDR6v8aK9PjdDZk2Ae9yPw/LNqtEfV
FGGm3QeRkZmsxjIb1aqk+UMuKWiPAioeLQq406roRJfTidA8uv7/R8pmEFSd
vyL/DxC2A8L+b76jTXfI3niXzOcL/pdGcqgy3Woq8bRb4fKquDmIiI+R4uOY
tMzQC+hu9iSe9hCWXCqe+6S0Ioy0hxBU1JKiPsjd80s87b75iqAjo8VtFko8
7SbQiqCjeNoFM+0+xoqgo4GPe6Qm8bQ7ECuC7vZVdxel3URX6hhEHLftMLFT
EpqTLXW0oV0QltlMLwoQrYFUVOsHLEQrPsrEXXypkJBP5BtlUTXCRM5WyA9B
lAh53VcETg4cRZCu+Sl6cWuU4TED2QNP0VJDw4CRozzypLlgYtqxQEWA5TAM
5B59xjCRoVUWTGsLLn5V1tM26Smy2LXgcrfh4rlmPd14ENfcRiqY3mWaiCU7
gzo8bdqKIGFPAy9wiLZ6TDWO8OV716REG7kaBtdwX9dDaiMka7taNpDcbWhf
ST2kHlIPqVZIfbPYgUrqIfWQSoCElia8x7l2aaJGpPblIetearSC0FQYmsZo
E74+TB+WWOoNl+R5qZfmoKM35ERXavmGHPN9kNBt0WcvsvQZjWq5roWi+Asl
fs5n6d/exq+AIh88\
\>"],ImageRangeCache->{{{0, 287}, {287, 0}} -> {-3.34304, -3.29279, \
0.0223347, 0.0223347}}],

Cell[BoxData[
 TagBox[
  RowBox[{"\[SkeletonIndicator]", "Graphics", "\[SkeletonIndicator]"}],
  False,
  Editable->False]], "Output"]
}, Open  ]]
},
WindowSize->{637, 740},
WindowMargins->{{173, Automatic}, {Automatic, 15}},
PrintingCopies->1,
PrintingPageRange->{Automatic, Automatic},
PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}},
CellLabelAutoDelete->True,
FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (November 20, 2012)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[557, 20, 1171, 30, 172, "Input"],
Cell[CellGroupData[{
Cell[1753, 54, 4520, 150, 612, "Input"],
Cell[6276, 206, 220, 7, 48, "Output"]
}, Open  ]],
Cell[6511, 216, 1113, 31, 112, "Input"],
Cell[7627, 249, 203, 6, 52, "Input"],
Cell[CellGroupData[{
Cell[7855, 259, 346, 10, 52, "Input"],
Cell[8204, 271, 321, 9, 31, "Output"]
}, Open  ]],
Cell[8540, 283, 121, 4, 31, "Input"],
Cell[CellGroupData[{
Cell[8686, 291, 57, 1, 31, "Input"],
Cell[8746, 294, 63, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[8846, 300, 94, 2, 31, "Input"],
Cell[8943, 304, 45, 0, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[9025, 309, 62, 1, 31, "Input"],
Cell[9090, 312, 46, 0, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[9173, 317, 208, 7, 31, "Input"],
Cell[9384, 326, 45, 0, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[9466, 331, 152, 5, 31, "Input"],
Cell[9621, 338, 46, 0, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[9704, 343, 4797, 142, 272, "Input"],
Cell[14504, 487, 134, 4, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[14675, 496, 633, 19, 52, "Input"],
Cell[15311, 517, 10816, 622, 296, 5621, 534, "GraphicsData", "PostScript", \
"Graphics"],
Cell[26130, 1141, 134, 4, 31, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)