(* 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[     29073,       1222]
NotebookOptionsPosition[     27489,       1166]
NotebookOutlinePosition[     28002,       1185]
CellTagsIndexPosition[     27959,       1182]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"CH4", " ", "molecule", " ", "with", " ", 
     RowBox[{"C", "[", 
      RowBox[{"1", "s2", " ", "4", "sp"}], "]"}]}], ",", " ", "overlap", ",", 
    " ", 
    RowBox[{"parametrized", " ", "ES", " ", "06.11", ".2012"}]}], " ", "*)"}],
   "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"Variables", ":", " ", "P"}], ",", "Q", ",", "R", ",", "S"}], " ",
    "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{"with", " ", 
    StyleBox["exchange",
     FontColor->RGBColor[1, 0, 0]]}], " ", "*)"}], "\[IndentingNewLine]", 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Clear", "[", 
    RowBox[{
    "k1", ",", "k2", ",", "s1", ",", "s2", ",", " ", "Z", ",", " ", "P", ",", 
     "Q", ",", "R", ",", "S", ",", "Vne", ",", "Vee", ",", "Vnn", ",", "Veec",
      ",", "Vnec", ",", "Veecb"}], "]"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Clear", "[", 
     RowBox[{
     "Vnecb", ",", "Vnncb", ",", "Veebb", ",", "Vnebb", ",", "Vnebb", ",", 
      "d2", ",", "d1", ",", "d3"}], "]"}], ";"}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"Z", "=", "6.0"}], ";", " ", 
    RowBox[{"w3", "=", 
     RowBox[{"Sqrt", "[", "3", "]"}]}], ";", " ", 
    RowBox[{"wd", "=", "0.94281"}], ";", " ", 
    RowBox[{"R", "=", "1.046"}], ";", " ", 
    RowBox[{"S", "=", " ", "1.012"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", " ", "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[{"Veec", "=", 
     RowBox[{"3.0", "*", 
      RowBox[{"s1", "/", "P"}]}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", " ", "bond", " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Tb", "=", 
     RowBox[{"9.0", "*", 
      RowBox[{"k2", "/", 
       RowBox[{"Q", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Vneb", "=", 
     RowBox[{
      RowBox[{"-", 
       RowBox[{"(", 
        RowBox[{"2", "/", "Q"}], ")"}]}], "*", 
      RowBox[{"(", 
       RowBox[{"1.5", "-", 
        RowBox[{"0.5", "*", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"S", "/", "Q"}], ")"}], "^", "2"}]}]}], ")"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Veeb", "=", 
     RowBox[{"3.0", "*", 
      RowBox[{"s2", "/", "Q"}]}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{"core", "-", "bond"}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Vnecb", "=", 
     RowBox[{
      RowBox[{
       RowBox[{"-", 
        RowBox[{"(", 
         RowBox[{"2.0", "*", 
          RowBox[{"Z", "/", "Q"}]}], ")"}]}], "*", 
       RowBox[{"(", 
        RowBox[{"1.5", "-", 
         RowBox[{"0.5", "*", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"R", "/", "Q"}], ")"}], "^", "2"}]}]}], ")"}]}], "-", 
      RowBox[{"2.0", "/", 
       RowBox[{"(", 
        RowBox[{"R", "+", "S"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Veecb", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"2", "/", "Q"}], ")"}], "*", 
      RowBox[{"(", 
       RowBox[{"3.0", "-", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"R", "/", "Q"}], ")"}], "^", "2"}], "-", 
        RowBox[{"0.6", "*", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"P", "/", "Q"}], ")"}], "^", "2"}]}]}], ")"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Vex", "  ", "=", " ", 
     RowBox[{
      RowBox[{"-", 
       StyleBox["1.2",
        FontColor->RGBColor[1, 0, 0]]}], 
      StyleBox["*",
       FontColor->RGBColor[1, 0, 0]], 
      StyleBox[
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{"P", "/", "Q"}], ")"}], "^", "3"}], "/", "Q"}],
       FontColor->RGBColor[1, 0, 0]]}]}], 
    StyleBox[";",
     FontColor->RGBColor[1, 0, 0]]}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Veecb", " ", "=", " ", 
     RowBox[{"Veecb", " ", "+", " ", "Vex"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Vnncb", "=", 
     RowBox[{"Z", "/", 
      RowBox[{"(", 
       RowBox[{"R", "+", "S"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{"bond", "-", "bond"}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"d1", "=", " ", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"R", "+", "S"}], ")"}], "*", "ad"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"d2", "=", " ", 
     RowBox[{"Sqrt", "[", 
      RowBox[{
       RowBox[{"R", "^", "2"}], "+", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"R", "+", "S"}], ")"}], "^", "2"}], "+", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"2", "/", "3"}], ")"}], "*", "R", "*", 
        RowBox[{"(", 
         RowBox[{"R", "+", "S"}], ")"}]}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"d3", "=", " ", 
     RowBox[{"R", "*", "ad"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Vnebb", " ", "=", " ", 
     RowBox[{
      RowBox[{"-", "4"}], "/", "d2"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Veebb", " ", "=", " ", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"4", "/", "Q"}], ")"}], "*", 
      RowBox[{"(", 
       RowBox[{"1.2", "-", 
        RowBox[{"0.5", "*", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"d3", "/", "Q"}], ")"}], "^", "2"}]}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"3", "/", "16"}], ")"}], "*", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"d3", "/", "Q"}], ")"}], "^", "3"}]}], "-", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"1", "/", "160"}], ")"}], "*", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"d3", "/", "Q"}], ")"}], "^", "5"}]}]}], ")"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Vnnbb", " ", "=", "  ", 
     RowBox[{"1", "/", "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[{
    RowBox[{"Vee", " ", "=", " ", 
     RowBox[{"Veec", "+", 
      RowBox[{"4", "*", "Veeb"}], "+", 
      RowBox[{"4", "*", "Veecb"}], "+", 
      RowBox[{"6", "*", "Veebb"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Vnn", " ", "=", " ", 
     RowBox[{
      RowBox[{"4", "*", "Vnncb"}], "+", 
      RowBox[{"6", "*", "Vnnbb"}]}]}], ";"}]}], 
  "\[IndentingNewLine]"}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"c", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"k1", " ", "->", " ", "1.0"}], ",", " ", 
       RowBox[{"k2", "->", " ", "0.625"}], ",", " ", 
       RowBox[{"s1", " ", "->", " ", "0.3"}], ",", " ", 
       RowBox[{"s2", " ", "->", " ", "0.3"}]}], "}"}]}], ";"}], "  ", "*)"}], 
  RowBox[{
   RowBox[{"c", "=", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"k1", " ", "->", " ", "1.08"}], ",", " ", 
      RowBox[{"k2", " ", "->", " ", "0.54"}], ",", " ", 
      RowBox[{"s1", " ", "->", " ", "0.25742249"}], ",", 
      RowBox[{"s2", " ", "->", " ", "0.22958899"}]}], "}"}]}], 
   ";"}]}]], "Input"],

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

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"t", "=", 
  RowBox[{"FindMinimum", "[", 
   RowBox[{"func", ",", 
    RowBox[{"{", 
     RowBox[{"P", ",", "0.28"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"Q", ",", "1.607"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"Method", " ", "\[Rule]", " ", "Automatic"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"MaxIterations", " ", "\[Rule]", " ", "500"}], "}"}]}], "]"}], 
  " "}]], "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "40.478782207821645`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"P", "\[Rule]", "0.28319422692484925`"}], ",", 
     RowBox[{"Q", "\[Rule]", "1.6074562617432329`"}]}], "}"}]}], 
  "}"}]], "Output"]
}, Open  ]],

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

