STEM Mathematics.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "[Oregon Curriculum Network](http://www.4dsolutions.net/ocn) <br />\n",
    "[Discovering Math with Python](Introduction.ipynb)\n",
    "\n",
    "\n",
    "Crystal Ball Sequence\n",
    "============\n",
    "The face-centered cubic (FCC) lattice is not always presented in this simplest form, ditto the cubic close packing (CCP), which amounts to the same thing.  A nuclear ball is surrounded by a layer of twelve, all touching it, and adjacent neighbors.  The shape so formed is not a cube, but a cuboctahedron, with eight triangular faces and six square.  This is where I can type stuff.\n",
    "\n",
    "As the cuboctahedral packing continues to expand outward, layer by layer, the cumulative number of balls or points forms the Crystal Ball Sequence.\n",
    "\n",
    "cubocta(), a generator, yields the number of balls in each successive layer of the cuboctahedron, according to a simple formula derived by R. Buckminster Fuller, a prolific inventor and philosopher [1].  cummulative( ) delegates to cubocta( ) while accumulating the number in each layer to provide a running total.\n",
    "\n",
    "![Image of Cubocta](http://www.4dsolutions.net/ocn/graphics/cubanim.gif)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "=== Crystal Ball Sequence ====\n",
      "  Layers     Points  \n",
      "         1          1\n",
      "         2         13\n",
      "         3         55\n",
      "         4        147\n",
      "         5        309\n",
      "         6        561\n",
      "         7        923\n",
      "         8       1415\n",
      "         9       2057\n",
      "        10       2869\n",
      "        11       3871\n",
      "        12       5083\n",
      "        13       6525\n",
      "        14       8217\n",
      "        15      10179\n",
      "        16      12431\n",
      "        17      14993\n",
      "        18      17885\n",
      "        19      21127\n",
      "        20      24739\n",
      "        21      28741\n",
      "        22      33153\n",
      "        23      37995\n",
      "        24      43287\n",
      "        25      49049\n",
      "        26      55301\n",
      "        27      62063\n",
      "        28      69355\n",
      "        29      77197\n",
      "        30      85609\n"
     ]
    }
   ],
   "source": [
    "from itertools import accumulate, islice\n",
    "\n",
    "def cubocta():\n",
    "    \"\"\"\n",
    "    Classic Generator:  Cuboctahedral / Icosahedral #s\n",
    "    https://oeis.org/A005901\n",
    "    \"\"\"\n",
    "    yield 1  # nuclear ball\n",
    "    f = 1\n",
    "    while True:\n",
    "        elem = 10 * f * f + 2  # f for frequency\n",
    "        yield elem    # <--- pause / resume here\n",
    "        f += 1\n",
    "        \n",
    "def cummulative(n):\n",
    "    \"\"\"\n",
    "    https://oeis.org/A005902 (crystal ball sequence)\n",
    "    \"\"\"\n",
    "    yield from islice(accumulate(cubocta()),0,n)\n",
    "\n",
    "print(\"{:=^30}\".format(\" Crystal Ball Sequence \"))\n",
    "print(\"{:^10} {:^10}\".format(\"Layers\", \"Points\"))\n",
    "for f, out in enumerate(cummulative(30),start=1):\n",
    "    print(\"{:>10} {:>10}\".format(f, out))\n",
    "   "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Octet Truss\n",
    "======\n",
    "\n",
    "When adjacent CCP ball centers interconnect, what do you get?  Why the octet truss of course, a well known space frame, used a lot in architecture.  Alexander Graham Bell was fascinated by this construction.[2]\n",
    "\n",
    "<a data-flickr-embed=\"true\"  href=\"https://www.flickr.com/photos/kirbyurner/23636692173/in/album-72157664250599655/\" title=\"Business Accelerator Building\"><img src=\"https://farm2.staticflickr.com/1584/23636692173_103b411737.jpg\" width=\"500\" height=\"375\" alt=\"Business Accelerator Building\"></a><script async src=\"//embedr.flickr.com/assets/client-code.js\" charset=\"utf-8\"></script>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "[1] Siobahn Roberts. *King of Infinite Space*. New York: Walker & Company (2006). pp 179-180.\n",
    "\n",
    "\"Coxeter sent back a letter saying that one equation would be 'a remarkable discovery, justifying Bucky's evident pride,' if only it weren't too good to be true.  The next day, Coxeter called: 'On further reflection, I see that it *is* true'. Coxeter told Fuller how impressed he was with his formula -- on the cubic close-packing of balls.\"\n",
    "\n",
    "[2]  http://worldgame.blogspot.com/2006/02/octet-truss.html (additional info on the octet truss)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Pascal's Triangle\n",
    "=========\n",
    "\n",
    "Pascal's Triangle connects to the Binomial Theorem (originally proved by Sir Isaac Newton) and to numerous topics in probability theory.  The triangular and tetrahedral number sequences may be discovered lurking in its columns.\n",
    "\n",
    "pascal(), a generator, yields successive rows of Pascal's Triangle.  By prepending and appending a zero element and adding vertically, a next row is obtained.  For example, from [1] we get [0, 1] + [1, 0] = [1, 1].  From [1, 1] we get [0, 1, 1] + [1, 1, 0] = [1, 2, 1] and so on.\n",
    "\n",
    "Notice the triangular numbers (1, 3, 6, 10...) and tetrahedral number sequences (1, 4, 10, 20...) appear in the slanted columns. [3]\n",
    "\n",
    "![Image of PT](http://www.4dsolutions.net/ocn/graphics/randtrianim.gif)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "==================== Pascal's Triangle =====================\n",
      "\n",
      "                               1                            \n",
      "                             1    1                         \n",
      "                          1    2    1                       \n",
      "                        1    3    3    1                    \n",
      "                     1    4    6    4    1                  \n",
      "                   1    5   10   10    5    1               \n",
      "                1    6   15   20   15    6    1             \n",
      "              1    7   21   35   35   21    7    1          \n",
      "           1    8   28   56   70   56   28    8    1        \n",
      "         1    9   36   84  126  126   84   36    9    1     \n",
      "      1   10   45  120  210  252  210  120   45   10    1   \n"
     ]
    }
   ],
   "source": [
    "from itertools import islice\n",
    "\n",
    "def pascal():\n",
    "    row = [1]\n",
    "    while True:\n",
    "        yield row\n",
    "        row = [i+j for i,j in zip([0]+row, row+[0])]\n",
    "\n",
    "print(\"{0:=^60}\".format(\" Pascal's Triangle \"))\n",
    "print()\n",
    "for r in islice(pascal(),0,11):\n",
    "    print(\"{:^60}\".format(\"\".join(map(lambda n: \"{:>5}\".format(n), r))))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": true
   },
   "source": [
    "Each number in Pascal's Triangle may be understood as the number of unique pathways to that position, were falling balls introduced through the top and allowed to fall left or right to the next row down.  This apparatus is sometimes called a Galton Board.\n",
    "\n",
    "For example, a ball could reach the 6 in the middle of the 5th row going 1,1,2,3,6 in four ways (counting left and right mirrors), or 1,1,1,3,6 in two ways.  The likely pattern when many balls fall through this maze will be a bell curve, as shown in the simulation below."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/jpeg": "/9j/4AAQSkZJRgABAQAAAQABAAD/2wCEABALDA4MChAODQ4SERATGCgaGBYWGDEjJR0oOjM9PDkzODdASFxOQERXRTc4UG1RV19iZ2hnPk1xeXBkeFxlZ2MBERISGBUYLxoaL2NCOEJjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY//AABEIAWgB4AMBIgACEQEDEQH/xAAbAAEAAgMBAQAAAAAAAAAAAAAABAUBAgMGB//EAEwQAAEDAgQACQgGCQIEBgMAAAEAAgMEEQUSITETIjIzQVFxcrEGFCM0QmGBwRUkc5GhsiU1UmJjgsLR8EPhFkRkoiZTdIOS8ZOz0v/EABgBAQEBAQEAAAAAAAAAAAAAAAABAgME/8QAIBEBAQADAAMAAwEBAAAAAAAAAAECETESIUEyQlETA//aAAwDAQACEQMRAD8A+foiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIt4onSus21/et/NX5st23U2urXFF2fTPYwuJFgbLm+MscQdwmzVaoiKoIiICIiAiIgIuogJiD2vYTrdoOoXMC5Ava/SUGEXWaEwutnY8ftMNwViCHhnlvCMZpe7zYIOaLLmlpsRYhdRTk03DCSM62yX43bbqQcURdqmnNM8NL433F7xuuEHFF2pad1VMImPYwnpebBcnDK4g9CDCLu2le6jdUhzMjXZSL6/cuCAi71VLJSOY2QtJe0PGU30KxS00lXO2GK2Z21zYIOKLLmljy07g2K7NpZX0j6kAcExwaTfW6DgiLvU0k1KWCZuXO0OGt9EHBF1pqeWqmbDC3M92wvZc3NLHFrhYg2IQYRdW00zqZ1Q1hMTDYu6iuSAi61FNNTOa2ZhYXNzAHpCxBBLUyiOFhe89AQc0WSC0kEWI3C3FPKacziN3BA5S7oug5oi6TQS07g2aNzHEXAcOhBzRdIYZZ35IWOe7ezRdc7ICLpwMvA8NwbuDvbNbS65oCLeSKSF2WVjmOtezhZaICIiAiIgIsgXWQwluboQaotmsLtlvDTvmlEbbXPWi6rkimnC5wWi7OMbDX3X+S51FFLTszvLbXtoVPKL41GREsqyItgwlpd0BMhy5uhB3ouePdUlvPff4KNRc/8CpQ5/wDzqXPLrthxpPzDu0KNUC0rlKnH1d3aFGqOdcriZOCIi25MIsoiCwsoAToBdAU/C4qaolihkY4zOlFj7OXpv+KjU8zImStfC2TO2wJ3b7wuebK/NEXNtsb6or0VTiGBCVzfMuEseUxoAP4rtQ/QmIT8FBREODcxz/8A2vOumhdSsi4ENe25MgOru1WXkr+tj9kfELNnpqX2s8RiwfDXwieiu2S+rei1vf71EZWYBwoaaMtB9oi4HbqunllYmkIFuXp9yoKWSGKQmoh4VpaQBe1j1pJuLbqp+Kw0Ta+qaDwZAbwTWDinTpVTqNApEUMtdWNijJe99gC7qA/svRzNo/J2haWRtlrH8l7h/lgm9M62pH4XI6mhfTQ1D3kekBjIA7Fw8zmjkYKqKWGK4DnFh0C7vxzEXvzGpc33NAAVvTeUsbqhscrTwDgAXO3aff7k9np5yoZHHO9sLy+MHiuta4W9PHDK2YzT8G5rbsFr5j1L1ldglFM7zhsRuASWRm2f/deMdbMbbKy7LNFyAR1rtLDEylilZOHSPvmjtyUaabzN4cH+cZuKRybLgqjdgEkjGvflBIGY9AXSeIQ1RjgkEwBs1zRykqhTZ2ea58uUZs/7XSrryToWS1D6uUXbFozvdf8AnWpbpZN+maTAYKemNVi0uQAX4MG3w7fcFwONUUPEpsLiydcmpK28q6oyV7acHiRNuR7z/tZVLW03mb3Oe8VAcMrQOKQpPfurbr1HpIqzCaiCJ1XT00bpSRZoHF7dNExDybZPaSjkLTpxXm4t7ivKL2eAVkD4/NIpnycGwOBeLEdY+Gn3qWa4su/VeVnifR1j4WSZnMOXMxdKfDaidkri0x8G3NZ4ILvcFd4/Q04rYaqV7oo33bI5g1uBof8AOpb+StRLL50x8r5GMLS0uN97/wBld+tpr3p53gazg8gimydWU2WZqSSlhinc9l3k8S/GbbrCv8XxWalqahsdWWvY5uSHgwQRYXufvU3DaqHHKBwqYWFzTle22naOpN/Txm9PHt4WqmY1zySbNBcdAP7Lu6nrKSZzIQ8lp5yIEj4FZxekjw+vdBBK54A1uLWv0KV5NyyfSULDUOZFcjJckHTqV+Jr3pDhwyrnZLIIy3g25jnBBd2daw5tcYuDMcwjO7cpAK9V5RVs9DRMfTvyOdJlJtfSx61UsxushpOFfWRyS5rcE6MbddwpLatkiqkppqOOGd2S0oOUXuR2rVoqK+djS8vebMDnnbq1XqKDFKPGfQVUDBMRyXC4d2FVuPYL5g0z0jyac8pv7P8AsmzX8V7W1tLK+KEOzRuILohf8QsRYbWSRyTCFwEQDiHCxPYOlbYaHiS/nvmrdXZr7ke5eyxCaSHD5ZYjaRrbtJ60t0SbeKd59LGYeDm4Mm+QNNrrMtLU0LInSRgcO3ig6n/NVaOxuuipRMaqB0mfKYsmtutQsTqfpKsbOxwbaNuY3tlPSqnpxNNiFfUN4SOV0juKHPFlxloaqKR0boJLtNjZpIXq8AL2+cwmr85awtLXg33GoUetrq1mLTRNrIoImNzNDwLHQafipv2vj628q5rmGzmlp6iLLC9th9TDjlA7zmFhc05Xt+Y6l5LEKbzOump73DHaH3bhWVLNe0ZZWFlVllu6dBRu/wAE9korpD090+CkYcL1rfiuEI0d3SpGGeujsKzeVvH4tnDjxd535SoeL6Uw76ngXkj/AJz/ANqgYz6u3vfJcsex1y5VOsrCLu87sDamPvKSNLYWHrQj6uD71tPzEV77f2UUouf+BUoc/wD51KLResDsKl/8wsZddMeNZvVn9oUWp50/DwUuX1aT4KLVc7f3DwCYmSOiIujkIiIgutLUyUk3CwkB1rai65IitnOzG9gOxdo6t8dM6ns0xvdmII1UdEGSvQeTtU6oxJrXtYODgLQWi2lwvPK58lf1q77I+IUvFnVh5UTGnnoZWhri0vIDhcdC8zNJw0jn5Q0uNyBsvQ+WH/Kfz/JebspjxcuvTeR8elTKWtOzQbaje/yUPyreXYm1p2bGLD71L8kJhaogO+jx4H5KP5WwFtbFPbivZl+I/wDtT9j9VTNUNlgijELGGMWzNGru1coniOVjywPDSCWnY+5aItsPeYLVeeYcyTKGWJblGwsdF5jEXQ0eJ10RgZIH8kn2Li9x96vvJdpbhDSfae4jw+S8xi8wqMUqJG6tL7A9dtPksTrpeRDXeSaJ9LFG2BrHsvmeDq5cEW3NtE5rZWue3O0EEtvv7l7fAXRPw/hIYhEx73HKDe3R8l4Zer8kqkPpZacnjMdmHYVnLjeHVbi8sMeL1wmg4QuADDmtlOUaqmV95WUxZWMqQOLI2xPvH+yoVZxMupEslOaWFkcOWVt8783K+Cn+S78uMsHQ5rgfH5KoV15KxF+JmS2kcZ1950S8J1c+VeU4Ve1vSC34qH5G5bVl9TxPhuseV1SMsNKDrfO73dA+ax5Hb1f8n9Sx+rX7Kzyi/XdR/L+UK8wDzakwuSqcHRA6Pe86Ot1D4rhi2LRU+JyQy0EEwYRx3DjHQFTMVo24thcb6R3JGaNoNgfdZX4Tt08zUVUdZiT6mcERudfKN8vV9y74WYj5QQmna5sWfih24FlWOaWuLXAgg2IPQp2Bfrmm73yWrxiX29J5TmnGGDhw8ku4mXodY7rxrQC4Amw616zyu/V8X2o8CvJKY8az6lSvp4m076QyCdmr3Ha/RZe5jMdfh7TIOJNGMw7Qvni95Su80wWJ8unBwAn7lMlweHnEYneIs3Bhxy5t7L3uIQmowueKMF0j2WaPevn5Nzde8xZzmYRUOaS0iPQhMvhj9eYj8ncQc6z4g0deYH5qPitJDRSsgY9zpWtHC32DrdC4trqtjgW1MwI/fKV9T55WST2Iz20PZZa9s+novI4N4Kou7ctuOrdVHlAB9OVAvYXbr/KFa+R/N1Xa35rXFcWhp8UkhloIJgwi73DjEWBWf2a/V18nWRUGHVNXPIQxxFiRoQOr715yuqTWVstQRbO64HUOherxeibiuGRyUjrlozRtBsHDqt1rxpBaSHCxGhBVx/qZfxhZWFkLTDZm57CnsI3p7Fr0IrvDyH90/JScKF6w90qNDzb+78wpeDi9Y7ulYy5XTH4thzrO68+Crsa5mPvFWbR6UfZv8WqsxvRkQ95+S549jplyqhERd3ndj6qNen+62qeai0Oyw63mrdenZZq2lojB/ZUUofWR2FS/+YHaolD6y3sPgpf+uO0LGXXTDhL6tL8FEq+c+A8Apkg+rTfBRKw8do6Mo8AmK5IyIi6OIiIgItm2s64vpp7tVqgIsnUDTVYQFeeSbHHEZH24rYyCe0hUjtHEBTI8SrGRZI6mRrWtFgCpVi78r43OhppADZrnAntt/ZeXUz6Wr9vO5T2uuo0sj5ZDJI4ue7clJNFu3WgrJKGrZPHqW7jrHSF65/meP4eWsfruP2mO7F4kalbMkfE/NE9zHDpabFLNkukyrwaupXkOgc9vQ6MZgVmiwatq5Q3gXRMvq94sB/dZZjmJMbYVTj2gH5LlPildUjLLUyFvUDYH7k9npe4pisGH0QoKFwdIG5C4eyP7ryyFAkmkt2IiKoKRQVklDVMni3G46x1KOiD3LZaLHaF0Yde41b7TCvM1uA1tK85YjNH0OjF/w3VdG98bg5jnNcNi02Knx47iUYsKkkfvAFZ1Zxu2XrjDhddO8MZSy363Nyj7yvQwupfJyhLZHtkqn6lrek9HYFRzY1iMrcrql4B/ZAb4KASXuJcSSdSSrrfU3Jx0qqiSrqHzTG73m5XovI9hEdU8jRxaAey/915k6FTBileGi1XKB1A2Sz1ol97SPKWMsxmVxGjw0j7gPkuvk9iwo5PN53egedCfYP8AZQJMTrZYyySpkc1wsQSoia9aN+9x6zH8F86HndIAZbcZo9v3j3qlwKN4xunaWkFriSCNtCuFLXVcLcsdVKxg6A7QfBbHFq+5+tydt7KavC2deo8o6Oasw9rYGF72PDso3IsR815LzCsvbzSe/wBmVI+lcQjt9clJIB1N1n6cxI6edO+4f2SSxbZUzDcAkLxPiFoYGalrjq7+wW3lBjLKpvmlKbxA8d37XuHuVNPVVFQfTzSSd5xK55DkD+gmyuv6m/kGNL3ta0XLjYBe/wAQgdUYdNDHy3sIF+teFpqmemceAldGXaEhSJMTxCORzTWzEtJHKKlmzG6YOD4iHZfNJL9iziGHuoKam4VuWeTMXC97DS3zXSPFMSdHI8Vj7RgE37bKHV1M9TKHTyukcBYEq+z09F5INIp6l1tC4AfcqrykYWYzMSDZwaR9wC5vxHEI2s+uSAObmADlqK2urDwDql7g4G4cehTXvZv1pN8nsXFJJ5tUOtA88Un2D/ZWGPYL5zerpG3k3ewe17x715Yts0G/wUuGsraemDoqp7WZsoaHHT4Jr7CX1qobmlri1wII3BRdqmsqKq3nErpMu11xtZaRsz2uxYHJKy32uxY9koO8XNP7vzCmYL62/uqHFzTuweKm4L6xKf3VjLldMexcDnCeqI/maqjHNoe13yV1G2/CE9EX9bVS45/o/wA3yXPHsdM+VUoiDdd3nSHtcIGX2Oy3r25XsBvfLv8AFZl1ggNj1LGIm847qy00ofWm9h8FMI9O3tCh0PrTPj4Ka7n29oWcut4cZcL0tR3R4qFVgBzPexvgpxF6afu/NQqzeP7NvgmK5IwcRe3ShFgFhbvHFZ2fNdHFjhHcGI78UG9lqiIOrpDMCXBoyt9loHT7lo1xY7MLX94utomFzJSPZbf8QuaDdw4rXdLiUDy1rmgCx6whHomHrJWiDZ4s8jqK6F+ePKWtGRvQLX23XOTnHdqyzkv7vzCDDHZHh1gbG9nC4KPdncXWAv0AWC1WT0IN4pMgc3Ix2bS7he3YuR3WzeUO1YO6Dd8gdE1nBsbl9oDU9q0aQ1wJAcAdj0p0LCI2ecziQAPcFtG8MDrxtfcW1votCsIouj5GOiY0RNa4buBN3LmiIywgPBIzAHUdazI4OeS1oaOoFarO5RXSJzAHB8WckaHNay1s1u+pTkj3rVF46ulYWNAiAI3N90a1hcDa46W3tdWGCUwqOGBIGWx1+Kn4ph8ZY+obdpZHrfQlZ3qtSbm3npix0rjGzI3oaTe3xW0ZjbG7hGFxIs2zrZT1o8Bwu3dczyQqzxiy3kdEWNDGOa4col17/wCarRYKqOkWXUvBLekBayFpkcWAht9ATewWzOaf2hc0HdxjEZDw4vLRlINgNt/hdcWi5W83KHdHgsRC7/gfBBmXguLwYcNNcxW8OTIwzhxizG4YbHoXA7rsR9UYf4jvAINDbhfR3Db6X3XSbJmlz5uFzHbbdco9ZWdoW9T6xL3z4oN4M/m9QRbLlGb7wtKjgs44HPa2ue267U3qNX/J4qK7UoJNoxJEKkvyZBzYBNvitKQONSBHbNY79hW1ZyovsmeC2w31wdx/5SoOJ4LgN38Lm6hlt/lloDoASct1gp0Kjabgs/oc+Ww5e97a/issDMzuFLhobZR0rmtnco9qA32rdS3AZwDiS7PcWFtFiMcWTu/MIB6En3hB2i5l3w8VPwIXllPuUGIehP8AL4lT8BHGmPuC55crrj2LqLUTe6Jv5/8AZUOOnjQj3HxXoIRxKn7OP8zv7Lz+Pc9D3CfxKxh2N5/jVUg3RZG67vOmS8xTN1326Ohc6+/D6jXKPiukou2mt02+S519vOTbqGvXosxpig9bZ8fBTXc83tHioeH+ux/HwUx/PN7R4rOXW8OOjRelqfdH8wodeNYT1xM8FPj9Wqvsj4hQcQ/0PsWeCmK5IK6ScmPu/MrTpXWYejh7nzK6uTiiId0RKpB6CrP8L+pqjsYX7KdhwDqTEBbaC/8A3NUWB4bcFRTgiY9NwSuKl5w2O52zOUQ6lCtn8s9qzHyZO78wsSc47tW0fJk7vzCqOaydgsLJ5IQG8sdqw7crZnLb2rV3KKB7KwtvZK1RGSsLJ6FgICdCJ0IC6AZRcrRouVl5ubI1GL3OqzY2vY2612oI2y1kTHi7SdV6Y4XTGnEZbdodo/o7Fm3RMdq/yZBvUWy7DlfFW2Jj9G1B1Pot3bqBgkbY6uujZlLWusL/ABU/FTbDajfm7cbdYvXSceOYbOW0o0BGytqXDMlMJZY8zJm6HcjRccRoW0gjYC45wTxlvbGvSrQrNrOsetYO60y3ZzEnaPmua6N5l/aPmuaI6z8pvcb4JTC8h7jvylZqOWO43wCzR8677N/5SiuBUhw+oxn+I7wao5Ul4/R0J/iv8GoOVOL1EY63DxWan1iXvnxSk9bh748UqNZ5D+8fFB2p/wBX1Z97PEqKVNprfRNaf34x+ZQjug71fOMv/wCWz8oXXCxesPuik/I5YxNuSpYB/wCTF+RqzhXrTz/Bl/I5QQyuz43No4nnZz3AfCy4lTZ/1TRj+JJ/SqIXSt5BaR/eK1HKC6Tj00vfKDEfIl7vzC3c0Ckv1uHgVrFzUvYPELtKPqDT++PBRSHmvu+asMBHFmKgQDifd4FWeAECGYkarGXK6Y9i5hb6CoPXwQ/F6pcToKmuq42UsRkLYrmxtbjFXjNKeY9box+Dyt8IP6Ukb/0zT/3lc969ulm48y3yZxV3/Lgdr2/3XVnkriZ9iIdsgXuyD1rXW6z/ALZH+UeDxShnw+WjjqMpda4ym/SqyrN6h1vd4L0vlhfz+iLTxspt9681Ui0xHVYfgu2F3JXLOarfDxeti7fkpr+db2jxUTDPX4u1TJOdb8PFMurhx1jH1eq+yPiFAruTB9k35qdH6vVfZHxCgVvJg+yb81MWskNdpz6OHufMriuspBji7vzK6uLksndYREWOFH6tiQ/6b+tqrlOw8gQ1wPTT/wBTVBUVu7mmdp+S0XR3Ms7T8lzVRvLzz+1bRcmXufMLE/Pv7StoBdk3c+YQcVseS1aro8WiYeu6DRnLb2rB3K2Zy29q1dyigz7JWq2HIK1RGT0LAWXdHYgQYGhUumySNIcxuYe5RF3pbiUHo2KLGrrsuDylZYdhBqQ18nIe3QnQXUKtjIkzgHKfcvT4SL0FMLXNtncnZZt1G5N1CfR09FVUnBtsXHdwV1Yuu4G/G5XR9yq8TNsTw8BxsCeUNFag6k3vxxxvZ+5YrcVWGNLa3EScmryQT8VIxS7sMnyXk4oBLhYhR6QB1fXOLA+zrhwO2pUzE2l2HVJbxjlGu1kvScbUDC2hpr8Q5BqTfoVL5SDg5qbQt4p3PYrrDxkoqcniHgxq7XoVL5Ucum0I4pSdS/iquCdKWvaL666rg5pDrHdTMPOrgojrsk94K6Rzv9dWwSCF/F9odPaufAS/sLLSTC8k+0Pmup9RHb81UKiCQyCzPZb0+4LFC28sg6opD/2lazvLZgQdmt8AumHH0sx/gSflKCIVLkH6Lpz1yyeDVFKly/qqm+1k8GoONHrWQD+I3xWs/PP7xW9H67B9o3xXOXnHdpQS6f8AVFZ9rF/WoRU6AfoSqP8AHiH4PUEoJmKm9WPsovyNWcK9Yl91PL+QrGK+uH7OP8jVnCxeSo91PJ+VQ+oZU2qFsJoT1ukP4hQSrCsJOE4cOgCQ/wDcqIUYvKwdZC3n56bvn5rFP6xF3x4rMpu+bvlQIuYmPuHipEwthjPtB4LhF6rN2tUuoH6GY7+Nb8EWOdOPRn/OgqxwEfV5T71X0/NHs/pVlgHqsneXPLldcexdgWoCekztH/Y5R46+mw+t4WWojje6IMyva43be99AelSP+Qd7qhv/AOsrzHlD67H9k1Zk3dNZXUeoPlLRX9dgA90Mh/sub/KSkFy2sYdNB5s7/wDpeGWw2Wv8sXP/AEyX2N1UlZiVPmc1zALss3KR1319ypasZah496sqq3n1Nb9g779Kr64Wq5BYcpaxmkvt0woXxCHvKZLpK0jXbxUTCBevj7VMIvIL9AHipl1rDjI5qq1BHBHUX/aCgVnJg+yHzU/aGo+zPiFArOTB9mPmmK5Ii2cbtZ7h81qtjyWro4tURCg70zsrJx1x2/ELgrHCImyMr84vlpXuHuNwq9B1e21JE7rc75LipcwthlMf4kn9KiIV0qBaoeP3it6fkT/Z/MLWq0qZO8V0pBdlT9ifEKL9Rl2lH1aE9ebxXFS6htsPpHdefxVRGj5xvaFq7lFbR843tC1dyigyOQe1arYcg9q1QZcsLJQIjCtcFpxJUgv5LdT0qsY25V55OO9JPxi3QbC6zlxvHqbX0oGEy5QRI0Zn6aKVhAvh9N7XF2dtssVrQ2lnuC28Y0GoK3w9pFBTsAzcXZ22yx8dPqBiJ/SmH8Y9PKGitgeNe9+PyvZ+5VOIG+K4eMx0vyhoFb+1mv7fK9nbqUvFnVXRM4StriADZ24NutScXcW4dUOAzWDdTp0rjhoJra4hrXcbf710xvj4XUOtn0aL7W1V+pOJVFY0cDuQTHudehVeJxOlxWkiAyksduVY0RyUkB5Ho93a30UStA+n6KzLHgzue1J0vG8VFe0Msd+D2eNAfivL1TAHuLdg4he4ABuCBueL7IXipDmmla79ohXGpYjt5h/eHzXY+ojt+a5EFsbwf2h811PqA/zpXRycp+c+A8AulCbSS/ZPH4Fc6jnT2DwW1LypPs3eCDgVMm/VVL9pJ/SoamTn9GUg/ekPgg40htVwn99viub9XHtW9N6xF3x4rm7coJ8RtgNQOuoj/K9QFMYf0NKP+oZ+Vyh9KCbiwtXO+zj/ACBZwveqP/TPWMXN8Qf3GflC0oX5BU++Ej8QgiqbVPLsPoW/stf+YqEu8zwaWnaN2h1//kUGkHrEfeHikh1f73FYg5+PvDxWHHftKDrH6rN2t+am1BH0CwdPnB/KFCjP1SX3ub81IqHfoqEfxXeAUUg0id3T+VWeCObHQOc4gXd0qsi5l/dP5QrHC42SYaA8Xs+4XPLjrj1c5x5nJcgfWBv9mP7rzOPn9ID3Rt8F6FzBNRv4W7vrF7km9+Db/dedx5xfiTiTc5G3PwCmHVz/ABVq2bqtVszcLs4raqBbiVPc6hn91XVvrculuMVZTZn4pAL65PiN91WVnrcu3LO3asxal4KL1QNtnN8VKHOHu/MKNgnrX8zfFSmc87ufNZy66YcanmJ/s/mFArORB9mFOJ9DN3PmFBqweDhNtMgCmJkiLJ2Cwsrq4gQ7lYCyd0FtgQBp8UPVRu8QqhXOAgea4qenzN3iFTKRalTG+G0w6pJP6VFUiT1GDvv/AKVHVSutUb1Mp/eK6UnIqfsT4hcZuef2ldKY2ZP74/mFD64KfV/qjD//AHPzKAp9Yf0Rh/8A7n5lRCj5xvaFq7craPnG9oWp3KDI5B7QtVsOQVqiMndbQxulkDG8p2y1K6U5LZWubuNkWJDKKZ7H8G3MGuDS5Wnk410UtTG4kObYG2qnYfTCnp2jc5hd1+LsuGEn6/iBze37I06Vz3vbrrVTsW9HQ1BaC30Y0bqCs4eeEoKc8slvJdoBotMWu3DqkWLfRjRuoKzhgzYdSjl8TZ2g2Wfi/UOvd+lsP4+1+UNArbNY5ib8flezt1XVRXu/S2H8bYHlDQK2BBJJN+PyvZS8J1V4U8vqa8lgdZ5GZug3KkYy9pwud1uFtlF+jdRcH1qK8kF3pTq3QDU/gu+O64XOSc3J4zdt1fp8TaQhtLCbZPR9Jv0KBUgf8Q0lmW9Edz2qZTXbTREHg/R7u1uoM7f/ABDSjIR6L2j2qTpeLa9uLb2jxegdi8LO76zIQLcY+K9z1NA2J4vQOwrzNPhoraOaVrQHiU8a+tuqy1imStl5u9tyLrThzweTK3KratoBDRscWZb6HXdUpFjYrUrGUdqz1h3w8Fim3k+zclUb1D1rE7Ln97SFplzUmY/UqYd7xUZSJj9Vpx7neKDWkF6uIfvjxXI7rtQ+uw98eK4ndBLZ+qJPt2/lconSpjf1Q7/1A/KVD6UEzFdMQkHUG/lC4QcmbufMLriZvXy/DwC5Q8ibuf1BByXSTmouw+JXNbvPEjHUPmUCDn4+8FqfmtodJmH3rRB2b6pJ32+BXWd16GAfvuP4BcR6q7vjwKPdeBjeok+CipMR9DJ3D4BW2D6YePeVURaQv7rv6Vb4P6gL9a558dcerIH6m7/1B/IxeYxo3xF/db+UL0h0pD9u/wDIxeZxc/pGT4eAUw/Jc/xQltHy232utVvFzrO8F1cVy5tsXiF+TFrfcKoqfWJNAOMdu1Wlz9KgHW0QB9yqZedftyjspFqdgptVtHW5qkxn0rr/ALHzUHCjavi7ylB1pvh81nLreHGxNopu58woFUbiLuBT3c1N3fmFX1O0fcCYrk5MjfJfI0mwuVoiLo4txG4x8Jbi3tdaFEQWmHmWjhrA+I2npXAXNtLjXVVrGl7g1ouSukE3BtlG+dmX8QuKg6yFwhawi2VzvktWRPexz2sJa3cgbLLj6Fg95+S53VG8xvK82tc7LdjXxxvJabPZp94Wk/Pv7xW0TjkluTyPmFBza0vcGtBJJsAOlSap580poS1wdHmvf3m6igkG4NiukjiYowTe1/FUYiY9xzNaXBupsNlzO63je5psHEA72O60O6DcxvawEtIDttN0bFI4gCNxJ6gsFznNALiQNgTspWZzaMOBIPWO1Bq+ikDBo7P0stqF3oqdwqo7NLm5hmdl0CsOFEeGxOt6RzeVuXLanhdTz0bOTwly650Oqxt01r2u2t3P72/s7Kpwxro6yuL87Q9xy2bxXaq4eNd/a3HJ2VVhkjn1de0uccryBlFxZYjdd8WH1CoaAWEtFgzW664aWuw+njPHcxurHaW0XDEYvN8OqS3MCACMnRsu9CGnDKfhbP4lyBvsnxPqvrncJitEWvzBoJNxo1WxdrmLva39nZVVYXDFqJokaWvBuBoD2q1c0Egacr+Xb/PxSkVeE3bNW52vGeQlpaLNOvguuOXOHyMAJe4tADNiuWElzpKwvkcRHIQ1o1HZ2LvjLSyhkLfRSXa4ZToNetX6fEuDi08fsHg7HNrdV8zb4/A7gzlEdiSbjpU+mY7gY/SHhHt4xcNCOmygvyfTcTGl+Ux3MZvqpOl4siQBa21+L0D4qswBhFA8uZb0l79I2Vi9hPF4j22JEZ0y/FQcAL5KAySEkh9g8uvYXT4XrGNhzqWzG5yZBY9LvgvMTt1DxsV6nHOJTCRj3NfwoAkboVTS4dIIzdwIva46FrH0zl7V9QC573gHLfdc2dI6wpWIRy00zoHu0NnWB0UeBpc8hpsQLhbc71o4FpsRYrq4F8cbWgnK03+8rSWR8r80ri51gLn3CwXUyPijiMby0lhBsfeVRpTODKiNx2DgVoWkAEg2K7ULC+qZ7jcrk6R5YGFxLWm4F9Ag7seDh/BDlma9vgo7gQ4gixCkA5cPDho7ht/guEkj5ZC+Rxc525PSg7VnpaqWRgu0HdaQuDY5gd3MAH/yC6TzPZUThji0P4rgOkLnDbgp775Rb/5BQcy0ttcEXTlWAC2fNJIxjHvJazkg9Cwx7o3Nexxa4agjoVGGmzgepYWS4vkLnG5JuT1rZsr2xuiDiGO1I60DMOALenMCjgcjDbe6z/ofzfJHyvfHG1xuGDQdSK7jSN3vB+SuMJ9RZ2nxVN/pnsPyVzhPqDO0+K5Z8dcepzyTC77Z/wCVi81i36xl+HgvSHmnfav8GrzWKfrGbtUw6v8A0/FEW0XOt7QtV0gAdPGDsXBdXFYB/wCkhf8AYA7NlWv5x22/QrBrc2IuDTazRa/Rsq47lIVKwz16LvKQ4+mb2hRcP9ei7ykv59vaFjLrphx1cfRS935hQKnaPuBTnc3L3fmFBqf9PuBMVyR0RF0cREREZCwiIrJPFA6lhEQbzc8/vFIzxZO78wsS6yu7SjdndnzQarJOgWEQZbuFqsjdYRGehST6kP8AOlSsDibNJI11htqRtupmI0UUVK8gZLGzQOnXdZtamO5trh7fOJ6fW0cTRc2vrZTqhtsSorNDdDqTus4LCIqOOS2UvJOffo6lmpFsUojkDdDubrH11+LMg3Pf39nZVOEtLKyvuSLvJs0XG5Vqdz3/AOVVOGSu89ro23uZCARsN1Jxb2JmJs/R9QC0tBDebOh23TDS3zON7SXBzbPzbNFlriYMeEVHBjMW84QdDst6AsjwyncYXBrmcYB176J8T6hV5/T9Gwt4rQbEiwKtsozAfvbezt/n4qurw5+L0Ly5rgL2HUNVZGwtty/5dv8APxSrOqnBgc9eA5oHDOu0dPuC6Y+bYRK1oDRmHFdq5cMFDzJWubGxw4YkuvqNehd/KR36OkABsXN1Ivr2q/WfiTT3fTMc5mfQNawm5b71Cnc53lJCGOJLYrHTUbqya13AscQQ5rLsybnbdVhDG+UzHNBZeLM6xvmKTq3i0kBc4scCWalrhyr9N1DwUDzCwbYNfyr7Ka4HJxczS4OIc3Vx7VEwRlqEvt/qcq6nxb1wx8fV4nWv6ccf/ZTKuITRX0NmjVuw7VE8odYIXWveccfo+5WZGZhvY2aNW7DtT4fXj8Xa7hWvduBl+5RqTQvedmtVrjsbeFdGLX3FlUxaU8p67BdMeOWU9lVHlfnGzlrLzcXc+ZXWI8NTmM8puy0DMzowdg25+8rTKRh4DJWA8p2vwUJ4s4jqKk0smatDu233LhOLTO7UHU/q5o/inwCjLuT9SYP4jvALj0oOlRz8neK1a6zXDrHzWZjeZ/eK060GFk7DsWERBZWEQdD6uO8fBa9Sz/ojvH5LCKkf6R7D4hXGGeoR9p8VTO5r4HxVzhulFH8fFcs+O2HU1pvG77R/yXmsSN8Qn75Xo2n0bvtH+K83iBvXz98qYdP+nEddKe/nEdt8wXNdIeeZY21C7OSaHlte8iwGUb9gVcd1OZx6uToFulQSpCpFB63H2qTJz7e0KLRetx9qkv55vaFnLrphx2cfRy935hQaj2O4FMefRS935qFP7HcCmK5OCIi6OIiwsogiIgIiINn8t3asDY9iO5R7UGx7EVhFhZRALCysILnydL2yylgBNhe/UpOJF1RBJwZLmRtzEnqUDCXiKOoeegBWLojFg88sznCSVuZtj71zvXXHixwoZcPpnWycTlb9C41QtilEQ0N0OpK74ZduH0zuRdnKGp2UOvJbi2HgAMBvudCpOtfFv0k/v7+yqfCHv89r2NBDDI67wNBurkakkn2t/ZVRgmY1uIWB4LhHZi09qk5VvY741GPoiVgdlEevVn2UjDeLh9K7kej5W/R1KPjuZ+FzPa9122BFtMt1Iw7i0FMdWejHG36E+H1FmaG4tSNyBl2E6nQqyzN4bg78a5P7tlUYiP0xQts5t2EcU3urZshdxAfSDi5SLadaUnVVhMIHn13xgCa/VfXZd8caZMLkAGWPOBodP8uuGBNAhrM7n5mzGwte50XTyju7DHuOUcduh0P3K/U/VYMjcyJgjYQeD0IO6qS7/wARMdC1o9DqR81axhzGguaB6PQtO3aq8PI8pmt4RpHA9SkLxYTHIx5c7KGtuZW9JKi4AeEwtrwLcfe+2qkzsa6mcy4DSHE5uS5QvJwWwlhtbjnjX96fFvWvlCDwEBte9QON/srUcZpJNwGjUaAdqq8fc4Q07h0zjjf7K0vdpN81gNW7DtT4fXlcaeTjEmoIDRsoMrMlO62xddWeI0rqzGKngjd7GtIsNFxnpeCjLJZHHiEtI6wty6Ys2qYZODkDujpUmrLWMGXd417FDW7zcNv0Bbc29IbVDfj4JVi0xPWtIiQ+43sfBYMjzoXEjtRHQ+qs77vALkN10PqzO+7wC5jdFbSc47tWi2kN5HH3rCDCIiIIiINv9MdpWFk8gdqwEVIcfRDsPirjDvUo/j4qldzY7PmrrD/U41yz47YdSmniO+0f+YrWl8nI8SY+pdOWZpHC1r7FIzdh77/zFW+APBw94vqJn9PvXO2zenTUvVUfI4X9e0+y/wB1szyPAdfz/Ufwv916UORt9Vn/AEy/p/ni+ezt80xKoiLs2RxYXW3t7lAKtMYbfGq3W3HJVWV6cePPk7UfrUfapUnPN7QotH60ztUp/Pt7Qpl1rDjd/NSdnzUKfdndCmyc1J2fNQp92d0KYtZOJREXRxYWUWERlERAREQDusjY9iwiKwsoiI6UrBLURsdyXGxV99Ft4PO0iNjXiwKpcP8AXoe8F6atqY2B0mrnDkuI022WMq6YxT09M+bEZIWnM0uu732VxiloMMnbazJG8TNqehQcBBkrqgvYHEN6TbdTsVBjwuVrRwgLLkvNi3sWb1qcSsMs3DqU2ycTlDU7KBXtAxjDtGs0Ou6nYWcuH0ptl4m7dTsqyeJ0GL0GZ9731Juk6t4vXaEm9+PyvZ+5U2ElokxIOFy6Qt0dYK5vc5r+3yuj7lT4M3PJiF2scwyuu73KTi3sScas3DakBwZYN4rRodQu+G3ZQUxHo7xjjHXoUbFy1uF1gbI612+z2KXhoy0FMbZPRDU69CfD6gVseXHqAZXNuw9N+tW1nEkFxLC/+ZVVYD/xBQ5Wlp4M7u7VZAODsob6TMbydfuv/myUnVXgYBgqyRJfhzbW2twuvlGD9GSFziHZ26WvoueANvTVZcDfhXXdfbbVPKJubCWuYQQ2QZnX5XQr+yfqsomh0b7EPcGgknQhVrXud5UnPYWhs27baKyY9vAcKIw06NdffdQWEHyqvnNuA3cPepFqfKcsRNwLtdq7Vp7FX+TothLHZQOOeMN9+pWU2kLjmDeI7jO1B+CrvJ1p+iYzlHLPGG+6fC9aeUA9HSusNajln+ytr3BN81gNRoB2qpx4aUZsPWOWd/uVt1uJvbLxtrfD/N0+H1VUZzeUNccxPFGrB7ljFKa9KXZWtsM1gfvW1IXPx+v4J4zZW67dC74gQcNqS1vBgRm7XbnREjxcrMjyFh2w7FIlbwkLX9IC4P6OxdY5WMM3+B8FqtmbnsPgtVWXU+rs7x+S5jdbu5hnePyWg5QQZfyj2rHQsv5Z7VqgIiICIiDY8gLCz7IWEV1PIHZ81eYf6pF2KjdyB2fNXdB6rF2Lnnx2w67RniHvO/MVDOLx00j2xUz2u5L3NmIzG+9rKTEeJ/M78xVBU+syd4+KzjN2md1Fx/xFJ/5T/wD8zkb5QSOLWiN4vueHf/dUa3h51vaunhHPzqRmc+omdK8ucb3cTe6iFSM1p5b2F77BRyrErtR+tM7VJfz7e0KNR+tMUl/PjtCzet4cbyc1J2fNQpt290KY/m5Oz5qJNym90JiuTiUQotuTCIiIyiBEBERFEREBEVrT4cDSNmtnz6EdXYpbok2g0NvPIri/GV069bViNrbQRkF3UFBrKdlPVRCnsHdv3K6oKZtODEHZn3Bceg6bLGVdMZpHw1gOLV2rDbpOikYkDLhc5kHCFrRZ21lFw8skxSuL2O3GjehTMUucOnzHO7g926AKXq/HXDNMPpjbLxN2anZQ6/i4vh54rdDqNVKwy/mFPcGPiaObqTootfpi+HkBjDY6hJ1fi3ze1+9yuj7lUYOGufX5iSTMQLaAq0zWOYn2tX9B+H+bKrwNoM9dKQ14Eh0PSpOUvYlY7xcLqQDk5NwOnULrhk5GHQuaODyxi+bW4t0KHjJc3CJBEHWLvSB3RqFKo3k0NNdokdkFmW203T4fUGpGTyjpRdzmiPTOPcVciQZG8U2ueIRoFT1LnjHqS0mdwjOw7VaMkJtI8PY4EtcTsT2JVir8n3HJOWXeOEN2DYhb+URZHhgYGg55MwcBtqdFpgNuAneXkemPRytlnH3P+jCI2BkIlFxtbf5q/sz+q1a9vBniiM5QNRy1XNIj8qjo2O8PbdTg4mFvCAODhxTvkVZS8Tynm4+hjuDl7FItXEjssDnXy3Y7jHUH4Ku8nXZcMju3K3MfSdWqnSvyU7nE5fRu429/gq3AXOGGxFsgPHI4M9OqTi3rXGS1zqJ4z5nT3zkaHXdXlwWEnW2XjdX+fNUmNSAz0MWYh7ZgS0bDbZW5N7u3sW8Y7j4f5ul4TqsoSXY/iBa0S8ka6KRiRP0bUlr+EaYtzuFFonNGNYk9/pNh1dC74wHHDZuFdZwj4oaNwn1JxTU1Gx9LnddriwZBuHFQa6NjWMcxmUkkEX6rL1FE0Ow6kDgXgMByga7KlxqkbHKBmBL7kdq1L7Z16U8e57D4LUrZoLXEHexWhXRzdDzLB7z8lo3ljtWx5tvaVqNwgO5RWFkrCIIiICIiDPshOlPZCDdFdXckdiuqD1aLsVI7YdiuqE+gi+C558dsOusOrB2u/MVQVHrEneKvYDxG9p8SqGfn394ph2p/040W8NhK0na60WzOWF0cnUi8srgbWvsFHXYHjSXXFCu9H6yz4+CkO58doUej9Zb8fBSHesDtCxeumHGz+af2KLOOSbjYaKTJzb+xRZuUO6ExXJyRbNDSHZnWI2Ft1qtuQi2yjg75he+y1Q0IukojGXg35tNdLWK1YGl1nGw67INUWSNdNQt2tYYiS6z76NtuiuaLNtV1fG2zchuba6WsoacVdw1ohooY4m3ltuNVAbQyBpfK0sYBmuR0diusLgjZSMkc0MzG4kO5Wcq1jNK+SmdDUUbpXHNK65I16l6FgYA+LNms4cVzbE6XVViIy1NGRZjhySOnb7lZB+ZvCFpe5g4+vuWa3Ffh7smM1u7GX1AXfFPSYdNwjXMcxvFFtwo+FlhxOoGj3u3J0AUjEzahlMx0cMode9re5PqTjvhrcuH0xHo+JyhqdlExDTF6AgNZvruu2HuMdDTmK+Z3FDzta3Uo1W9hxanGjZI9BpcFPp8WxdY5i72tX+z9yq8Cbnmq5HXdGJSRl6+tWU7yIS5ou7Nyuj7lW4RJEGztu52c8Z2WzQ4lScW9dMeJGGva8jhCWlzgdCpdM4upqYCbLaPpGmwUTFiPM5WPktFnGYtbfqUilAfRwAkvayO5aBa2gV+H1Bkc0+UFMGNy+jN7ntVq0EF9iJHlxuy9x2qpkLn4zCZyRI0cRuTQjt+9WUxkMD3tY0OebNc3cBKRBwDKKSS7Tm4XTXbboW3lCC6gL3Buko33WMEMQgfHG67b3Mj25Tda4y5nmrGSHiPdmztbcg9SfU/VaxEvij4LKA9trDTRVUL2/wDEUhgLj6K2mtjorNoBja4gEtZoWdHaqoOyY1I9jrShou1rdCOnXr2SLVrK4Clc83jLojc7qFgV/ouJpaeDc43IGxupNRmbDM9hysdFoHa9GqjYa6OPDYREZHNLwGuOmut7/ip8PrnjDrVNBESLtl6uNuFbvOthyrjUjUKnxUtfWUscjiyUOu14F7no/GysgXtvZwM7GjM87OFtkvCdQKGRzMZr3Z7tLgHO+9dsXjvhUzbZeC2JNi4KPgwZw9UW2zknO145PVr09P3KRirwaCbW7LAOLtx2K/ScdqAn6PpzEHueyK9hbTQKrxiJkldSs4zMwcXG6taQg0ETYMpaRlaTxSRbQlQJnRPxqAch8bSACLgn5aJOpeIDYmQNfw0DJQ61nbEKqmaBM+wsLm3YvZTwR2keWts69z0FeWnfDcNuXE3zDLbKbq41MoiHkN7StQOMFMihiL2EuPB31Nr6LiYCJQG6gnQre2PFxRSOCZwDrkiYO0aGnUdq4W6wqmmtlldahsTXjgX52loJ0tY9IWYGxOeRM8MaW6OsTYoacLLNllblrOBDg/j3sW/NBodgsDdZ6At5Wsa4cG/OLC5tbXqQHckdiuKHmo/gqhws0e8BW9FpHH2Bc8+O2HXSDm2/HxKoZudf2lXsB9G3/Okqik513amHan/TkarLOWN/gtVszlD5Lo5Og1MhJXNgzSNHWQFtc8dbUgvVwjrkb4oVmk59vx8F2cQJwT7lEjeY3ZhutuGcX5rC6zZ7axykibM0hhv0j5qE83WfOH2Is2x30XPMUk0tylbsjL9it+APWFyDiEzlVNx04A33CzwJ6wuWcpnKezcdOCIvqE4Ow1cFzzlM3uCezcTYojUBsLGxtO+Y6LRtMPOxTue0caxcNQozZXNN26Jwrs2b2r3umjcXRwqmjjkc6pzZDYADdaYHwfnMjJA3XQOPQqs1UpBBO+6xHO+M3ZYFTVXyj0mJObHh01pWSknIOu3Ws4YYoqGnL5GvL+LkOzfevOPq5XtLXEEHoRtXK1oaCLDYKeN0eU2ucTEcmIUlPmADDYvbt0K1klhM8cQmAydNtCLLyDqmRzw91i4bFbmvnJvcX7E8avlF1QuhOJVxlls2+lhva6kYjPFLhUsjS0Etyhg8V5tlXKxznNsC7c9ay6slczISMvVZPFPKL7Cqlop2tly5coaGnX4rniMrW4tS8G4OEe7ujVUbKqRnJsPgjqqR7g51i4bFPH2eU09m+aESP9MORm069dPwVZg8kfA1ZklyuMmYC26o/P6j9ofctGVcrAQywB3U8avnHo8anifh1QBJd+ZoAA0Oy7NkpzSxFs8jXCLjW6dl5d1ZK9pa4gg7hZ8+nsBcWGyvjTyi8rKhgx2nfFISBHu7o3Vm+ogaZJGTEuyZgLbleOdVyueHusXDYrf6QnO5HVsp4nnF5gUsIoXGSQtfwug6Fvjs0UuHFzZSX8LyQNCF59lbKwEMygHoAWH1sr2ZXZS297EK+N3s8prT2Mc0ckMBMzGueLOA6FW0szHeUczi5sbQ3LcC4NulUXn89gLiw2WorJWyF4sHHpU8aXOPUy10L6eZrCWObGRffMuWDzxR4RCHS8YOPEI06V5w10xBBIsd0bWzMZkbYN6rJ4084vsVfG6soHtl4rpMzgBydlZRTQecebsmGRhzZnN3NtvxXjnVkry0uykt29y2+kJ73uLp4084u8OmYzEKyocQ28mjbaFScemiOHuhZJc6O4o39y8y2skaSWhoJ3WZK6aQWflcFfG7PKaeuo5IG0NO5z2vIi5J0tooVS6P6eo3CcFuQ3OXbQrz4xCYNDRlsBYCyw6tldIJDlzDQHqU8aeU09g2op3vfBwmVzS4l5GhVNhtJSyRzVNQWPMbzxCOUqr6QnBvdt+xc2VcjAQ0AA7jXVJjYXKVfYs2nZRsjp2sAccwA3b7vxUahoGy07pJKtrXNcBlsqqSrfJbOGm2yMqpGNytAA6hdXV0nlFhUQMiqnxMcHtaL59lzMMY9q1/eoLqhz3ZnNaT16oJyBYMb+KeJ5RNMMZGUSsA942WfM2acdmllC84d+y1Z86f1NTVPKJslDmPEcwaLQ4Y82tIz71F87f1NWfO5OpqayN4pJwuQgWfH96jVNM6meGvINxfRPPZepq0kndIbua0qzy+pbj8dHS3pmsc3bklWlFxmxsHKyg2VIHkW0GikCvlBByR3G3FUyx2uOUi1pbyNa1mrspdb3alUUnLd2ru2te23EjJGxIUcm5umM1UyyljC2YeOFqsg2N1thtuXFdsOF8Rph/Eb4qPc6+9bwSugnZKwAuYQ4X2QrmiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiD/2Q==\n",
      "text/html": [
       "\n",
       "        <iframe\n",
       "            width=\"400\"\n",
       "            height=\"300\"\n",
       "            src=\"https://www.youtube.com/embed/9xUBhhM4vbM\"\n",
       "            frameborder=\"0\"\n",
       "            allowfullscreen\n",
       "        ></iframe>\n",
       "        "
      ],
      "text/plain": [
       "<IPython.lib.display.YouTubeVideo at 0x112141860>"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from IPython.display import YouTubeVideo\n",
    "YouTubeVideo(\"9xUBhhM4vbM\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "[3] for more on both Pascal's Triangle and Cuboctahedral Numbers see:\n",
    "http://www.4dsolutions.net/ocn/numeracy0.html\n"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.9"
  },
  "widgets": {
   "state": {},
   "version": "1.1.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}