(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 6.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     30639,        606]
NotebookOptionsPosition[     29594,        569]
NotebookOutlinePosition[     30129,        589]
CellTagsIndexPosition[     30086,        586]
WindowFrame->Normal
ContainsDynamic->False*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[{
 RowBox[{"%", " ", "Ellip", " ", "comparisons"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"%", " ", "MSSG"}], ";", " ", 
   RowBox[{"Jan", " ", "26"}]}], ",", " ", "2008"}]}], "Input",
 CellChangeTimes->{{3.410366080260127*^9, 3.410366086026516*^9}, {
  3.410366572257091*^9, 3.410366580872253*^9}}],

Cell[CellGroupData[{

Cell[BoxData["ClearAll"], "Input",
 CellChangeTimes->{{3.410375088383728*^9, 3.410375105431608*^9}}],

Cell[BoxData["ClearAll"], "Output",
 CellChangeTimes->{3.410375105769997*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"Remove", "[", "\"\<Global`*\>\"", "]"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData["e1"], "Input",
 CellChangeTimes->{{3.410375100892053*^9, 3.410375101051083*^9}}],

Cell[BoxData["e1"], "Output",
 CellChangeTimes->{{3.410375101587072*^9, 3.410375106829569*^9}, 
   3.410375301905866*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"e1", "[", "e_", "]"}], " ", ":=", 
  RowBox[{"1", "-", " ", "e"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"e2", "[", "e_", "]"}], ":=", " ", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"1", "-", "e"}], ")"}], "/", 
   RowBox[{"(", 
    RowBox[{"1", "+", "e"}], ")"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"e3", "[", "e_", "]"}], ":=", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"1", "-", 
     RowBox[{"e", "^", "2"}]}], ")"}], "/", 
   RowBox[{"(", 
    RowBox[{"1", "+", 
     RowBox[{"e", "^", "2"}]}], ")"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"e4", "[", "e_", "]"}], " ", ":=", " ", 
  RowBox[{
   RowBox[{"e2", "[", "e", "]"}], " ", "/", " ", 
   RowBox[{"e3", "[", "e", "]"}]}]}]}], "Input",
 CellChangeTimes->{{3.410365752678855*^9, 3.410365976455923*^9}, {
  3.410366021493353*^9, 3.410366023139984*^9}, {3.410366114169287*^9, 
  3.410366183125592*^9}, {3.410366323765868*^9, 3.410366416123093*^9}, {
  3.410366451527303*^9, 3.410366458468592*^9}, {3.410366533214636*^9, 
  3.41036655798258*^9}, {3.410373834640908*^9, 3.410373886712046*^9}, {
  3.410374084387968*^9, 3.410374092336384*^9}, {3.410374274707078*^9, 
  3.410374297517956*^9}, {3.410374423181376*^9, 3.410374435147748*^9}, {
  3.410375072684499*^9, 3.410375085727875*^9}, {3.410375305435964*^9, 
  3.410375346907455*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"p1", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"e1", "[", "e", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"e", ",", "0", ",", "1"}], "}"}]}], "]"}]}], ";", " ", 
  RowBox[{"(*", " ", "Greyish", " ", "*)"}], " ", "\n", 
  RowBox[{"p2", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"e2", "[", "e", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"e", ",", "0", ",", "1"}], "}"}], ",", " ", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"RGBColor", "[", 
        RowBox[{"1", ",", "0", ",", "0"}], "]"}], "}"}]}]}], " ", "]"}]}], 
  ";", 
  RowBox[{"(*", " ", "Red", " ", "*)"}], "\n", 
  RowBox[{"p3", "=", " ", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"e3", "[", "e", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"e", ",", "0", ",", "1"}], "}"}], ",", " ", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"RGBColor", "[", 
        RowBox[{"0", ",", "1", ",", "0"}], "]"}], "}"}]}]}], " ", "]"}]}], 
  " ", ";", 
  RowBox[{"(*", " ", "Green", " ", "*)"}], "\n", 
  RowBox[{"p4", "=", " ", 
   RowBox[{"Plot", "[", " ", 
    RowBox[{
     RowBox[{"e4", "[", "e", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"e", ",", "0", ",", "1"}], "}"}], ",", " ", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"RGBColor", "[", 
        RowBox[{"0", ",", "0", ",", "1"}], "]"}], "}"}]}]}], " ", "]"}]}], 
  " ", ";", 
  RowBox[{"(*", " ", "Blue", " ", "*)"}], "\n"}]], "Input",
 CellChangeTimes->{{3.410365752678855*^9, 3.410365976455923*^9}, {
  3.410366021493353*^9, 3.410366023139984*^9}, {3.410366114169287*^9, 
  3.410366183125592*^9}, {3.410366323765868*^9, 3.410366416123093*^9}, {
  3.410366451527303*^9, 3.410366458468592*^9}, {3.410366533214636*^9, 
  3.41036655798258*^9}, {3.410373834640908*^9, 3.410373886712046*^9}, {
  3.410374084387968*^9, 3.410374092336384*^9}, {3.410374274707078*^9, 
  3.410374297517956*^9}, {3.410374423181376*^9, 3.410374435147748*^9}, {
  3.410374506025864*^9, 3.410374506455927*^9}, {3.410374574054235*^9, 
  3.410374597343533*^9}, {3.410374647530468*^9, 3.410374650137961*^9}, {
  3.410374688304855*^9, 3.410374690864213*^9}, {3.410375317480659*^9, 
  3.410375328476749*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"p1", ",", "p2", ",", "p3", ",", "p4"}], "]"}]], "Input",
 CellChangeTimes->{{3.410365930298932*^9, 3.410365991352042*^9}, {
   3.41036642056925*^9, 3.410366432272588*^9}, {3.410366474255439*^9, 
   3.410366474700268*^9}, {3.410373869854572*^9, 3.410373870540362*^9}, {
   3.410374071653224*^9, 3.410374072017476*^9}, 3.410374389553131*^9, {
   3.41037460164993*^9, 3.410374602174541*^9}}],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJw91G8w1AkcBvAlLkeD1RVLklK5uBRKoR6tvRSq829ibo/qlKNsmC26uWqZ
