(* 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[     18449,        831]
NotebookOptionsPosition[     16806,        774]
NotebookOutlinePosition[     17798,        810]
CellTagsIndexPosition[     17755,        807]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["\<\
(* Fe+++ ion  + 6 F- in Oh 08.11.12 *)
Clear[k1,k2,k3,k4,k5,k6,sig1,sig2,sig3,sig4,sig5,c,z,R1,R2,R3,R4,R5,S2,S3,S4];
c = {k1 -> 0.965, k2 -> 0.96, k3 -> 0.96, k4 -> 1.26, k5 -> 0.98, k6 -> \
2.108,sig1 -> 0.3, 
sig2 -> 0.302, sig3 -> 0.302, sig4 -> 0.302, sig5 -> 0.302}; z=26.; z1=9.0;

(* He+Ne shell *)
T = 2.25*k1/R1^2+9.*k2/R2^2 /. c;
ad = Sqrt[3./8.]; sq2=Sqrt[2];
Vee=3.0*sig1/R1+12.*sig2/R2+16/(R1+R2)+24*ad/(R1+R2) /. c;
Vne=-3.0*z/R1-8.*z/(R1+R2);
S2 = R2*4^(1/3);

(* Ar shell *)
T = T + 9.*k3/R3^2 /. c;
Vee = Vee+12.*sig3/R3+80./(S2+R3)+24.*ad/(S2+R3) /. c;
Vne = Vne-8.45*z/(S2+R3);
S3 = R3*4^(1/3);

(* Fe+++ shell octahedron, full with 5/6 charge each*)
T = T + (5.*9./8)*k4/R4^2 /. c;
Vee = Vee+(90.+12*(5/6)^2/sq2+1.5*(5/6)^2)/(S3+R4);
Vne = Vne - 6.0*z/(S3+R4);
S4 = R4*6^(1/3);

(* F- He+Ne shell *)
Tf = 2.25*k5/R5^2+9.*k6/R6^2 /. c;
Veef=3.0*sig4/R5+12*sig5/R6+16/(R5+R6)+24*ad/(R5+R6) /. c;
Vnef=-3.0*z1/R5-7.5*z1/(R5+R6) /. c;
S5 = R6*4^(1/3);

T=T+6*Tf;
Vee=Vee+6*Veef;
Vne=Vne+6*Vnef;

(* Fe+++/6 F- *)
Vnn = (6*z*z1 + 12*z1^2/sq2+1.5*z1^2)/(S4+S5);
Vee = Vee+(6*23*10 + 12*10*10/sq2 + 1.5*10*10)/(S4+S5);
Vne = Vne -(6*z*10+6*z1*23+24*z1*10/sq2+3*z1*10)/(S4+S5);

func = T + Vnn + Vne + Vee;

t = FindMinimum[func, {R1,0.0496435}, {R2,0.155502}, {R3,0.360538}, \
{R4,0.741828},
    {R5,0.166514},{R6,1.4630045},{Method -> Automatic}, {MaxIterations -> 500}]
    
N[Vne /. c /. t[[2]],10]
N[Vnn /. c /. t[[2]],10]
N[Vee /. c /. t[[2]],10]
N[-(Vee+Vne+Vnn)/T /. c /. t[[2]],10]
N[(S2+R3) /. c /. t[[2]],10]
N[(S3+R4) /. c /. t[[2]],10]
N[(S4) /. c /. t[[2]],10]
N[(S5) /. c /. t[[2]],10]
N[0.529177*(S4+S5) /. c /. t[[2]],10]

\
\>", "Input",
 PageWidth->Infinity,
 InitializationCell->True,
 ShowSpecialCharacters->False],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "1849.8230189339763`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"R1", "\[Rule]", "0.04964775057971061`"}], ",", 
     RowBox[{"R2", "\[Rule]", "0.15571785940071023`"}], ",", 
     RowBox[{"R3", "\[Rule]", "0.3624692802825318`"}], ",", 
     RowBox[{"R4", "\[Rule]", "0.7405145740990833`"}], ",", 
     RowBox[{"R5", "\[Rule]", "0.16652471410381428`"}], ",", 
     RowBox[{"R6", "\[Rule]", "1.4592815576785754`"}]}], "}"}]}], 
  "}"}]], "Output"],

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

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

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

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

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

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

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

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

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

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"d", "=", 
    RowBox[{"S4", "+", "S5"}]}], ";"}]}]], "Input"],

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

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"xy", " ", "plane"}], "*)"}], 
  RowBox[{
   RowBox[{
    RowBox[{"plot2", "=", 
     RowBox[{
      RowBox[{
       RowBox[{"Graphics", "[", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Circle", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"0", ",", "0"}], "}"}], ",", "R1"}], "]"}], ",", 
          RowBox[{"Circle", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"0", ",", "0"}], "}"}], ",", "S2"}], "]"}], ",", 
          RowBox[{"Circle", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"0", ",", "0"}], "}"}], ",", "S3"}], "]"}], ",", 
          RowBox[{"Circle", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"0", ",", "0"}], "}"}], ",", "S4"}], "]"}], ",", 
          RowBox[{"Circle", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"d", ",", "0"}], "}"}], ",", "R5"}], "]"}], ",", 
          RowBox[{"Circle", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"d", ",", "0"}], "}"}], ",", "S5"}], "]"}], ",", 
          RowBox[{"Circle", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"-", "d"}], ",", "0"}], "}"}], ",", "R5"}], "]"}], ",", 
          RowBox[{"Circle", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"-", "d"}], ",", "0"}], "}"}], ",", "S5"}], "]"}], ",", 
          RowBox[{"Circle", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"0", ",", "d"}], "}"}], ",", "R5"}], "]"}], ",", 
          RowBox[{"Circle", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"0", ",", "d"}], "}"}], ",", "S5"}], "]"}], ",", 
          RowBox[{"Circle", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"0", ",", 
              RowBox[{"-", "d"}]}], "}"}], ",", "R5"}], "]"}], ",", 
          RowBox[{"Circle", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"0", ",", 
              RowBox[{"-", "d"}]}], "}"}], ",", "S5"}], "]"}]}], "}"}], "]"}],
        "/.", 
       RowBox[{"t", "[", 
        RowBox[{"[", "2", "]"}], "]"}]}], "/.", "c"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"Show", "[", 
    RowBox[{"plot2", ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", 
       RowBox[{"Axes", "\[Rule]", "True"}], ",", 
       RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",", 
       RowBox[{"PlotRange", "\[Rule]", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"-", "6"}], ",", "6"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"-", "6"}], ",", "6"}], "}"}]}], "}"}]}], ",", 
       RowBox[{"Frame", "\[Rule]", "True"}]}], "}"}]}], "]"}]}]}]], "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.0833333 0.5 0.0833333 [
