(* 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[     28301,       1615]
NotebookOptionsPosition[     27035,       1574]
NotebookOutlinePosition[     27574,       1594]
CellTagsIndexPosition[     27531,       1591]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["\<\
(* H2C=CH2 with spherical bananas: Two fused methanes with CH4 data \
CH4Mt(min).nb  17.06.2012 *)
Clear[z,sig1,sig2,sig4,k1,k2,k4,nc,R1,R2,R3,R4,w,p,vee,vne,vnn,
xc,yc,zc,xn,yn,zn,oc,ch,rr,cs,ss,d1,d2,d3,d4,d5,pi,i,j,d,t];
z=6.0;
pi=0.0; (* distance of bananas from C-C axis, pi =0 is regular tetraedric *)
nc=8;   (* number of clouds *)
sig1=0.3; sig2=0.3; sig4=sig2; (* screening const. from e-e interaction in \
doubly occ. clouds *)
k1=1.027; k2=1.23; k4=k2; (* parameters for kinetic energy of clouds; k=1.0 \
Kimball's lowest value *)
bohr=0.529177; rad=57.29578;

(* C He-shells *)
Ekin = 2*(2.25*k1/R1^2);
vee=2*(3.0*sig1/R1);
vne=-2*(3.0*z/R1);
(* this is the common edge assumption *)
R4=R2;

(* bonding pairs *)
Ekin = Ekin + 2.25*(4*k2/R2^2+2*k4/R4^2);
vee=vee+3.0*(4*sig2/R2+2*sig4/R4);
(* cloud occupation *)
oc={-2,-2,-2,-2,-2,-2,-2,-2};
(* nuclear charges for C1,C2,H3,H4,H5,H6, banana7, banana8 *)
ch={6,6,1,1,1,1,0,0};
(* cloud radii in the same order *)
rr={R1,R1,R2,R2,R2,R2,R4,R4};
(* w is half angle between two C-H of CH4, i.e. 109.47\[Degree]/2 *)
w=ArcCos[-1/3]/2;
cs=Cos[w];  ss=Sin[w];
(* edge length of tetrahedron of 4 equal clouds *)
a=4*(R1+R2)/Sqrt[6]; (* 4/Sqrt[6] is also Sqrt[8/3] *)
(* x is C-C bond axis, xy plane of molecule *)
(* nuclear coordinates in terms of radii; C nucleus assumed in center of \
C(1s) cloud *)
d1=(R1+R2)*cs;  (* R1+R2 is radius of outer sphere for tetrahedron of equal \
clouds *)
d2=d1+(R1+R2+p)*cs; d3=(R1+R2+p)*ss; 
xn={-d1,d1,-d2,-d2,d2,d2,0,0};
yn={0,0,d3,-d3,-d3,d3,0,0};
zn={0,0,0,0,0,0,0,0};
(* cloud coordinates in terms of radii *)        
d4=d1+(R1+R2)*cs;  d5=(R1+R2)*ss;
xc={-d1,d1,-d4,-d4,d4,d4,0,0};
yc={0,0,d5,-d5,-d5,d5,0,0};
zc={0,0,0,0,0,0,pi+a/2,-pi-a/2};
(* potential energy of protons in CH-clouds with eccentricity p *)
vne=vne-4*(3-(p/R2)^2)/R2;

(* cc: sum of cloud-cloud potential energies *)
For[i = 1, i < nc, i++,
  For[j = i+1, j < nc+1, j++,
    vee = vee + \
oc[[i]]*oc[[j]]/Sqrt[(xc[[i]]-xc[[j]])^2+(yc[[i]]-yc[[j]])^2+(zc[[i]]-zc[[j]])\
^2]]]

(* nn: sum of nuclei-nuclei potential energies *)
vnn = 0.0;
For[i = 1, i < nc-2, i++,
  For[j = i+1, j < nc-1, j++,
    vnn = vnn + ch[[i]]*ch[[j]]/Sqrt[(xn[[i]]-xn[[j]])^2+(yn[[i]]-yn[[j]])^2]]]

(* cn: sum of cloud-nuclei potential energies *)
For[i = 1, i < nc+1, i++,
  For[j = 1, j < nc+1, j++,
    If[i != j,
      vne = vne + \
oc[[i]]*ch[[j]]/Sqrt[(xc[[i]]-xn[[j]])^2+(yc[[i]]-yn[[j]])^2+(zc[[i]]-zn[[j]])\
^2]]]]

Epot=vne+vee+vnn;
func=Ekin+Epot;
(* results of CH4 computation; if this is not available, decomment the \
minimize function *)
(* R1=0.2623610; R2=1.2461360; p=0.53986226; *)
(* minimization function for R1, R2, p *)
t = FindMinimum[func,{R1,0.26},{R2,1.24},{p,0.54},{Method -> Automatic}, \
{MaxIterations -> 500}]

(* func *)
vne  /. t[[2]] 
vee  /. t[[2]] 
vnn  /. t[[2]]
-Epot/Ekin  /. t[[2]]
2*d1*bohr  /. t[[2]]
(R1+R2+p)*bohr  /. t[[2]]
(R1+R2)*bohr  /. t[[2]]
2*w*rad  /. t[[2]]
4*(R1*(R1+2*R4))^(3/2)/(R1+R4)^3 /. t[[2]]\
\>", "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "77.46830376342251`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"R1", "\[Rule]", "0.2627360697917163`"}], ",", 
     RowBox[{"R2", "\[Rule]", "1.2563232439286869`"}], ",", 
     RowBox[{"p", "\[Rule]", "0.5811519153807172`"}]}], "}"}]}], 
  "}"}]], "Output"],

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

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

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

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

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

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

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

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

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