Cell[CellGroupData[{

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

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

Cell[CellGroupData[{

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

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

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"4", "*", "Vex"}], " ", "/.", " ", "u"}]], "Input"],

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

Cell[CellGroupData[{

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

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

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"-", 
     RowBox[{"(", 
      RowBox[{"Vne", "+", "Vee", "+", "Vnn"}], ")"}]}], "/", "T"}], " ", "/.",
    "u"}], "/.", "c"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"(", 
    RowBox[{"R", "+", "S"}], ")"}], "*", "0.529177"}], " ", "/.", 
  "u"}]}], "Input"],

Cell[BoxData["2.0702130539999075`"], "Output"],

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

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"plot1", "=", 
  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[{"R", "/", "2"}]}], ",", 
          RowBox[{"wd", "*", "w3", "*", 
           RowBox[{"R", "/", "2"}]}]}], "}"}], ",", "Q"}], "]"}], ",", 
      RowBox[{"Circle", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{
           RowBox[{"-", "wd"}], "*", 
           RowBox[{"R", "/", "2"}]}], ",", 
          RowBox[{
           RowBox[{"-", "wd"}], "*", "w3", "*", 
           RowBox[{"R", "/", "2"}]}]}], "}"}], ",", "Q"}], "]"}], ",", 
      RowBox[{"Circle", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{"wd", "*", "R"}], ",", "0"}], "}"}], ",", "Q"}], "]"}], ",", 
      RowBox[{"Disk", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{
           RowBox[{"-", "wd"}], "*", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"R", "+", "S"}], ")"}], "/", "2"}]}], ",", 
          RowBox[{"wd", "*", "w3", "*", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"R", "+", "S"}], ")"}], "/", "2"}]}]}], "}"}], ",", 
        "0.1"}], "]"}], ",", 
      RowBox[{"Disk", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{
           RowBox[{"-", "wd"}], "*", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"R", "+", "S"}], ")"}], "/", "2"}]}], ",", 
          RowBox[{
           RowBox[{"-", "wd"}], "*", "w3", "*", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"R", "+", "S"}], ")"}], "/", "2"}]}]}], "}"}], ",", 
        "0.1"}], "]"}], ",", 
      RowBox[{"Disk", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{"wd", "*", 
           RowBox[{"(", 
            RowBox[{"R", "+", "S"}], ")"}]}], ",", "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[{"R", "+", "S"}], ")"}]}], ",", "0"}], "}"}]}], "}"}], 
         "]"}], ",", 
        RowBox[{"Line", "[", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"0", ",", "0"}], "}"}], ",", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{
              RowBox[{"-", "wd"}], "*", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{"R", "+", "S"}], ")"}], "/", "2"}]}], ",", 
             RowBox[{"wd", "*", "w3", "*", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{"R", "+", "S"}], ")"}], "/", "2"}]}]}], "}"}]}], 
          "}"}], "]"}], ",", 
        RowBox[{"Line", "[", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"0", ",", "0"}], "}"}], ",", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{
              RowBox[{"-", "wd"}], "*", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{"R", "+", "S"}], ")"}], "/", "2"}]}], ",", 
             RowBox[{
              RowBox[{"-", "wd"}], "*", "w3", "*", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{"R", "+", "S"}], ")"}], "/", "2"}]}]}], "}"}]}], 
          "}"}], "]"}]}], "}"}]}], "}"}], "]"}], "/.", " ", "u"}]}]], "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 .0472 0 365.73 arc