[.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
[(-4)] .16667 -0.0125 0 1 Mshowa
.33333 0 m
.33333 .00625 L
s
[(-2)] .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
[(2)] .66667 -0.0125 0 1 Mshowa
.83333 0 m
.83333 .00625 L
s
[(4)] .83333 -0.0125 0 1 Mshowa
1 0 m
1 .00625 L
s
[(6)] 1 -0.0125 0 1 Mshowa
.125 Mabswid
.04167 0 m
.04167 .00375 L
s
.08333 0 m
.08333 .00375 L
s
.125 0 m
.125 .00375 L
s
.20833 0 m
.20833 .00375 L
s
.25 0 m
.25 .00375 L
s
.29167 0 m
.29167 .00375 L
s
.375 0 m
.375 .00375 L
s
.41667 0 m
.41667 .00375 L
s
.45833 0 m
.45833 .00375 L
s
.54167 0 m
.54167 .00375 L
s
.58333 0 m
.58333 .00375 L
s
.625 0 m
.625 .00375 L
s
.70833 0 m
.70833 .00375 L
s
.75 0 m
.75 .00375 L
s
.79167 0 m
.79167 .00375 L
s
.875 0 m
.875 .00375 L
s
.91667 0 m
.91667 .00375 L
s
.95833 0 m
.95833 .00375 L
s
.25 Mabswid
0 0 m
1 0 L
s
0 .16667 m
.00625 .16667 L
s
[(-4)] -0.0125 .16667 1 0 Mshowa
0 .33333 m
.00625 .33333 L
s
[(-2)] -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
[(2)] -0.0125 .66667 1 0 Mshowa
0 .83333 m
.00625 .83333 L
s
[(4)] -0.0125 .83333 1 0 Mshowa
0 1 m
.00625 1 L
s
[(6)] -0.0125 1 1 0 Mshowa
.125 Mabswid
0 .04167 m
.00375 .04167 L
s
0 .08333 m
.00375 .08333 L
s
0 .125 m
.00375 .125 L
s
0 .20833 m
.00375 .20833 L
s
0 .25 m
.00375 .25 L
s
0 .29167 m
.00375 .29167 L
s
0 .375 m
.00375 .375 L
s
0 .41667 m
.00375 .41667 L
s
0 .45833 m
.00375 .45833 L
s
0 .54167 m
.00375 .54167 L
s
0 .58333 m
.00375 .58333 L
s
0 .625 m
.00375 .625 L
s
0 .70833 m
.00375 .70833 L
s
0 .75 m
.00375 .75 L
s
0 .79167 m
.00375 .79167 L
s
0 .875 m
.00375 .875 L
s
0 .91667 m
.00375 .91667 L
s
0 .95833 m
.00375 .95833 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
.04167 .99625 m
.04167 1 L
s
.08333 .99625 m
.08333 1 L
s
.125 .99625 m
.125 1 L
s
.20833 .99625 m
.20833 1 L
s
.25 .99625 m
.25 1 L
s
.29167 .99625 m
.29167 1 L
s
.375 .99625 m
.375 1 L
s
.41667 .99625 m
.41667 1 L
s
.45833 .99625 m
.45833 1 L
s
.54167 .99625 m
.54167 1 L
s
.58333 .99625 m
.58333 1 L
s
.625 .99625 m
.625 1 L
s
.70833 .99625 m
.70833 1 L
s
.75 .99625 m
.75 1 L
s
.79167 .99625 m
.79167 1 L
s
.875 .99625 m
.875 1 L
s
.91667 .99625 m
.91667 1 L
s
.95833 .99625 m
.95833 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 .04167 m
1 .04167 L
s
.99625 .08333 m
1 .08333 L
s
.99625 .125 m
1 .125 L
s
.99625 .20833 m
1 .20833 L
s
.99625 .25 m
1 .25 L
s
.99625 .29167 m
1 .29167 L
s
.99625 .375 m
1 .375 L
s
.99625 .41667 m
1 .41667 L
s
.99625 .45833 m
1 .45833 L
s
.99625 .54167 m
1 .54167 L
s
.99625 .58333 m
1 .58333 L
s
.99625 .625 m
1 .625 L
s
.99625 .70833 m
1 .70833 L
s
.99625 .75 m
1 .75 L
s
.99625 .79167 m
1 .79167 L
s
.99625 .875 m
1 .875 L
s
.99625 .91667 m
1 .91667 L
s
.99625 .95833 m
1 .95833 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 .00414 0 365.73 arc
s
newpath
.5 .5 .0206 0 365.73 arc
s
newpath
.5 .5 .04795 0 365.73 arc
s
newpath
.5 .5 .11213 0 365.73 arc
s
newpath
.80517 .5 .01388 0 365.73 arc
s
newpath
.80517 .5 .19304 0 365.73 arc
s
newpath
.19483 .5 .01388 0 365.73 arc
s
newpath
.19483 .5 .19304 0 365.73 arc
s
newpath
.5 .80517 .01388 0 365.73 arc
s
newpath
.5 .80517 .19304 0 365.73 arc
s
newpath
.5 .19483 .01388 0 365.73 arc
s
newpath
.5 .19483 .19304 0 365.73 arc
s
% End of Graphics
MathPictureEnd
\
\>"], "Graphics",
 ImageSize->{288, 288},
 ImageMargins->{{43, 0}, {0, 0}},
 ImageRegion->{{0, 1}, {0, 
  1}},ImageCache->GraphicsData["CompressedBitmap", "\<\
eJztXc1y2zgSpi3Ldmwnjp3EsZN4LGcm8c/uYWYPe9qf264vWzOvkJraqpnD
1G7N6u4X0bvoUfIoWgEku0nwQ6spgn8yWZEiN4jG1wAINIDu5r++TH/5929f
pr/+/GXyz9+//PeXX3/+3+Qf//l9SRptRdHWxHwmkfm9iKL0a7H8t7g2X4bC
v/9m/9uJvs4eovh6iGZfF58teTdazKdR/ppG88XiUyHXMt/s662c6w6XFeca
L7/nEeezSdFWjPLv9r+RZbwsx/51k9BiZraAhDaOZtOZ/WU+s4c47Rbc/0fA
N4VjZJjOFwmfeTSfE5y/pHLmEDsCjagaDJdP2SwPjG4+tVnuQJZbuZQbJznF
b2p3mT2B+meHsbm27fcuSUeSxrVwrcjCH1sxK6S7Tzkm1cjVvOMWnvYR02oP
1E1m27E0f01SY5mXvGc2E1XF19ky04yQTbO9ddkrHuJOlUmKVmS9xQXmmivJ
Nco+UKOYYdIsKRfncfmcMk9uTRq++JjYu+8B43twK7e975kv+/vY/GeHkWPz
NcpQl3+YK23v/Qq0Z4H5VSljz3yPDWUPpI4p9aBDmEOXseVv5nuiPRLtDtA+
mG87dBwBLrbyLOkS5EVl3NcD84KAvABcnpMQ79uF+ZaAnAAuxyTEebsw3xDM
V4DLKaWetQvzlOrrDeDymmC+bhfmSwJyDricUeppuzBfEJB3gMs5pb5sF+Zz
kBcBRrTn7UJHI+h7H83WdSjAltkflDAPKMc3gO0lpR4qwVnVyTIR1A5tHaJp
DsG8ArR9JeCamh/N2gjmBND22oWOVCgEs4PQd0HeayVtt13oY5D3o5I2bhf6
Dsj7rZK20y70Ecj7nZI2ahf6Nsj7SUnbbhr6OEezDJbXLeDxme6PnPsfafun
ZYg3Sc/N0pZ510G3Yq7eydGW1101OMFmZ3Mh2dzqqr3VYlp6b08gZVsxzY86
QVBI6D4trWVIqDaDQsr2kWKDtALpxrkv3239efNP5qZDSrdks/flxj43uzCi
bgAw73SSbqlnM65uQhWmteeU9TCtmqaDtuLGQiyjSj5luEjvl5XnAfqTho7W
uPJCUeYnL1BV0OOtqMIxPuLcE4nSfcGsGUCPxYl3nryGA7YEm/Eny5HVVPP5
Xkn70Xxnp+molJEB2uPpSfXuJf3fsQloSyJ5t6zEE503LeixONqxFW19yvuH
vYNulUSrHfUEMN58gjPuzSBROxL5t1TTySDP4/Mg1tD/Qg9qo0y7pFf5wyQ9
fu9WSKlZxMtFC3Pt3ZEez9MDdJB3bfV1gD5AH6DnWfR4g2iAHhx64FOP9UA2
fAzSWifY2CORHkDc7MPLOoBxdnMhlgFQdfBcfB2Dhg5CSuk1QWraEkVhma0F
oqjrdSajTlswddLOq/jUNAPJuvbFjlhVrAQviVc1THZCGwM+R5QQO8RUNrM8
J4bH9AtlF8B+Sw8jspw+oYfnTUnEku0qu/4gD6VtKhOZ6wqiTCgjMr1nt6hX
vkK9WqNMQ3o/VxzyFRsTFrTwFET8hjIeALZnVCjyVQpt383D+AVI3SMsE7+I
dNqLne4RYuSXyaVeUOoLkBraOF92rzygVOSnkzuQzXv2I/cH5OfJRbHjFKqU
0O4UB4QFOb8eEZYPfrHp1DYTM2CH2L41DC2LH+jXK1sIz7Lmc62knQB+5+ZX
Vi3IoDHXiqgBqCsh71GulytKRe5YoR2Ndqm0jyBV9sx0H04cXgA9nK8BDdUF
8jlCPU3rOYZoqJblJ4HvOwX3IVdbx6rCiXqA+gjyK0aT4EdqIVQvzygVSY6G
m2eAi/Ys4RWVJtQBmq1OKONbwPY7SkWzyyHghwYUeT7YplR04im7ewuCvQQ0
JKJnjY5c9b3+sA7tKOnI2g2aM8AD+eAL0iIactRHa3+Ut5zXsmo/5S3IjhpJ
kBI1ClJskJSoqAslzZ3PzYUWWmikRnIIIsou5FwU2p6wIloSAoLiO8j7A3JV
qMSR3cq5KLQERLEc5JlcXgbLXvcqceQRD4nDNYxiPjjCjtxUaa/CP/L4xXF8
4/kmpCWibQNevqABinUWt8V8WwlounhGZQi7QbKvvStEYYDikBHehzBRc9Ps
udkBiXIJGB0AmtC9ZF3KFQpPMLwKQZ3hkOS+IB7FxT8ST9ZYahAPaAloqN4z
FKQ02s6zbyjvS4p6pcQuiCqvFcqLavMeGQrSy7X+0QEEk0MWyIKh5YgdlMe+
NuT1Jnq6tB7t8oJMJbZ2VaIVm8cfNETwxgyagrUxCFoWW15+eyN2OWuR8mLL
wTJUYmtXXdrAGEhDQRXA+zJIRdCGzZD3iGqvAFTXrBho1XqkXmgrQN4bVFWA
1jhIWwEcwMobXM0JYIU6b+1iaw0ztHF4UCdHFYBoKK82Xo+8Py7oxuimKhWA
dG20c4VoKG+5CsjrrQukVWN7dLRLoQ0bhZRCNAAimhyMUhZZPvbxdP/w4iOa
vOfaovgFs30kZWjJEQ09O2g/tH3JURuhER/twyIayquVfLVLSR3SIl0C7dAg
GsrbHWmRLo5mBKSeyipreWk9fjVPU1DmbC6NjYtPSjQwoN1SREN50Rl/N6W0
j2PGsmWTZOO9VZRqkbFxTHH0QMdSawpaKprpIebHq2T0XCHaKWakDdGK7e9C
iMLnkmieQ7OBqwmXlsPJW1GiPDOktCHBEM1d43RKMLT+QDoZorlr1k4JhlaW
3pD8qMVKxToOKg5aAKP9gX3ijLS1l6WFuKMctUO3z8koW03J5e1nASRx2KKM
aIxF+55ITv+uuWfiWD0myierjwAi4hegXthWw2mvKpUzTr6ztPzDE66GJj4B
VNVUwohIctJBAlwpKxDtlSJ+SA9CNJQXlYEq8wrchxoMHU956ha/W8kfLKnJ
CpZ20zpaqblQTcUXJ9m7bJUdl8wEG6Nwyo5GcXRI38V2QOfd/pHqMZ4aYYcG
r4XqRofeiBGj8H4rbWQpnsv4FzrJrl6neGpDWqi2HicAOzofraYAOd4P3Z7Z
2RoXnZfVpNN6db09XwtzuU2r497TR+1efisriSprujBCDGtvvWCFFQ6Sracb
Jpu3t4U3HhXTc2NbkI68obeNbQvuZKsySWT/xoa2jWsW1J7ZWCHQ6U2/9/67
cnrzZA7jBkFXCfq0Ts2flkVE7y1fBouvCuKHtvFDQ6JW/N4ZOdZj5akwclxd
R1rbXq+Rq/OWXP8KzXzQNo/WynedVg9YAWhwls2Xte6qDZo5hzZvRx0beRZ2
0NC9FlcHZ38FPSotOziEdmxhv0g0onfGsSW0G9M7EsxrU2FJLbsxhXZa+0CC
ddppLbQToj0SsyeFCFyDToi1OJJaxRBV2ZhSPyR/Z1MtyKgjPrP49JCfUq/X
XeYp1YsXwOM5hBu3zY2mVo79dE7Zs2u+YG7cld3qUWwGR1uypNrd6qsEOJAj
MnA0h+yaP+UZKsCB4vZS4SdGLuA0VY7mUFP4iSrBQdiYEW2Ocri/0wy/9I6a
goNUDudiAaN6ZfdqN8qeuWoK3fI0YvCguGRcwCdKReqA/Dg6tMwiIi9bmoZe
J7N2uCg5nBgXwLE4kf4qxxO8pNRDkFe7Ktc6EwrClopwplDb5VkcKS5IwXkE
ZciLUlmOQlWsEcoSDfNc/oRS0ZaVNjKerOkzPxSRlFNRILdCDWhDWcoRyR5F
cRD00KETtSsbOQhcoXq0IS//ZAvh+dV8TgCtVPBLh/YDoJ2Z76wKR28WjSNe
0l9wVnXiL3AtcXw+NC6FDk3KeqjXFNSmoolh/bCX7pQClJmRLzV0TFoUflU7
h5QPZ3lJsqHTbl55NBGP+AVhQXoeO7KgcWBVVGlv3NbXVOgpSK0lqrQ3oOaY
sKDHRRDR3j7mRz/PlkOLh4wB/oa4ooqzCSXDf98Qx0P69a4kYk+c9XeAdcXY
+fEzE6URXKvErbccDg1R+waiHsb4FyCl9/YEUqj3RSiM8hBgLY35BXiDhr7Q
Ys0FbaQOvnmlg6/M6SCkzXkPU1PAuv1erY19MVkPIJZRUgPMR7149WBXX+g4
QO8o9B6/MnaAPkAfoK8HvccvYA8F3TvRa490QgNee9LXArYFOLs+rXQJLWCP
BwqaHtFJ5yBRAxKtdifJ80BvMBrEGvpfpUFtlGmX9NLaJXdh2vPORQEA7yeM
+BxsRYyazk/fdKQ7n5oTPQdS78TJBROKshc3lM34k+U4ShLj8+7vlbQfzXdW
uUiPR805vn07ZFyZ2tAwHdSlKdiNrcdl3c36LhEZWyThe0xn77E4T2Cd1uOd
oAF6G9C7tPU9wO3N2Uc4iIHPO1QHb2vvdHT2lLI0MM5uLsQyAKoOHoCvY7nQ
QUgpvSZIHTQ5CW2OExRSM6ZKASG1ZNBVfJLCQ1ph0FXFRBDBqWjWVdmostZW
C2Ok2lDHMp/QpsINjBwpLbQhd4PQQ1vfNwg9tOtEg9DlkAkyTRvJpSbooZ2N
GoSudRlDNG3YqJqg231rOz0icHIkhGqAKwXgNR/kLc8+UdoonwFMsJlWylPW
oWkB19QROOoCAid7ojUIU/bWlmktQ0cu2MxF9vFuBTDHfkCOpchnuxWYckwH
FA+iFZgcawoFC5H98xuEiZ505vKchEABLErCTF+1kT9kjF+7oT07RSLYAcwu
Ow4Al31KRWFe1hSBTixtcuzVi0TQhiHVxuuswq9yGbuGgpTlEaWiebtVzAHL
8ByVp6f/xrZh2SVsXSzi5ykqnHSbVyU9kA3EzDI+y96SeLWn3tLmO54pgA2F
7YL2lnioi2+/Lv6Otv4P7dbm/A==\
\>"],ImageRangeCache->{{{0, 287}, {287, 0}} -> {-6.6861, -6.58559, 0.0446696, \
0.0446696}}],

Cell[BoxData[
 TagBox[
  RowBox[{"\[SkeletonIndicator]", "Graphics", "\[SkeletonIndicator]"}],
  False,
  Editable->False]], "Output"]
}, Open  ]]
},
AutoGeneratedPackage->Automatic,
WindowSize->{844, 740},
WindowMargins->{{433, Automatic}, {Automatic, 64}},
PrintingCopies->1,
PrintingPageRange->{Automatic, Automatic},
PageHeaders->{{
   Cell[
    TextData[{
      CounterBox["Page"]}], "PageNumber"], Inherited, 
   Cell[
    TextData[{
      ValueBox["FileName"]}], "Header"]}, {
   Cell[
    TextData[{
      ValueBox["FileName"]}], "Header"], Inherited, 
   Cell[
    TextData[{
      CounterBox["Page"]}], "PageNumber"]}},