Cell[CellGroupData[{

Cell["\<\
(* projection on xy-plane of molecule *)
plot1=Graphics[{Circle[{xc[[1]],yc[[1]]},R1],      \
Circle[{xc[[2]],yc[[2]]},R1],Circle[{xc[[3]],yc[[3]]},R2],Circle[{xc[[4]],yc[[\
4]]},R2],Circle[{xc[[5]],yc[[5]]},R2],Circle[{xc[[6]],yc[[6]]},R2],Circle[{xc[\
[7]],yc[[7]]},R4],Circle[{xc[[8]],yc[[8]]},R4],Disk[{xn[[1]],yn[[1]]},0.08], \
Disk[{xn[[2]],yn[[2]]},0.08], Disk[{xn[[3]],yn[[3]]},0.08], \
Disk[{xn[[4]],yn[[4]]},0.08],Disk[{xn[[5]],yn[[5]]},0.08],Disk[{xn[[6]],yn[[6]\
]},0.08]} ]  /. t[[2]]

Show[plot1,{AspectRatio \[Rule] Automatic,Axes -> True,GridLines -> \
Automatic, PlotRange \[Rule] {{-4,4},{-3,3}}, Frame -> True}]

(* projection on xz-plane, perpendicular to molecular plane *)
plot2=Graphics[{Circle[{xc[[1]],zc[[1]]},R1],      \
Circle[{xc[[2]],zc[[2]]},R1],Circle[{xc[[3]],zc[[3]]},R2],Circle[{xc[[4]],zc[[\
4]]},R2],Circle[{xc[[5]],zc[[5]]},R2],Circle[{xc[[6]],zc[[6]]},R2],Circle[{xc[\
[7]],zc[[7]]},R4],Circle[{xc[[8]],zc[[8]]},R4],Disk[{xn[[1]],zn[[1]]},0.08], \
Disk[{xn[[2]],zn[[2]]},0.08], Disk[{xn[[3]],zn[[3]]},0.08], \
Disk[{xn[[4]],zn[[4]]},0.08],Disk[{xn[[5]],zn[[5]]},0.08],Disk[{xn[[6]],zn[[6]\
]},0.08]} ]  /. t[[2]]

Show[plot2,{AspectRatio \[Rule] Automatic,Axes -> True,GridLines -> \
Automatic, PlotRange \[Rule] {{-4,4},{-3,3}}, Frame -> True}]

\
\>", "Input"],

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

Cell[GraphicsData["PostScript", "\<\
%!
%%Creator: Mathematica
%%AspectRatio: .75 
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.125 0.375 0.125 [
[.125 -0.0125 -6 -9 ]
[.125 -0.0125 6 0 ]
[.25 -0.0125 -6 -9 ]
[.25 -0.0125 6 0 ]
[.375 -0.0125 -6 -9 ]
[.375 -0.0125 6 0 ]
[.5 -0.0125 -3 -9 ]
[.5 -0.0125 3 0 ]
[.625 -0.0125 -3 -9 ]
[.625 -0.0125 3 0 ]
[.75 -0.0125 -3 -9 ]
[.75 -0.0125 3 0 ]
[.875 -0.0125 -3 -9 ]
[.875 -0.0125 3 0 ]
[1 -0.0125 -3 -9 ]
[1 -0.0125 3 0 ]
[ 0 0 -0.125 0 ]
[-0.0125 .125 -12 -4.5 ]
[-0.0125 .125 0 4.5 ]
[-0.0125 .25 -12 -4.5 ]
[-0.0125 .25 0 4.5 ]
[-0.0125 .375 -6 -4.5 ]
[-0.0125 .375 0 4.5 ]
[-0.0125 .5 -6 -4.5 ]
[-0.0125 .5 0 4.5 ]
[-0.0125 .625 -6 -4.5 ]
[-0.0125 .625 0 4.5 ]
[-0.0125 .75 -6 -4.5 ]
[-0.0125 .75 0 4.5 ]
[ 0 0 -0.125 0 ]
[ 0 .75 .125 0 ]
[ 1 0 .125 0 ]
[ 0 0 0 0 ]
[ 1 .75 0 0 ]
] MathScale
% Start of Graphics
1 setlinecap
1 setlinejoin
newpath
0 0 .5 r
.25 Mabswid
[ ] 0 setdash
.125 0 m
.125 .75 L
s
.25 0 m
.25 .75 L
s
.375 0 m
.375 .75 L
s
.5 0 m
.5 .75 L
s
.625 0 m
.625 .75 L
s
.75 0 m
.75 .75 L
s
.875 0 m
.875 .75 L
s
0 .125 m
1 .125 L
s
0 .25 m
1 .25 L
s
0 .375 m
1 .375 L
s
0 .5 m
1 .5 L
s
0 .625 m
1 .625 L
s
0 g
.125 0 m
.125 .00625 L
s
[(-3)] .125 -0.0125 0 1 Mshowa
.25 0 m
.25 .00625 L
s
[(-2)] .25 -0.0125 0 1 Mshowa
.375 0 m
.375 .00625 L
s
[(-1)] .375 -0.0125 0 1 Mshowa
.5 0 m
.5 .00625 L
s
[(0)] .5 -0.0125 0 1 Mshowa
.625 0 m
.625 .00625 L
s
[(1)] .625 -0.0125 0 1 Mshowa
.75 0 m
.75 .00625 L
s
[(2)] .75 -0.0125 0 1 Mshowa
.875 0 m
.875 .00625 L
s
[(3)] .875 -0.0125 0 1 Mshowa
1 0 m
1 .00625 L
s
[(4)] 1 -0.0125 0 1 Mshowa
.125 Mabswid
.025 0 m
.025 .00375 L
s
.05 0 m
.05 .00375 L
s
.075 0 m
.075 .00375 L
s
.1 0 m
.1 .00375 L
s
.15 0 m
.15 .00375 L
s
.175 0 m
.175 .00375 L
s
.2 0 m
.2 .00375 L
s
.225 0 m
.225 .00375 L
s
.275 0 m
.275 .00375 L
s
.3 0 m
.3 .00375 L
s
.325 0 m
.325 .00375 L
s
.35 0 m
.35 .00375 L
s
.4 0 m
.4 .00375 L
s
.425 0 m
.425 .00375 L
s
.45 0 m
.45 .00375 L
s
.475 0 m
.475 .00375 L
s
.525 0 m
.525 .00375 L
s
.55 0 m
.55 .00375 L
s
.575 0 m
.575 .00375 L
s
.6 0 m
.6 .00375 L
s
.65 0 m
.65 .00375 L
s
.675 0 m
.675 .00375 L
s
.7 0 m
.7 .00375 L
s
.725 0 m
.725 .00375 L
s
.775 0 m
.775 .00375 L
s
.8 0 m
.8 .00375 L
s
.825 0 m
.825 .00375 L
s
.85 0 m
.85 .00375 L
s
.9 0 m
.9 .00375 L
s
.925 0 m
.925 .00375 L
s
.95 0 m
.95 .00375 L
s
.975 0 m
.975 .00375 L
s
.25 Mabswid
0 0 m
1 0 L
s
0 .125 m
.00625 .125 L
s
[(-2)] -0.0125 .125 1 0 Mshowa
0 .25 m
.00625 .25 L
s
[(-1)] -0.0125 .25 1 0 Mshowa
0 .375 m
.00625 .375 L
s
[(0)] -0.0125 .375 1 0 Mshowa
0 .5 m
.00625 .5 L
s
[(1)] -0.0125 .5 1 0 Mshowa
0 .625 m
.00625 .625 L
s
[(2)] -0.0125 .625 1 0 Mshowa
0 .75 m
.00625 .75 L
s
[(3)] -0.0125 .75 1 0 Mshowa
.125 Mabswid
0 .025 m
.00375 .025 L
s
0 .05 m
.00375 .05 L
s
0 .075 m
.00375 .075 L
s
0 .1 m
.00375 .1 L
s
0 .15 m
.00375 .15 L
s
0 .175 m
.00375 .175 L
s
0 .2 m
.00375 .2 L
s
0 .225 m
.00375 .225 L
s
0 .275 m
.00375 .275 L
s
0 .3 m
.00375 .3 L
s
0 .325 m
.00375 .325 L
s
0 .35 m
.00375 .35 L
s
0 .4 m
.00375 .4 L
s
0 .425 m
.00375 .425 L
s
0 .45 m
.00375 .45 L
s
0 .475 m
.00375 .475 L
s
0 .525 m
.00375 .525 L
s
0 .55 m
.00375 .55 L
s
0 .575 m
.00375 .575 L
s
0 .6 m
.00375 .6 L
s
0 .65 m
.00375 .65 L
s
0 .675 m
.00375 .675 L
s
0 .7 m
.00375 .7 L
s
0 .725 m
.00375 .725 L
s
.25 Mabswid
0 0 m
0 .75 L
s
0 .74375 m
0 .75 L
s
.125 .74375 m
.125 .75 L
s
.25 .74375 m
.25 .75 L
s
.375 .74375 m
.375 .75 L
s
.5 .74375 m
.5 .75 L
s
.625 .74375 m
.625 .75 L
s
.75 .74375 m
.75 .75 L
s
.875 .74375 m
.875 .75 L
s
.125 Mabswid
.025 .74625 m
.025 .75 L
s
.05 .74625 m
.05 .75 L
s
.075 .74625 m
.075 .75 L
s
.1 .74625 m
.1 .75 L
s
.15 .74625 m
.15 .75 L
s
.175 .74625 m
.175 .75 L
s
.2 .74625 m
.2 .75 L
s
.225 .74625 m
.225 .75 L
s
.275 .74625 m
.275 .75 L
s
.3 .74625 m
.3 .75 L
s
.325 .74625 m
.325 .75 L
s
.35 .74625 m
.35 .75 L
s
.4 .74625 m
.4 .75 L
s
.425 .74625 m
.425 .75 L
s
.45 .74625 m
.45 .75 L
s
.475 .74625 m
.475 .75 L
s
.525 .74625 m
.525 .75 L
s
.55 .74625 m
.55 .75 L
s
.575 .74625 m
.575 .75 L
s
.6 .74625 m
.6 .75 L
s
.65 .74625 m
.65 .75 L
s
.675 .74625 m
.675 .75 L
s
.7 .74625 m
.7 .75 L
s
.725 .74625 m
.725 .75 L
s
.775 .74625 m
.775 .75 L
s
.8 .74625 m
.8 .75 L
s
.825 .74625 m
.825 .75 L
s
.85 .74625 m
.85 .75 L
s
.9 .74625 m
.9 .75 L
s
.925 .74625 m
.925 .75 L
s
.95 .74625 m
.95 .75 L
s
.975 .74625 m
.975 .75 L
s
.25 Mabswid
0 .75 m
1 .75 L
s
.99375 0 m
1 0 L
s
.99375 .125 m
1 .125 L
s
.99375 .25 m
1 .25 L
s
.99375 .375 m
1 .375 L
s
.99375 .5 m
1 .5 L
s
.99375 .625 m
1 .625 L
s
.125 Mabswid
.99625 .025 m
1 .025 L
s
.99625 .05 m
1 .05 L
s
.99625 .075 m
1 .075 L
s
.99625 .1 m
1 .1 L
s
.99625 .15 m
1 .15 L
s
.99625 .175 m
1 .175 L
s
.99625 .2 m
1 .2 L
s
.99625 .225 m
1 .225 L
s
.99625 .275 m
1 .275 L
s
.99625 .3 m
1 .3 L
s
.99625 .325 m
1 .325 L
s
.99625 .35 m
1 .35 L
s
.99625 .4 m
1 .4 L
s
.99625 .425 m
1 .425 L
s
.99625 .45 m
1 .45 L
s
.99625 .475 m
1 .475 L
s
.99625 .525 m
1 .525 L
s
.99625 .55 m
1 .55 L
s
.99625 .575 m
1 .575 L
s
.99625 .6 m
1 .6 L
s
.99625 .65 m
1 .65 L
s
.99625 .675 m
1 .675 L
s
.99625 .7 m
1 .7 L
s
.99625 .725 m
1 .725 L
s
.25 Mabswid
1 0 m
1 .75 L
s
0 .375 m
1 .375 L
s
.5 0 m
.5 .75 L
s
0 0 m
1 0 L
1 .75 L
0 .75 L
closepath
clip
newpath
.5 Mabswid
newpath
.39037 .375 .03284 0 365.73 arc
s
newpath
.60963 .375 .03284 0 365.73 arc
s
newpath
.28074 .53004 .15704 0 365.73 arc
s
newpath
.28074 .21996 .15704 0 365.73 arc
s
newpath
.71926 .21996 .15704 0 365.73 arc
s
newpath
.71926 .53004 .15704 0 365.73 arc
s
newpath
.5 .375 .15704 0 365.73 arc
s
newpath
.5 .375 .15704 0 365.73 arc
s
.39037 .375 m
.39037 .375 .01 0 365.73 arc
F
.60963 .375 m
.60963 .375 .01 0 365.73 arc
F
.2388 .58935 m
.2388 .58935 .01 0 365.73 arc
F
.2388 .16065 m
.2388 .16065 .01 0 365.73 arc
F
.7612 .16065 m
.7612 .16065 .01 0 365.73 arc
F
.7612 .58935 m
.7612 .58935 .01 0 365.73 arc
F
% End of Graphics
MathPictureEnd
\
\>"], "Graphics",
 ImageSize->{288, 216},
 ImageMargins->{{43, 0}, {0, 0}},
 ImageRegion->{{0, 1}, {0, 
  1}},ImageCache->GraphicsData["CompressedBitmap", "\<\
eJztXc1y2zYQRvRjO07ixHbs1H+Rc0hOPbW959RElz6EJ9OZ5NBpJ9XdL6J3
0Wv0lkdRTZBYkOAHcEGAIiyRE8kKAC72WwALENhd/nG3+PrnX3eLb1/ubj9/
v/vn67cv/95++vv7Q9L4iRBPboUQ/92K7Pf64WfxtX74t36XfQnxq/wzET+W
c0HXYnUtk/fEerUQ1WshVuv1lcper0S1wFwsf6wvANH58se1+64rd53XmNO3
tWRJbJLD+0XR/LEU8/lS/i+vXpbKaY5Ffq+sRsyKtIyNB6bl/y7l91QsF4rE
+uGevPwNKH8D6N64WXkLyBQgfiuylFwWq7UoF0dpisRFWeZzzf1qIWtFzLtI
l1loaM6C+Z8Natk1kt97kopiiMQxX5573rJare23TKu3FDgv/cVC/SKvsSqq
EhnVTzNOZNU5janRIXXLF7+rA6vWSUpjYJz3vYK3K9U35zVqApW390MwtGsk
FTrKXpbRKW0iadMQLcRZ0wdlKRdJBRoLhWrR+sgvUa1RoKYqkjMKBddlfRjy
+zj7I1XsYd5ZdepDxdmlutkEpI0jl0u9jtT5G+rouw44dN5S2n1AWnf0BtYH
1luzLpcss4QYviGuHLPZY5LwwPrA+sB68qwPinDoHLFZlyTkrW8CQJwRlb3N
wZG97yArNaZfLwNAvAD03KOtNetXNHDyuijzFWW8CEDyjKgcF9yp3OK66QaW
7EITjaqae0y5TwOw7ROVUwNb8ZF5Fy0AWtXrGeXue+dyYU2JChqH7tzWaviU
bkQt8oZypwHAxkQFtck+5Z510yV1yz4DuReUizqsHWKt00kaV4DGU8o97Qbg
K1AOQbWqnJm1oipKi9owdU10gEdUwZGt7xdUmhHl96rLOgkY96J6IwB7TsBe
2SowVJKmp4f8rbfIWuvBZ+BGK+ultKqs9b3XRlpWHsHhdvHWetAPmGNY+GND
LB52098OQTn7mJ1W0lQzmrCyjzqdKKc93IvQnoDbuYsBT7RPQbkTZ1U6bQTu
vQBp6gS0nOZq7FNA46Ab+AegnHsG0mlo8/knpkhGzDpeg3u566dCFDlE69E4
QoE4kbpCao1zJu9ow/6NLY2xgNJpZ4CK+4mqJhXV78un7FNQ3LoINcrFFso5
SNugeJQSLJ/j62fgD5KMXk1kn48g7ZyZhu59n32XZ6mRwVJ2NZzsI2H7N4pb
2HvgXtQA/o2i0/w7g10L4JN7rhYIkdQ+uBfptxBJoaE0cd5r0Qw10wBuX3LP
Cm7u3bORW2ruWeG+yoscN55y4U6p7lVBHAH4rxDuAS/c6bgXUbgXZnFEwVuY
NT7+NK2B4wD1XxTzgM4ITePjUBhQ94NF90CbWhRCrj3FSf2RP8XYNkj8Ifs/
3untiRFx1SluObPsFbhjNXUI7rHG7XiobKW7skvqqWlkvNwtir7wcvu1eyck
Dl7Zry8NAq37cxUs99iLi/MlSHvOrMO9BeSx9+YD7DkTRAgwtOeqtyuZe28q
PXVgM0CP2SkHiK0hHoE094loMhDdJwhxILoP3zwhupf6A8SkIFpPhFIH5nk4
4w/MzfBWdMoBYmSISLv3AlHOizIJzYaPGJisSppSCfX0kWfKNLRz2iXYYO2J
II6IMtrx0RYu6HQlBOxmNKp8UkaPSxJs9blSGDJKure+IPFZh5yRm7Ti1HY4
1uYwtKdMeiTTgjajQlvP3PWL1Q4sAmI6h8KOv+5zKH+BoAOFMNGsa0osqnzK
B69VV+54QtFH+UjXckWBuAoxa24jFFl6JkmPi26Unx7/nn0X1qb13CNm2pSo
zEHuDeU2nlBjR/ImewOVzm1XNHjQ8TR3JwEtO9wGoDrNc3dBm2NoZ/OUxIPS
kHi4dj0txVN2YY8tHqQWEUSueHLVy7f6irwflVfl2hHGp3noKIMLWZkOlNME
NvHbONragRYCyt3bN++1telG9/bVFXI4UwZWbxY3sBlx72FX2wzMLvXswz1l
4wJrtFL16Josw/gQoO7jU/fCaONA0YhHXTU20OgmAXwriBDITVYQKj0xyEgx
IMh+Fi4sz7ToqJvXB8iYyTKwuEKw2zOxRlkkk6ao0NFI9bNq4/c/rnVbaycQ
rn2jqyG4ArHbOcYXSOtdmnujeLUcMo51ieaA2gUZ/XJNY7kd1t93sSORcFvJ
bXrfnaAiWGBzRcb1LogtMi7ERywy5J2DNjRiO2lwRRbBcYMrMrfbQnci4xro
c6eZWCKzzv1u1+s4cIw0q3U+Vyhc/4DWiwCu1yBXPDE8CVGPRlPNBr0Lua6k
djEFu5dypbJBl1OuG7FdKq1di7nS2KC7sdvUjz+dqkuWuAKF5XgaF2VUT/IR
iuTKK+iKpyi4QQWsohgTOPkLKdhzUA7tcqPlg78xjedcw60ACUAHEEEzvfWJ
4jXllnWgurgLK27LtZ5w3H0PAdMntaJqz4/mLr+n8dpuDFr8IDFxx3qEscTV
YOeEqDKjGjYEunxsDz0jqoxbfofEmjtEQWupnVAFdlmVZxRh8q8Lh7u11aLn
mMo7+4sUtd5DeA1yI4gJK42KmEodCC/XwjdF8dpIj3u0MuZGMaiJJtArWNYm
Rxla/kqMU924e1BiigKar7mSGxMV3TZ1GU6oFOpAI8p1+zVbZIjfhMB1QL+k
yhFXmbFHPmmVZI36X33M5uUt4uWeyUyoIY1JZ9zEyYQ43zd7ibpkt3U7m9ZE
XjG3qr+rQFZdnAYTdNZNjKZiRO+bdtlA1v7v20AT4jdOnEAaB+CtC1zhcmMI
GtPFqGjaOkytX9By2Eu+Jf0Cxmh9sSGlHjFModLVxlscuLqaG8fQ3mW6m930
Kt8//mGEid8dCVE/W6GGFjEOwmpP/iOqE01GPcdNxM9plsERe2VtVOReWScT
f9HrUQ6JDO192vdDWz/KhcVzdO8AWEMSIuHoyUl3drRzZmg1mZnM4707xiNj
v8eInGpsZGl1PSYBMGCHxYdkwObuhHLOmXKIwrH3h4IboB6A5g7/QAbD9vBw
hNB71MptPn+yziHcmJaxzy65BgnJnGK6N4M6Ofm22tps+rybIZ7BqmKw3Wmg
l6CFWOx4nAmKLIYJJppB/R6d2okpQoDObbVUDYvMORgvM7rE9tuwQyGk476Q
Juqu/VSyayN+KrvliLSN/mUtvR/9sXXgJDgjAFEdIO/p/r48O1v1Ti7E2hZ0
qo65QSi3ydlapW+pZ31QiBZ/6TzisAz+wVr8xYPooW0BblAPFPxlM0KRpWeS
hCaWfT6BtClI8wrbYqR9Bmm1AC4bC74DJdfcdj2H4fEP0uIvHn2IGSIUpFO4
L8hoKRR+aBZ/oWiLEFPhC7gyOCoYNMpbDVsjiCF2JDNt6otmT/uBtAZdXNyI
zz3jdUe25J5EJwNHsjmpyt9YpU6rYLlw0AzKepRqCaI+jGQCOmTbo1yku9KP
YYlWGcc2sI842OoOBModIAZA7E9jJghxiKSeGrANdspHB9H/+XoHIO7A20SG
d8L0BnF33sXkCyzVt2cxIO7U29Ck4UXsV/sl9ao7/BbH2G/3S+pthjvw4spo
r2Dt+82kM0LT8StY+37XbKCyyj7b/E7hjuFv2dulU3/TeARR0FH0qnoQLTb3
1np3+Jzgd9W3Cb6pz+FIEo1vA4ghitjBN1EdEcyoy6e4UkCr8gH0B0lmXHSj
/CD4I0g7Z6ahe99n3+X5TC38MpOK1co+fP2l7pZmbJcLVEcE9woa5pnDPgWc
oKAIXElxY1qGSMrffj2qpALD0HD94bmTjJHGMF7Xaf5TS+uZkusu6L/g0Gm8
dY+mm2AYSt4KW9NoNHHVZBIML2l5TOT6CiOUsfyCrb4c/gHkkAG6vte933df
rcMafq+1OwY3rgB3/5x71BErkGPEqrg7r9wd5BDRekJ0dw7uBqU5HIsLmfQn
GEbxkmTAfcLx3c7oOfahv3OadTI1tnJ1bljcQqvSlMSMyGE6V0cr5C7FrMtH
WQdamrsdzVurT3ebnBBL3GWFdaPA2iatQwFyIVrD/SEDM390h0SlvEJQV5uA
e1xgMu3YrNmX/yOigpYjsdj0SuuO3sB636xbVWxfDLdWnalKeGB9YH1XWB+0
ydA5BtZTYR3tLIvFyt/9rf+0xlMFqx8rQhu7XOp1pM7fUEffddAGkTqVpxNW
uX1VDiVfPvikUPLZd378Yfia18Khl4taonrLIvneRE75Xf5bPPkf0EGusQ==
\
\>"],ImageRangeCache->{{{0, 287}, {215, 0}} -> {-4.53044, -3.37249, \
0.0302733, 0.0302733}}],

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

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

Cell[GraphicsData["PostScript", "\<\
%!
%%Creator: Mathematica
%%AspectRatio: .75 
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.125 0.375 0.125 [
[.125 -0.0125 -6 -9 ]
[.125 -0.0125 6 0 ]
[.25 -0.0125 -6 -9 ]
[.25 -0.0125 6 0 ]
[.375 -0.0125 -6 -9 ]
[.375 -0.0125 6 0 ]
[.5 -0.0125 -3 -9 ]
[.5 -0.0125 3 0 ]
[.625 -0.0125 -3 -9 ]
[.625 -0.0125 3 0 ]
[.75 -0.0125 -3 -9 ]
[.75 -0.0125 3 0 ]
[.875 -0.0125 -3 -9 ]
[.875 -0.0125 3 0 ]
[1 -0.0125 -3 -9 ]
[1 -0.0125 3 0 ]
[ 0 0 -0.125 0 ]
[-0.0125 .125 -12 -4.5 ]
[-0.0125 .125 0 4.5 ]
[-0.0125 .25 -12 -4.5 ]
[-0.0125 .25 0 4.5 ]
[-0.0125 .375 -6 -4.5 ]
[-0.0125 .375 0 4.5 ]
[-0.0125 .5 -6 -4.5 ]
[-0.0125 .5 0 4.5 ]
[-0.0125 .625 -6 -4.5 ]
[-0.0125 .625 0 4.5 ]
[-0.0125 .75 -6 -4.5 ]
[-0.0125 .75 0 4.5 ]
[ 0 0 -0.125 0 ]
[ 0 .75 .125 0 ]
[ 1 0 .125 0 ]
[ 0 0 0 0 ]
[ 1 .75 0 0 ]
] MathScale
% Start of Graphics
1 setlinecap
1 setlinejoin
newpath
0 0 .5 r
.25 Mabswid
[ ] 0 setdash
.125 0 m
.125 .75 L
s
.25 0 m
.25 .75 L
s
.375 0 m
.375 .75 L
s
.5 0 m
.5 .75 L
s
.625 0 m
.625 .75 L
s
.75 0 m
.75 .75 L
s
.875 0 m
.875 .75 L
s
0 .125 m
1 .125 L
s
0 .25 m
1 .25 L
s
0 .375 m
1 .375 L
s
0 .5 m
1 .5 L
s
0 .625 m
1 .625 L
s
0 g
.125 0 m
.125 .00625 L
s
[(-3)] .125 -0.0125 0 1 Mshowa
.25 0 m
.25 .00625 L
s
[(-2)] .25 -0.0125 0 1 Mshowa
.375 0 m
.375 .00625 L
s
[(-1)] .375 -0.0125 0 1 Mshowa
.5 0 m
.5 .00625 L
s
[(0)] .5 -0.0125 0 1 Mshowa
.625 0 m
.625 .00625 L
s
[(1)] .625 -0.0125 0 1 Mshowa
.75 0 m
.75 .00625 L
s
[(2)] .75 -0.0125 0 1 Mshowa
.875 0 m
.875 .00625 L
s
[(3)] .875 -0.0125 0 1 Mshowa
1 0 m
1 .00625 L
s
[(4)] 1 -0.0125 0 1 Mshowa
.125 Mabswid
.025 0 m
.025 .00375 L
s
.05 0 m
.05 .00375 L
s
.075 0 m
.075 .00375 L
s
.1 0 m
.1 .00375 L
s
.15 0 m
.15 .00375 L
s
.175 0 m
.175 .00375 L
s
.2 0 m
.2 .00375 L
s
.225 0 m
.225 .00375 L
s
.275 0 m
.275 .00375 L
s
.3 0 m
.3 .00375 L
s
.325 0 m
.325 .00375 L
s
.35 0 m
.35 .00375 L
s
.4 0 m
.4 .00375 L
s
.425 0 m
.425 .00375 L
s
.45 0 m
.45 .00375 L
s
.475 0 m
.475 .00375 L
s
.525 0 m
.525 .00375 L
s
.55 0 m
.55 .00375 L
s
.575 0 m
.575 .00375 L
s
.6 0 m
.6 .00375 L
s
.65 0 m
.65 .00375 L
s
.675 0 m
.675 .00375 L
s
.7 0 m
.7 .00375 L
s
.725 0 m
.725 .00375 L
s
.775 0 m
.775 .00375 L
s
.8 0 m
.8 .00375 L
s
.825 0 m
.825 .00375 L
s
.85 0 m
.85 .00375 L
s
.9 0 m
.9 .00375 L
s
.925 0 m
.925 .00375 L
s
.95 0 m
.95 .00375 L
s
.975 0 m
.975 .00375 L
s
.25 Mabswid
0 0 m
1 0 L
s
0 .125 m
.00625 .125 L
s
[(-2)] -0.0125 .125 1 0 Mshowa
0 .25 m
.00625 .25 L
s
[(-1)] -0.0125 .25 1 0 Mshowa
0 .375 m
.00625 .375 L
s
[(0)] -0.0125 .375 1 0 Mshowa
0 .5 m
.00625 .5 L
s
[(1)] -0.0125 .5 1 0 Mshowa
0 .625 m
.00625 .625 L
s
[(2)] -0.0125 .625 1 0 Mshowa
0 .75 m
.00625 .75 L
s
[(3)] -0.0125 .75 1 0 Mshowa
.125 Mabswid
0 .025 m
.00375 .025 L
s
0 .05 m
.00375 .05 L
s
0 .075 m
.00375 .075 L
s
0 .1 m
.00375 .1 L
s
0 .15 m
.00375 .15 L
s
0 .175 m
.00375 .175 L
s
0 .2 m
.00375 .2 L
s
0 .225 m
.00375 .225 L
s
0 .275 m
.00375 .275 L
s
0 .3 m
.00375 .3 L
s
0 .325 m
.00375 .325 L
s
0 .35 m
.00375 .35 L
s
0 .4 m
.00375 .4 L
s
0 .425 m
.00375 .425 L
s
0 .45 m
.00375 .45 L
s
0 .475 m
.00375 .475 L
s
0 .525 m
.00375 .525 L
s
0 .55 m
.00375 .55 L
s
0 .575 m
.00375 .575 L
s
0 .6 m
.00375 .6 L
s
0 .65 m
.00375 .65 L
s
0 .675 m
.00375 .675 L
s
0 .7 m
.00375 .7 L
s
0 .725 m
.00375 .725 L
s
.25 Mabswid
0 0 m
0 .75 L
s
0 .74375 m
0 .75 L
s
.125 .74375 m
.125 .75 L
s
.25 .74375 m
.25 .75 L
s
.375 .74375 m
.375 .75 L
s
.5 .74375 m
.5 .75 L
s
.625 .74375 m
.625 .75 L
s
.75 .74375 m
.75 .75 L
s
.875 .74375 m
.875 .75 L
s
.125 Mabswid
.025 .74625 m
.025 .75 L
s
.05 .74625 m
.05 .75 L
s
.075 .74625 m
.075 .75 L
s
.1 .74625 m
.1 .75 L
s
.15 .74625 m
.15 .75 L
s
.175 .74625 m
.175 .75 L
s
.2 .74625 m
.2 .75 L
s
.225 .74625 m
.225 .75 L
s
.275 .74625 m
.275 .75 L
s
.3 .74625 m
.3 .75 L
s
.325 .74625 m
.325 .75 L
s
.35 .74625 m
.35 .75 L
s
.4 .74625 m
.4 .75 L
s
.425 .74625 m
.425 .75 L
s
.45 .74625 m
.45 .75 L
s
.475 .74625 m
.475 .75 L
s
.525 .74625 m
.525 .75 L
s
.55 .74625 m
.55 .75 L
s
.575 .74625 m
.575 .75 L
s
.6 .74625 m
.6 .75 L
s
.65 .74625 m
.65 .75 L
s
.675 .74625 m
.675 .75 L
s
.7 .74625 m
.7 .75 L
s
.725 .74625 m
.725 .75 L
s
.775 .74625 m
.775 .75 L
s
.8 .74625 m
.8 .75 L
s
.825 .74625 m
.825 .75 L
s
.85 .74625 m
.85 .75 L
s
.9 .74625 m
.9 .75 L
s
.925 .74625 m
.925 .75 L
s
.95 .74625 m
.95 .75 L
s
.975 .74625 m
.975 .75 L
s
.25 Mabswid
0 .75 m
1 .75 L
s
.99375 0 m
1 0 L
s
.99375 .125 m
1 .125 L
s
.99375 .25 m
1 .25 L
s
.99375 .375 m
1 .375 L
s
.99375 .5 m
1 .5 L
s
.99375 .625 m
1 .625 L
s
.125 Mabswid
.99625 .025 m
1 .025 L
s
.99625 .05 m
1 .05 L
s
.99625 .075 m
1 .075 L
s
.99625 .1 m
1 .1 L
s
.99625 .15 m
1 .15 L
s
.99625 .175 m
1 .175 L
s
.99625 .2 m
1 .2 L
s
.99625 .225 m
1 .225 L
s
.99625 .275 m
1 .275 L
s
.99625 .3 m
1 .3 L
s
.99625 .325 m
1 .325 L
s
.99625 .35 m
1 .35 L
s
.99625 .4 m
1 .4 L
s
.99625 .425 m
1 .425 L
s
.99625 .45 m
1 .45 L
s
.99625 .475 m
1 .475 L
s
.99625 .525 m
1 .525 L
s
.99625 .55 m
1 .55 L
s
.99625 .575 m
1 .575 L
s
.99625 .6 m
1 .6 L
s
.99625 .65 m
1 .65 L
s
.99625 .675 m
1 .675 L
s
.99625 .7 m
1 .7 L
s
.99625 .725 m
1 .725 L
s
.25 Mabswid
1 0 m
1 .75 L
s
0 .375 m
1 .375 L
s
.5 0 m
.5 .75 L
s
0 0 m
1 0 L
1 .75 L
0 .75 L
closepath
clip
newpath
.5 Mabswid
newpath
.39037 .375 .03284 0 365.73 arc
s
newpath
.60963 .375 .03284 0 365.73 arc
s
newpath
.28074 .375 .15704 0 365.73 arc
s
newpath
.28074 .375 .15704 0 365.73 arc
s
newpath
.71926 .375 .15704 0 365.73 arc
s
newpath
.71926 .375 .15704 0 365.73 arc
s
newpath
.5 .53004 .15704 0 365.73 arc
s
newpath
.5 .21996 .15704 0 365.73 arc
s
.39037 .375 m
.39037 .375 .01 0 365.73 arc
F
.60963 .375 m
.60963 .375 .01 0 365.73 arc
F
.2388 .375 m
.2388 .375 .01 0 365.73 arc
F
.2388 .375 m
.2388 .375 .01 0 365.73 arc
F
.7612 .375 m
.7612 .375 .01 0 365.73 arc
F
.7612 .375 m
.7612 .375 .01 0 365.73 arc
F
% End of Graphics
MathPictureEnd
\
\>"], "Graphics",
 ImageSize->{288, 216},
 ImageMargins->{{43, 0}, {0, 0}},
 ImageRegion->{{0, 1}, {0, 
  1}},ImageCache->GraphicsData["CompressedBitmap", "\<\
eJztXUtyGzcQhfmVKP8kW3JsiRK9cqUqySLJGRJucoZE5UqVvUgl5XCvi/Au
vEZ2OgrDwQw+M/MANgYYDsYEyqQpANPo18A0Pt0Afrtfffrzr/vV54/3i1+/
3P/z6fPHfxe//P1lFzV8wtiTBWPsvwXLfm93P4uv7e7f9n32xdhP/L8Re1wv
mQyrzQ2PnrDtZsXKYcU22+21SN5uWDnDkq0ft28B0eX68cb+1LW9zBvM6W0t
mhMb5fB+FDQf12y5XPO/8uJ5rpzmkOXP8mLYXRGXsbFjmv/1jn+P2XolSGx3
z+T55yD/HNCd21m5BWQKED8XSUIuq82W6dlRnCDxVpf5UnG/WfFSEfM20joL
e6qzYP67CrUsDPj3hFMRDElxLNdXjo9sNlvzI+PyIwXOd+5ike0iL7EsKo2M
aKcZJ7zonMa40iBVzRe/yy9WrZFo78Awb3sFb9eibS5r1BjKb26H4NWukRTo
ZPJaRye0CactX9FCnDV9oEu5iCrQGCiUs9bffI1qjYKsqiI6o1BwretDn9/n
2X9cxc7yxqpidwVnQTSzEYgbBs4Xexmx85fK6LoM+OrcyrgHj7j26CXWE+uN
WedDlruIGJ5Lriy9WZ8knFhPrCfWo2c9KcLUOA7BOq/WYRuNrSWGL7PvE04u
HzgLwjzuTT9AvJIv0ymgPJWpl/2A8yL7HmXUzgDlU5l63iIcb4X5TFJ5AYpC
qREq0TNZAKL8QqY+60fDOpUMvwKUVb2f9QPOCaCHgCH1EDWwiWT4ClC+BHGT
fgAbSWDfAMpXMrVNON6aDS0wITgobuhRbuvajgrsLYgb9KMBDohwhIFOj+MM
7WgtYoA1LsXlrDFpn9HpCguIHrd7NgoUo1LcLswBySqoQHWQ9x5G8zB6GZq0
PBEQJzpa8by7ZqohE72cbi0eB4GDnm0dzqxI0u3HPInn/oOTUYVmn1ti3O/Z
t65RB5WisrDHUoxMEE4a3Sgwu7ADKFf5BmDLrc8b4I7sQcYHQHZa/FEx5frU
lV0LVekpxVaEuT8oH3W7txMwAyhA9A3Avg7SD4XTIJI2mAgLIOhgsRkAl5FX
J23KfRybIMYH0X2adfQQ7RO9TiC+MZKND6KT8g8NDI1hqMv0QTuFmID1pFFG
CBEtMI6JxR8pRDSVviLykiAeA0SnDiImYEE7iJiApUYZC8QJyOdulw8KUS0m
IlN1aLBPZWnuVq7GsPmLzV06hvLX9ABgR6Bcu2J01KjXMpUJ/5U88aVMQC9V
aJwDWdp5wbt4ogho+ayxtuUDqpHCXE49l6ms7NJjb28E0HwF6FUFYPHhBdrn
yY1b76UU1NQ51Un1cipouGpPDQDR7qPAlfCw3JjdIQ4lFTQTprpJNYZo9y9B
/czA1MTNiGvtktNAazqoX0UOMiQZSFMK3rtINaW8BPnowtpiTXhnLLy2yI9m
b0hOyC+KJCfdllfelbpfOM9loc+JItEp75dGTk8E6lDsGciH+GsiFJ77+4wM
h/0tJ6jgZJ+xTL0DqXOZOgWpH2TqD/LXXkMo3v9KNS8jh8aXIB9l1JkFJFVk
qUN9NtLk1BfQUpvK+q72yIYWj3lhqrFHBVVCSNvPQD6ChPTNt1QJzUA+pIPM
EvL2nKFK6gLkszsgNu5akdvzBchHXdAM4SZFFRPy2DwB+QKIye4w6i6m0A5z
VKPFa5CPOvYlzaqQA+oDYIkB1h8qRZlaQCWuMlVzFwrVLbbxrIs6IggtHuoK
DNVw6e4/3PiFo4qM+oLY6bUnMp9FMEeRUQ0oVJHZ19vscahV9Fhktt4K9aBI
zYYWk4+J0VFMVEVqExMafqEO3NxD0d2aOG9cUbckEGodGcaBVFmYBzh0Wfj4
8ZB6fE9ZoIkLGv+aZUFyOLQ1zTuJb2+/Tpc7milRsZrH/95YF00aPx018t5F
qFFZ+2Y9IqAXMjhkkqB90Npnw3b10w1apCHsk+owaB0ceB3UuKkl+aBEa0D6
s/WVxO5RKu95+xJSGJRUX/i9PY8dWM2IhmoQ4UBm3OqzZZ7t2By6Fio2rIwR
QAQGxamuxdD07BrvDpAkNlERT7W3o7EGMkZTkavBlE2RRgMfjcPRLI0KH4mO
uh2zc1FkwWSltsN+DuKUlTpWAQwlYP4L1R0ChuLGgB5q7nbzMElb27v2G/lE
JUEtI6KFGyrOoaSiLwuLQO39HZf9Gg1mlPkYLUlQ8SpDrsNsoUklhwU/k2wj
BWcGXyYkaMyihG4bsp/JlxD1b3YzsRhqFEbWM1lMfYhg3wDWvQRqLy/DbnXK
mo/EddquBBrtD1VxfGGk7GH1oL0aEygd8YjRfYynPgWp/F0Ym96soXx2QBSa
CAOUQJIcPrnc7uVSEd+4EF+9cRnFN1biQw2QaaLg/c1U/uJxykCu6y1tzFCp
TARHXxITgee4dhNjyV+ofl44L5pzp5/8THmoS/Gr013sx3CNLTTGXdSGbNTg
yHMnqTbVCVPQNIUc7GeAVQbPSJdWFcT+yRiSa3NlMaoenU5Vs836HzyesB9C
hiq55S6orRFYrWepNxLjWMw4enFYoo9EDEYo9FkHi3UUipao2pl6qHKzQJ1q
BZ1cqnnbVNarz2xyVKZH2KrVeA5JXw8RoZWFATQzD1WVAbArsoxBDy3q4g8X
mbGfi2HxpxKXFgIBm3Y4PVkSzov6qtb/vew3Rd6vyH5T1tuxGt28IZrln31i
tp4ejTH8uHwc/Fx30OpaUyeWnkAO6a1UhNDeSskpLTj0PvomEqTSNxfWA4om
hBP0iWytPueOO7k+G7d9HFAkoXcf9NiVPm3YIHZNx7YZqPGuKep5L61sKoM1
26xLOeCuKftJDKE3J9q3WRx6SyJBPGmzK0lMx7VTOtABHMexAd/nFA40hUZ3
OO3rqEQw2lmKOEETdbhIPNRLmehC4bn5CRkaS9nnA4gbgzjqdQToWVSGOhFk
4l+fyJKCJBfyfBVTfbZ0vorPiRz25oSEUlU1phWQ6updEahH8FDfQ4JcmpzD
QdVCVxJlqdOsOP7Z2bUtgyHFjsSlvFPRUCFAd3shC6CIoVbZFAmAYQty70Ba
Uc3sXxOLdIT/WhaAhjPUm70I8/nK7WEq1X5mRQCInKzmSQjqnw8OfU6hm8gy
UD0pN1H7UeaOc9l3kizqGyqI3efT5llhTS3wCH28JULN+0w85Hdi7Xm5wYbE
Bujx0hZEyuno1gSxFYhHcMZwOtsbxKWzvRPEBNEXItWmdECIR3BLSZ8gptuB
2gAW4ZVVR3Ar1xHcrXYEN+T1CeIBLtYMByaCSzbbvea0IZTILms9FIrY7sxt
emVxOO5bvGO96Q3TWaCWYUEW+o71jq8CP+Ad68LEXtyfLv8KiQrFtXQ7fYuX
qFeb7oGRtXCJelXTFIHqImFfO29dn9L6KkV3bwfSIQqXQZCKo41RooEV+nz+
AIrWfZDrPs2vxBFcrAl9U6hhLooL7SUcqp4CQgx9DHuEEJHbqqLs7iAaDTC7
w2klzugCEiEw5BSjKKsbFJ+2CMdbF3LmeCdrPB/BeGlYNBoQyVpRnkmIyBc1
woZ1IRlGzhV2J5cI4fB+h3tvuK/HdsJwmFliYj1u1r1VZ2iGW1eTqXEk1hPr
SZv0U8KJ9cQ6lXW0+MpWG/ejM7uP27tQbtyhh9CGzhd7GbHzl8rougy52iSM
rNJ4x00T+gnSuu1NniCdfef7kiu7aGuHIetZDWf78iz5tD6n/D7/zZ78D70/
q4Q=\
\>"],ImageRangeCache->{{{0, 287}, {215, 0}} -> {-4.53044, -3.37249, \
0.0302733, 0.0302733}}],

Cell[BoxData[
 TagBox[
  RowBox[{"\[SkeletonIndicator]", "Graphics", "\[SkeletonIndicator]"}],
  False,
  Editable->False]], "Output"]
}, Open  ]]
},
WindowToolbars->"EditBar",
WindowSize->{882, 964},
WindowMargins->{{Automatic, 412}, {Automatic, 0}},
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[CellGroupData[{
Cell[579, 22, 3021, 92, 1374, "Input"],
Cell[3603, 116, 323, 9, 31, "Output"],
Cell[3929, 127, 63, 1, 31, "Output"],
Cell[3995, 130, 45, 0, 31, "Output"],
Cell[4043, 132, 45, 0, 31, "Output"],
Cell[4091, 134, 46, 0, 31, "Output"],
Cell[4140, 136, 46, 0, 31, "Output"],
Cell[4189, 138, 46, 0, 31, "Output"],
Cell[4238, 140, 46, 0, 31, "Output"],
Cell[4287, 142, 45, 0, 31, "Output"],
Cell[4335, 144, 46, 0, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[4418, 149, 1311, 26, 201, "Input"],
Cell[5732, 177, 134, 4, 31, "Output"],
Cell[5869, 183, 10620, 688, 224, 6112, 611, "GraphicsData", "PostScript", \
"Graphics"],
Cell[16492, 873, 134, 4, 31, "Output"],
Cell[16629, 879, 134, 4, 31, "Output"],
Cell[16766, 885, 10116, 680, 224, 6092, 611, "GraphicsData", "PostScript", \
"Graphics"],
Cell[26885, 1567, 134, 4, 31, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)