HGEvXchFtSlzIrby5xRFbFxuixDK2mWttfav5mwo5fbXi3vxne98Zp63z2N/
6HjwYX0SibRXd8T3Pyx72TwRtv2ynYl2fl4DA1/rBBHFHfVzZN4vHzVwu15Y
XE2hYeC11aXZ9xoEPO+/f5USCkremo3/ajXYHz396HdKNIqNqPFShQZDLhm+
bAoTnLcpY/x+DThL9v2aRTmH8lbpqz+4GqQaVa89T8lDY0xLnVuUBuERN7iT
b28hTm6oHb6rhsl8Kt1vORdmwx5ZNqZqrA5Sy1qzazC4SZyWwVAh+cVWFE49
wMY2q5ljT5SwrCEL9fSacPrZNp7ITok6sS+/07sFwQs1i/gJCrTEXbWJUPLw
XanVgwmeHNo128wq49vhW7ehMc9BDv3ZCnOb1R1o1nvl9eDkBDbPH6gLqeLD
tTnG4mivDPxz7i/a1nXiuTApcO8mGexyFbO31V2Y9DBeUJE5DsW0SPhqQzeM
DXLN88eluCVv6DD9rQdVxiVXZrZIYWr7p0Yi6cWpROdMNnsMObG1RzgOffAM
KOVzVRK886mQ/HS2H+lR5ewAHwmYrryapa0DaJNeLWVeG4Woipf90v4Ncsts
PS01YnyI6y36+udBpHUscVVQxXDasfDTkQYB4qLmWIElI/ArsY+5ZymEYP3x
wLVTw0i2jVo5QheBltjFuBA0DGdlmEn/nmGUrQ2dCIoWYbHZpBHv/TDejNip
k42FMGgq4q4qHsFM6t+c+nYBjro3vP4yQIx1l88HKmMHUbH7EOfxOzHKj3k1
pVm8gVZhbS4tGMVSxw6RW/cAakP8FuTRJCgIT5ztP9mPs7vSr9dIJdA/rm0u
cO2D/ZD/wM30MTz+Nm5f/VQvPD3O5VFcpci5ON3HL+oBfXnnrF6PFLtaDy6a
i+iGtoMtTkgdx9I4iefyj10oNIhkBjvIwFFlrK9jdGKxsLZEnydDTYbjCa2G
D8PTqaFbEyegNHlfae3VAdLtSK7YQg4SPZ16oqYdSVe4bIsWObYU7HhWO83D
I6dMIzJDgUVBoqEGvxaQ/ZOc00yVSOHGF6uzmyD4pqvnVKMSKy6eObhx2UPc
vQYHSaQKLOMFQqFRLX7gbibbflLBKt3/pKXZXbDSb6QKytUglY2PeQSXolJm
W2QaqMEBObXY+UwBgsimW9p1PSI1Pmy3G6Uhy1DfRkB4YGL60gUaeLPaj5OE
J60cjbxp2CwearX+oLNjcuZkni+sqyv3MQgXuu1+spMKceiemCVzRD48/2ko
wCzMKfhR11vSAe+igFF3VOawTqUQLrzywc3bHVIWk84m3DXz/bJ8N4TH0lf+
Rdindplqpyu2eTpXGX/SeQPd2zbMBV8I/nl6jzAH1JujDkBXU1kbYfM5fWb+
KqTwqrMHCbPqW2l+KyG/UxRkqNsRUvOhLH6YHV6cPjoUTthnZ26w5CssTIpq
jie8ook85WcBnyMhN9MIkxLWOO03w/29XrF3CLNCnO+NGUBBdQl88tlZkbKI
+e0OHqtc+j7nGdEM2dT2SCdLC8Vnl/TKVJL/d0w3bd3Erv0Hh/Nydg==
      
      "]]}}, {{}, {}, 
    {RGBColor[1, 0, 0], LineBox[CompressedData["
1:eJwV1Gs01AkYBnBZUmwYQsNOFjOhsbFykDKPLa1Zl0q6UFKrct3K7RzZatHG
Ftntgp0NNdSezbqMkLYUk5l0jMlt3O/GmMYM/rqodGH/fXjOc37f3/exCD2x
86i6mpraNjKf2/uoop0/uZvlAmPl4iIBjS2mMSNUJ+w8Zj49/57A+huc/Cqq
J9Jn391/PUfA51lPZQF1FwxXHKyYUhHYe+TtoyvUI9iX3KQt7yEwZP/blixq
AjxGun0XeAS4RtvTM6jnUL60o05wiEDqsirrC9RsdB9N8OcIZxAYVMibfXEL
RVnjjtTT09BZTA32Ws1DPV6pPHdPgeE/oxBkVsNQ2ZHJZ6mQ2LIBnNf3Ee95
vUFopoRJNWV4yZJ67Nf2i6YaT6JGukXcuqkBoRHVff9pKNAQVWAWNCWE60m9
Z1T155hb465XduwpirTOx1lQ5VCfL9E3Y4gQZTNdrM+YgPPioZqAcjHKi3XG
+9kyiM85tTSubQVTIBlqShyH+SXVfPFMGy71Gcw450uhejsy3OXQgaicQnpX
1xhuKWtFur9LcKKXvzJTcwy6tH8ImawTjsMDn65gFBcj74Zx6d1wXZP/4OOq
EbzxKJFFJPfgwGELDyPuEBIchdXGgl6sbJbEl3gPYqRcmNlu0Y8XruHbL4z3
40NUZ57tzwPwNA5iRWT1gfmd1kJY7SB2nJkuW7e3F15FFuF3TIaR2l8/uk6/
B4m0g5ZjwSNINnEL9H3eBbup3To9fqPwqAxN6mjuhKHe7DLh+1EY5Q6fjuVK
oFGfx7PKHwN37Tat6EsdiHaq7VvuI8Xi1l8sfWLbUfJDKLfujRSrrDatl4e0
YU5lqi/PHQffQce6WdSCuwFeX2R7yvAislHAqBUjmZ12o1oug5/fy7TWfhEs
hrx7b6ZN4M2NpQVvdZvg5nIum+ooR27sQKPC+imCV7fOL5HI4Soc99j60xPM
ibKkManP4Rmb2lKRKwBHIyRhJ10B7VoFLUX5GIbDd4vUhQqoHDhVmjv40DyT
umtD7CSa3P64XLPhEdSKQ3hSAyW+z+hfe9/8AeKu8bIMGpRopzUkXemuwSPm
+WWU4yqIZYzHHF4VKN5xdmd1p7DR8mEYtbUCg9+0SZIeTqGJYmb9468lqLgO
uixkGuH+9feuEn/jAM+ZQluYxrvlcZuds7lISStMHfx3Bq7xwfHsxhyUKWh5
ur4EbL9kf726Lh3+FF3Xp+RfWTW3fNspc0KGprrZIOlyzfQ+t8tOEM7PfZol
zfTtuFfo7gRn6ZDA9AMBm422q2Jy18O0qmz7cdKnbovC9diOkO7yCzf6SJB3
GJRWscceCZyLuYc/EVje9lLEkNFRdjEl6STpAJvuOjV3OuQpCcFZpCM8u80H
c6wQGBlseY90wVdM22wvS7i72ZVrLxAIZrpUGuwxx9LB5id3SNvcSdivL1sJ
tNXfbiSdVHVCLskxxElhVeYA6Vlprg/HywDK0jx/TXJHyibEbK89emg5Ez0U
SNrw1YT5BZkGtOIO8o+RnteL6aWx1eERFnDzLGnK/oyMiL2LrMptGyNLSb88
dTVRPPGapdps7/uYtENUj5HWPoJFd7Gy7yZdavfXitMKGSuEaWKgIn2NISol
pjtZf5rrzH3esTUh/C7+5FbW/5gVdvI=
      "]]}}, {{}, {}, 
    {RGBColor[0, 1, 0], LineBox[CompressedData["
1:eJwV1nc819sfB3BJJcrM7/Y1y2pJblYpXoSI1DUyyhUqK6HSNVq4FEWX4kvI
Voi+ZsrITteWkcqIvr62DzcqGf1O5/E4j/N4/nneZ7zeW+3dTc6xs7GxJZD5
azU4N9peOXZC4/vPX4MCh7awxwBNCcMJS7G/rJgYE19A08HE2YXzK8RMUWl/
L9lj8HH6emmJ2LD5Xf4jmhlYLV8YC8TSoWLpPyWsENg0K/mVeMngTLSdrA3C
0mYaZ4nHAztr7eh2UMicyZoktjj7rfw+7SwiDvxXwSKu1RYurol3QGXjPNcn
YgUp9dw5CWe86lm820Mcz26bKZN6HgEzHNptxJxDASnmsm6YZvAr1BOXfHrb
XGPmgZ5qSdNyYo8/zhpY0y9iDXNfdj5x357b2mE0T/x3y4stjjgycXOlzKkr
mH6TvBhGbMCbdfBV/F/IfNgp50f8fLpJmZLwgfgGy9/tiV1tbPJv2/niwXAW
hwmxZAslvyX1KtTLuXi0iO/lCGw3lr0Blu/XalFiXbG0tDHHm6DxBbmtJ77x
OpavyswPrFWyf8yvUFBSkQuzovtjYjmtqIE4Sej4rTu0QLA/q13vSny1ObrW
vCQQXUKCDibE5kGf2KVOBYF+xv+7KvGG+Ys3S+NvQTyiroKN2Lsz0ntSIgRL
MbI9wcsUTEP7nr+oDMEjuWfDTsTyOrLzgXZ3cHrcYYc+MbOg+KJY6l0UJF8x
5SA2fvDB5ZjsPWzwOL738hIFOUOpTOH6e+iyfjJ7lJhztesIy/EfrLqkOyJD
fPrlSl6FWTgYjNNXuhYpCO0Il7ekR0Bx1lZsF7E/Z8G2EFokIredT8xdoLBT
LyI43yMS6uo2nj7Eb4Pcxz7WR2Jez8dai1h6tdxT+b+isNSk5tDyncLr5TT5
zrd0JHL42PZ9I/v9QleWCH2I68zwuZZ5CkW/X6Hrf34IxSfeK8HENh6m3y6q
xSJjbEDkEDFjivdl7Wgs3sbrRuTNUTAZDT7oohuPrGVxx6AvFGL6fLWLVhIQ
mt9dyzZLQUvUKq3/RCIS1erWZ8+Q+3pSdQ1nTiKsNRLtzYnFJ33O5OxOgqy/
4h8ZFAV3ZxUZtt4kmAzmiWhNUxA4y8hI35cCPckZ38MTFCytkhkzs2koez8R
lcSkQH0040oRTseBq9fYVIhv2XCeM9VOh9zmnQkNnykUnHEXfh6ZDoUXS5dm
hyjwuKkHXVV9jB2CDn8rD1KoC3h/cu31J4g/xTgc3Udhbzb/WlHOLBjFukZV
dVH4V77OtlkhC4wsT091Yts879IbVln4YFVy+0UnhbDnAx6DWVnYniRrmN1B
YbQy5+Njo6codU44fKedQkLXkbzfH2TDS5lG29JMgfunv7WeOANjm6q8CmpI
vZwtbOQOM/Dfgy/cNGLpTjlbfjcG1A886rteTUEj892Zj+UMyLwJ2K9bReFl
Jr2oWD4XJzQ7OhteUZgynSwzbMmFSsrD1LyXFE5kRDde3piPcNW181PPyPv+
M+Cow9Z8dI22HDtAPCRwodlSOR9cA/dGbueQel8/1Kr+Zz5ubdveLZ5NIcRk
6u26nHz4ze3K1M6kUL506EPc0QIY1LHvdU6lIGM8PVpztxD7ORZ3S0dTEL6p
p+aYXAiuFzsM9OkU+HKS7nIVF2Lk4r6g81EUfnCayhsPFSJVT9OS8YBCW2Xx
5f59RchNWeMpF07BV8F/5ftwEZZyVz9YCqbQwrtpk7xmMRpb7b/N+VDwatmP
mLmXoMlV9BhYkfOTLxr+l7ME9zWWC49YUkgLUwhdFC2Bu31Wu54Fhfaj297b
6JbgKT2+SesEhV2Nmy7L0EugHT7dIGdMYaCeepyvWgoWbWihX5+CXlU6T7Nv
GfiNLmcUqlL4rZC/f9WqCrDWSubsECT/1ZB2U+vBatjzbvFnL5lGtcsjEavJ
WhSxBv+OPzGNeVl13pwL9ZBcftTNmpwC+8JTPhGZBrjuCg1ud5+Cyk/b56bP
mtC56ORRMjOJpkClltc7W7F5zGP7osskJMInFjKn26DSJaOp0jaBiW8D/V0K
b9HkNEg3UZ1A2nhpA8+9DrgwWofWhoyDR+wJxWR2IlImXWymawyhzkUOSdLd
GJUJyHHePoavmk+ZTjffYank+eM350bhube28H81PZCNNbYoyhrBwLPau+1b
P2BOaEXwzQcWFl0643b4fsSWey7+67lY2KW1bsWhtBfxO/tz67SHoZey1THv
t34MvObmjnJmwkvstOSg9QCSLuyg2cV9htzkCe53Rp/wzmrkmH3rEAR5Zzhr
f3wCy8goO2RmEBwVcQyp+EEUCH7/Hi0xiPNKpe/XGw5h3QGF6LCDn/D0iH3S
q69DaB4pVuct7cf8hDAfi/4ZHuxxd8JF+lBkqrc6UocJbsnP1lYnP+KmflBi
IYsJxUdzOfMv3mNrn0FPatAwkkSTNPdI9kBNNTCStpeFCoWAwzyu3bAWb11Y
1cGCKl+gWVd3J+YbwoY8/EeQ76q7kZLqQAyHjaeJ9ChuHLsmfvp+OwT7i1LY
a0eR/JfRBv2NbVhz3d9s/8UxcK1Ou37lfRPYMm0YQwLjuGWi49j04V9cimWE
CVSPo0/fNdQp4zXKdwVz8rtNoCzhfez3iBrwG1ySC+CZhJebemhSWSV6d7d1
+JRN4pqnrjvNvRS5CZBm2kxBXCux7gh5B38yVPjFVqawar3o/tp/cuEXlOzf
mzWNDmWLGzU16cgZFYvjOUrBSjTUyqGcDmN+nn31P8g99dvN08nUwZ017CK9
xN+qQ4rVInRQuzC/PEPsNnYhOlldBypDfTXCJEfoE80/3OnaEC7IOe5GbFA5
eZxX/xCGzIwchUgORXgbOuaaA54xofQzJMdETnFXyzCVkBPq5+NNPNBW/pxN
XQksP0/rMOJrQn6be6MUYelsLVlM3CyiJBWptxfqanLPuEiOKrIVZgmY78Ha
3sa6POLbfWpWfExpoK0i4zVxnWwssyNKCt61BXc/EvdtoRnE6EliPDvOeA3J
+c18k4f1zCXQcv18nyXxK18LiRDmJqy7dLryAvHjcL4eMX0BaDqYpgYQ+wgK
3HGy4EX+sQPO2cRVje5eTcMcmDi052gVscvZl0LrTv7UkFaV2tNNbIfwjddG
5zRsdv0mMEGs3F6fTU0xNaIluOd/9W2GZX5dv/q6/wNJexjv
      "]]}}, {{}, {}, 
    {RGBColor[0, 0, 1], LineBox[CompressedData["
1:eJwd1nk41dkfB/Br33cxl5GUTEKLSstMvTNJoYhKlLKkaJOkyHSLFvVTtJeh
BSU1QlmKqCTLIyQhzQhf13XDvd97KmNJYs79/XGe87z+OJ/nnM/nOc/zNgvY
67FNlsPhPKRLujtv62ko7V2/BDDoGx8nkF9mFNrBnYvhfl2PkRGCObcSr+dx
HWBcH716cIDApa4l9wZ3HbzSD45+FRFsCBx6dpEbiMI4R6vxFoK2maeWxXPD
kdpg4B7wkCBlgltsHPcErtwIkbgFEMQo5/3yP+5lPD3zNPev1xJ4eafmfP5y
B6qedRFBV1iojcf4rJiYg7w3Db0aZ8SY6i7peXUmH84Wq37M/FOEiDcLkfhv
EY49ielLvtsHw3yddhmZF7im4cVzae7FY/6y2vrfyhC7K3I0Tq4XZTtvGHuL
y2HhqXrj2PIeDFgs1sraU4WQqM8f16R+guy3TG3jqa9xcH/t9uBRIezG/R6v
za7FrDKPZ617hKg9MfdN5fR6FPl2aXPE3TA9L/p2X/IWUZv4s312dEM01NHe
POsdGo8NPdg6LMCdvuLXmgmNcCFk89dYATRNMohA0AQzj9vpypMEOLujYHuK
+XuEuyaka1V1YXBppiD4aAuaeFzRvT1dCLctzzd49QH1E3tW7prUhY7s8jMN
Zv8gP4SdJNvIx/edTcmWUa2Ye8o6UnyWDyt7pbHtxR/Rae+77+AaPlakmQU9
MmxHtcWKsF49PiJMfCd3+nTA9uSGmFcdnbAWr1drWc2g9Fvo9D8edkJP67Ny
+QiD73uvbZ8X0wn5F8k5U653QnUg+WSKVyd2zS3+W8WFD47juMmSWZ3IdApI
eT7IR4a5w7/jap0YEBlpC692QWZzorYqy6Bg7Qq5yw4CuOKlw5Y3DI6uPHkr
XyjAhcC1Li75DMzanD/cPtkNwzErt+/XGSyaf+Iy11aIS6KJx9VjGfhMrP8m
0yhE0uNqp5IwBgOv4/mhMZ9Q2XS5hPVlkCi/JdzDvAeabPK/S9Yw0GsvSJMt
70HI39Ehsb8zUODFrFu4rxex10YVAuYx4NzfksPX7QOsHnllTGcQlpQTr1vW
hzqXgDV/TmLwzOq0sk6ICIdZ9sKCnxjYNRnWXQwVoYLz7ocx9cPDGRf094vQ
aCzHGzdkkFZXafRTpAiWxosKKqlP71WwMT0mwr1Nf6h7UnvkHXO3uSqCN8/N
I9KAwadFvCSn5yLYvLthXKPPQMc5zPqYphht/JDMezoMUl5qqS3SESONc6g8
gXrWwqzeL3pi+BV3fw2ndp3Wc9efK8Zjw4Zz9tRxilvM7M3FKLT2rG7VZiD/
ysmAs0iMsUvr7PSpB3814xzZJsbqvQ7vLmsy+GjztvFQiRhBzpmfGtUY3BLp
8OReiBE+JzehlDrg/lqL+JdiHGqduyKLute8JTKlUgy9vOviWOpB43aTqrdi
jJfHKP9KraMqDtLvFiO1fyz6rioDR6HSaLYGC+aLcEOcCgOVO07pC7RZHLD1
tzhEXet/xrVMl8WJyEaFYOq1bZqpzYYsVnWeG1pO7dc0wXHUjMXC/sF0Weo/
yqacd7JjgYqh6UeVaf9vwlywhUXbjIqiE0oM+qP6lfz9af1xD9EBarsNGaK2
rSzOKhRYB1EXa2nltQSzmCsQtjtRV0a329eEsXgflcTRpm4N4PnmxrIwCFW2
u6XIQPGX4qTobBZpD4LdqhUYOMntPfLjIYsRxYDwEuqzHZMDovJYbGt8ei+H
WicxzjK8kMWUwMO4Sv2zysbCoDIWAzKt5oHUtqLhZtf3LKw6u2tkqTfn2OmY
jLHw3OgT5iLP4PHGYc23HAk4frkJoNZUfKp+XE6CXKPkkjnUpT6LlXuVJfBn
dFf/TD1Z1WE8X0+CTY6mSyVyDISB7uwqSwlUVWp3XqEOMdpdfXidBOXRBXmf
Zel7K2yqZm6QIPhSyzchtek+Us73liBIm+PeRv22Kqx0pa8EsWprbV9T2x44
9ER/pwTaF0+X3qEeqo9Nf3BUgvaOfUqbqaNPpsZ8/EuCn/UjzP6RYaARoeTs
nyXBJO/Fso3UScF7dIU5EkSd9/tSQ53nsvA2yZdA6PBw5Bm1UKehXPaFBK7D
E4rSqF1ucpQtGyVwrFfkh1B/OLf9bVazBOltihnB1IExtYm2HyRIG7SJCqA+
Epho+VubBFf6upw9qR9On+Xi9kmCFQcvXFlCbfDEN+HAdwkSElXrdKhv36vw
HPkhQcGQKEydelaSlelRDoGn4wsLJWon3lDOaQUCjVMzCkc5DA4vO9eQrEUw
+XJjTA81v/65ftkUguEtG7PKqPe+NG9ztCA4vKJh1nPq0dy49JppBFl1CmVF
1BOueto12xCQoiSNR9QrfSQbeuYT2O/mOqVQZ/WYJGuuIvh64ZprDPWa+xd4
1a4Ev8QNbOFRf92h4HfCnWCmx69HD1EvEInNRzwJdquFkH3UZWxxdrc/wavQ
sfit1IHZMy+kBNL7zLA08aNW2nt7/6YgggRn4xIf6lWf4xY07Ca4E7nc0pP6
/VfvspIIggDB0lIn6si8ujsRUQTC8nl1jtTG4fanbHm0flND7zJqv4Fpq+4d
I6g8QtyXUPcNDTVfiidoMXJvm0N9tnBXoet5gt54raDZ1DMPdSSpXCI4sGw9
ZyZ1+Eil79FEgqKL2/2nU4+PXundeZtgX+XsNjPqtGcqtVPvEnx0KuabUjsc
4WUz9whMAx8MmVDHjQfu98wmqEkqcDOiti79sF7nEcG1DxUXf6Kuj161oDaP
4A2x6Tag1ped+8O+iGDJ8MATPerHZRkdo8UEzbquC3SpvY8blz15TnDKW1Kl
TX1TXu6UTTntR/WQvib10oqDO3oqab0p3u/UqbtO9rncribgavTfVKOepvRO
m1tP4DNSuk2F+nXV8v7GBoLH81lfZeo9p4uaE5oIHh38fYcStbaTTaETzUHh
70tiFKlzVVKT5P+h5yO87itQr3utz3vxkeB6lA4jTz0Ud9o3qoNgqqxgqtRJ
LqP28/gEhrNro+Sof1MPNf8soP23qmyXpW6v7VLM/EQAwzfuUsfEb+jd1kdQ
ymUaZajNXWtqJrEERhuHA6Wu1ER2KyFYw9FTkHpHfe75q18JxPOs8znU6uct
9rvTXGc5dWmY1NlrktarDxNEDDlDancdzQVVNAdmtq40lrr828CPz9T9nGSu
1Hb8tldG3+n8GkYMpTbKy3ILoRYEdelLPXac1/qSOvK6rK7U/HWrgyaMErCh
XlpSV0416Q+mftEhUJc6c1B8pIR6fm+GqtTnq0pUtX8QvD2Tqix1eOLZq1up
NSobFKVevMg6W3WM/tf+ZjmpFT/WVDyi7tmVLyP1G96uNi+aa6cY7+FIfc1U
bUCacy0Ozvi//wNWDWjY
      "]]}}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesOrigin->{0, 0},
  PlotRange->{{0, 1}, {0., 0.9999999795918367}},
  PlotRangeClipping->True,
  PlotRangePadding->{
    Scaled[0.02], 
    Scaled[0.02]}]], "Output",
 CellChangeTimes->{
  3.410374395243453*^9, 3.410374603011682*^9, {3.410374652953019*^9, 
   3.41037466113528*^9}, 3.410374692921341*^9, {3.410375332240746*^9, 
   3.410375349055377*^9}},
 ImageCache->GraphicsData["CompressedBitmap", "\<\
eJztfXl4Vdd1r4xAYAcPxUmc2o5NqV+a9nmK46RpmrbpS18c5z2/1HbapKFO
jB0XY8CAmYwNBk+AjW3AzGaewQLMPJlJzEhCjBIIJAGS0IAQmkBoOmt3/dY+
+9x15XvuldPve1//yN6fpnPPWWvttdde01776Ikeg158oW+PQb2e79H1sVd6
9H+x1/MDu/603yt8KfmGpKQbXkxKSirtmoTfDf/qf5PWC9/8P77E70/Ij69Q
QUGBefbZZ2nWrFnyMz8/H5+Q53mmW7du8vucOXPMXXfdldTO/OAHP0hKNg88
8IDxW9J35XunWHAIzYfTTm5rTw0NDTRt2jRzww03UGpq6vfk8q1touLUqVPm
hRdeEGr69etH6enpQOBT0DYYS5cuNf3796d58+aZgQMH0pIlSwTGowLjFiov
L6eZM2cKeTNmzKDS0tIvwIiBp53maQfz9NNP08aNG2X8W7ZsMU8++WSSvaUj
ZWZmCpSmpibTs2dPqqqqkvvq6uqSvhfveaEz2T0vd1qYyUItN/nr0YQ0dApg
PPbYY3T48GG5p7Cw0Dz00EP49btfAsStt94aDKG2ttbcfPPNYNb35bbOQhlG
6je5zT0aA3wrPnbp0oUuX74sd1y5csXcdtttSe19QUKrqakxw4YNwwxYhCmx
HsGvPqdkHu+44w5ZPsn2u4zccfSRhJhTTOfOnenq1atu1ky7du0A97ttePTe
e++FeNH169dp7ty55pvf/GaEW8maW/4jscFFsSlFCMC6QmtubhYwLS0tHfwx
NzY2mu7du9PWrVt9CQt9xAJOkcuYmry8PBEsX8ASPJ2sPp48eTKYZEaNGmXu
u+++iFglelpEoWPHjrhsUlJSaMOGDbgznEehEG/4ojRVVFTEmJgOTlPJz9zc
XHP33XcHvIrxHH51gjh16lQhKIpH8bClyLQvWrTI3HnnncCEIet1F+9RLNHp
06fLcPmnefDBB+NzJza4Vgvtl7/8ZdRSfOqppwLWHD9+3Dz++OPBguMVH7Am
xnOaNdXV1aIlW7EmHFuK6du3L61atUo+u3btGpStZk28R7Uq4p/yN7dHlcLA
8vdZ1CEOxGju0Llz58xzzz0nQ+nRowfUTcAdiBtoHjt2rBk6dChsm+NOrOc0
d9DYhomQR7gTB1uKDGvAgAHy2fPPP0+TJk1S3In76IoVK8TiOcv36aefau7A
8L355pvm9ttvV9yJCdHnjihM35noYsUv4pTwX/L9l/9/Pv+KvcTtGXy7Ifyp
m2Ncu6WN1/747H/t2Rvle4pdiIcOed6ECV7a293p6XUP4b5O8vEfZ+2/2bMd
5XsnYkNANGaM5w0b1lwxaqj3q32/9f4ja7jX2Nz0x9n7b/us9ao6E/tDRHv2
ECt5Y157jWjxYlhJej97Fv1o/69ajpSf9mAdjVWjcWbzH2Jc+3Ebr8V69juO
QjT2i9nhM2bECGOGDydat46ovl40xtGK0+afjnSncTlzQGlyDFAp4VTL9x/a
u+BZEHsZYIRhD5Vo9244jn6MCommnlkj6KmsXl7FtapbEkC1vlcn8gMCouXL
DYclYDXRoUNQdpEQpD1o99bmpZkfH/lXb23hztsSQH/U0VxTQ+w02sl7+20i
G7f55l2gVtRW079k9KWeJ173GloacSkR9IcddDT21mjBAkv7u+8SHTvWyad5
4uFl5ifHf+1lVWUDapcEUB90UAFx/nxL85gxRCdOAI+DmlNcaH6W3oPeKpgI
qF9NAPWvNK3g85IlFvLo0UQnT97kQ71+vYkG75hET+Q86xU1Xvy6ghoB+m0t
CewZBVw9cgTQHSy0bZnHzU8zn6HZlYu/kYDC+zSFzc1EmzYZ8/rrRJCxQ4c6
q48vXbpKPTe9R/9e2IsquXNLBL2bexzytH+/MSNH2kWyfTuWt4be3NzizVm3
i35+4re0umEdhnFnAuhd3ePZ2UTvvGOlYPVquySJbvE5kp1dRr/b9Cb1uTLI
q/aq704A9W7NkeJiovHjjXn1VWI3FBJ9mw+1tvY6jV+6mZ7Mf5a2edvvSQD1
TgcV3v7s2ZbWiROJSkqA5zaFdOeOPOqxfRQNqX+N6rhzSwT965pmjP6zzywG
yNmpU13Uxxcv1tA7c1bTb8p+T2ncuXVNAP2rehbT0qyqGzmSaN8+PK6hNzW1
eMuXHaUX971NI5pHele9q90SQA8eP3+eOCaxdK9fT9TYCIhf9fl97OhFenfm
Ovp9zUtempd2XwKot2iO1NYSzZljONohjvKJLl92DKupuW5mzcyg146MpxHc
fX5r6BHgWlxp1y6sE8uHY8dwRU/Cli25NH7BVurT0I+2kITK30pA8U3uccjc
Bx9YmVuzxq7KaOgXzl+hjyfuo4kFc81QM5SqqOrbCaCntJYQaHxgmDCBqKzs
G+pj1kbeooVZ3tyVe2moN5RW0+q/SgC9vYZ+9Khd5+BMRgauOOE/sP+CmTrl
AC2uXEmDuV8mSRIlgh6keJBTmjzZ0r1wodWEEeiXyuvMtKkHvS37TtBw7qnc
ud0fDt3TZIPN69db8Rs7FuLoNAGrJlq18oQ3b24mbWzYYgaZQVRGZQ+GA27W
gDGf771nAe/ahVWkVczJE6WYS+9kQZEZaUbSIlr0cDjgRvckfCOoarBi/Hgo
ZwCzgK1dd6K9eVOul+6l08vci7lzi4PgukNQVkY0bpyl+vPPnUegEeDKls25
QOJdrqqjD7lP5u6R90g4gjrNGoDYscOO4v33iUpLnU4/f/6KmfTxfu/YsRIv
38s3L5uX6QSdeFQBFrg1GhzkY+JEu86Rgmtp6erTCwvLMk2pnx6X6VzA/W3u
10lMxaPh9FZpBAcPWuMIx6CwMKI7LUOceJeV1cJAmsFmMO3gTi7ZFRNBpUOA
/OysWcYMGUKUmur0n0YACZ865SDt3XPOWgraKUgY2Q/CEVTI9/YwBJGFeepU
NwW4pcWDhNP8eYe9a9cavQavgd7lPo/m/TAccKlmDQQSVgfUz5qF0dznIziV
U46ppPw8mxI+Q2dMX9MXP/FnHAQXNYIrVyKmOC0toqU7MapGM3dOJnveOb6Y
EiiXETRS44/CERRpBCdPWt/knXfc5DoEVt/nmxnTD1Hl5WuWHLoivOc5+Ptw
BBccArAH2mXwYJgiFyHc56aGW2lJrZky+QBlpBe5nafdtNu8Yl5x3paPJ/YO
mL+e2rYD1i4GqeccLfBTp061M7l9uzVbN/kqsoU+XX6MFi86IisKGBu9RhrA
fTZ3bj/2aYy3IxR/R6l9DNryHW3nztlwECJ89qw1eskc5peKhLHSCDadCrj3
pb7eXm9v0v/yiQrda3H5TbU7E7kcTcpZRwrcWSgaOHHV1d8OIFity4oRsuhm
cin3N7n7+iZCUMg+gE64qs0QudwB9GRryW1qsvZ80CCrOZqbHTlKRzOR1x05
tVRrhplhzjX5pyj+xNk4aLXfEIs/UYRhWt2aRVhJ9G01sAsXqmTajh8raacu
Z1GW6Wf6wcT+NCAsJG1vQ+rYGf9Y1J3Q1F24YANefOF3R53dW2loaBZBh8Cz
W+sLVgeYN5rKfQJ3/I6rAZ1hCXRHZ6zceyw6j2k64VA5O4/w1HlM7aErRGeU
lNTgsuNiBVWYAWYAHeTO7WeWulj66YjGs3Wr1U8IdKxIOjys9ET5sRLUeNA4
hhA9WE3VPw/Hc9jhgR5EGmPgQGuirUA5PLBD69bmiDZnra7xtHAfx531zP8N
x5PhnsCSgJLAkkhLu18Bys29ZCZO2IefTqxLqEQs0nHu3OLAP+QAwYWeMQNK
0vMOHLhfsTGwRmtzZDyt+LWVtsrSu0pX/184ngPuCYQucGkQ7GdnP+hf5hBR
GyMNv4n7WO4LaeE/h8Pf556ALkTQBV4VFDj4xcXVZvKk/ayni/TyL6Ii4VMO
d24avoDf46BevGhXFdITpaUOKq+iwHw0NkbpuI200Qw3w6me6p8KpzrNPcHm
Tyh+7z2oqQcVINYlolNYt2j4zdzBlcW0OM7+zE73RF6e9fMQHNXWavhKv6co
6YF76nMlDvztDtDp03ZGP/4Y8hLknkg71elQ2Jr7aNtpu5Oefw3H87l7Ijvb
2qipU+FxPOxf9n1r8Ko1fCg0f5X9Ohz+FvcEPEoECNOnQ2M4+K2cbccnp5VY
vwNVHPib9DwA/pQpWHEPK0L37TsvrjD72inq8n7aL+EZa6PfhMPfoOUI84x5
uHZNwy+3wSTt23tOwwd/EONMo2n/Hg5/nXuiqAjJEqu56+p05tKPKeFxO1dK
8wn2L52kWCgOnjUOD6I/rDesh9raR/zL2dllBmmCgvxKDd/xCTjYi/2tgi/g
P3P3IX2LSBiQKyoc1NraBjNnNkfCG09DWv3MqJjEj7n7JvG34VSvdPARvkO3
gUPFxY8o8rZtO2tmfpLO/kO9gw+XZaAZCMpxx7Ph8FMdINgApNpgO/PyNPzC
Qut5HDlysZO6zCvY9DF9qJzKnwuHv9w9gbW6cqV1u7KyHlGz29gY8R1Y0+GS
w7OM+xju0Ejc4uBZqiijnTuNGTAA0bLb/na2Hwk9BZ8DHskjrOPO7ffh8Bcr
gsXHwDhWrHA22fGrouKqtv2aX7AHvU1vyqf8/wjHs1DjuXDBBlmffGLnJ4IH
tnL16mxCloVjUD+fHWg8tgy4Ow6e+RoPokVsCNgVofGgnT5VLivj7JkKhwdy
O5H7HO7cXgzHM1fjgZ80ZYq1RrB6Do/9uK6uIYhMfV9A84+1lMR3vNJfSoSv
vZU1BGVId0fKIHZst77AlSvCIDeek3RSVjhrE1zW8AX8bAf1MPtmr7wCydJQ
i4qq7QrJuqihoikJ7hNO9Uz3xMWLdsYR9ba0OPhYFUuXHKXly2SFaPjwYBAr
raE1L4fDn+GegBb54AOb6q6t1aUzWYeLZQzFRdWd/culVCoeTC53bnHgT3OA
sBqwbQMtcv68hs8cF87zDHRWl/fRPvGI66iufzj8KeoJmYOXXiLauzfYOyOL
mSVHJIglCZcsnmSawZ01rvw1IBzJJC2q0ROhJzs/77Iw6lROeQRJeyw586p5
lXZx5/ZKOJ6JrZcEknnIo9il5/DU1zfR++/topUrTsAAaqZxuC6ql1XKoER4
2ttUCtTVmTMavnB/7zkx3JfK6/wNINpMm80b5g2oRtzhw/+DC4BjxWvjHQXw
KEAYNjCiCXPJO/ZcHGHQOR9wZ18ddwzxCfuvVdPG2mX+0NGBtBPmBl/19VrY
tEOCeVJyoFVVJmW+6pP5B1SudtCVqz6YpAh1aG5mCwq+71/Oic4fOuZlU7a4
3b6Ci1D1hxWKRlIrYzQ5sFeIVdetc+TA/8GK3LD+FFboLeruFdyRdGTN+FoU
k9peixlLuMbo1QWjCXcVxs3uf2ghc8aAHSdNGIJnEMYEDg8I+3KVkOGEtYfe
MqZfv9aWCU3ZEX/nVDKnfjofdwQEfcniw1gEvetuvXzZ6roFC7CatZzDNVuy
+IgzPFqi0JgoIY6JfMMSFksPva1FBPUHwFVZ6UQkIyPKNXMDd+KB2eAWB/6b
Gj6CoP79UfHh4F++bF0y7BQr+GqWcXlUOPyRrSXqo4+sq2SXjZtAuCxrVmeL
uLNLdpuiCvsbcC+YX28lwtPebpy8/DICxtYCAlcMAsKuWBf/8jE6Zvqb/pgD
3BEH/ggHqKbGumCwbdHzHcMF0+NYzh2JAebdO+F4XtfzsXGj1VClpW4+EK4o
F0zDh5uPjabVtHq0gi/gh2mokFikBFatclBVkIIbHHcQHqp0w+hwqgP4qak2
jK6s1NyH9KuV0EVRc4SOyAyw1h8bDn+oe6KkxHLE1oTr2c3MLJJEUnFxtb+R
L+7j+9z9QDEO/MGaPwsX2hmOdvEqK6+5leDgg+NvcWf3EXe8Hw5/oIaPUzNw
g3fubO2CrV93Srtgmk9uLk7RqQ/agGfZMmPeeMMG1qHelxsH5gDpEp4DXI4D
f4CDX15uE7hbt2r4SC8sXHDYeV1fVfSrufgoHH4/zScEIUiZVFfreUg/VCjz
XFZWq+ErV35COPy+Gn5VlU1ZbdwYNc9+AjptV/7X/cscpIsrf5o7t7bA37bN
yql1WXSg1irxrOdBxJgyxfuppdqPFR5B09vdB98KVVMoU4xwX6WbcdlRD5fv
I+7+KnjRpwW6GKuUv7yyMs8rLfVQP9COjQB8v9OnbUkFf7E/iO/J3rZtnrd/
v4diHAbAU8/XUlM9MX3Ll+M+fO/ozZzp0cyZck2+xo2Tn7ger3ayl54cuJ5Q
4xcuuOG1yna74bksxCbu5HYg7hQY2MFA0RPcL4wFe7kQWGQeQNeMGTafiD12
dnI47jLy8403ZPcaYYW7Rm+9JXEfCnHk/t/8hh3BV5ld/AczA4+5j3hm8BhA
JMvHHbxXX/UIogxs+PrkE+zvgFvJ3pQpHi1dKpRprno7d3oR7rfz0tI8Nw8b
Nsgn3q5dcfnZU0md1JwgdkVRUnTWBU1l+b+uLpdRmQh+HuVZU9FVLqN8oKjI
1jyDpyiEgyhiZGzazahRBpk8qQAYMcIgN/+F35nH7h5xB/DsokUwSnaODhyw
c4awC/MHUb1yxW36hg/XD4wgOidPth4mh+dus+Eb/uV0Spdkn1+OZbflOsmm
C4Z49KitIIS4IDWMeebhdQTdSTcK5ZAiHBdB6QooPnvW5mfs/hklxcmev6Al
HqNH/vP69db5PZjOZUuPOidSz9Ai7qiG4cV9j1rYIB0HQNeutSUSTHY7kNwO
wim/8Rgh1kiOrVxpq4VQf4rl4gc0cfJfz2u6KyutKdi+vXW+kE1xa3phMqGm
N9AGe/k2kaacHEsCDpyOHy8sFsmAhEBqoOzAG/AahWoYGwonWBHGieEDMrdv
t9rYeiUh6UYnEbC0Stt39sm0k2HrAlARjekaO9YK+zPPWEEGmUhjo2IMUgMJ
tDjjueA9HJmYOUwKQIdmKx2ZaEfpqHOebMr9FrmMCczNRdWWrHdZbyATisnn
Ks2ebZ0ocJ6Duzi+0e/0ZJ8/b0yfPgAfleQ8fcklOe/0FxC83fe4+xGIpGey
suxKYT2aLDLYQbbtwC3Ud5aW2gUTx74GtKD+Bz6arb/V2x3MotabolbSbBE9
P+AdOdZM/5a6gp4ZnU3sJAX6HqkK6B8sepbK6eF0POPoAKsx6/ZQnuOJq+Dh
6AmXuyj8cMgOHiTYR5Fy6EHsSgx/9xr9fPlcb2VWnh9s0xyFX9B31zOB89uY
ifx8h7WkpMaV9eAGG2t2CAQXo2JLDaxJVgUki/5C/qiggJJcDgpCBW+wiqRS
Zl44FwJ6MKnggq0jc5tbuoanoaHZJRThEkFXslfg09FBJBJAILU2YUhRZ5mV
NC0Ip+ff3BOgIzIrerMt2LQ9f0Wnx7GrtXmzLN52mBK+jDh1506b3qfoZLqK
tBaFU/Or1uumb1/MlpbV6uqoEh609v76wRBgAydMIL+U3cgeMU6p8FpyiQQX
HGeT5EPaQg7SXggN7NKJvSOcgd1nXHIn2mEX2E4EpGDdQhGXSyARFYcir+HH
1EvCyfkX9wQcE+hM68s+7C+Vc+ciRZlNTZ6cPIAqs3XQRtwYXLMVy4EHvZK7
ynvEwf+0np2yMusx5ORY/J28+vpGb97cw97YMXmsFjyMWfDCD2Dj6pQw8GCD
wi8AXtYGfPv2WXNp89X3+/wVsWom78ABj/r2rTSP/6yIBg5sEB6zBN6pqOV1
Kfw9Tsc/tfhY7belTC/GKyJiZbae1JzBdj78HkdpCgHwtGkk7uygQY3U68XT
3uZN5+/xn3Kz7xdcOwrbVKQX65ULMWuWfqFJRBFh9+4eO8i2tucmqQXjpS96
HcLKutYrKPhC0vluBQUKBhPJE7rCJzk06exO8oe86kA+lvTuEw5Bbq5d/dbC
fcuXcTSY3w8/jMQOcGx8BRjkq/JsNlqTi2jXzzysCsiN++qADrFePBCTuU9o
5n72mXW8mpp0hSqQsLoUBoPy+fMDTeCsSGCGbTLMiQdCaaQxfOdqVRSvY+Su
HfExXgoQk/j/o4mHPoNqSU11Z48gK6mp5B9tM6Lg2evSpiYtUn17j7p8iA5J
nMBx+eqA6PilgrGO68cnmoVIMkP2nFVw6IhsHmHuXCsn0L0caOgiC2f60w8V
dvUvY3cZu8x+2B8QnahuMMYp+phEP6457dzCsjJHtD3cRmbYMBvf2XQs6cMh
CGeUf4BLjuPYD8RW/Hbu3NZa4mPp1YCMPXts5tGaNX3uAe4UJh2kwL1mr9xt
Ah/JuihLrLCwSqNHO0fn3Dbh+nD0j2ku8ASJYYqghwRyXCPTBofEvpQnQO9y
iju250Um3O4ijjAjXAlvHPT/W6OHLcRa3bTJoUcoiiQGRg43z7eWDr2koiK7
vg59q73CjW1Aj1QadmVyJFGsj/jYA1E2j4LojXWz3uw7e7ZCJoB1XVd1GRs3
sCJZlLU5HP1P9OjPnrU2/PJljZ7XibjYYL4ljvTZKlgCWITPVp10wYIjw7nB
fkW8JkOoCJAjuMSmjK0xiRz3kmjVDB1qJLtjxTLw2/bu8feOL11VSHEUkWZy
R6aUTdHW8LH/ox47km1gbiTWxrhh+pC2OnMm8M/Kymrl3MiB/Rf0WNEgbCNo
hDfJm7QtHO0/RAucDcI+/zw4t+eLIZbbkCFGnDr2ZB16RAQrUo9Ljvj69aZu
/uVCKpSIG4X03NqCvrbWClyuFFd0Uehh/OE0Q+DsqcEo71QV8HdTl7cx1tfM
a1A8O8LR/50ePbYieveG4dPo8/JsjA3G20rvqBOXMTx+PQmu9Hg+zd/VBjJw
jgSlWjZg1Hgw/+AAXoFSV/cFMoB686YoT99xAzWmqLTyK5TjkPG3mhs4pAhZ
aG7WeJCB6N7dhvh2fQSbGgUFlZIzwKk0hR4NhVg4bMi2dnci9MnY6yAep/x1
qw8FuBDrDh5s4Ki48AD7E1jucARZFjVOxfm94Tj/Rg8ZpZapqXpVw3NbtIjM
oEGGfdqovL4rFikvr3NHceALIXbzOR0H7V9rtFh1yBjs2RMUd3iS02bXwSC3
oNGWR0pB9AkglNL7HN7fBrTXrtl9YWvEtf5et45YAg3SwFoHqFNeruZUMxsW
Zjz3WTTroEIv2L+nx1pRYddYSYlGunu3RbplS6DvsPehfGU3Vo5CRLMcJgkq
DoaPNUCLPXDYERuf3aQGdeIEmYEDjSQq/LXrWA07gvzcmtXZGLtmNRLM2NBh
k5Yejv67etTp6bYGrLFRoz9/XvxWyZT7a0nP9P79tkyZLcq3/MvTuPvHJAVs
G9DDS4DL7KkKQoseeRAoNWwi+BpFo0cyFfLNhs2hB/Mxcr+ONjMc/Xf06JED
hkGLRo+lBXsK3WpPoH5B4GDE58+TJKmWASfvKJTjSchKREay3YxYaotib/Sh
sIDL+vKjUyd3QU3B2YrWOJElQs5qNs0+Go7zIT10bJTZcMKlm9iEiB4HU3hO
HNpWpQaO4+fpvN6wj4P2AY0WbzaBg7R3r0OLDUZX6ciTrTP4wVZR5bVvqcur
OJbDNhwb0ONtQIu3sGHv1UY6OpQBo/v3N8gzaLSuuudwZrFGe5Wuiqu8mlaf
DEf7P/Vo6+ut63DmjEbL4i07Day+9fl3KK8PxqXpGh6HHrvmUN/+SWeNXrAH
SOGrIPovLtZjhcKGX+7HvG6sBw/Yk72lJbUaGZbwJO6I6NhU5YSP9S/1WOGX
A/WlSw5tUZGVY1tAEaD1T7S7876axXAI4ItnUMbpNqBFdQlUtt330JE1sgQQ
s6tXA7SRZStZGXdkbx2tM6+b192mXBy0f6FHizQeUFdXa7SlpXbEbDW0QLmz
wmwxHFpGJ0eDkEXkdqYNaIuLLYNtGsihZQshUUeUnvxikZBvRDSzUdELoSql
0rxw9P9DjxqxN+T52jWHnvWjeN+t5lit3NZoG8iewmaDXNAGtLCR8Dttskmf
CVuzxoY87IU7tG7TNSNDtiUcs9nrEV2VRxL7xkH753q0SKWhSrKpSaPl6Eos
w6pVAVpb9n1ESr/5d4cWfh4WEdwP/v28RdvGDGqMF+TGSo0EBCMmBbE27nYE
nz5Npk8fI28hoWgr6sr5iouDY74iiSyL4NUxOlboExwvoZrwBbSxCnf/THMZ
E4xEVPDaJ7mKfUkkUNj70EQHNYI784M3BxEFp1dR8+bXQzviE6ZWQ14FG0mt
3mtHYLfBQGljY1C2zJ6baFaUUrT2EdauCWoMtTDmU772E4sDOuPnVGO9jjVm
suxezVzUJ4DkhgbHXLaCvISN2wF1zD2TWyEFPfxTEwvGoip/NHdeuCVRTI2T
Q43xdtQ4xCZbzWIJDXiLgP7DD8UdS3Jkgpdwfpm34LEW3LN0VmzHITpUFpCZ
IGsa402lMcn8pnz33yyG4mkowWCqbCpsw4YoSd2dVmCmTzsIZ1VLKrLpOGMN
bw2akFtAbaJ0aYw3h8anFhYD8YVPKa8BmfrTkg+PKvlzr2lIP1ToTi3jwDoU
gb/dd8lSGUtp3qW5g2o9WCl/Y/n4cRs/VVdHZU5ivHPBzeU17u9wR0TBWvNy
G/A6P8DfIEEaHBsc/uaGG6dKx2h8aDtoh4QP7NReCcf3p3qccHeQhmbpgW2F
B2sD9ABfq5cmaHwIx3EGZTktr06EL8V6dX36eF5lpVdf74m9zcgI8LR6G0Jr
vch2Vrxlf489Dr479PiuXrXzWF4uctO3b2BEdIpLHbp1cgNjC7eC1yHuqE2E
LwXyCVyeV1Li5eZ6gqumJiqJ1PpNN2r+2C2XHXJEn7y6rip8gu5relRIaMCb
KCyUyB5euB/fulG5192Ulta2lpJUSjVDzBCk6evDR/U1RzZCncGDPR6SZGXx
Rg2Fh2Nn7QVrPO6QG0tlQyI87YN3dSFFxJIflM2rKDVwd/3wmjsKo+DkQx9x
i4PndvUkR41ey569si20S45PRR0zQGZEanPOVmh88IPm0Bzxd1nyRQ0mwtfe
82bM8JrXb5TUcno66VefwQrMnZPhrIDmH5KrSKUf4A7U4Xj+RI9ryRKvYeEy
ggNlT3sE49qz55zT4xoPpBx7dByLyu2J8PB4Vq/2rk2dIxqY9bqDz3Immvfg
wQsaPg6hgl/TuTdZdpk4b1O+VQ/mwAGv5r0pBI1o48AgjeaS4lhKvKQcMqiG
0X6vtSvWJCdCxiM6eZKqho5m1St5Ep02O3miVCSB9YPW6x/RR6KLLpK8BMl0
CEdysx4RL9eq3w+kl3p55Ftrna3C2SX/nWCt4xyYWGahbPAVWMffdFRIBefN
bjwcmFf97mXq3bMJqlajiH6BQb1GAeleQAtEM5y0obm5MXxcX9HjYsNc3Wso
9X6+3q/eDvbLkLIWxXq8VCPDcl1IC0X00m2WzXwlETIeGZvCmgFvUO/f1XqX
LnkOiduS8nPUDgm2Zz6hTySLxV6qgLo5MZLaUeOo968rXMpMVzHt3l3gtqAc
Eqjtj+ljmRkcekO7NRzJjW4k3Oomz6GXflGEUnKNRPkwDkkZ9zE0RtI28LvR
/iQxkusr1lLfx067OgO9kexWkHJcXCIOKgFLFkEIDDza7eHIOqoRNWVkUb+/
OehduBCMCOcDsYLOnRNvpZuiAVYdrxbBKyHmcuf5EpBfawMydl5H/mibd/JE
C/C6TI3LFPin/zQy6IYpNEWE4XPuvq4zd4Qj66AmBRHS6O+v8vbvbtQpMJUf
0MjwAgewEP6fv0MnIP80ITLmYRPN/Oly2ri0CiPTWRmXz8vMLHLI3KiwMYEU
jB80Csi7FDLB5f5/x7J/XkJLPrIxULKiOsbhu3vUx0jOrqf1sueHxIfTQ/eG
Dyo5Ihu07ZX1NOXlU3pQKHpGmsM/iaeRYSt5FHektHDCmWVFQP5ZImTJ3oH3
d9HoXx3RBEiqHdYvzdYtO0xIpHxIH4oPu5g7dLtj359rTJ00JjQR94JNOdT/
hwdZK3kuLPUP5bnKY11TgTersGUS/rFakpMCaH+hEXX+AqIrZ8pN/+/s8q7X
i6gLKBXluuJhNGhVvHUQjAMWHCCqoZpgSH8ZE5NdudV1puf9e70rF23s4WqR
lGvi1wh4GV6GaDxWRmIEeQnDPgmGBzSGmyIYWlpkueamXZQxpEUX+7rsMd45
OJtmY9Jl+YzijlQq2ndiTgdNemq7l/ZJjrgfECVWmoDYRXEFR7RQFDOOxglk
HDaHusmlXIFst7GS/vB/1uS/XxU00drhB2nBS3tbUj895rwT68kmQ3F7E72J
0KtJ7UAFX0PRAmhjGu2s291Sv7D4S/+7I39hCSVD/nqzaAtWvNZOdhDtgCoy
6D9kQ96it3iF8Byy4HJ45a3z1rn8sWD4O01M2/+bkL8Qm4+fOJPkmizVpFtp
JHfkEVE4yNMrAsRRj4gsrNtz3OG+4RWgu7lDe/ovY7EM+kkimhIcMW/fnH3q
LGhp/qz5M6yQdsjSs9rI9DKlApNVqsX0i7aNPvY/1olWTfLdf/WU289t23+O
aQO0L/GfVhJBi6VW9efye4z//ZZ0w38CfpJMDQ==\
\>"]]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.410374935788305*^9, 3.41037493843894*^9}, {
  3.410375356473483*^9, 3.410375357222987*^9}}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.410374745779085*^9, 3.410374811260517*^9}, {
  3.410374898308358*^9, 3.410374913945377*^9}, {3.410375031625885*^9, 
  3.410375037506136*^9}, {3.410375353179807*^9, 3.410375354661168*^9}}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.410374814721533*^9, 3.410374816112363*^9}, {
  3.410374947353578*^9, 3.410374996550804*^9}, {3.410375359681493*^9, 
  3.410375360344409*^9}}]
},
WindowSize->{1125, 1123},
WindowMargins->{{107, Automatic}, {Automatic, 0}},
PrintingPageRange->{Automatic, Automatic},
PrintingOptions->{"Magnification"->1.,
"PaperOrientation"->"Portrait",
"PaperSize"->{611.25, 789.5625},
"PostScriptOutputFile"->"/n/22/msgill/ellips.pdf"},
FrontEndVersion->"6.0 for Linux x86 (32-bit) (April 20, 2007)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[568, 21, 330, 7, 55, "Input"],
Cell[CellGroupData[{
Cell[923, 32, 100, 1, 32, "Input"],
Cell[1026, 35, 77, 1, 31, "Output"]
}, Open  ]],
Cell[1118, 39, 74, 1, 32, "Input"],
Cell[CellGroupData[{
Cell[1217, 44, 94, 1, 32, "Input"],
Cell[1314, 47, 121, 2, 31, "Output"]
}, Open  ]],
Cell[1450, 52, 1355, 33, 99, "Input"],
Cell[2808, 87, 2258, 55, 121, "Input"],
Cell[CellGroupData[{
Cell[5091, 146, 436, 7, 32, "Input"],
Cell[5530, 155, 23473, 397, 246, 10559, 183, "CachedBoxData", "BoxData", \
"Output"]
}, Open  ]],
Cell[29018, 555, 140, 2, 32, "Input"],
Cell[29161, 559, 236, 3, 32, "Input"],
Cell[29400, 564, 190, 3, 32, "Input"]
}
]
*)

(* End of internal cache information *)