PrintingOptions->{"PrintCellBrackets"->False,
"PrintMultipleHorizontalPages"->False,
"PrintRegistrationMarks"->True,
"PrintingMargins"->{{36, 36}, {57.5625, 36}}},
PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}},
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[CellGroupData[{
Cell[579, 22, 1758, 61, 949, "Input",
 InitializationCell->True],
Cell[2340, 85, 511, 12, 52, "Output"],
Cell[2854, 99, 62, 1, 31, "Output"],
Cell[2919, 102, 45, 0, 31, "Output"],
Cell[2967, 104, 46, 0, 31, "Output"],
Cell[3016, 106, 46, 0, 31, "Output"],
Cell[3065, 108, 46, 0, 31, "Output"],
Cell[3114, 110, 46, 0, 31, "Output"],
Cell[3163, 112, 46, 0, 31, "Output"],
Cell[3212, 114, 45, 0, 31, "Output"],
Cell[3260, 116, 46, 0, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[3343, 121, 131, 4, 52, "Input"],
Cell[3477, 127, 46, 0, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[3560, 132, 2937, 84, 172, "Input"],
Cell[6500, 218, 10153, 547, 296, 5069, 461, "GraphicsData", "PostScript", \
"Graphics"],
Cell[16656, 767, 134, 4, 31, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)