(* Content-type: application/vnd.wolfram.mathematica *)

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

(* CreatedBy='Mathematica 8.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[     53547,       1133]
NotebookOptionsPosition[     51779,       1066]
NotebookOutlinePosition[     52142,       1082]
CellTagsIndexPosition[     52099,       1079]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["\<\
Least squares solutions of linear systems and orthogonal projections\
\>", "Section",
 CellChangeTimes->{{3.625568347222041*^9, 3.625568382004538*^9}}],

Cell[CellGroupData[{

Cell["\<\
Let A be an m x n matrix and V = Col(A) be the associated subspace of R^n.  \
For instance:\
\>", "Subsubsection",
 CellChangeTimes->{{3.625568392836306*^9, 3.625568444035939*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"m", "=", "3"}], "\[IndentingNewLine]", 
 RowBox[{"n", "=", "2"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", 
  RowBox[{"A", "=", 
   RowBox[{"Transpose", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"1", ",", "0", ",", "2"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"0", ",", "1", ",", 
        RowBox[{"-", "1"}]}], "}"}]}], "}"}], "]"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.624355479294818*^9, 3.624355525906296*^9}, {
  3.62556827720107*^9, 3.6255683166081333`*^9}, {3.625568456229638*^9, 
  3.625568466750351*^9}}],

Cell[BoxData["3"], "Output",
 CellChangeTimes->{3.6255760390238166`*^9, 3.6255779420094204`*^9, 
  3.625578092952118*^9}],

Cell[BoxData["2"], "Output",
 CellChangeTimes->{3.6255760390238166`*^9, 3.6255779420094204`*^9, 
  3.625578092967717*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"1", "0"},
     {"0", "1"},
     {"2", 
      RowBox[{"-", "1"}]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.6255760390238166`*^9, 3.6255779420094204`*^9, 
  3.625578092967717*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 "If b is a vector in R^n, then the system Ax = b has a solution if and only \
if b in Col(A).  Regardless, there is always a vector x in R^n such that \
minimizes ||Ax-b||.  This is called a ",
 StyleBox["least squares solution",
  FontSlant->"Italic"],
 " of Ax = b.  "
}], "Subsubsection",
 CellChangeTimes->{{3.625568510579257*^9, 3.6255687137853737`*^9}, {
   3.625568769544742*^9, 3.6255687811442337`*^9}, {3.625568818024311*^9, 
   3.6255689231521683`*^9}, 3.625568971816285*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"b", "=", 
  RowBox[{"{", 
   RowBox[{"1", ",", 
    RowBox[{"-", "1"}], ",", "1"}], "}"}]}]], "Input",
 CellChangeTimes->{{3.625568984478387*^9, 3.6255689972390823`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"1", ",", 
   RowBox[{"-", "1"}], ",", "1"}], "}"}]], "Output",
 CellChangeTimes->{3.6255780957132645`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 "It turns out that x is a least squares solution if and only if it satisfies \
the so-called `normal equation' ",
 Cell[BoxData[
  FormBox[
   SuperscriptBox["A", "T"], TraditionalForm]],
  FormatType->"TraditionalForm"],
 "A x = ",
 Cell[BoxData[
  FormBox[
   SuperscriptBox["A", "T"], TraditionalForm]],
  FormatType->"TraditionalForm"],
 " b associated to the given system."
}], "Subsubsection",
 CellChangeTimes->{{3.625568510579257*^9, 3.6255687137853737`*^9}, {
  3.625568769544742*^9, 3.6255687811442337`*^9}, {3.625568818024311*^9, 
  3.6255689231521683`*^9}, {3.625568971816285*^9, 3.625568977982808*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"MatrixForm", "[", 
  RowBox[{"ATA", " ", "=", 
   RowBox[{
    RowBox[{"Transpose", "[", "A", "]"}], ".", "A"}]}], 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{" ", 
  RowBox[{"MatrixForm", "[", 
   RowBox[{"ATb", "=", 
    RowBox[{
     RowBox[{"Transpose", "[", "A", "]"}], ".", "b"}]}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"x", " ", "=", " ", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"x1", ",", "x2"}], "}"}], "/.", 
   RowBox[{"First", "[", 
    RowBox[{"Solve", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"ATA", ".", 
        RowBox[{"{", 
         RowBox[{"x1", ",", "x2"}], "}"}]}], "==", "ATb"}], ",", 
      RowBox[{"{", 
       RowBox[{"x1", ",", "x2"}], "}"}]}], "]"}], 
    "]"}]}]}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.6243555391945963`*^9, 3.6243555583343363`*^9}, {
   3.625568323030569*^9, 3.6255683391349907`*^9}, {3.6255689386346397`*^9, 
   3.625568961778758*^9}, {3.6255691615526323`*^9, 3.625569237212492*^9}, {
   3.625569889922791*^9, 3.625569911608953*^9}, 3.625570570668743*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"5", 
      RowBox[{"-", "2"}]},
     {
      RowBox[{"-", "2"}], "2"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.625578103341518*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", 
   TagBox[GridBox[{
      {"3"},
      {
       RowBox[{"-", "2"}]}
     },
     GridBoxAlignment->{
      "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}},
        "RowsIndexed" -> {}},
     GridBoxSpacings->{"Columns" -> {
         Offset[0.27999999999999997`], {
          Offset[0.5599999999999999]}, 
         Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
         Offset[0.2], {
          Offset[0.4]}, 
         Offset[0.2]}, "RowsIndexed" -> {}}],
    Column], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.6255781033571177`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   FractionBox["1", "3"], ",", 
   RowBox[{"-", 
    FractionBox["2", "3"]}]}], "}"}]], "Output",
 CellChangeTimes->{3.625578103513115*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
The underlying reason this works is that the vector x that we just found has \
the property that Ax is the orthogonal projection of b onto Col(A).\
\>", "Subsubsection",
 CellChangeTimes->{{3.625569268083712*^9, 3.6255692697396803`*^9}, {
  3.625569864213821*^9, 3.625569951037139*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Show", "[", 
   RowBox[{
    RowBox[{"Graphics3D", "[", 
     RowBox[{"Line", "[", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"A", ".", "x"}], ",", "b"}], "}"}], "]"}], "]"}], ",", 
    RowBox[{"Graphics3D", "[", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"PointSize", "[", "Large", "]"}], ",", 
       RowBox[{"Point", "[", 
        RowBox[{"A", ".", "x"}], "]"}]}], "}"}], "]"}], ",", 
    RowBox[{"Graphics3D", "[", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"PointSize", "[", "Large", "]"}], ",", 
       RowBox[{"Point", "[", "b", "]"}]}], "}"}], "]"}], ",", 
    "\[IndentingNewLine]", 
    RowBox[{"ParametricPlot3D", "[", 
     RowBox[{
      RowBox[{"A", ".", 
       RowBox[{"{", 
        RowBox[{"s", ",", "t"}], "}"}]}], ",", 
      RowBox[{"{", 
       RowBox[{"s", ",", 
        RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"t", ",", 
        RowBox[{"-", "2"}], ",", "2"}], "}"}]}], "]"}], ",", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "2"}], ",", "2"}], "}"}]}], "}"}]}]}], "]"}], 
  "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.625569953890942*^9, 3.625570054650062*^9}, {
  3.62557008600945*^9, 3.625570125177227*^9}, {3.625570162920803*^9, 
  3.6255701920170717`*^9}, {3.625570499543186*^9, 3.625570504906995*^9}, {
  3.625570586053319*^9, 3.625570590444787*^9}, {3.6255706462128887`*^9, 
  3.625570655381824*^9}, {3.6255707079559927`*^9, 3.6255707161159887`*^9}, {
  3.625570750434341*^9, 3.6255708652754993`*^9}}],

Cell[BoxData[
 Graphics3DBox[{
   Line3DBox[
    NCache[{{Rational[1, 3], Rational[-2, 3], Rational[4, 3]}, {1, -1, 1}}, {{
      0.3333333333333333, -0.6666666666666666, 1.3333333333333333`}, {1, -1, 
      1}}]], 
   {PointSize[Large], 
    Point3DBox[
     NCache[{Rational[1, 3], Rational[-2, 3], Rational[4, 3]}, {
      0.3333333333333333, -0.6666666666666666, 1.3333333333333333`}]]}, 
   {PointSize[Large], Point3DBox[{1, -1, 1}]}, 
   GraphicsComplex3DBox[CompressedData["
1:eJx1nHuMXeV1xY0fMzYD+MWrgDGQURKprRopBhoV9btJmtQGTFD6R9Iq6QPs
gBSmSdMk5ZUURVXVSlVahodU9aGIJM0YO4lDo2mkjnovDqGR0kd0SwnmQoAb
8sTg8di8YuCW3rN+6+hbtv3Pldece+4+6+yzv73X3t+58JqPvHfn8mXLlp23
atmyFa9/Xvj49n8ejUbd/PzQx/f+7fzMUf9/99TG22YGL3Tv1Wd93DPd5+75
xIemty0ZPzD+f7+7U38HX2zw0hw/zPMI7xvX75Vl+hd2cnwBv/J9v3F46+zK
js5vfHXvqVNmp9d0ZI/xb775T984mD9Vxy8Z72497YLZ6fUdfh9806uP7Zmf
2dgRP8Zv+drutw3mz+hgF/j+nTc/OBqd1Ulem99b3uP/+QnP9fFL/gQ/IB7h
Gbxha+9l8Ay+KF7gGfygrh+ewx7zHPab59ruCfNcn2fKPNf+tdY81+ffYJ7B
m/t7unmuf/fMTvgN/tBLXi+K4/ITnvn/DvELz+D4OTyD137ej/v5jHmO3zXP
8ZzxvJS4LvMMjp/DMzh+zu+A9+Tn8Ax+vvyc3wcPPw+eJnrprwfl58nrifya
T3iu8WF8LnXh8bnqPvj4Uv/d5zfP9e8dNc/18cvNc1xvpz4f55nq1PdzyX5e
P09L9vP6Pi+ln5sP3V/zDL5dfg7P4BfFfUr/3Rl+XfM+jP/3/Rm4eY7zmOf4
XfMMfoGun++B1/F8yfik/JzzgePn8AyOn19Y/47jee3PxMup8Oeh/bz2Z/xx
eS/9dWfcp/Tf4/v1MPge2s9rf/bx4c8+f/iz7Sn1/bH94c++3vBn8xP+PJS/
rQ1/5vwbguc+9yt47tvPa577juc1f3385xh/Pf59OsZ/j/k8flzpB8/98PN2
XQ2/LGFnqXn2dQXP5iF47ncfkJ9HHuS8BZ7rvGltr84PWOenzHO9Lk2Y55qP
5T34q5/79jk5ftxNXo/16+A94rf9OOO6eQ57zHPYb57jes1z8GOeg8/05yL+
05+L7lf6c8HPg6eCn2ecvSDiVOYL6a95n9J/T+DX6c8ZlzNPSX/GzpJ5DX6e
eSV+nnn9N2s/j3i3IeNzUfzJ+Fzw88wP8POIz6WO522c3RlxKvOFE/hrN3k9
kV/Hupj5SHzanpLrKn4e8Znrzfhcws8d17Q+Zr5RtJ5mvsH9Ms/gWq/NM3gd
z51vhJ+3cTbjVOYL6a95n06Ql5zAv9u4Enmz8YvqOsS48rHMN9LPvU4r38v8
2X4e+XMhntf+jL1T4c9LpY7nmR/g5/04z9GM45G39tMf83nIvPAY/z2+Xy8d
E1eiDjSu+iLz5/Rzr8eqX7IeLKp3zDM4fs51g+PnmdfW8dz1YPh5mx/kOp31
WcbZWm9xPVgWj/+clANxn8B3nMCvow73//HzqAcdz0Pf0N/PTH3Dfh76hv0c
nmt9Za15rv1gyjzXvzthnuvzL+/V6/rRkus0+MFYJ8AXI07V+cLehfTXA3Gf
0n/Tr0NXiuPt59YFpC+ZZ3D8PHUE4jk8g5O3hF4X8dx6Xfi59Tr7eZ0ntXoY
PKeuAM/1dbfrBDh+XsfBkf08nivOcwyvu+O4/GzytVe7+cl18P+Gt5N6+r+/
d8fvPvGB6W3Wo4zv+8a7Pz8/c1Ynzlt0fC/OX3R8nv+YT37nOl1n2iu8G8db
N4vj/fyG/fKTkzpLn/nRDTOD17on+Cz8/5LhnvWz02fl8UV4L473J/bxvbm/
e/SVrbPLevn3+8Unx2GvcNvP93SePF51aHv9/H37Xz+5ZTDf6recR3jao+fl
1fKpC59eGI1aPffO/zrrya2zr0gPftG/28TVg8ZrnfNI9+J71313NDrs46fH
n08Y5/jmefypz1evK3sXwGsduN99dPy7P/bx5zXnL+Ac//vjz+8Ufp/jz2/O
a5zj58b+86PypjH+3z7+e811Kl7OLYBTL3A85zlvHK8O+3jwHfLzc5vj7Q9b
xrys6IBz/MPPv/2dg/mTOufr+jj+kbE9k8Y5/uZx/FzVaez9sY+/dXxfp4xz
/PkzZ7xlML+ms1l81Ne11jjH71Nd09h72Mdf8x+Dh0aj9cY5fsvHPv7ZmcG6
ziPir77ejcY5fuLqXX85M9jQaext+xQTf/Gu74xGZxjn+Ed++ft/NjM4vYPf
Epd57vDPWmd+2X5Y68mH7Ie1Hj60vzmP//z/4/fZr8EbfxvYH3z+MW/P2h9s
p+I59x38kJ5P7i9441crfB99/nHcmPT9qvX/Kd8X8DVvOedv5mfWmn/wxt/W
m2fb+T+f/PXB/EbzzDrc3K/l5hmcuAzP4LeO/fN58wx+vdZNeAbfO/79h80z
+JPjz33mGXzrOP48bZ7BD2jdhWfwi8/Z9lvT235unsH/TesFPINfs/LGM2en
V5ln86D6CJ7Bd99yxYrZ6VPNM3jjn+vMM/jNN619aevsBvNMXrOj+V3zXOvD
J/XgGfxXm+syz+CL+h14Bn9Pw5t5Bn+l4bkLz+Cfb+6LeQa/RnkLPIPf1Nx3
8wy+T+sOPIOvavzKPIPvUn0Ez+4j/Putn5wZnGyewZvn6DTzDH7F5rmvzc+s
D38eKs+ZDH8eqq+0olf787B7WOtt7c9D951rfx52H2ziQPjzsHtDEze6tT8P
uw+NP+8Lfx52n1HeWfvzUH57KPx5KP98udT+POx+SXlF7c/D7hWqj2p/HnYX
N3zlp1tnV4c/D5UnnBL+POye963Oh6e3rQue+93NzboTPPfVp1sVPPe732vW
wV7Nc9sHrnnu634dDp773S8363vw3O9+rMkTujXP/e5PlJ/XPPe7K5u8JXju
dy9v8p/gud+dbfKoUvPc7z6r+qjm+fXjx3n1RPDc714yfq5PDp5fz4f+975t
09tOM8/UMc3/p8wzePP/SfMM3vx/hXkGpx8Hz+CP6rzwDD6t78Mz+Juq4+YW
llX/vrEAzyDnil94Bid/g2dw8jR4Bt+sPA2ewcm74Bmc/Aqewcmjwp/LoOE/
/blc2tyv9GfVTRPpz+hCvfDncrfy8PDncpXy7fDncrLy6vDnckT+Ev5cfkd5
cvhz+ZLy4fDnQn4b/lzIb8OfC3ls+HMhXw1/LuSlEZ+L4knG56L4k/G5KF4F
z0PX+xGfqbN6EZ/LGuWNEZ/L/cobIz6XQ4p/EZ/LSuWNEZ8LeWPE50LeGPG5
kDdGfC7kjRGfC3ljxOdC3hj5RtH6mPmG6svTMt8oVzXrby/yDetakW8U8sbI
N1z/R75RyBsj3yjkjZFvlK8qb4x8ozytvDHyjbJNeWPkG4W8MfKNconyxsg3
rDNEvlGuVd4Y+XNRvpf5s/xwXfjzUtnT5JPhz0vWtSJ/LuSNkT8X8sbInwt5
Y+TPhbwx8mfFn6fDn5fKS8obI38uX1DeGPlzIW+M/Fnx4fnw56VC3hj5cyFv
jHqw7Grql6wHS5NXrO9FPVhObuqjXtSD1m+jHizkjVEPFvLGqAfLEeWN8Ay+
R3lj1IOFvDHqwULeGPVg+a7yxqgHy8/03Ec9WMgbox4s5I1RDxbyxtCLiurx
XuhCet439EL/KdIHeqHzuO8cek65QHlj6DaKD6t6ocOUR5Q3wjP49cobQ28p
m5U3hq5SvqK8MfSTcqPyxtBJFGf6JfSQQt4Yukchbwx9o5A3wjN63WSjh/RC
lysNvxt7ob+Vaxt9phc6m/vOoaeVTytvDN1M+tpkL/SxQt4Iz/X5R93Quwp5
Y+hahbwx9KvyZtkROpXyst5C6FEFPTB0p0LeCJ/5CZ93Kd/aJH+Azy+P4++r
0tl+6LqPdWtRcQ8+Hxvnha9Jd9pvvZG6cZ/Wd/hs4s5Iek6rN6KD4xfw+WnZ
3ehvrU6I/dQR8In9+xaHT2+dPdl1H/ZrfXc9MpD9t5/8iZ9vnV3v/A37d4Uu
1OjXL6rua/W3/ITPQ+P87CX97nPmk3xrh+IbfDbP/8vS2f7VfLJuXal1HD7f
M85Lfi49Ycl8xly4+dwn+xt9bOR8GH+lXoBP7N9x2envn9622nxiP3UxfGL/
NXtP/fr8TKuLYj/6D3Hgatk/+471lw3mW51T16PjW/0Nf+X+w3N+wvN21Q/3
Kr7Bc5M3PC/97UHXHeRhc1rH4bmpZ15QnXvAPNf9yFZ/m5D9jW7T6m/U7dQL
8IzdjQ7c6m/YT10Mz9j/oT1XDrbOnmaesR/9B57vlP3X/mTX81tnW3++o7ke
6cwrzfM3tR5TL2DfYfnxovwCnvMTnpvvL3a/PT7f4+YZPfpJrePw3OQrh7rn
jHn6QfQXnuleG/rbXbL/9ub6zPMDsp96AZ6xv9F7V5hn7KYuhmfsX/3Sp66a
3naKecZ+9B94fkT2b/nsO746P9PqnE198gP13Vr9Tc+zdWZ4vlXxd4/iGfbt
lh/Th4Ln/ITnNaqTb9A6Ds/N8/izblMXPGSe6Qc8G/rb+bK/8dtWf9su+6kX
4Bn7Lx/74UnmGfupi+F5P/nmO796/2i0xjxjP/oPPL9N9q/aMrh+elurczZ5
xkPdO+++5R/nZ9ZEPvywdWZ4bnSE/erjLHP+cL3iMv0UeH1QfoxfwHN+wnNT
Pz/VfWTM5z3m+ZB0imdCf7tK9m9prs8875H91AvwjP0NP685H8Z+6mJ4/pbs
b/qkk+bZ8U/6Dzzvk/0P7//ji2anW52zeW7v6S5dfvkvDuanor67zzozPP9h
s551Z9/4znNnp1v97YjyCdZneCYuE+ewDz/GL+A5P2v9rf2E50anuEd68lPm
+SXZT70Azx+V/Y1u87J5xn7qYnj+iuw/+Oj6ufmZVeZ5hnxa+g88/xLr+ejs
b4xGp5rnadkrPPSK7xR0Znj+YvO7Rb9rnn8iXYC6GJ7frzyDdRueicvEOXh9
SH6MX8Mrfoxfw3N+wjN96wnVC/D8Bdnf9KGeN8/YT10Mzx+U/bc3fmWevy/7
75f+Y3+W/YcavzXPPfmB/Nw8XyB70ZnhudGtHte8waR5vkY61ybVxfC8Uvkx
+SY8/5PyDNZteN6quEychmfmlg6G/oYfr5ZfO9+IT3jG/qbOXTTP2E9dDM+r
ZP9bmzhpnrEf/Qee3y3772jisHm+VHFNcds8T8oP0JnhGXt3N+uCeWaOgX4K
PDd9uB9K3zvJ9d3lypsvVR4Kz+cqzyDvgGfm8K4L/e1i2f8exWl4xY/xa+zO
T3jG/itUF8Oz6q7S9FVfdX2H/coTzDN5kvIK86zrKcpDzHNXcQ2dGZ7xA+U5
5hl76afA82KTR0l3WmGe71e9x5wOPJM3k0fD87PKm+kPw7PzDOUd8ExcXlSc
hlfsx6+xOz/h+ZDsVx5rnrGfeSR4vl32K09u82flncqrzfMtWqfRmeH5WsU1
5e3mGT9Qnm+esbfRQ1qdk/4EfUN4pg6kLoTnXTG/Cs+zsv9u5dHwTJ5B3gHP
2P85xWl4xX78GrvzE56xX/NX1oVWyn7VfeZZ11NUJ5pnPc/WmeF5s9Zp1aHm
eVJxjX6K/Vl+sLOpc80z9qL/OK+TfjElPaPW7dv5Hed1sn9SdSE8kzejMztO
yP7blHfAM/bTT4FX7Oe5xO78hGf0F+kY5ln6TLn/72/8zMyg1d+Yp0Nnhmfy
Tukn5pl5O/op8Exckz5jnvED+obwLLuk77U6J3Npqb+ha6BzwDP2ozPDM/af
ojwanrGfORx4Vj1WblSchlfsx6+DT9fL8I/fiGfHGX12N/3C7v8cjVo9BB7O
b/BeHH/Mp85fnpVfxLxZoY8cuusx847x/Pk4cOlynmeO9cdzS+DKUzxn435+
k7d4rgJc8d9zAOATld7FHEDLE31WcOIJfUFw1k36WODkh+g64NRB6Grg1Pvo
l+DoWvAJv4djzhVc+pv5jPXEfIIr7zCf4MpDPGcDrnhuPsEfqOb3mJ+wnmk+
wVkf4RNc9bX5BKfegU9w6nr4BEe/gk9w9Nh6zv71/FHzk/BMfLpSeSA8h85m
nmPdMM+RH5nnqAPMM7jqRfMcuo55Dt3SPIMTb+AZXHW0eQanrodncPQreAZH
p4VnX9eOxzqD+Xb/O7jij3km3tMHzPlp6WnmOdYH8xx5UOyTd75vnqOu9ZxT
6DfmOfRJ8wwuncc8g6teNs/gqjvMMzg6bb0/56D7DvAMrv5sp+4nHZQ/+H0m
1kUWVdfDs/u0qo/g2fmjdC14Bj+o9QieI683z1G/mufQacxz6JDmOXR38wyu
utjzBOCqL8yzz6N+BDy7f/LFr/3RzOAU8+x5gT//0rsG8+s6dX/ux90rf+2T
/dHIcdt1OfpVzvdLHzPP4Oi38Ax+h/oU8Bz5u3mOOtU8hx4T82fWG81z6Ovm
GZx8GZ7BVUeYZ3D6a/AMrnkY8wyuuYJO3e98oqv+rHkmv0Onrfe/o/dPmX/3
u6SD1fso6Me2+4s5nr5GvT+VeYR2Hz3H09cD5/jnNLfGfeR45dG+v7V+0M7d
Rh1sPOpg3/fQgYyHDhRzh9Y//TyG/nkZfhJ9APsPx9fvIxpmH8w4x/OeIvwq
+r3xHoB2TgF/M5+33fsH09va94oE//ZDjr/1p795y8xgdbw3wPfX/unrfe8P
75yfOTXeJ2B/s99y/M0ffPv3RqN18Z4B+1sn6zyd3/4JflNznl7W37Lf/gYu
Huxv4KtVN+A/4Op/2U9Cr7KfhB7recroP/j+Rp/N9zH6yb6P4NSV2X9GP8k+
v+KAeQZHD/d+HsU75pfq/U5PFO3TMc/gzO/BM7jim3kGV3/W+yujXm/nsGtd
yjyH7mqeo89gnqOfZp6jb2yewak7PY8lHJ0QnsG1rnXqeRHr2+bZc23Numme
wa9o1k3zDK711zyDa702z+Dqw5rnqMvb+Xjh6guY59BXzXP0EzzPGn0z8xz9
YfMMjr7iuTrhyq/Ms49X36eut61jm2f0Bub04LnWK9v97+DoKvAMjn4Iz+Dq
t5pncPbJwjO49H/zHDqqeY6+gXmO/ph5jj6weQZXX9s8x/xOO8cmXHWHeQ69
2jyj36iuMc/g6CrwDI5+CM/g6OTwDH6l+kHwDK7+lHkGlx5mnkMvNc/RHzDP
4JoX8jwx+APSw+E55hrMc8zvmGdfl/r48By6tHlGD0MnrPdzvli0z9E8g6OT
wzO49AHzDK75efMMrj6UeQZnjgWeQxc1z+CaWzPPPo/0cHgGv159H3gG1/5f
8xxzOuYZ/EjsWw792Tyj56GHwyc4fR/4BKe/CZ/gi+rjwyf4xZpXgU9w5rLg
M/VL+ARnngs+wZmL8Ry8cOYL4NO6tfqc8BlzN+Yz5svMJ7j05A7zfDl3xie4
+myel3LfWHZyHeCamyjMo4Bj5xuU34BjD9dn/fPim57bOjvRoW8K/ug/nLl6
dnqqQ58P/O6Jr184mF/bIU8H3/Uvc5tmpzd06KOAT97wJ/fOz5zeIU8H33HB
uy8dzJ/ZYd4OnY3ngk9wzb17ngmc54L5G3Cei4uUf9Rzkys69IPB2ZdH/9J6
pvaf0m+zjqp92dS74I1utpH3DBhv5iXO6DCnid5FvGWeKfQ9zaHsN74v+urG
NY9K/wz8079922u3T6/t0O8B3/TUaX81M1jfoT9he+56/OzZ6dM7e6SHoRcx
J8/cT8590J8DZ26TfpJ1Nu3bpf8Bzr515p9yTof8yPPsmkuk3wPOnC15q/UQ
zdExB0yddZ10gN1Rv+9W/c78ZV2H8v7eth7nvUnMt0Xd2mVuJepfz2dF3e33
x4DTx2ZuIH638Pf6+OWdXVG3sv+X79Xn8bxKCT2kQz+M43m/AX2wqLN6zKdG
XdZjnjLquB59eHDiJvNqOY/JPErOTfCZ/Wf6/dnH4+9RF3fop0Yd3eF74Jeq
3j8/9j2g7z0Q+07Q8bbHvh/eE0I/POd5mf/IOYXcr7yFfnXsF2ce+JLYr4/+
thTvS0B/Y56V+oK5L+JV1Ck95gWjrvH8R9RZ7vdnv5R4BY5Oy9+j7usw90l+
znwU8Srn9IlXUY+4/x31i+MVOPq88l/n28wFEa+yr8snOH0W4lXkk6ybXt/p
N/Gceh1X34Q8BHzxhSe+PRqd2aEfTH6i9df5s/O05vge/V3w7crPH1M+knMT
d8V7bsiHD0hXBicv3aF1O/v/7K+IPlSPfQK5j+QGrZPg1B3MKbqvJH9i/QTX
/fJ7kFiP2A/FPoTo1/SYm8/9Fi+oXwTO/l/WyZzHYR4OnPlD5lzAmbNlP4zX
QekezFe6XyPdg/nW6IN0d6nf6H2lOo65UvAJ6RjvU78x57OYO4y+T2FODpy5
Pr7nvol0DJ6H6Dd16GNHX6DHfg/Pd0hnY87R8zXS2ZgzjT6F53xz38CH1d8G
f1p6GnOrOffHnGLOT/H36Nd4DgycfZd8z/ZITyMOhM7cY27C8yzSe5mD8zyR
9F7mJsC1nnleP/ejMIebc/3M6eY8KZ85l8ccZPQ1Cn/3vKR0XebMvG9Dui5z
UbZTui77uFiXr1YfhznsnH/nfQCsd5s0r5j76tjHd0RxDJw65PcUx8Bf0XHE
J9Yp5peIY+DXaU7vY4pj4Ee07jKvnPP4zCHmXDNxzOua9o0Snzwfqv3RxDHw
uzS3tqQ4FrqZ8/3cb3Gt4hg4+/iIYzn/Sxyy7tTUHY5j4KpTHMfAVdf02CcB
ztwXcSz3jX1AcSz33xDnQtfyvC/41fo+8505V0u8Yh1vdNozHK/AVSc6XoGz
n5c4BM6+dd5LAX5I72f4iOJV6Dxd5tRz/xbxCpzzEq/AP6rzEq/A2YdOvPK+
JR3HvjvyB+kAPfaJgUs36LGvCVw6Q499OODSJXrsGwHf3OgYvV9RXHCeo/kx
5vJzXyBz5OD7Y38Y+AP6PeZ0wT+n8zJXCn62vk+cy/lTdCxwzeN5HpL8Kvcd
Wp/ROkl+AP5mXT9/Bz9b8Y91xfNmih8xD9l9q/yauBnzeJ26r3PU+wqJR35/
gNZxnn/wq1RX8vyAq4/p+X9w9YM68ABPuS8QnOuEh/r9PXsX+Ds4+53hYX+1
Dhx0fI85Ol8vOPv76j7Bi973VO9He9HvCeR6wXkfoP2y+t2NnXqO5/Ax+/PA
8Yd8fyT+wPWC4w/Y7/c66D5iPzj3sdaRD/s+Wucl39Z9TP/N/W31+5P2LmAn
+CbqXdnjvF08Yw84PGOP99OI5/QXeON3wc9Bz6B/IBweHKeEw8Ob4j6go3A8
OHbW76GaWzj+e6jmFo7/Hqq5heO/h6p972b9HirOk++hmvP7Q3PfG8fV75ua
W8j9QLVelO+bmluAx/p9U+38df2+qbmF479vivO0+7CwN+9LzftUh+vmOtgH
BJ/g7I+AD3DmxuEz42XyxHmd7wrn+/AGnvt7wOEF3qxT0Z+J+E3chjfvI1fc
5nt+b4qed9fPVfxr96OB4295HczPw0O+DwAe6uPb91eA5z4bcOdN9FmF4z/u
51Xxon2PU32eyU7+Pved6/I+Fd13z7cJpw5x3SLc87/xHOT+Fe9vpu8k+62b
6T7yd3DuI3bW8bvdP1XvR2/ryVyv0r9ynwc4PDsPEw7P/C7ngzd+1/k8/au4
n7nvwH1z8eA5vCret3oWeOrXPJ8cl+8vyfm/en1vn596HW91N88/63rzvR0Z
38HzufVcd/hrnN/rBf7LnH/mSTmnnnlTxv2ML/Xcdftc1etsO48S9nQyf2VO
3vP6wuEt85u8zsx3sD+f14yzGS/A87kKe5zXhf0n1K9yHWBePddb7hfXke/b
zjwp+Tp+3tSu07mu1PlGG09zXxo81Me3czBhv3mI6zUPPK+ss1wHOPE210ni
RuYj6U+Zz+b7EjLvyXWYz4zf8JNxi797XjTiIDj+w/fA8R/49nvOwn9iXrGH
/Wmv+wGyl7nlup5p95tn3pb5SuZx3n+l62bulPODp535PGV+88bIA8DPi3zK
OpLmML3/RzjXm+ttXm++lwV7wN8Qz02dn37O15P5iPdrMdekOcZcVzX328u4
D2913DzsvLTOgw+XSyJPzThU5z2HC3VN5h9nx/WDUw9+6rjrc7sfGZz7wnXV
7w9r3wdY1/ft+wDrOnjS11vXW22+nnVe5oWZn9V1w96FzG/OC77q/KHVL7we
h24Dzn3P9Zv7XterR33fcx3ivtf7Q446L824k/tYwVMfz7oc3rJ+zXwU/4E3
613hf3X+1fKeeVzmX+xL/T9UDRY/
    "], {{
      {EdgeForm[None], GraphicsGroup3DBox[{Polygon3DBox[CompressedData["
1:eJwtmXngFtMXxufOvS1Ei6KoJCFLZGlTya4ooRBZirSphPBTWSJLURIpUZGo
kOxLi30ru7K0a0OrtColv8/jmT/O937mvPedd7537pzznDM1r+7d5ro8y7In
+RMZX2J8CPgePiTg4/gouBacw0fCB8MBPgKuCWfw4Qxb+d67Kct24zsUS/jr
4P8L//v4S3BcHf8ufY7/RI7357Mb4G2M7zHnXz4ri23E9se/Ef80/Bv0m9hi
rAT+Uozrsf3gD5hTC7sFLoHvF6wUHBmXYCXhMozLsT3gs/ndY5l/F7yD8UPO
XwrfHny+DCuNvxrjTuwQuCrj31gtuBLjFqwavA/jZqyq1o+xH+d4Bj6B8+3C
7uT4eMad2B36jDl3Mb7AnHr4/sEGcPwE/r6M4/HXwbcVu57juozbsf/pehl3
YLfCxzH+jd0O7+Q7C8J/tyAby3gb8Cx8AJ+vwJpwvB/+rVh1/Hsx/olVgcsz
bsIOgDezDjOY/5d88HR4mz7DdmAHM2d/fMuxxpzzUnyNGK+SH96O1YTHMPbH
PwE+irlbsN4cv885a2I34/8XW6T9wbgff5Zg9ZlTmfEXrAG8D+MCrC68L+Ni
rB5cgXE+dgz8Huergd3Eed5lrI71gcvx+c9YHeaUZ5yHHQ1XYlyEnaD9ioXo
66+IbyF2PL4TuK59GBvjr8K4EbsA/oe5nzKnnO4pcyownoj/SO0p+bQn+fxb
7CCO2+A/TuuE/wL4WPgS+HStJXwOfCnjKZx3OLwn35uDHap9yPg9dghcmvE7
rBZcke/+hu3N/IxxDVYR3I2t1nUzJsZ12L56LBnXYpXgCoy/YnvBezOuwPaE
P+T3D8X+p3XDtxIrA+/CVun/1P+OLQyOD1dxHWcDT8AbsF+xlXqmgtdUa1WG
652LHcbxadoX+o72Hr4fsNocl2X8CTtKa8GcI7RmzNnOuT/AX1LxRc+i9rvu
Bb6lWEOONzDnbfh3rTFWHl8j3Qt8mzEFtEfw35D7Ov8HXwzfCd8BXwHfCx/K
3D+wDtpL+B9i/Bz/LfBF8B3wnfCV8H3wzfCF8O1wf7g9fLfiJOdYj12ia2bc
hHWDH2BO99z39z64s2IrfBPcFr4N7gO3gfsrrvK9tdgFHN+DvxPjEPwD4avh
B+G7g+/BA/DhzN2IddX/g/8axqH4D1O8xDpy3C94nynW3Q5frnPDt8GXwQPh
j1nP2titcDPFG8VcuAncB34LbgrfrHWHT4Jvgd9RDNd6w6/AzeEKBOdxnPdR
fS/3PnwEfjP3PjybOfdxPBtuDN+kz+BOXM9pzOkBn4V/IP7P4JPhN7BH+OxD
5gxn/AL/+/AweFbmOLRf8Syfytz+ykPwafBtyit67hQvlYfgM+G79TzDp8D9
4JnwGYrT2n9aE87/KOf9Ep4A3w2/CF8Dnw73hGswfxV2KsevFdc2Hf90eLSu
F54BPwF/BL8e/H/MKNa5FGv1sPIB/sGMrxZrW7KE/7dGcIDvz73OOTwo93Xu
DT+e+//aCx6V+/8tA4+Ej1Es133SveD8NzI+yflHwH3gMfBj8E3wWPghuBf8
GPwsPBCeUqznnpzzsdxrtQc8Am6o4A3fBzdQ3sbuhefz3Vdzx40X4Qe0NsW9
jswfwnEX/GcyXoe/K3wW3BvuDbfOncc7w2fAveDqfPc3rBnH1Rh/xU6CD1QM
wE6GazKuwc6HqzKuxJrC13Oe83LHt6vgU+BucDe4uT6He8GtcucR7fnEdQ7N
vedLlHBcmMKcBxlfZ85QuCc8Ah4C99Aaww/C1yr2wD3hlrnzUQ/4HPhG+Npg
/SHNcx18bu6cuIt99wm/V5bjkozfYDXgoxSnlZf0XEjvwOfBFwbnhyvgi4Lz
xpXwxcFxt4P2s2I/fC5cP1h/NYOPCdZoDXSP+N3eyTqvQXDcPRl/u+A43RG+
pIjter7aBueiy+EWwTlSMb85XBk+X8+X8g/cGn5NuYXz76/PgnOj8uDr+Psm
54bLgvNAJ/y/4H8heb++AfdLzpdT4OuTNeUL8HXJevRV6T64Ilw3WD82zPw8
XR/9DC5lfDF53+fSUNgB8GL8k5Pv5SJ4UvK9XAhPTL6Xw4t4qzwyD/8zybF4
PjwhOUYvgZ9Pfn7GF3lHGu8k5Wr4TMWxYO3QHG4YrDdPgZ/nu72S9fSxIftP
yCmXNSpy36nwS9LIyTrjZemdZP1xSpH7WjBnKv4bk/fKK9JZyTqmWfCeOos5
E4PjwVT4XLg6fBF8XvBeawfXDtaAx8LPcJ7OyRpxItw9WQseHqz3lJcPC9aJ
dRU3mNM1WR8/B3dL1p1PwFcka5cJcJdk3TweviZZ6z8Nd0rWpjWC9X5tzjkW
f4dkrdMhWA92xf+U4n2y1hkHd0zWOmPgK5P10Gj48mTdUyVYz9Tku5WDtcNB
8LhgvTAJfpz5lyVroMlwj+QaoLX2L3Fgop497R14UpHjLkquSUbAFyfXG08F
a97JRe67MFnvDofbJmupSfC1yXVFa+xA5l/M/MuDdes18NPBdcDz8MPMb5Os
9Vth1fBfiH8k/kuTa5vX8K3iuDz+lvj24zon6PmHq8DPwa3gyvCz8Cjmtk/W
gi20X/A/lTtHV4SfhocxZz02h3Oej78q/sn4z4ErweOVM/j8kuTaaQH8bLIW
GRWsTZ7iu6ODtdB4+EvmjEiOubPgh5NzQN8iNkqrzMb/SHLc/BoemZwbqvGb
q/U8w18xPpYcQ7+BRyXH+gHBeWAw5/kC/6PJsfiuYG0l7fQt/seTtdegYO32
MP6ynP8C/FM4nsOcJ5O14OBCR0nLzcU/Nln/fS8tnKwLv4NHJ2u7H+BxyRrr
R/ipZO2l+lrruSrzepwfvSbSPOX57bG59Uk5eAz8L59/nhxPE+NXyc9tCcav
k/fNbmml5Fg/A64cndfqBcfFJnAduAxcD97J5x8nx6bp8L7R+XcmY9XoPPU5
47DkPD0Nrhidr99hrBCdl99iLBOt2aSntG7SVG/D5aI1UsA3Ozk3nFjkC8W0
xsE54bTMtVHp6Lwcmftl8v5uEhynT8evYnFWcv77WWuWrMt1/Xcl14z63+9O
rhHeDNYj0oEzg3WHrk2x+dLo+LxSujVZ262AX06un5fDU5Pr2JeC62xp4MnB
9be07jLmvJRc878brBekV18Orvuljd8K1kfSn6r7u0frw0nBNf3Leq6Ddfpo
xZ/guuFp5YXgnoA04dvBGlDa9Z1gDSiNKr3ZN1oP/8r4WrLemhac/5T7fsP/
erJ2UYx6O3fPZGpwT0J6/pXgPsS04p7enlwLaz/ck1zvq76+N7mf8CZ8W3J9
2zH4/unefYb/oWTdpn1yZ3LNrv0zILnG1964I7kG/wQekqz7Xw1eI63Pp/iH
JuvsG4L1db/M69Exek2GBddtj8NXBvcDusA/8fnTyTWeeg73F1pRPaJBhYb8
CH4gWQuq/h1caEhp9gHRvS/1Fu4rNOQf8Fvw0uAa557oOvGK4L5C5yKPz8/d
E5CmWlg8d9IqC3L3E37XXi1iqeqa0iVcY6gHUj5av6kHsk+0ftM9vSG6dlDf
o0R0jpAeWF48O+qNlI3Wgepj7KUcVtz3G/WcZ+6BpOhcr75KyehcL91yIfZc
Zn21iHM2zVzPNtM9y6y3D4z+H6XN2kfXQdJUl0XXCNLb1aJzlp6Lq6I1f338
u7Veuev3StF6sntw708a+Gj827CbOa7FOCa5VpfmPyj6npYOXvvKme9Fu+Qe
n/LpG7l7IKoFDo7WBsoFDXTvM/d89ozWq+oLxei8+SfjOxz/EVzHHaF8nFk7
rcndn6ktPYJ1ya0Hfijip/ThBdG9Pl3zOqxdbh3ySe6eofRkm+g+mHLBadE1
gnLB6dG1w/OKJbl7d9KTLaPrMuXQVtF1mTRq2+ie3hbGmfzWP8F6b0URG5Ur
j4vuRSgfNY3uAyh/nRTdc1AOrR/dx1CeahLdW1AebBzdW1A+rRfdG1GOPiG6
V6Acenx0r0B54ejo3ohydd3o+3iQNBfWooifV0fXSsrvx0T3XvRcHxXd91Cs
qBPdS/lbe1x7N3duPTW6hlJ8vja6XlZ87qT7rTUM1jKL9X8F58yfdP3BOfNH
3V89b3oO9LvB9evXig/BOm4J/ElwLfuV9kxw3+Jn3a9gfbRIaxWsg+Znjqm9
i/gj/XN/EXsVh28tYvvcYM21EP4rc4+sXBH/u0XnnbWFf6XWIbhe/0b3t7ju
dZl1/arc2l51wLIib27FNmLrsW16HrL/WoHZ5uKcOvfPwfr0l+LzNdjy4rx/
FmsyK7h/8B2sDTgvd93xeXAv4dsiXr2Z3L+YHRwfFRvnBWvApdo/wZpxHry9
OP/c4vp0nasK3/bifuwo/MuK/dwz+X2A9Jz61OpRHx28x+sXc/Wd1Zl//87o
vPZH8f+vKf5vrbV6kepRV9H+K353S3E9G4o1Kld8d1OxbhuLNdX1qMZVk1P1
7JfB70R+yNwP/xaukblP/k3u2kEx5Ork9xBfBGs7rc+BwXHlsOL8q4t1mKqc
X6zbvnz+Alw2c13yYe7aZFNxfSv0XGTup6iXspbfWqf4Aj+pGIUdmLnPf0B0
7Sy/tLH2leqz73PXaOqXfpC7LtD7lM9yx1LVUJ/mfhej2vac6B6Raq6Pc7/v
0LuMLtG1jN5lXBNdK+n9Rdfo2kRx/tXknpHeTXTGPzGzLl2P76TMuf7I6D6n
3gtcFJ271bO9GB6VuW95VrRuV2+2OTwocw98F+dplbnv1DO6hn0jOLeod6fe
+1/F86L1fD/3uj0XnFteytxz2JG7Jy+dcEZ0v0t9hr/xt8+syc+L1vl6n/Vj
7vuo2nBm7nc90pCXRPcE9F7mymjNpnc37aK1nPLalOT+mt7jXBFdB6kvdzk8
LrOe6RCt//WObAnHx2R+F/ZT7vpX7wJ25u7h6B3BP7n7KqpJH4yOG6rpBkfH
MdV6D0THOtUO63L3fFQvbM3dB1Ctd29073QO6/MM/gVwS7gq3DZzXTk99zs4
1XfTcr8zUh36eu53DdLqm+EzMtebg6L7peqXtIius/QuYHfunpXqx6GxiDnK
vbl7TcqzzyX3x9X72pa7T/JBcH9UveKRwXlY66b3CP/m7jWpJtqSu8eiWmND
7r6K+vCto2s61bkzcr/nUk+1V7T2Vj13cnS9qR5mj+h9or76ufCwzPX+kOh4
qHcTi3P3OtTTeC93vFpRPLOKseqBfJX7vaf60suStY76IV8zHp25r/JF7nep
eg+SR6+5eiNf5n7fujf+H7Ejc+eC34v4rP7P6tzvU9T7Wgu3zNzreDf3e0nF
859zv5PVbw9MvhbFhO+KOKPz909+7yI9vwfWHf+twTl8AHx/sK7ROii/D0/u
meqd6S+5e0Ttg69R8VY9ojm534Wpflb8VK5Rb/PX3P1GxUDFScVD9eK+SN5z
6t82Yl7fzLpIcU/rqT7J3Nz9T9V9Z2IjM7+P+y13r09z1ffQfNW2K3PnDumK
Psl1r2qUUtG9x6bBWkP7Vu9cGkbXF6pN/yzyrPpRs4o1VJ9qdu733epnbsr9
rkp5eWkRY9V3+ij3e9L/A8oL5aY=
          "]], Polygon3DBox[CompressedData["
1:eJxFmnf8l1Mfxu/7nPOrhNDQRhGpaBoltNCwykxWQ5FRohAaCkVlr1Syd5rI
Lqkke5eRmZG9H8lzvZ/rvF7PH9++9/t77t933Pc513V9PqdGA4b1OSsURdFf
/yQ9t4tFcaMOxpRFcam4mXiq+BzxheIfdV4dPbqLDxfvo/HpGp8oniLuKL5V
fKl4qriD+BbxBPGV4j3EN4nHii8T/6T3Ok/cQ9xb3F7jN4svEV8h/kvjo8V9
xH3Fhcab6rVjxScLP9XxCI13EHcRrxEPE7cVdxC/Lz5L3EbcXvym+HRxc3Fb
8VrxcHE78T7ir8UjxV3EPcTvic8UtxbvLe6nz39IfJf4IfG7XDw9Won30uEE
jbUXnyY+VY/1Ot5aj046Pkjjn+j4bM4RdxZ/JT5X3FncXXyw3v9u8Q3imeJv
ND5K3FXcU9xK49eKR4vHi5uLp4nPFV8k/lPnXyDuLT5W3ELjV4lHii8WtxFf
L75IPEG8m/hq8SjxGHFT8RTxCPFocUtxD95TPE5cRXyRxvvz+8TfiqtV6Fro
nC/FV2hsp+ixTnptmo4f1vMcPdbptcriVjo+WX8/RNxafJ3+5kLxJeK24hvE
F4snitfp3Cp67C3uJH6ce6/xquLa4gt0/kLxUvEq8R8afyX5enM9fxe/LN4/
X//fxKvF+4kPFP8rflN8UJ7PzK+3kud373w9X833h/mwSfyG+EDxYeKTdP4F
eu1+8dy8Pl5M/r7c304anym+QnyN+B+Nvy4+QHyo+G/xa+Ju4kPEfXX+A+I7
xA+If9D4SvFe+feXef6zXvqIG4g/0/jx4kHi7cSfi08QnyLeVrwueb30F9cT
fyruJx4oriP+RHyceAA/X/yO+GDxkUxt8dviXuIj8npl/bfP6+17Ha/Q+J7i
/cW/iFeJO4q7iTeIl4v3EO8n/lX8knhf8QHixnr/L8Un5flwnPhB8Z3iB9El
8Rzxg+L54rHoCddY/Jr4MPFp4lvEs8WHi+/R+dPFt4vHiReJV4lfF58vvkbn
LxG/KD5bvEDjT4qXioeL54ufEC8Rz+Xe6vEpa1GP3hq/V+O3avwO8SN5bn+F
dugxTOOTxYs1/pz4DuYLn833Yd6zHvT3D4jnafwI8f3imeK7xL3Ed4mvF88Q
nyoer7+bJ14sPk08V+PzxU+IB4gfET8kXiA+Vny2zr9dfL+4e7S2XyO+RdxD
fIfOv1Y8XdxNfLt4mvgGcU/xneLrxLeK+4jvE88Q3ynuKp4tniq+XtxZPEt8
pfhacRfxbeIp4uvEO4i/EJ8oHiyuKf5YfEzW7+rij8RHi08Sr9J3HSzeUdxM
/JJ4CJoibi5eLT5V3ETcQvyU+HhxdXF98ZPifuKtxfXQV73/+zr+j4631mOe
xvtofKOOt9FjofhIjE/n1NDTi+JTxI3Fu4qfEB8nriauK14gPkK8ie+ux2Jx
X/EWGq8jflR8tDiJa4lfxwvEu4pbi1eKB4kbiZuKd2c+is8TjxW/ofGh4mbi
NuLa0Z6ClwwX7yueIb5MPE38PHNKXF+8o3iZuL+4gXgn8VLxSeJ64sbiFeKB
4h3Eu4hfE//Jd9GjlXiJ+ESN1xU3Ej8tPkFcU9xA/Cxznd+qxw7iReKjNB7E
NcX76/v11WuTxVeLP9TxsxrfWdxS/Ir4CXHt/Pd8/pPJ79coXy++T938favq
/dZq/Ajx8eJq4g/ER4pPyH6En+B3fbO/rNH44dn/thZ/KD5KfGL+/Yuzf3C/
KjE/xIeJjxEn8bviQ8RHiSvE74kPFR+drw+/f2O+/1yPx5LnF/NpufhxcZXs
T+/lfNAoz6ePdfycxnfJ15v88Ey+H8yHt8VPibcT7yx+R/y0ePt8v4bq+8wT
LxA/WdgflyTfP+YXeeOF5DxCXvlcvFS8m3iPwvkDP901z68vdPy8xncX71k4
j5BXyDvkpS91vEzjLUvnG/z3UXElcZXCeeEr8QDx6eJnNN5AeeBBPTcsrG0f
kWkKaxze9nG0tuNxZB0yD3P5lGCt+TBaC9EctPaDaO1Bc9H2tdHai8bjLXgM
3lkzr1XWbFXxZ6WzGpmNLHVu1go0Ywt+e+m5zRzfkt9aWjvQkM25dqW9aF10
1sOT0OY1+f3QaLTy/WitRjPR0nfJ0IU1FW18K1oL0Ui0851orUVDX9BnDdBr
DXXcpHBWfi06S5OZWYvMOeba16W19e1o7UVj8ZJPor0RT0Hb34vWbjSe7PFm
tDaTQdAONAStOTo4u78endXJ8GT3N/L9IMOTtbjn3OtPSmsdmreZ+NPS2oLG
oJ3fiKfqb+vp/q/M/kz+q19hv8Z/J4lriZ8XvyS+TFxd/Fz248vFNSrsz+Q5
tBbN1cvFWr32Xc6T+PkbzAnyhI43FM6f1BJkFLJJ3az1aD7LcE1p71id5+Zv
zEcd19Lzt9xvrnH0/OV1GK/9PHqt4bk1K/z98XX8fLLGthUv09jqwt6Bh+A9
P5fWWjQHrfmxtPbjAXjJT6W9Ds/Du34p7XV4Hlr7q/gKvX9tvf8LOn65sNaj
+Wjtt6W9Bs+pJP5IPEXn19X5K3T8amEvwBPwog167UqN19H4ch2/Uli70fAK
tLq0F+KJWOIHpbUITUI7/imtRb9FaxOahHb9Hq1NaNh41p94mHhw6VrvFfHl
hWs+apEV0bUGNQnehwfibQcF1yIr8/ykJqHWeTG6NqDmOUf8mc5rJ948uJZ6
OPp+k0mp1e6P1i40jNrqkejsToYfx/wVn8W9Kl07zhF3LVxDUgs+JO5YuCak
dnwgupakhqRWWRVdO1GzUKusjq6VqVkuEz8jPld8hvgS9CT/viHiSeLnxKPE
Z3E/uJ/i88XDxZeLnxWPFJ9ZurZbHl27UeORBcgE5Pq99fsvRb90fI7GTy9d
270UXTtR401E78Qj8I7StfwL0bUiNT2fzXegblsaXfsti671qQGphamJ+f3V
g2vrB6NrZWpsatN54l6Fa1Sy9Jfi5wtnarwQT2S+fF86W38RnbXJ2GTXz8QL
C2dYvG1o1tsJwbXdBXn9XZu9+6dob8XD8b6/or0UDyRr/CHevXDmIGv8LN6+
cObAe/+MziJ4MN6IR3YTb5GzyC/R2aRO/v1zo2tvegLUEuujtQcNejTrI99n
nfh8Pi94PtXT88Wsz2hvxCPHcr74TPGg0no1P7qWo6YbI35MfIZ4YOls8UN0
diVjUHsvEp9WuAanN7JQfFzhHgle/nf0esXTL0yuGXn/XYN7K69G92LosdAr
WBCdnchQ9H5eju7NjMn3lx4E97tG8OfRk+F6NA7+e2pk1lej4LnUJLoHwJxi
Pe+W9eDA4DHOYa6fh4ZEz3nmOmumYfQaYu0wpxtFrxnWCmsC/WiR1+8BwfO5
Wb5eXYPXZ+u8HnoFr5/m+Xp3Q3ej1wRr4ezC871pvj9d0OHoNcVaGpGvd5V8
fVprfEcdT9Zrw0qvYeYLGZb590eef2Reat2KnPXIfNTGlXLWI/ORpf7N6+V/
Naj4ypz1yHxksU15fQzIeeiinPXIfJ3ElXM2I6PRy6iWszUZm/n9e+ksSCak
97FVXj/4EXr9Q54vZEBq8S1zliXTMv/X5/U6POeLK8hJ0XOauTxUr9WK1ji0
jXtMfumQ9eWI4DzSI8+3AXn+tc/60ifnj32yvh0ZnD86Zn84Klhf22Y9PCR4
frbL+nlocL7plP3huOC81Tnrc79g/6HnyFw7LDg/dcn6eHxwXuuW5/uJwfms
a9bzE8iB0WtyYM635L2eOZ8NDM6zfXPeGhHsV62y3vYM1uM22Q8OFm8T7UF4
D3Nyq2hNQAuYwzWiPRTvZI6TP+mJkC/PCe6F0BO5E/8Jzsf0JMiXpwfnQXoM
5OPRwfmbHgN588zg/Nor571BOW9T45MnzwjO1/RQyItDg+f6knz9mPPkIXoQ
5PdLND5Ix//Ra7PI6oX7MT8n13/0uw4V/ySeVLoe5L3p15BF+Qzqxe+S+4vj
ucfiH8WXkxXQDPH34nHiyYXvNfUlWYB7zrWmPzk0r9e9xN+wpkr3L48S/yK+
Snxz4bn3bbK3MQfpb36d3K9FD/YTbxCfX7rfeYz4V3Ss/H//hH4K2smcniDe
SuHpafHywrU7NfyO6HXpXjM9553Ef5fuBVDjNhb/lfMJNRd5ZZvg/ELNSK1Y
BvdvftBr40v3w8kj9LTpjYTgvEIPm15I0vNEvELf5xmdv6KwX9MTp3cSyQsa
30bjz2p8ZeH8Qw+dXnjVnKeoEakNNwv2u3HZr6eLeyf3e59CW/TaeRqrJJ5b
uj91gvgPnXOT+DbWmPh38Y3iWYX7pX8l9+tuL9yv+i25P0U/6kTxn+KbS/fz
zhD/K75HPKfwWqO/RW3DmjtdvEnjd4sfFg8W/y2+TXxf4bVAf5ZakTUxRLxR
47NL98eYy/TXqEWZ0/Tb/knun9F/HSWu0O97hCzD9ReX4vtK93svElcVP0oW
Ltz/CxXuB9PfQzvpB1KboaFnxuJ/YftetJ01LY7ih8kS4pHiJJ6D97PGxZuJ
F5G9C/cvK1e4H0j/b7S4inghWV18sXhz8WN4UeHsRL+T2oUMNUa8hcYfL92f
5N7STyWLcY/pr25Z4f7lEq6/rsVbwbXpxtL682n2K3qy1N+bxPsW7vFSb10d
Xa/SI6Feu4o1XrhHSm1Ijfixju8NrlWoWT4S3xNcK1IzrhPfJ/4GbY2uxb7Q
Y0P0ngO13vrCWoPmMFdODe73zdfzZ4V7tvjjxmg/pYdB7UMN9KH47uC1wpp5
RzwzuNai5vpAfFdwbUeNtxbtD66Pro2uv+gZspZYU+8KZwXXrtSw74tnB9e2
1LjviW8Lrh2pIdmvoAZGO9CQtzU+I7jevSa6nqbHSC+KntQgHbcI7lXRs6KX
untwr4ueF73zlsF7Cewp0CvbMXivgT0HemU7Be9dsIdBL6xJ8N4HeyD0XncO
1mI0Ga0+Nlh70CC07+Rg7URD0aLewVqL5qLNhwdrNZqN1p4UvJdC8X6ynnYJ
3kshL9Kb2yF4r4Y9G3pv2wfvPbAHMVDcPLj3Rg+OvYzdgvcqqHv6lt6zYC+B
PQX2KhoE7+2wx0PvrmHw3gR7FOy11AreC2JPiF7edsF7NezZsPdTO3ivhT0X
9na2Da6vr9fnb1u4p0uvhp4N/F3p+n52dO+O+pf6/vZo/abnhzagEdTuNwSv
ZdY0vYLrg9cia5Ja/abgtcsapja/MVh70CC0fFrw2kcD6JVcHay9aDBaflWw
tqFxaNXYYC1D09DmccHagsbgVVODtRWNRVsuDdYiNAltnxJcy1DT9CxdQ6AF
aEKL0jX4Y9w75mDpa0KtQc3B/elfuFanZt+ndE1EXbA+ef+PPUK0EE1EOy8L
1k40FG2aFKy1aC5adnmwV+AZaPdw8S7RNSy1687Z6/F8sggb23gBnoCWjgnu
n8yK9mN65vSvZkb7NT1y+m0zov2ZHjT9rhuj+zX0+OnP3BXtt+xJ0F+5LVof
6bnTX7pbvFvhPQr6QXdE+zV7CvR77oz2Y/Y46K/dEN0fop7lWnJNK4s/Lt2f
uS66f0OPHe/EQ/Ge84O9Fs/FG0cFZzEyGV5ycbCX4ql4zXnB2o7G05u7Obj/
OT26H0aPn/7pTdH9K/Y86JfeEt2PYo+F/ujN0f0q9ljIRmQkvHZYcC+Yeors
1Tk4+5EByQJDgntx7MewpxYK62uv5F4gPTTmxnbJesYcYW+wYbLesUeIvjQh
TxbOhOjTzsn6Q0bCPw5P/n30HPGbw5LXFz1J/ObQ5PVFj4zvu31yLUNNQ++P
Hna1wj1A9q4aJOsPe1jspdZJ1jP2VNnLqp2sn+xpsfdaP1mv2YNlb7Zesp6x
R8teKnuqA8TNguvR9lnX2ItFf3dKzp5kYvzo4GS9oAfI3m3d5L9nD5e9s22T
9Zw9NPxw12T9JoOzV1Iz2Q+OzP53SLIeLc/59L7otcz/iSD7koHJ3scEZ3Ey
OdltcPaXKnocXHjPGj2vnLzXTn8NPS2T61P2wPGfzZL78+x5o+dbJ+szeyr4
2TbJesweDPpfLXmvmT1o/G+rZH1nTwb/qpqs5+zBo9c0K6mv2SNH7/+Nrr/Z
k8Yfqyf7C3tC7E2zR8371w/O/tQAXO/u2c8a6XFB4RoCf9w82Q/4PwLsJdVK
9lP2lPDDxsnznZqCvaoaevQrvGeF31Qk1+f8HwT8KSXX2/tl/4rJ/Sf2xPHD
Ssn1Ov+HgfXTNPl+TMp+0jlZr8mc+OsWyf7F/1nAT7om6yUZj6zRM7mXTuZA
P1sk+zmZm7zSI9lP6LmjD92S9XdxXt/Nk/2cGgi/2TfZX8j8+FPHZH0lI5PF
yeRk0YnZf9ol6zWZmjzUPdmv6PlTD90TXW+wJ4W/7ZfsV2RcagdqCLL9yOD6
6t7oeoU9Uvob/0T3U9jzIl+xx7CmcI8eP2ybrE+zc147KNlPqYHwvz2T358a
A/9ohgYWrhnxww7iRYUzOnq7e/J6oMYhe5PB0a/rsl/tnay31Aj4217J+kzm
x9/aJ+s3NQJ+vE+yXs/JefCAZP+lpiSfHpjs19R06Htrrlnhmht9b5msr9RM
5ItOydefmgH/a5Psh9RE1G7k52WFazj8pFVyfqYmI390SfZjaoj/Aj92vX4=

          "]], Polygon3DBox[CompressedData["
1:eJwt12eYVOUVAODZnQGkiYhAEFCRYqQZqYYauqFLTWhGqlGKSjNSlS5NBYFI
01gTUWOoKqChiQVQEVBpiu3BAoiUKC3vee78OM93znvP7O7ce7+y5foO6zQ0
N5VKvS3yiE2ZVOrXdCpVEFbNSaWqiS/lL7MCxirqWvoeVH8rnmBvsM3sDvUX
4gG2mD3NjshfYvmNlVlN9q78H+yk+gpRmm1ly9kP6owowd5hT7Gf1JeLq9lG
9gT7Rl0snK3zN59gP6rzicvE+/qeYafkV4my+o6qF4tr5CXYi3q2iyfZCXUh
Ucq1HexZdjon6Yv+/uoafs/4uBdsDdsg/zv/Wl1UFGbr2SL2lbpI/Ez2BlvI
voweUZB9xJ5jv8T3EuXZx+wFdk5dRlRk0/zOa+NvFmXkDfVsFkvV3+lJxfcz
bGHL2PfqXFGcjdFfWFR2rQf7s5iuLqfvOlFW3si1TWKJ+qjrF32umJjsWin2
G1FKXk/Pe+Jp9c/6roy/Ud+N6mtdb+LavexvrC87IsbEz2XPsPvimYmZbAVb
xfqof+ez98f3Zq+w3uwmNpo9z15mvVh1NiruGXuJdWdV2T3x/NgLrD27kd0Z
z4AtZ11ZFTYs3i/2PPte/iovZqzFGrB26t/qGxTPii1jK9TVxMh4H+J7iGfV
N4iBbIH68bhn6tpiUm62RzTy8yqwP7Hp6rl+XktWif2FPcIWsVasIrudzWUL
2RJ1edGdTWNTxS36yrGObGI8R323qR/iH8sXx30Vd8W7Jx6KZ6VvhWt11dfr
vY09yGawe+O9FQ+zF9lK1k1dWV8sCMvYc2yKujS/OkJe37XLLBSN5Y/JV8az
FBPVRfVcIUrE/XDtTPb+jchN7lPcr56ufy6Gs/nqpX5HD/VhcR+bx5awzuqD
YjCbHfeZ9Ys5JMbG+x/3mT3q518nbmX3s9FivfqUvvMx72P+xj1kZWMOidLy
BvpnxH2OexPzS96YdZIfEHfLZ/nc/LjP6v3iLjaTzWNd1IfEEDYnvhsb72cU
jHVSXCWv6dpDxuJRx1yS141nYCymvlKUlNdh42KuqgvFeymvwTJxP+Meynv7
Hf39jjQrwqqyXqwfyxs/i93EbmcDYy1nRVl11ocNYLnsclaF9WR92QR1kXBR
XF7Ltd0x99S/xtosKsRaom7i+uOubY31lX0q/yfPMV7PKrOf5Wl9ZYxtWWd2
Kr5L3H9jO9aFnZPnZxWMXVgvdr86N9431oa1Fj/KL8b9MzZXt9X3Q3Z+dMtN
nnc89wZ6PhRd2Ej1eH111DtFB3ZP9LEL8gI+W9HYlfVmx2IDZCWMLVn7+B3y
S3FPjC1YO3ZcnqOvpLEV68B2yn/Wd1Z9g7iJnWZ59F1jbM+6soZ6PhJd2Sg2
gQ3OvqcXsvts7Ldfu34+noexkbqFvu/k/45nbKzJ6rO9Mb/Ypfidog7bw06z
i+rqojb7la2Kd9zYjLVh38gvxHtkbMxasrfkx2NPV5cU17LzMdfjXhs7x9xl
Z+X5Ys4ZO7Lu7ER2TRyQm6xVsWadkeeNvcvYQd1N3y+xXsQ6Zmwaz5ftkp9k
Z9SV4u9m9dUfiM6ujWDjYk/J7vHHsvM55vU2139ix+Pvjc+nkvn8VjrZ177J
JHO7lnhftGND9Y1KJXvgm+lkP/06k+yHR7NrSSt2p76/xjvNt4imrJ96cJwj
5K/E/DLWYPVYRfV/08ka8m281/Ka4j3Rlg3RN1JfFfVW0Yz1Z0NYDfW7og0b
zEaweupdohMbzsbGnufnHov1Myc548RZJ84GG9PJevZVJjknNBN7RK9YV/RM
8dmmsU+Inmwsm8yaqHeLHmwMm8TaqD8V/dhkNpu1VX8m+rMpbA5rrt4rerPx
sb6y1upPRF82ic1itdU7RHs2LNZnVkm9SfyB3RH3mVVXbxct2SB2D6umflu0
YAPZMNZCvU/0YRPiXWPDs2fWQqxBzA/xlfwcK2xsqG6u75D8LMtjrM0asiHq
W3x+GlvLNrJXs/v55Nxkb4s97m59ddgUtjr2mJi/rC6bytawDey12Ov4F+r8
Eexm9TuitWt3s+FsXDwnn13I3mY72Sr5o7HXqC8kS1LqgbhPkvmubeHvsZVx
fohnoj4Xe4KYFeucvudd28cPsdmsI3uBfcIOs9Xyx2KPU1/KnvkHxruqb0K8
13wtG8p+z6azdexNNoDdzMaxl9hqtj/ODzxtrBBrGDsgX8EyxoqxtrHX5QvY
EXVBUYCtyb7Ph9V5RN5499X1+Cz969Vb4vmy+mw228C2xr3KJGeQeWwz2yQG
xVxnE2O9VK/TNzrWAvYIe4ttZyNjvrA5bGOsI2xUrAVsLnszJ/l/aGa8t+w5
tpcdZGvl82PvV6djDsb3YPPYQXVOfI9Usr+/lk7OQYcyyV4f56j16eQsfSST
nKnizLQuney7BzPJ+amkeCOdnBU+j30pnZwhXk8nZ4XDmeQ8EeeZDenkTPZl
JjnbnJT/x1jKWM/f0szfclj+C8tnrMsas/zq1fE8WSfWgx2U/4/lNdZhjdi7
2T356txkXsX82pP9v+LG3GSPjr06zj1r0sm5ZX8mOQPFuWJtOjkHHcgkZ4xp
MQ/kT7EPfe7TODex5mwR2852sbGx5rIFbBvbwabHGsmeZrvZZ+zhWNfZs2wP
O8BmxJrLnmEfs/3sM/m/eK6xPKuSSs7OLfQtYe/nJOfoSbEGsaVsR3xXdlFe
kFUydmN92GR9rdgytpPtZZfkhWIvNHZnt7Op+v7InmQfsE/iTM1uZcvZLraP
VWXb4j6wAWwo+z9nGNI8
          "]]}]}, {}, {}, {}, {}}, {{}, {}, {}, {
       Line3DBox[{251, 474, 1106, 252, 488, 280, 1211, 503, 295, 1216, 518, 
        310, 1223, 533, 325, 1233, 548, 340, 1244, 563, 1150, 355, 578, 1164, 
        370, 1259, 593, 385, 1263, 608, 400, 1268, 623, 415, 1275, 638, 430, 
        1284, 653, 445, 1295, 668, 460}], 
       Line3DBox[{253, 475, 1107, 254, 489, 1118, 281, 504, 296, 1217, 519, 
        311, 1224, 534, 326, 1234, 549, 341, 1245, 564, 1151, 356, 579, 1165, 
        371, 594, 1176, 386, 609, 401, 1269, 624, 416, 1276, 639, 431, 1285, 
        654, 446, 1296, 669, 461}], 
       Line3DBox[{255, 476, 1108, 256, 490, 1119, 282, 505, 1129, 297, 520, 
        312, 1225, 535, 327, 1235, 550, 342, 1246, 565, 1152, 357, 580, 1166, 
        372, 595, 1177, 387, 610, 1186, 402, 625, 417, 1277, 640, 432, 1286, 
        655, 447, 1297, 670, 462}], 
       Line3DBox[{257, 477, 1109, 258, 491, 1120, 283, 506, 1130, 298, 521, 
        1137, 313, 536, 328, 1236, 551, 343, 1247, 566, 1153, 358, 581, 1167, 
        373, 596, 1178, 388, 611, 1187, 403, 626, 1194, 418, 641, 433, 1287, 
        656, 448, 1298, 671, 463}], 
       Line3DBox[{259, 478, 1110, 260, 492, 1121, 284, 507, 1131, 299, 522, 
        1138, 314, 537, 1143, 329, 552, 344, 1248, 567, 1154, 359, 582, 1168, 
        374, 597, 1179, 389, 612, 1188, 404, 627, 1195, 419, 642, 1200, 434, 
        657, 449, 1299, 672, 464}], 
       Line3DBox[{261, 479, 1111, 263, 493, 1122, 285, 508, 1132, 300, 523, 
        1139, 315, 538, 1144, 330, 553, 1147, 345, 568, 1155, 360, 583, 1169, 
        375, 598, 1180, 390, 613, 1189, 405, 628, 1196, 420, 643, 1201, 435, 
        658, 1204, 450, 1300, 673, 465}], 
       Line3DBox[{265, 481, 1112, 266, 1209, 495, 287, 1213, 510, 302, 1219, 
        525, 317, 1227, 540, 332, 1238, 555, 347, 1250, 570, 1156, 362, 585, 
        377, 1261, 600, 392, 1265, 615, 407, 1271, 630, 422, 1279, 645, 437, 
        1289, 660, 452, 1302, 675, 467}], 
       Line3DBox[{267, 482, 1113, 268, 496, 1123, 288, 1214, 511, 303, 1220, 
        526, 318, 1228, 541, 333, 1239, 556, 348, 1251, 571, 1157, 363, 586, 
        1170, 378, 601, 393, 1266, 616, 408, 1272, 631, 423, 1280, 646, 438, 
        1290, 661, 453, 1303, 676, 468}], 
       Line3DBox[{269, 483, 1114, 270, 497, 1124, 289, 512, 304, 1221, 527, 
        319, 1229, 542, 334, 1240, 557, 349, 1252, 572, 1158, 364, 587, 1171, 
        379, 602, 1181, 394, 617, 409, 1273, 632, 424, 1281, 647, 439, 1291, 
        662, 454, 1304, 677, 469}], 
       Line3DBox[{271, 484, 1115, 272, 498, 1125, 290, 513, 1133, 305, 528, 
        320, 1230, 543, 335, 1241, 558, 350, 1253, 573, 1159, 365, 588, 1172, 
        380, 603, 1182, 395, 618, 1190, 410, 633, 425, 1282, 648, 440, 1292, 
        663, 455, 1305, 678, 470}], 
       Line3DBox[{273, 485, 1116, 274, 499, 1126, 291, 514, 1134, 306, 529, 
        1140, 321, 1231, 544, 336, 1242, 559, 351, 1254, 574, 1160, 366, 589, 
        1173, 381, 604, 1183, 396, 619, 1191, 411, 634, 1197, 426, 649, 441, 
        1293, 664, 456, 1306, 679, 471}], 
       Line3DBox[{275, 486, 1117, 276, 500, 1127, 292, 515, 1135, 307, 530, 
        1141, 322, 545, 1145, 337, 560, 352, 1255, 575, 1161, 367, 590, 1174, 
        382, 605, 1184, 397, 620, 1192, 412, 635, 1198, 427, 650, 1202, 442, 
        665, 457, 1307, 680, 472}], 
       Line3DBox[{277, 682, 1308, 683, 278, 501, 1128, 293, 516, 1136, 308, 
        531, 1142, 323, 546, 1146, 338, 561, 1148, 353, 576, 1162, 1163, 368, 
        591, 1175, 383, 606, 1185, 398, 621, 1193, 413, 636, 1199, 428, 651, 
        1203, 443, 666, 1205, 458, 685, 686, 687}], 
       Line3DBox[{459, 667, 684, 1294, 444, 652, 1283, 429, 637, 1274, 414, 
        622, 1267, 399, 607, 1262, 384, 592, 1258, 369, 577, 1256, 354, 1149, 
        562, 1243, 339, 547, 1232, 324, 532, 1222, 309, 517, 1215, 294, 502, 
        1210, 279, 487, 1207, 250, 473, 681, 688}], 
       Line3DBox[{466, 674, 1301, 451, 659, 1288, 436, 644, 1278, 421, 629, 
        1270, 406, 614, 1264, 391, 599, 1260, 376, 584, 1257, 361, 569, 1249, 
        346, 554, 1237, 331, 539, 1226, 316, 524, 1218, 301, 509, 1212, 286, 
        494, 1208, 264, 480, 1206, 262}]}, {
       Line3DBox[{690, 1000, 473, 689, 1106, 912, 691, 1107, 913, 692, 1108, 
        914, 693, 1109, 915, 694, 1110, 916, 695, 1111, 1006, 1206, 696, 1112,
         917, 697, 1113, 918, 698, 1114, 919, 699, 1115, 920, 700, 1116, 921, 
        701, 1117, 922, 702, 1308, 1001, 923, 1002}], 
       Line3DBox[{704, 1007, 1207, 703, 488, 705, 1118, 924, 706, 1119, 925, 
        707, 1120, 926, 708, 1121, 927, 709, 1122, 1008, 1208, 710, 1009, 
        1209, 711, 1123, 928, 712, 1124, 929, 713, 1125, 930, 714, 1126, 931, 
        715, 1127, 932, 716, 1128, 933, 717}], 
       Line3DBox[{719, 1010, 1210, 718, 1011, 1211, 720, 504, 721, 1129, 934, 
        722, 1130, 935, 723, 1131, 936, 724, 1132, 1012, 1212, 725, 1013, 
        1213, 726, 1014, 1214, 727, 512, 728, 1133, 937, 729, 1134, 938, 730, 
        1135, 939, 731, 1136, 940, 732}], 
       Line3DBox[{734, 1015, 1215, 733, 1016, 1216, 735, 1017, 1217, 736, 520,
         737, 1137, 941, 738, 1138, 942, 739, 1139, 1018, 1218, 740, 1019, 
        1219, 741, 1020, 1220, 742, 1021, 1221, 743, 528, 744, 1140, 943, 745,
         1141, 944, 746, 1142, 945, 747}], 
       Line3DBox[{749, 1022, 1222, 748, 1023, 1223, 750, 1024, 1224, 751, 
        1025, 1225, 752, 536, 753, 1143, 946, 754, 1144, 1026, 1226, 755, 
        1027, 1227, 756, 1028, 1228, 757, 1029, 1229, 758, 1030, 1230, 759, 
        1031, 1231, 760, 1145, 947, 761, 1146, 948, 762}], 
       Line3DBox[{764, 1032, 1232, 763, 1033, 1233, 765, 1034, 1234, 766, 
        1035, 1235, 767, 1036, 1236, 768, 552, 769, 1147, 1037, 1237, 770, 
        1038, 1238, 771, 1039, 1239, 772, 1040, 1240, 773, 1041, 1241, 774, 
        1042, 1242, 775, 560, 776, 1148, 949, 777}], 
       Line3DBox[{781, 1149, 950, 779, 1150, 951, 783, 1151, 952, 785, 1152, 
        953, 787, 1153, 954, 789, 1154, 955, 791, 1155, 569, 793, 1156, 956, 
        795, 1157, 957, 797, 1158, 958, 799, 1159, 959, 801, 1160, 960, 803, 
        1161, 961, 805, 1163, 963, 807}], 
       Line3DBox[{806, 962, 1162, 804, 1255, 1054, 802, 1254, 1053, 800, 1253,
         1052, 798, 1252, 1051, 796, 1251, 1050, 794, 1250, 1049, 792, 1249, 
        568, 790, 1248, 1048, 788, 1247, 1047, 786, 1246, 1046, 784, 1245, 
        1045, 782, 1244, 1044, 778, 1243, 1043, 780}], 
       Line3DBox[{809, 1055, 1256, 808, 1164, 964, 810, 1165, 965, 811, 1166, 
        966, 812, 1167, 967, 813, 1168, 968, 814, 1169, 1056, 1257, 815, 585, 
        816, 1170, 969, 817, 1171, 970, 818, 1172, 971, 819, 1173, 972, 820, 
        1174, 973, 821, 1175, 974, 822}], 
       Line3DBox[{824, 1057, 1258, 823, 1058, 1259, 825, 1176, 975, 826, 1177,
         976, 827, 1178, 977, 828, 1179, 978, 829, 1180, 1059, 1260, 830, 
        1060, 1261, 831, 601, 832, 1181, 979, 833, 1182, 980, 834, 1183, 981, 
        835, 1184, 982, 836, 1185, 983, 837}], 
       Line3DBox[{839, 1061, 1262, 838, 1062, 1263, 840, 609, 841, 1186, 984, 
        842, 1187, 985, 843, 1188, 986, 844, 1189, 1063, 1264, 845, 1064, 
        1265, 846, 1065, 1266, 847, 617, 848, 1190, 987, 849, 1191, 988, 850, 
        1192, 989, 851, 1193, 990, 852}], 
       Line3DBox[{854, 1066, 1267, 853, 1067, 1268, 855, 1068, 1269, 856, 625,
         857, 1194, 991, 858, 1195, 992, 859, 1196, 1069, 1270, 860, 1070, 
        1271, 861, 1071, 1272, 862, 1072, 1273, 863, 633, 864, 1197, 993, 865,
         1198, 994, 866, 1199, 995, 867}], 
       Line3DBox[{869, 1073, 1274, 868, 1074, 1275, 870, 1075, 1276, 871, 
        1076, 1277, 872, 641, 873, 1200, 996, 874, 1201, 1077, 1278, 875, 
        1078, 1279, 876, 1079, 1280, 877, 1080, 1281, 878, 1081, 1282, 879, 
        649, 880, 1202, 997, 881, 1203, 998, 882}], 
       Line3DBox[{884, 1082, 1283, 883, 1083, 1284, 885, 1084, 1285, 886, 
        1085, 1286, 887, 1086, 1287, 888, 657, 889, 1204, 1087, 1288, 890, 
        1088, 1289, 891, 1089, 1290, 892, 1090, 1291, 893, 1091, 1292, 894, 
        1092, 1293, 895, 665, 896, 1205, 999, 897}], 
       Line3DBox[{911, 1005, 685, 910, 1307, 1104, 909, 1306, 1103, 908, 1305,
         1102, 907, 1304, 1101, 906, 1303, 1100, 905, 1302, 1099, 904, 1301, 
        1300, 1098, 903, 1299, 1097, 902, 1298, 1096, 901, 1297, 1095, 900, 
        1296, 1094, 899, 1295, 1093, 898, 1105, 1294, 1003, 1004}]}}},
    VertexNormals->CompressedData["
1:eJztXXtQVccdvlhtEQwxCaIYnWBtrFqNjE0sGvSogRobo4NgEmvUVkEeF7i8
LghG0EQEJDERohaUpFWjIRqDqI1K9NpMojhjVR4Zy0mKj6Byb9pJawxkEiP1
fOdsF5hZ4OC59+6F89fO9/lj7+5vf/v6dvc4YplpfkQfg8EwrJ/B8JO76dRx
AxKO+9kscioKJJ2m8NPaYDavVT5q+ZzsreCr40LjKBaFb6vXmCR+r+2TJIrZ
PCuf3HZ8rsIT+yrG75J8ye/Y2w/25htDQhIl3moqj6aYbf/mT0/HS/ymxohU
itXn/2TJ58kSP7OuJFzCQkVgrJS6hZXFSfzhZwtXSvhD4y2H+EFQeKEN5o9X
Wy/r5bXwp/VjPzOwkl6o/znaZejpesQ5SbXq71qV3978ieRfIJ6PnipIp1h9
+e/3z0yReN8fEv4o4YfmXkI8j5v4T9g3NJVHSfjq7BEpjvAnb3GrVflZ9qx4
ZvG6n7X1sx7Pup9d0c8s/v30bzEvLNs0Ges8kup+1v3Ms59Z/tT97Jh4Zq33
okcXGyV+xeI9mRKO8pic3J3yHG/4BvPCyU3ZKyjmb12tlT9Z+xHW/FiQtiVJ
4rdsbo6RMEm1Kj9vvL3bkbUfn4r9us0yJStjlYQDsZ/XTh9g2avN31V434D9
8MN+91FGCS/3i0G9WHpUZ/pVV3Uzlj2Lt7cf7K1bsvoRyz69/2XEm7fYaKZY
FG48dwb9IkDwQb/YHn9fYkf2hW7zMC4F+fbBvFBv9TB1pzwsntVechzZLOsi
J61sHVdq/VO26/fI54aiH509eBz1/exIPeoV9dL8CAlXlK7qsF5H8e82Sy3+
jmBReGPDEXleW/4A4j911QKUM3NNKuwb0jbDrvj1uqSO8mH97pwXC8CXhx4A
3ho0N64j3lXWIbzNR1rxS4auB//n7PfiKBaFX78/GnGSsmd1DMV34yLrFOLB
nB/VCotCqXkx8hlz3BxPMXs8PHRrOvibM/qmUEz9L7Tzf3RRNOzJek7G7P5I
6vVWu3qxxvMTynpv+u/iolr7hzUvDHrcB7p6Vc6wSAnPO1ATy2P7ugpP2teq
rEtIutvQhPHf9EoA4oSkvJXf3vzXL3yOeeHsHEMqxaJwZ94Q+O2l9PMpFLPz
If26aWlEbOt+FzbQA7p64ZSFyyS88bXTDhmfexv/Rp8J4Mm4SVIWz1v5XYXX
/an7uTfw/1pgxfh/9R8hwCTlrZyuzut+dow/dT87hv9L3kzoilsuNEJXfP7B
57Cvmr33a6wzc78rB65cfrHDfDLdDVhPGi49baT47vqxJRB8meeQGIr584O9
+Yy+UdBtTq0Kg39I+pBlP+yJrkpS3srvKvzzjxkRz55FwxHP932Yifi1KPeO
Ar32prW2f2RaAcaZ+uq4JIppu3gr94pkLApn0A9slg0+3rEUi0Lt1WPIJ3Tb
jkSKaT/62YFiM8X8+U0rXq1OzsqHjP/iomOtcOc6z9vtdB5n+cFZ+qGzfpel
E5J9wc7bi1vtb0Uh4OAAtO+ivP5mikWhqJ2uWKTsI8g573avZiPFoiDb0XOc
JkVX5K1faMWz9FWWPUsvZdm/q+Rf2ncjxrfRjwVjXIzPXged88emT3HucOx2
Qq/UOVn6W3f7S3v97YFn/g2dM6p/obzuqjGYefSDs3i1OieLt376IHTObYfu
QOecNLWlV8azvXld53Sun3krp87r/L3wui7kGH/qfnYMv6mgH/SK6zt9V0t4
1sO+yd3Jh6W/8VZfZ/HX5i6F/nbS7QPsTx85OwzrPbU6J4vXdU45PZBVh3gO
cD+N930ZfYbBr0Qfqw7dABw2tqbDewUs/Y1lz9LfePOPVvy63WPAX61caaKY
3rNqq4+w9TRyjyul3T2uoHc8Ec9+572MFLP1nJ6qc7L8POfIdownaUX/wT2x
cxuvIU7/jtRmuTnIJ4Pa8Vcv3vjFQ6Yg3va8PiqGYlEY9LfvEZ8DgpYnUCwK
ftPrEW+1lf81UywKkVuHI/9E78Mmitn5Bw9Nhq5Vm38Iulaud2WP0AGITji2
izohGT/Dujh+svIn+kyWsv64c/IG4n944G9h/3FjPuxm+u/q1rtv3vysltdK
P1xWcAnx7H69GuNPQ10z2ss2vxn234XvRJxffKGyw3VIT/WzVjqhrH/bLPsK
T8DPRHcn6+2W+q6tt5/4IAV8yKgnUQ6S8uY33njWuzOWPTknIu8uSMpbvXjj
1fpZ53WeB56Mw+Q7Mp19T4bFtyjzBXnfS1Le6uvqPOtdPG/l5I2X1ws2i/xu
lbxf1f2sNf906c42+glJZb3KZqn+ZCHujxGdiegqI4NSMV/uiBuv67f3wPue
CkDcpicFw8/FmX7grSr3NSy+t+m0rP3darJ/UfZFWVv9Eb+sfbpafbW7+kNP
029lPclmub/kMPTwWdCbRCHJxw31nVz0lJlitn7lMXMw2muGpTaGYrY9K3/e
/OMqPOt+I8v/LH2SpWeq1Sf/MOM25qNzZY0YJ8eEmhziB7X3D9X2aznfu/kr
9w5Lr4+DfmP0qoJ/duVvxPxI5ltnxcOO+BroQnnjv8J3gSJGDIQ/Dn5TjPqu
GPkiyhc8eK1D3ikn5lxAeR6+sgfv0P/ksbRb/mHpVF/i/N9mIfdre/o9W53X
+Xvh1e7XdN65POs8t+LNv2KeHeyRgu/5JO77Hvas+Tey5SLmr20GE8b9pzwz
uKyvzvdOnpyPhCvnIuH6+YjOa8CTe0Rtv2enXf5El/uVoseR1N718qyahXE+
PC0B4/+1Z49yqaNO3Pcy9keR5a/iu47Bu+Xvq9VZg1D+5rHjoT8M+Y07l+Vn
8f/XVw8HY/1E7p3J+0mbZYFyj/CzxK+MHeWjVo/ijc9RdCSvLupIRGfw7KLO
wNIxXOVen6vcu2PpSER3vf6FNZ5iUcjPm4B2PLlwUQrFnX9ftP136ibu/QK/
+2P0FaxvD9685dT9u711pBrFntiRv6swvQq+35nBsFu3tsUh+7IJT5zH75Lz
Bhl3Hift34e+k2QEv/bRCozzS3wvOzWeyXct3G+sRFxlD/TXdHxQ+65Nfwen
8zrvvHVaT3svNm/kbPDknq6M+SunVvyGnC/RjuS7NjK2/+/euvA4eDLPyZg/
/ziLV7t+6Km8Pr/rvCN5e89rzpo3yf8zkpb4Ns4vmmc1cnk+FfLLK9AfSq+t
x/43b/1U8Grf7fI2j0dsrgLvP8mM8ats9WvYr5N9MdkPk/0xKx9X17VYuk3J
K0vQvg0ncvFOiughRAceWOGTRDG9txBb4hZPsSjcPNcI/qM7HxkpFoVy6B82
yw/PrGmFXUfvYvGsc8z/AfApDqw=
     "]]},
  AutomaticImageSize->True,
  ImageSize->{305.3318805704282, 300.52790471774347`},
  PlotRange->{{-2, 2}, {-2, 2}, {-2, 2}},
  ViewPoint->{-0.3368369794216185, -3.3170204128826195`, 0.5778550249103456},
  ViewVertical->{-0.11840380202157647`, -0.43546560625736913`, 
   0.8923845838167181}]], "Output",
 CellChangeTimes->{3.6255781404667263`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 "It turns out that if the columns of A are independent (in which case they \
form a basis for Col(A)), then ",
 Cell[BoxData[
  FormBox[
   SuperscriptBox["A", "T"], TraditionalForm]],
  FormatType->"TraditionalForm"],
 "A is invertible.  In this case we have the formula\n\tP = ",
 Cell[BoxData[
  FormBox[
   RowBox[{"A", " ", 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{
       SuperscriptBox["A", "T"], " ", "A"}], ")"}], 
     RowBox[{"-", "1"}]], 
    SuperscriptBox["A", "T"]}], TraditionalForm]],
  FormatType->"TraditionalForm"],
 "\nfor orthogonal projection onto Col(A).  This is one way to see that \
orthogonal projection onto a subspace is a linear transformation.  "
}], "Subsubsection",
 CellChangeTimes->{{3.6255690327439528`*^9, 3.625569052210061*^9}, {
  3.625569097414259*^9, 3.625569154676292*^9}, {3.6255708997960787`*^9, 
  3.625571042415386*^9}, {3.6255711279461317`*^9, 3.625571146298271*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"MatrixForm", "[", 
  RowBox[{"P", "=", 
   RowBox[{"A", ".", 
    RowBox[{"Inverse", "[", "ATA", "]"}], ".", 
    RowBox[{"Transpose", "[", "A", "]"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.625571049973566*^9, 3.625571073956952*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      FractionBox["1", "3"], 
      FractionBox["1", "3"], 
      FractionBox["1", "3"]},
     {
      FractionBox["1", "3"], 
      FractionBox["5", "6"], 
      RowBox[{"-", 
       FractionBox["1", "6"]}]},
     {
      FractionBox["1", "3"], 
      RowBox[{"-", 
       FractionBox["1", "6"]}], 
      FractionBox["5", "6"]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.625578740470354*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"P", ".", 
  RowBox[{"{", 
   RowBox[{"1", ",", "1", ",", "1"}], "}"}]}]], "Input",
 CellChangeTimes->{{3.625571080722232*^9, 3.625571084068117*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"1", ",", "1", ",", "1"}], "}"}]], "Output",
 CellChangeTimes->{3.625578763822506*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 "Projection matrices P are distinguished by two special properties.  They \
are `symmetric\[CloseCurlyQuote]: P = ",
 Cell[BoxData[
  FormBox[
   SuperscriptBox["P", "T"], TraditionalForm]],
  FormatType->"TraditionalForm"],
 " and `idempotent\[CloseCurlyQuote]: ",
 Cell[BoxData[
  FormBox[
   SuperscriptBox["P", 
    RowBox[{"2", " "}]], TraditionalForm]],
  FormatType->"TraditionalForm"],
 "= P."
}], "Subsubsection",
 CellChangeTimes->{{3.6255711578256197`*^9, 3.62557122797753*^9}, {
  3.625571259312985*^9, 3.625571288160906*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"MatrixForm", "[", 
  RowBox[{"P", ".", "P"}], "]"}]], "Input",
 CellChangeTimes->{{3.6255710856944323`*^9, 3.6255711037163687`*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      FractionBox["1", "3"], 
      FractionBox["1", "3"], 
      FractionBox["1", "3"]},
     {
      FractionBox["1", "3"], 
      FractionBox["5", "6"], 
      RowBox[{"-", 
       FractionBox["1", "6"]}]},
     {
      FractionBox["1", "3"], 
      RowBox[{"-", 
       FractionBox["1", "6"]}], 
      FractionBox["5", "6"]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.6255788489166875`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"MatrixForm", "[", 
  RowBox[{"Transpose", "[", "P", "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.6255711100943747`*^9, 3.625571121589094*^9}}]
}, Open  ]]
}, Open  ]]
},
WindowSize->{1276, 774},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
Magnification->1.5,
FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (January 25, 2013)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 161, 3, 130, "Section"],
Cell[CellGroupData[{
Cell[765, 29, 190, 4, 58, "Subsubsection"],
Cell[CellGroupData[{
Cell[980, 37, 595, 15, 104, "Input"],
Cell[1578, 54, 121, 2, 44, "Output"],
Cell[1702, 58, 121, 2, 44, "Output"],
Cell[1826, 62, 717, 21, 110, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[2592, 89, 502, 10, 141, "Subsubsection"],
Cell[CellGroupData[{
Cell[3119, 103, 194, 5, 45, "Input"],
Cell[3316, 110, 146, 4, 44, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[3511, 120, 632, 16, 100, "Subsubsection"],
Cell[CellGroupData[{
Cell[4168, 140, 1064, 29, 133, "Input"],
Cell[5235, 171, 671, 20, 85, "Output"],
Cell[5909, 193, 694, 21, 85, "Output"],
Cell[6606, 216, 179, 6, 65, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[6834, 228, 296, 5, 100, "Subsubsection"],
Cell[CellGroupData[{
Cell[7155, 237, 1815, 50, 133, "Input"],
Cell[8973, 289, 38325, 621, 473, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[47347, 916, 945, 23, 229, "Subsubsection"],
Cell[CellGroupData[{
Cell[48317, 943, 262, 6, 45, "Input"],
Cell[48582, 951, 930, 30, 146, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[49549, 986, 172, 4, 45, "Input"],
Cell[49724, 992, 125, 3, 44, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[49898, 1001, 555, 16, 100, "Subsubsection"],
Cell[CellGroupData[{
Cell[50478, 1021, 156, 3, 45, "Input"],
Cell[50637, 1026, 932, 30, 202, "Output"]
}, Open  ]],
Cell[51584, 1059, 167, 3, 45, "Input"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)