s
.5 .5 m
.5 .5 .01667 0 365.73 arc
F
newpath
.5 .5 .26791 0 365.73 arc
s
newpath
.41782 .64234 .26791 0 365.73 arc
s
newpath
.41782 .35766 .26791 0 365.73 arc
s
newpath
.66436 .5 .26791 0 365.73 arc
s
.33831 .78006 m
.33831 .78006 .01667 0 365.73 arc
F
.33831 .21994 m
.33831 .21994 .01667 0 365.73 arc
F
.82338 .5 m
.82338 .5 .01667 0 365.73 arc
F
.012 w
.5 .5 m
.82338 .5 L
s
.5 .5 m
.33831 .78006 L
s
.5 .5 m
.33831 .21994 L
s
% End of Graphics
MathPictureEnd
\
\>"], "Graphics",
 ImageSize->{288, 288},
 ImageMargins->{{43, 0}, {0, 0}},
 ImageRegion->{{0, 1}, {0, 
  1}},ImageCache->GraphicsData["CompressedBitmap", "\<\
eJztndtSHMkRhhuG4SCQBEJIQisJ0AoB0kb4zhc+3XnNjZ/BVmysY/fCYcea
e70I78Kj6FEwnUNn9uGrInv6MDPQHQI11dVVf2ZlV2VlZlX9/cvFLz//+8vF
rz99Ofzbb1/++8uvP/3v8Mf//HaTNFpKkqXD9OcwSe+vkyT7dX3z7/oo/ZWm
2P2f5L/V5PrqIileF8nV9fXH7PH1VVLMcJ5cfrs+lscrybfLc3tw+e00/tZJ
vM5TKDS5uFqeAP6z/DdKJg/lheTTbVpa4E318tcE+Ti5vLiUu/Tn8nyS/wzy
n0G5ZxnQb5fJ+XlWTkqT0LE0QfSXXGkTqFzDcZ4r54bq6kJKo1c+QdFZWsa8
NO2WN78vZU+vZfm9KtmyGjM6rq5uxYBeGRdfuUX1sT4Rn7O0SY1FwnLFLJUk
ssJ0aZQTbaebP29rOcna+lxA5LMnp5A/3K65KjLxrRS5lBfEgOwbQw3fyqQ+
FegCs7Ksn7jUMy5hufpBN7l/mv4n/cXT9Ncol3rzR3plzb06R2njOcIypDVL
WwqL3CdN+6ppZ5D2yfluPG2ANEAaIM0ekugin3sAIuO2vBoZ9uaBIwOkAdIA
6YFDGrrFBWikAdLDgrSZ/t5N39yaP3Ay2xqnb+7q3VqHMO/ooVYKaRmgt+nd
apLrbm5Lp4La765SULkivJW20GKnuXxZ+mdlVKK8uLniH0OrqD4oU7Q5GNrt
j5T4sT9830O+kbOqFRXSD/0BPtJKV50VkJHkqD/A77RZ150VrKnMHPQH843y
ddNZwSakvekP8HdRSPGqtpTYtgAH++p9ffrEWewTSNt3wpy6836pMLedMI1s
L7gWGv0F5PMC3oG0F/1Bf6782nVW8Ezf6BHmLuTzAn5e+91WoT+LQopXtRfi
fz/QSTb3nFXRN7HTHPpk2Ay6mEdOyrajiONIXkJa/Gt3UZbp7Hk/NXnAqFhq
0JfO6l9B2tPm5GRqRd7DLaqZ5P6jFGNqcfrzV0jbh7R/pb/zo0lSyzG+4uTp
Y8gXH1Ms32tIexx9t57ws3vbK/w0iX3tREL6zVZzyibCX/GOe9uKlL7vnAje
QFpcY2u1j30UhRSv6i2kPeoP+kYUUryqd5C20R/09SikeFUHkBafRA3QJ2lk
ljpwVnUIaWstQQ9Ok2iSfuis4AjS4sYBS5t6wtQE8HtIGzsBtyAcpOsIdDEg
nsFTshKt9AeYBidpdLFznsLTD5AWN2t1Dlg4vJ6mnMwfYFJoBPBGCPAxpC33
B3gZ8r2PcngAfBfgcaEIyXxzHUdBftQ3k9KbX5MfpgCWGyCK/o3k1iwuSEQ0
JcBStKBtxTl9716pTd5dS+9e6R0Nq2N9+lTvxlpkkivy65pm8KoRNYQsvQKO
IXlRLCECiQyMW/p0Xd8oz8vSUp/oU6/W3w4JAm1bKydTjz1NzPOVXiWzphGz
oyV7p/BTEQMdgfFX7sj6JOBWy8RmJT0vtmt70LP0IHTh1yP9VoNGjZKPhkwe
Ro6Zn+NGt1bJGSkP1/SO7Alb+nRDGfAYiNjXfHGrTk0iLC1IhEAyd1xwtopf
QvoTnNYvF4mdmpy8VpkVE9Sciv7SyXUAhdK0eEdRE0315/GRYeqTk6ax0kQz
lTUt74Xe5efnWd53+pQAu0Y6byPEAY8VCH0sqwC9/oyyVZkSjqyGyLH++K3e
jYGII33qdTxHoJfnLGlRpPVknLavnHgpWWiOVlYX06v+TKNVijL+jYiY9Pp4
Kx75wrMCszz1VeOGNBRZWHLQC/uzBXv58nLv+PRkF8psKVSppjzQ7NG6AiUr
4KmCL9F3kn+QJ+Rzy0hLXbl9hVL5D/h6JNCnKTDq0B8pC8sG+SLEr0nuOlMa
7gKWk50KnLIlJq3EdJHXzJ08olxyWXXQau8aMwKsMhGh786YZopd2GfC885T
KHe5+MDV1rEv7yPkJ3o2lVPkbA1Txp35CVJWLYN6uSkJpaKo09/UhiNffthJ
tIZNmP1P8WzENBqapiSYiiKLpc17KfAhTDDLLLXi91Brx8STKkD2ZZuA0nw6
7BJn4omo91ArMZUQT0k8GddJY3usLU8TIpq4xogPmqhJtQ4yqUK8+qF5GwBi
eRu8MSsPcd8b87wcRJLgV1GTN/m4h2yhPDWDV1k8iJJaYk7FCJA9Jaeh1x5M
aYfwbjyi1c8oyX0sRZvSM+lAqmmv09+StAVPNyDtd5AmhNzOnrO7Ay05FmbC
mymQzHgbnIwB5Bl/qg1OIQSU5m1wQk9+fTJ0kEZfCXCyDRuaMIoCN8gIsa2M
oq+A7DFeRpHEU+wGNQY1c4lR+e0j+mDUjn59C8IoGga8TPEGaps1nojwrgci
ptCARe6LeER7q0wJ2plL+Z6ppJBc0gDkZQpJHsWmEeOpSftmitm87zdT9iEf
GeR3VVJIKySfn5cpNPpQfCYxnpq0BaaQ14nWV+yqpCTmTkj/J51/Hvghuo7X
LOHlgszlRQ+nOHMv2TRkvIJ3SdReloSzUeNTWDh53Pa0Upre9kE2fbstSL93
FYzUsZnmaiLuXgYQRSShPTJAunOL3BiVOoL0Eot9k0kbKQ0UTtAje4LhBESO
MEC6qrGyYmx8elBMSa/SEq1EH8hfZi+k6CZiBakKtACGPlHq6bpkRWmxmT2V
v+IRvEQ7jYC0rCm+NM9Fe8zsRUVJGVmHUBksqwuzMXJK5YWdFQvFkMfauuQY
IHeAec5p7FvW8gzJ1xISVkgooolGcBLjii51opVX2UAFmKuHLObkJDDfPLWu
Gdq+74gNZcW/tjQQGzaUgbRqIegrGYXYYFEuRzOShoYMIRWhFmtEQmgwMFYf
RFgTbHlH/fTx1veOUynWrGTRIcyrygvSEcIWmEXiCtmvvPwxRZV4Ud64oOLn
9y7fji/VbsIQKsVCxmh6GmRD0OYenmstCENGyhCathBDLCybGFLejKHCBiKP
Qm37ZQOpSMQQGhcsFjcYBnqH1p1e3nX/8XwuvlB02TPgVdaSJKVE1iOVDGJE
eBOJrsmXSimGiYguYHHQG/QyefagKHIjPkcmym2lQUTjCjb3iBqjKPmlDl+I
JSfGXlBYmpLYmmzHP8EQ2RaZQmTftS9Elu5lQDyfiwEUi7ytrUPmsiDZQRfY
XftHLCzZFm7ShGzvbiEdkR3fD4lKsbgC6UbGxXfj2k6PxFIQtemhSdGv4jUM
RxZrVTObMJV0yLnhx0RX51AnrzhY9AR5Nc2cQt9DXifM6qF8xJI462qyxD5H
UlOpcb1m4W39RsgwbYufiJy8YtQHeypG0q6YQmk5xaT4gDpXUpeIW9Qvkc+N
uGXdu3OPhSyd9FkizMutJ7eikk/Lvjlv3EjQClLSfsjR62VYg69vUj2rxtSt
eFlHNj1bP+llHmmaZamcU+aR3HkHuyDzKDKBhJ6ELrDGvHfW8fIuh+g14Z4Z
mx3cI6lbmRPu+aPJKf7LG4dALJTyxPlLgxLxLBzPzV+suVK0kTxD70JwT4Rc
Iguos3iB0oULUIgXG9o0zcSu4V5vy0wFEdxYEGVbD/p4ywp7HT6ua8mkvtzJ
OD750BsWYSZeIqsxw/4xBVmF1STVkwcFsXA4f3ad56XZs+OfaUrNfkVbGc48
rBf80vk3Yt1N29/IhpY81TdSOU+x7pK53jrrO4a6mXfW8zjKda8jmDG/Kx2h
c/2UfPrCONJPiY/hpbaz10+HiVGHrOt3Qu71ls3hhHwwBDHrMBgzvQY7o0e8
svR+TdU0SAdD8ErsmYlN39K6d3NQGjFnbhweJT2gJZ/YmhaQXfkCpQFst7AZ
+8PadoaabY+6zjkk9oE6vB9AeMODD2u530FMQ4RahXKO0oS2idIuSfMVjVh6
NOdBqgsSvzuEMctPfGuGOEPuUaB70BgT3ES4/gKI8Jr6BWEIjYC0RcyUa2dy
+z/PfMVM3SVWtpiIGFJ/YVV4j5uFYw2NHrQ3Wf2FeevK9MMIa+ZumaapGbR/
Xf1lmmawf98RGxov0yR90GwO5K6qv2jXBvTjGUlDQ4bYzDop+pPIy3cdWNat
f3l3au1xWXeTNf5meqOnJhSSxbvVw4LQTjMnyVEW0onoJAmFL01yeHc/WKSd
HyqcKXwbtjNI0HC7wDuDeHeTkQFXDsSSMWdLmXJ/WOHdWmlbWdHHjlJzyID5
31tq2EnsbrJr7a0X9A7fs731vLswmkL7AHZhHLamhKq8O7uax+AB7Oxaaw/g
oFPknu0BfG83Rp5ig/62tye3WR0NJnO4PXnbG/dTVbQ1k4UYkCAREV5GEXo6
4pOE+mAKRknuGWzcX0o7g7R3Wtuom9Ymft2LYxq8u/Y/mAM+vLvzP5jTYOoa
SB/0OUBynpJMgonkac99WkWSrwtG82JNHZ/1VOsIs4jtN0Rwr0eYxQ6bC9Bv
skumG8J096Fzdx1I14SWwFl/Jo+vgsSG9hxO8icIfs2/UI5tKkNe1NMHZ3BQ
ow/Y3B5tWd6MJL1Kp4NaNFyXZ4LWPSfWjlktteiE0D4OXK17WK+p6IdQmh3W
i5stp9eHYlGlD76DQ20nf1XfjZw4XE0rN1KxUAtfpTG+BTprnTAuckG6uDnA
aVa5AuSQvkYaR0fkOM7qLqXtK4lrQM6hPqUpdI/noq9qeeScIOgkgESEa+Qh
SsvDmIoELVeMo7G0XaUzaHoKbrBcU6a8NJn3N/452M7VOwDTGp8sWa2SYzIQ
JEdetZPHg6brkRIhd4+BnH1lD7nIOiJiWSs1rgc3C7SvuxSca0R4F0t1SY68
ajpl0A8/SVrRz+v2egZkPVc2UeReA2KyuifEVAZ+qdIMjDSPlm9E9JWy7zQt
y3wkFALQHXB50dy6ZCvbVOCrmm9TSytoX3aODLmmpqKiohbSWRsr+jlYqAqN
J+a6His9VSrWlOCDKajITScq2EVpXdOSd9K77/SOzAc1Biuetkqh61p56U2f
mtyaIF5fs+74PgrSaxqm2joCLBrLxpwCJlX+aJ45TPa2wyhgsqDRuTV9A5Zv
myxFMwZMEwERCUmiCupP01oFTCb3Q2dV9adkU/Xu7QEOzn8cgGvOZuxF8vJ4
xzuyCZATriPhoPGd/K1U1QGkxae6rUInNzP51Kmq4BysH+g0P6UQCqqKpoAU
APJwoGvkyFUxbiSZ7EDn3YCOQhkoAoqQUFROfF7uoswc2UpNxIlYLJaCwbxT
CZrTU4xeTXLyAQxC0FU+CuMPUszotukmMRA/QtprSPs5/Z3vurOPOY0murq6
1r88nKMoVAqqJUrjpw40F/F0i0DdOlI3WfRSRh+5d5J816YSjUR86p1VKcaA
VoYQAu+uHx31pHROaLyC8P4OHcL0Lv6lqjwbNHQI3bvGnaqKL03tHLp3eRFV
FV8s2Dl07yIgqspz2Ni00IPzjZf61AszvilTq3ML78IKghlf7NK5ILxWmuO6
Q0llkqS4StkqTHMMxXVR0lm9M49WAR9AvviUsTQ/yy1S7wXwoVbqhWkuS4p3
6wjme5U+rxPd/GEf+oN5rLwJ2xX5HESKNesI5IlyRmzfgoe8V6125PkKsssW
7xmaLvvruIvuQJkRjqfrpX0IpswsxYGxp3fheN2ZwZTedy9907uKpUdwVOwA
aYA0b5BcfWgbQFrtSB9YIw2QBkgDpKFbHBppgDRAahsSeR6Siyuv56FHpGHv
T3BpPM3bZpVGHst5wjek+dNKntzM466eXImUzB96l3fE6qF36e+9kmybPCc5
z2Y+a+CcNMkycRJMSj6q3idL/wfYZioH\
\>"],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->{723, 740},
WindowMargins->{{Automatic, 536}, {Automatic, 48}},
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, 1120, 29, 132, "Input"],
Cell[1680, 51, 6139, 189, 692, "Input"],
Cell[7822, 242, 684, 19, 52, "Input"],
Cell[8509, 263, 161, 5, 31, "Input"],
Cell[CellGroupData[{
Cell[8695, 272, 422, 12, 52, "Input"],
Cell[9120, 286, 263, 8, 31, "Output"]
}, Open  ]],
Cell[9398, 297, 121, 4, 31, "Input"],
Cell[CellGroupData[{
Cell[9544, 305, 57, 1, 31, "Input"],
Cell[9604, 308, 63, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[9704, 314, 92, 2, 31, "Input"],
Cell[9799, 318, 46, 0, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[9882, 323, 85, 2, 31, "Input"],
Cell[9970, 327, 65, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[10072, 333, 62, 1, 31, "Input"],
Cell[10137, 336, 46, 0, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[10220, 341, 336, 12, 52, "Input"],
Cell[10559, 355, 46, 0, 31, "Output"],
Cell[10608, 357, 39, 0, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[10684, 362, 4157, 127, 172, "Input"],
Cell[14844, 491, 134, 4, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[15015, 500, 633, 19, 52, "Input"],
Cell[15651, 521, 11685, 636, 296, 5620, 534, "GraphicsData", "PostScript", \
"Graphics"],
Cell[27339, 1159, 134, 4, 31, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)