Probability_concepts.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Probability concepts using Python\n",
    "### Dr. Tirthajyoti Sarkar, Fremont, CA 94536\n",
    "\n",
    "---\n",
    "\n",
    "This notebook illustrates the concept of probability (frequntist definition) using simple scripts and functions."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Set theory basics\n",
    "\n",
    "Set theory is a branch of mathematical logic that studies sets, which informally are collections of objects. Although any type of object can be collected into a set, set theory is applied most often to objects that are relevant to mathematics. The language of set theory can be used in the definitions of nearly all mathematical objects.\n",
    "\n",
    "**Set theory is commonly employed as a foundational system for modern mathematics**.\n",
    "\n",
    "Python offers a **native data structure called set**, which can be used as a proxy for a mathematical set for almost all purposes."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "{1, 2}\n",
      "<class 'set'>\n"
     ]
    }
   ],
   "source": [
    "# Directly with curly braces\n",
    "Set1 = {1,2}\n",
    "print (Set1)\n",
    "print(type(Set1))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "{1, 2, 3, 4}\n"
     ]
    }
   ],
   "source": [
    "my_list=[1,2,3,4]\n",
    "my_set_from_list = set(my_list)\n",
    "print(my_set_from_list)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Membership testing with `in` and `not in`"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Here is my set: {1, 3, 5}\n",
      "1 is in the set: True\n",
      "2 is in the set: False\n",
      "4 is NOT in the set: True\n"
     ]
    }
   ],
   "source": [
    "my_set = set([1,3,5])\n",
    "print(\"Here is my set:\",my_set)\n",
    "print(\"1 is in the set:\",1 in my_set)\n",
    "print(\"2 is in the set:\",2 in my_set)\n",
    "print(\"4 is NOT in the set:\",4 not in my_set)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Set relations\n",
    "* **Subset**\n",
    "* **Superset**\n",
    "* **Disjoint**\n",
    "* **Universal set**\n",
    "* **Null set**"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "Univ = set([x for x in range(11)])\n",
    "Super = set([x for x in range(11) if x%2==0])\n",
    "disj = set([x for x in range(11) if x%2==1])\n",
    "Sub = set([4,6])\n",
    "Null = set([x for x in range(11) if x>10])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Universal set (all the positive integers up to 10): {0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10}\n",
      "All the even positive integers up to 10: {0, 2, 4, 6, 8, 10}\n",
      "All the odd positive integers up to 10: {1, 3, 5, 7, 9}\n",
      "Set of 2 elements, 4 and 6: {4, 6}\n",
      "A null set: set()\n"
     ]
    }
   ],
   "source": [
    "print(\"Universal set (all the positive integers up to 10):\",Univ)\n",
    "print(\"All the even positive integers up to 10:\",Super)\n",
    "print(\"All the odd positive integers up to 10:\",disj)\n",
    "print(\"Set of 2 elements, 4 and 6:\",Sub)\n",
    "print(\"A null set:\", Null)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Is \"Super\" a superset of \"Sub\"? True\n",
      "Is \"Super\" a subset of \"Univ\"? True\n",
      "Is \"Sub\" a superset of \"Super\"? False\n",
      "Is \"Super\" disjoint with \"disj\"? True\n"
     ]
    }
   ],
   "source": [
    "print('Is \"Super\" a superset of \"Sub\"?',Super.issuperset(Sub))\n",
    "print('Is \"Super\" a subset of \"Univ\"?',Super.issubset(Univ))\n",
    "print('Is \"Sub\" a superset of \"Super\"?',Sub.issuperset(Super))\n",
    "print('Is \"Super\" disjoint with \"disj\"?',Sub.isdisjoint(disj))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Set algebra/Operations\n",
    "* **Equality**\n",
    "* **Intersection**\n",
    "* **Union**\n",
    "* **Complement**\n",
    "* **Difference**\n",
    "* **Cartesian product**"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "S1 and S2 are equal because order or repetition of elements do not matter for sets\n",
      "S1==S2: True\n"
     ]
    }
   ],
   "source": [
    "S1 = {1,2}\n",
    "S2 = {2,2,1,1,2}\n",
    "print (\"S1 and S2 are equal because order or repetition of elements do not matter for sets\\nS1==S2:\", S1==S2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "S1 and S2 are NOT equal because at least one element is different\n",
      "S1==S2: False\n"
     ]
    }
   ],
   "source": [
    "S1 = {1,2,3,4,5,6}\n",
    "S2 = {1,2,3,4,0,6}\n",
    "print (\"S1 and S2 are NOT equal because at least one element is different\\nS1==S2:\", S1==S2)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In mathematics, the intersection A ∩ B of two sets A and B is the set that contains all elements of A that also belong to B (or equivalently, all elements of B that also belong to A), but no other elements. Formally,\n",
    "\n",
    "$$  {\\displaystyle A\\cap B=\\{x:x\\in A{\\text{ and }}x\\in B\\}.} $$\n",
    "\n",
    "![3 sets intersection](https://upload.wikimedia.org/wikipedia/commons/3/3e/Venn_0000_0001.svg)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "S1: {9, 3, 6}\n"
     ]
    }
   ],
   "source": [
    "# Define a set using list comprehension\n",
    "S1 = set([x for x in range(1,11) if x%3==0])\n",
    "print(\"S1:\", S1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "S2: {1, 2, 3, 4, 5, 6}\n"
     ]
    }
   ],
   "source": [
    "S2 = set([x for x in range(1,7)])\n",
    "print(\"S2:\", S2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Intersection of S1 and S2: {3, 6}\n",
      "Intersection of S1 and S2: {3, 6}\n"
     ]
    }
   ],
   "source": [
    "# Both intersection method or & can be used\n",
    "S_intersection = S1.intersection(S2)\n",
    "print(\"Intersection of S1 and S2:\", S_intersection)\n",
    "\n",
    "S_intersection = S1 & S2\n",
    "print(\"Intersection of S1 and S2:\", S_intersection)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "S3: {8, 9, 6, 7}\n",
      "Intersection of S1, S2, and S3: {6}\n"
     ]
    }
   ],
   "source": [
    "S3 = set([x for x in range(6,10)])\n",
    "print(\"S3:\", S3)\n",
    "S1_S2_S3 = S1.intersection(S2).intersection(S3)\n",
    "print(\"Intersection of S1, S2, and S3:\", S1_S2_S3)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In set theory, the union (denoted by ∪) of a collection of sets is the set of all elements in the collection. It is one of the fundamental operations through which sets can be combined and related to each other. Formally,\n",
    "\n",
    "$$ {A\\cup B=\\{x:x\\in A{\\text{ or }}x\\in B\\}} $$\n",
    "\n",
    "![Union of sets](https://upload.wikimedia.org/wikipedia/commons/e/ee/Venn_0111_1111.svg)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "S1: {9, 3, 6}\n",
      "S2: {1, 2, 3, 4}\n",
      "Union of S1 and S2: {1, 2, 3, 4, 6, 9}\n",
      "Union of S1 and S2: {1, 2, 3, 4, 6, 9}\n"
     ]
    }
   ],
   "source": [
    "# Both union method or | can be used\n",
    "S1 = set([x for x in range(1,11) if x%3==0])\n",
    "print(\"S1:\", S1)\n",
    "S2 = set([x for x in range(1,5)])\n",
    "print(\"S2:\", S2)\n",
    "\n",
    "S_union = S1.union(S2)\n",
    "print(\"Union of S1 and S2:\", S_union)\n",
    "S_union = S1 | S2\n",
    "print(\"Union of S1 and S2:\", S_union)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Set algebra laws\n",
    "\n",
    "**Commutative law:** \n",
    "$$ {\\displaystyle A\\cap B=B\\cap A} $$\n",
    "$$ {\\displaystyle A\\cup (B\\cup C)=(A\\cup B)\\cup C} $$\n",
    "\n",
    "**Associative law:**\n",
    "$$ {\\displaystyle (A\\cap B)\\cap C=A\\cap (B\\cap C)} $$\n",
    "$$ {\\displaystyle A\\cap (B\\cup C)=(A\\cap B)\\cup (A\\cap C)} $$\n",
    "\n",
    "**Distributive law:**\n",
    "$$ {\\displaystyle A\\cap (B\\cup C)=(A\\cap B)\\cup (A\\cap C)} $$\n",
    "$$ {\\displaystyle A\\cup (B\\cap C)=(A\\cup B)\\cap (A\\cup C)} $$ "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Complement\n",
    "\n",
    "If A is a set, then the absolute complement of A (or simply the complement of A) is the set of elements not in A. In other words, if U is the universe that contains all the elements under study, and there is no need to mention it because it is obvious and unique, then the absolute complement of A is the relative complement of A in U. Formally,\n",
    "\n",
    "$$ {\\displaystyle A^{\\complement }=\\{x\\in U\\mid x\\notin A\\}.} $$\n",
    "\n",
    "You can take the union of two sets and if that is equal to the universal set (in the context of your problem), then you have found the right complement"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "S is the set of even numbers between 0 and 20: {0, 2, 4, 6, 8, 10, 12, 14, 16, 18, 20}\n"
     ]
    }
   ],
   "source": [
    "S=set([x for x in range (21) if x%2==0])\n",
    "print (\"S is the set of even numbers between 0 and 20:\", S)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "S_complement is the set of odd numbers between 0 and 20: {1, 3, 5, 7, 9, 11, 13, 15, 17, 19}\n"
     ]
    }
   ],
   "source": [
    "S_complement = set([x for x in range (21) if x%2!=0])\n",
    "print (\"S_complement is the set of odd numbers between 0 and 20:\", S_complement)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Is the union of S and S_complement equal to all numbers between 0 and 20? True\n"
     ]
    }
   ],
   "source": [
    "print (\"Is the union of S and S_complement equal to all numbers between 0 and 20?\", \n",
    "       S.union(S_complement)==set([x for x in range (21)]))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "**De Morgan's laws**\n",
    "\n",
    "$$ {\\displaystyle \\left(A\\cup B\\right)^{\\complement }=A^{\\complement }\\cap B^{\\complement }.} $$\n",
    "\n",
    "$$ {\\displaystyle \\left(A\\cap B\\right)^{\\complement }=A^{\\complement }\\cup B^{\\complement }.} $$\n",
    "\n",
    "**Complement laws**\n",
    "\n",
    "$$ {\\displaystyle A\\cup A^{\\complement }=U.} $$\n",
    "\n",
    "$$ {\\displaystyle A\\cap A^{\\complement }=\\varnothing .} $$\n",
    "\n",
    "$$ {\\displaystyle \\varnothing ^{\\complement }=U.} $$\n",
    "\n",
    "$$ {\\displaystyle U^{\\complement }=\\varnothing .} $$\n",
    "\n",
    "$$ {\\displaystyle {\\text{If }}A\\subset B{\\text{, then }}B^{\\complement }\\subset A^{\\complement }.} $$"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Difference between sets\n",
    "\n",
    "If A and B are sets, then the relative complement of A in B, also termed the set-theoretic difference of B and A, is the **set of elements in B but not in A**.\n",
    "\n",
    "$$ {\\displaystyle B\\setminus A=\\{x\\in B\\mid x\\notin A\\}.} $$\n",
    "\n",
    "![Set difference](https://upload.wikimedia.org/wikipedia/commons/5/5a/Venn0010.svg)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Set S1: {0, 3, 6, 9, 12, 15, 18, 21, 24, 27, 30}\n"
     ]
    }
   ],
   "source": [
    "S1 = set([x for x in range(31) if x%3==0])\n",
    "print (\"Set S1:\", S1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Set S2: {0, 5, 10, 15, 20, 25, 30}\n"
     ]
    }
   ],
   "source": [
    "S2 = set([x for x in range(31) if x%5==0])\n",
    "print (\"Set S2:\", S2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Difference of S1 and S2 i.e. S2\\S1: {25, 10, 20, 5}\n",
      "Difference of S2 and S1 i.e. S1\\S2: {3, 6, 9, 12, 18, 21, 24, 27}\n"
     ]
    }
   ],
   "source": [
    "S_difference = S2-S1\n",
    "print(\"Difference of S1 and S2 i.e. S2\\S1:\", S_difference)\n",
    "\n",
    "S_difference = S1.difference(S2)\n",
    "print(\"Difference of S2 and S1 i.e. S1\\S2:\", S_difference)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "**Following identities can be obtained with algebraic manipulation: **\n",
    "\n",
    "$$ {\\displaystyle C\\setminus (A\\cap B)=(C\\setminus A)\\cup (C\\setminus B)} $$\n",
    "$$ {\\displaystyle C\\setminus (A\\cup B)=(C\\setminus A)\\cap (C\\setminus B)} $$\n",
    "$$ {\\displaystyle C\\setminus (B\\setminus A)=(C\\cap A)\\cup (C\\setminus B)} $$\n",
    "$$ {\\displaystyle C\\setminus (C\\setminus A)=(C\\cap A)} $$\n",
    "$$ {\\displaystyle (B\\setminus A)\\cap C=(B\\cap C)\\setminus A=B\\cap (C\\setminus A)} $$\n",
    "$$ {\\displaystyle (B\\setminus A)\\cup C=(B\\cup C)\\setminus (A\\setminus C)} $$      \n",
    "$$ {\\displaystyle A\\setminus A=\\emptyset} $$\n",
    "$$ {\\displaystyle \\emptyset \\setminus A=\\emptyset } $$\n",
    "$$ {\\displaystyle A\\setminus \\emptyset =A} $$\n",
    "$$ {\\displaystyle A\\setminus U=\\emptyset } $$"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Symmetric difference\n",
    "\n",
    "In set theory, the ***symmetric difference***, also known as the ***disjunctive union***, of two sets is the set of elements which are in either of the sets and not in their intersection.\n",
    "$$ {\\displaystyle A\\,\\triangle \\,B=\\{x:(x\\in A)\\oplus (x\\in B)\\}}$$ \n",
    "\n",
    "$$ {\\displaystyle A\\,\\triangle \\,B=(A\\smallsetminus B)\\cup (B\\smallsetminus A)} $$\n",
    "\n",
    "$${\\displaystyle A\\,\\triangle \\,B=(A\\cup B)\\smallsetminus (A\\cap B)} $$\n",
    "\n",
    "**Some properties,**\n",
    "\n",
    "$$ {\\displaystyle A\\,\\triangle \\,B=B\\,\\triangle \\,A,} $$\n",
    "$$ {\\displaystyle (A\\,\\triangle \\,B)\\,\\triangle \\,C=A\\,\\triangle \\,(B\\,\\triangle \\,C).} $$\n",
    "\n",
    "**The empty set is neutral, and every set is its own inverse:**\n",
    "\n",
    "$$ {\\displaystyle A\\,\\triangle \\,\\varnothing =A,} $$\n",
    "$$ {\\displaystyle A\\,\\triangle \\,A=\\varnothing .} $$"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "S1 {0, 3, 6, 9, 12, 15, 18, 21, 24, 27, 30}\n",
      "S2 {0, 5, 10, 15, 20, 25, 30}\n",
      "Symmetric difference {3, 5, 6, 9, 10, 12, 18, 20, 21, 24, 25, 27}\n",
      "Symmetric difference {3, 5, 6, 9, 10, 12, 18, 20, 21, 24, 25, 27}\n"
     ]
    }
   ],
   "source": [
    "print(\"S1\",S1)\n",
    "print(\"S2\",S2)\n",
    "print(\"Symmetric difference\", S1^S2)\n",
    "print(\"Symmetric difference\", S2.symmetric_difference(S1))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Cartesian product\n",
    "\n",
    "In set theory (and, usually, in other parts of mathematics), a Cartesian product is a mathematical operation that returns a set (or product set or simply product) from multiple sets. That is, for sets A and B, the Cartesian product A × B is the set of all ordered pairs (a, b) where a ∈ A and b ∈ B.\n",
    "\n",
    "$$ {\\displaystyle A\\times B=\\{\\,(a,b)\\mid a\\in A\\ {\\mbox{ and }}\\ b\\in B\\,\\}.} $$\n",
    "\n",
    "More generally, a Cartesian product of n sets, also known as an n-fold Cartesian product, can be represented by an array of n dimensions, where each element is an *n-tuple*. An ordered pair is a *2-tuple* or couple. The Cartesian product is named after [René Descartes](https://en.wikipedia.org/wiki/Ren%C3%A9_Descartes) whose formulation of analytic geometry gave rise to the concept."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [],
   "source": [
    "A = set(['a','b','c'])\n",
    "S = {1,2,3}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [],
   "source": [
    "def cartesian_product(S1,S2):\n",
    "    result = set()\n",
    "    for i in S1:\n",
    "        for j in S2:\n",
    "            result.add(tuple([i,j]))\n",
    "    return (result)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Cartesian product of A and S\n",
      "{'b', 'c', 'a'} X {1, 2, 3}:{('a', 1), ('a', 2), ('b', 1), ('a', 3), ('b', 2), ('b', 3), ('c', 1), ('c', 2), ('c', 3)}\n"
     ]
    }
   ],
   "source": [
    "C = cartesian_product(A,S)\n",
    "print(\"Cartesian product of A and S\\n{} X {}:{}\".format(A,S,C))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Length of the Cartesian product set: 9\n"
     ]
    }
   ],
   "source": [
    "print(\"Length of the Cartesian product set:\",len(C))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Note that because these are ordered pairs, **same element can be repeated inside the pair** i.e. even if two sets contain some identical elements, they can be paired up in the Cartesian product.\n",
    "\n",
    "Instead of writing functions ourselves, we could use the **`itertools`** library of Python. Remember to **turn the resulting product object** into a list for viewing and subsequent processing."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [],
   "source": [
    "from itertools import product as prod"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "A is set of all possible throws of a dice: {1, 2, 3, 4, 5, 6}\n",
      "B is set of all possible throws of a dice: {1, 2, 3, 4, 5, 6}\n",
      "\n",
      "Product of A and B is the all possible combinations of A and B thrown together:\n",
      " [(1, 1), (1, 2), (1, 3), (1, 4), (1, 5), (1, 6), (2, 1), (2, 2), (2, 3), (2, 4), (2, 5), (2, 6), (3, 1), (3, 2), (3, 3), (3, 4), (3, 5), (3, 6), (4, 1), (4, 2), (4, 3), (4, 4), (4, 5), (4, 6), (5, 1), (5, 2), (5, 3), (5, 4), (5, 5), (5, 6), (6, 1), (6, 2), (6, 3), (6, 4), (6, 5), (6, 6)]\n"
     ]
    }
   ],
   "source": [
    "A = set([x for x in range(1,7)])\n",
    "B = set([x for x in range(1,7)])\n",
    "p=list(prod(A,B))\n",
    "\n",
    "print(\"A is set of all possible throws of a dice:\",A)\n",
    "print(\"B is set of all possible throws of a dice:\",B)\n",
    "print (\"\\nProduct of A and B is the all possible combinations of A and B thrown together:\\n\",p)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Cartesian Power\n",
    "\n",
    "The Cartesian square (or binary Cartesian product) of a set X is the Cartesian product $X^2 = X × X$. An example is the 2-dimensional plane $R^2 = R × R$ where _R_ is the set of real numbers: $R^2$ is the set of all points (_x_,_y_) where _x_ and _y_ are real numbers (see the [Cartesian coordinate system](https://en.wikipedia.org/wiki/Cartesian_coordinate_system)).\n",
    "\n",
    "The cartesian power of a set X can be defined as:\n",
    "\n",
    "${\\displaystyle X^{n}=\\underbrace {X\\times X\\times \\cdots \\times X} _{n}=\\{(x_{1},\\ldots ,x_{n})\\ |\\ x_{i}\\in X{\\text{ for all }}i=1,\\ldots ,n\\}.} $\n",
    "\n",
    "The [cardinality of a set](https://en.wikipedia.org/wiki/Cardinality) is the number of elements of the set. Cardinality of a Cartesian power set is $|S|^{n}$ where |S| is the cardinality of the set _S_ and _n_ is the power.\n",
    "\n",
    "__We can easily use itertools again for calculating Cartesian power__. The `repeat` parameter is used as power."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Cartesian power 2 with length 4: [('Tail', 'Tail'), ('Tail', 'Head'), ('Head', 'Tail'), ('Head', 'Head')]\n",
      "\n",
      "Cartesian power 3 with length 8: [('Tail', 'Tail', 'Tail'), ('Tail', 'Tail', 'Head'), ('Tail', 'Head', 'Tail'), ('Tail', 'Head', 'Head'), ('Head', 'Tail', 'Tail'), ('Head', 'Tail', 'Head'), ('Head', 'Head', 'Tail'), ('Head', 'Head', 'Head')]\n"
     ]
    }
   ],
   "source": [
    "A = {'Head','Tail'} # 2 element set\n",
    "p2=list(prod(A,repeat=2)) # Power set of power 2\n",
    "print(\"Cartesian power 2 with length {}: {}\".format(len(p2),p2))\n",
    "print()\n",
    "p3=list(prod(A,repeat=3)) # Power set of power 3\n",
    "print(\"Cartesian power 3 with length {}: {}\".format(len(p3),p3))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "---"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Permutations\n",
    "\n",
    "In mathematics, the notion of permutation relates to the **act of arranging all the members of a set into some sequence or order**, or if the set is already ordered, rearranging (reordering) its elements, a process called __permuting__. The study of permutations of finite sets is a topic in the field of [combinatorics](https://en.wikipedia.org/wiki/Combinatorics). \n",
    "\n",
    "We find the number of $k$-permutations of $A$, first by determining the set of permutations and then by calculating $\\frac{|A|!}{(|A|-k)!}$. We first consider the special case of $k=|A|$, which is equivalent to finding the number of ways of ordering the elements of $A$. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [],
   "source": [
    "import itertools"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [],
   "source": [
    "A = {'Red','Green','Blue'}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Permutations of {'Blue', 'Red', 'Green'}\n",
      "--------------------------------------------------\n",
      "('Red', 'Green', 'Blue')\n",
      "('Green', 'Red', 'Blue')\n",
      "('Blue', 'Red', 'Green')\n",
      "('Green', 'Blue', 'Red')\n",
      "('Red', 'Blue', 'Green')\n",
      "('Blue', 'Green', 'Red')\n",
      "--------------------------------------------------\n",
      "Number of permutations:  6\n"
     ]
    }
   ],
   "source": [
    "# Find all permutations of A\n",
    "permute_all = set(itertools.permutations(A))\n",
    "print(\"Permutations of {}\".format(A))\n",
    "print(\"-\"*50)\n",
    "for i in permute_all:\n",
    "    print(i)\n",
    "print(\"-\"*50)\n",
    "print;print (\"Number of permutations: \", len(permute_all))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.patches as mpatches\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {},
   "outputs": [],
   "source": [
    "from math import factorial"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Factorial of 3: 6\n"
     ]
    }
   ],
   "source": [
    "print(\"Factorial of 3:\", factorial(3))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Selecting _k_ items out of a set containing _n_ items and permuting "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {},
   "outputs": [],
   "source": [
    "A = {'Red','Green','Blue','Violet'}\n",
    "k=2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "2-permutations of {'Blue', 'Red', 'Green', 'Violet'}:  \n",
      "--------------------------------------------------\n",
      "('Blue', 'Red')\n",
      "('Blue', 'Green')\n",
      "('Blue', 'Violet')\n",
      "('Red', 'Blue')\n",
      "('Red', 'Green')\n",
      "('Red', 'Violet')\n",
      "('Green', 'Blue')\n",
      "('Green', 'Red')\n",
      "('Green', 'Violet')\n",
      "('Violet', 'Blue')\n",
      "('Violet', 'Red')\n",
      "('Violet', 'Green')\n",
      "--------------------------------------------------\n",
      "Size of the permutation set = 4!/(4-2)! = 12\n"
     ]
    }
   ],
   "source": [
    "n = len(A)\n",
    "permute_k = list(itertools.permutations(A, k))\n",
    "print(\"{}-permutations of {}:  \".format(k,A))\n",
    "print(\"-\"*50)\n",
    "for i in permute_k:\n",
    "    print(i)\n",
    "print(\"-\"*50)\n",
    "print (\"Size of the permutation set = {}!/({}-{})! = {}\".format(n,n,k, len(permute_k)))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "12.0"
      ]
     },
     "execution_count": 36,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "factorial(4)/(factorial(4-2))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Combinations\n",
    "\n",
    "Combinatorics is an area of mathematics primarily concerned with counting, both as a means and an end in obtaining results, and certain properties of finite structures. It is closely related to many other areas of mathematics and has many applications ranging from logic to statistical physics, from evolutionary biology to computer science, etc.\n",
    "\n",
    "Combinatorics is well known for the breadth of the problems it tackles. Combinatorial problems arise in many areas of pure mathematics, notably in algebra, [probability theory](https://en.wikipedia.org/wiki/Probability_theory), [topology](https://en.wikipedia.org/wiki/Topology), and geometry, as well as in its many application areas. Many combinatorial questions have historically been considered in isolation, giving an _ad hoc_ solution to a problem arising in some mathematical context. In the later twentieth century, however, powerful and general theoretical methods were developed, making combinatorics into an independent branch of mathematics in its own right. One of the oldest and most accessible parts of combinatorics is [graph theory](https://en.wikipedia.org/wiki/Graph_theory), which by itself has numerous natural connections to other areas. Combinatorics is used frequently in computer science to obtain formulas and estimates in the [analysis of algorithms](https://en.wikipedia.org/wiki/Analysis_of_algorithms).\n",
    "\n",
    "We find the number of $k$-combinations of $A$, first by determining the set of combinations and then by simply calculating:\n",
    "$$\\frac{|A|!}{k!\\times(|A|-k)!}$$\n",
    "\n",
    "**In combination, order matters, unlike permutations.**"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "2-combinations of {'Blue', 'Red', 'Green', 'Violet'}:  \n",
      "('Blue', 'Red')\n",
      "('Blue', 'Green')\n",
      "('Blue', 'Violet')\n",
      "('Red', 'Green')\n",
      "('Red', 'Violet')\n",
      "('Green', 'Violet')\n",
      "Number of combinations = 4!/(2!(4-2)!) = 6\n"
     ]
    }
   ],
   "source": [
    "# Print all the k-combinations of A\n",
    "choose_k = list(itertools.combinations(A,k))\n",
    "print(\"%i-combinations of %s:  \" %(k,A))\n",
    "for i in choose_k:\n",
    "    print(i)\n",
    "print;print(\"Number of combinations = %i!/(%i!(%i-%i)!) = %i\" %(n,k,n,k,len(choose_k)  ))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Putting it all together - some probability calculation examples"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Problem 1: Two dice\n",
    "Two dice are rolled together. What is the probability of getting a total which is a multiple of 3?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {},
   "outputs": [],
   "source": [
    "n_dice = 2\n",
    "dice_faces = {1,2,3,4,5,6}\n",
    "\n",
    "# Construct the event space i.e. set of ALL POSSIBLE events\n",
    "event_space = set(prod(dice_faces,repeat=n_dice))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(1, 3), (6, 6), (5, 6), (2, 1), (6, 2), (1, 6), (5, 1), (2, 5), (1, 2), (3, 3), (5, 5), (4, 4), (6, 3), (1, 5), (3, 6), (2, 2), (4, 1), (1, 1), (6, 4), (3, 2), (2, 6), (5, 4), (4, 5), (5, 2), (1, 4), (2, 3), (4, 2), (6, 5), (3, 5), (5, 3), (4, 6), (6, 1), (3, 1), (4, 3), (3, 4), (2, 4), "
     ]
    }
   ],
   "source": [
    "for outcome in event_space:\n",
    "    print(outcome,end=', ')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "metadata": {},
   "outputs": [],
   "source": [
    "# What is the set we are interested in?\n",
    "favorable_outcome = []\n",
    "for outcome in event_space:\n",
    "    x,y = outcome\n",
    "    if (x+y)%3==0:\n",
    "        favorable_outcome.append(outcome)\n",
    "favorable_outcome = set(favorable_outcome)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(1, 2), (5, 4), (3, 3), (6, 6), (4, 5), (2, 1), (6, 3), (1, 5), (3, 6), (4, 2), (5, 1), (2, 4), "
     ]
    }
   ],
   "source": [
    "for f_outcome in favorable_outcome:\n",
    "    print(f_outcome,end=', ')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {},
   "outputs": [],
   "source": [
    "prob = len(favorable_outcome)/len(event_space)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 43,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "The probability of getting a sum which is a multiple of 3 is:  0.3333333333333333\n"
     ]
    }
   ],
   "source": [
    "print(\"The probability of getting a sum which is a multiple of 3 is: \", prob)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Problem 2: Five dice!\n",
    "Five dice are rolled together. What is the probability of getting a total which is a multiple of 5 but not a multiple of 3?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 44,
   "metadata": {},
   "outputs": [],
   "source": [
    "n_dice = 5\n",
    "dice_faces = {1,2,3,4,5,6}\n",
    "\n",
    "# Construct the event space i.e. set of ALL POSSIBLE events\n",
    "event_space = set(prod(dice_faces,repeat=n_dice))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 45,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "7776"
      ]
     },
     "execution_count": 45,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "6**5"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 46,
   "metadata": {},
   "outputs": [],
   "source": [
    "# What is the set we are interested in?\n",
    "favorable_outcome = []\n",
    "for outcome in event_space:\n",
    "    d1,d2,d3,d4,d5 = outcome\n",
    "    if (d1+d2+d3+d4+d5)%5==0 and (d1+d2+d3+d4+d5)%3!=0 :\n",
    "        favorable_outcome.append(outcome)\n",
    "favorable_outcome = set(favorable_outcome)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 47,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "The probability of getting a sum, which is a multiple of 5 but not a multiple of 3, is:  0.11625514403292181\n"
     ]
    }
   ],
   "source": [
    "prob = len(favorable_outcome)/len(event_space)\n",
    "print(\"The probability of getting a sum, which is a multiple of 5 but not a multiple of 3, is: \", prob)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Problem 2 solved using set difference"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 48,
   "metadata": {},
   "outputs": [],
   "source": [
    "multiple_of_5 = []\n",
    "multiple_of_3 = []\n",
    "\n",
    "for outcome in event_space:\n",
    "    d1,d2,d3,d4,d5 = outcome\n",
    "    if (d1+d2+d3+d4+d5)%5==0:\n",
    "        multiple_of_5.append(outcome)\n",
    "    if (d1+d2+d3+d4+d5)%3==0:\n",
    "        multiple_of_3.append(outcome)\n",
    "\n",
    "favorable_outcome = set(multiple_of_5).difference(set(multiple_of_3))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 49,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(3, 1, 1, 1, 4), SUM: 10\n",
      "(4, 1, 6, 3, 6), SUM: 20\n",
      "(6, 2, 2, 4, 6), SUM: 20\n",
      "(1, 4, 6, 6, 3), SUM: 20\n",
      "(5, 6, 4, 5, 5), SUM: 25\n"
     ]
    }
   ],
   "source": [
    "for i in list(favorable_outcome)[:5]:\n",
    "    a1,a2,a3,a4,a5=i\n",
    "    print(\"{}, SUM: {}\".format(i,a1+a2+a3+a4+a5))    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 50,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "The probability of getting a sum, which is a multiple of 5 but not a multiple of 3, is:  0.11625514403292181\n"
     ]
    }
   ],
   "source": [
    "prob = len(favorable_outcome)/len(event_space)\n",
    "print(\"The probability of getting a sum, which is a multiple of 5 but not a multiple of 3, is: \", prob)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Computing _pi_ ($\\pi$) with a random dart throwing game and using probability concept\n",
    "The number $\\pi$ is a mathematical constant. Originally defined as the ratio of a circle's circumference to its diameter, it now has various equivalent definitions and appears in many formulas in all areas of mathematics and physics. It is approximately equal to 3.14159. It has been represented by the Greek letter __\"$\\pi$\"__ since the mid-18th century, though it is also sometimes spelled out as \"pi\". It is also called Archimedes' constant.\n",
    "\n",
    "Being an irrational number, $\\pi$ cannot be expressed as a common fraction (equivalently, its decimal representation never ends and never settles into a permanently repeating pattern).\n",
    "\n",
    "### What is the logic behind computing $\\pi$ by throwing dart randomly?\n",
    "Imagine a square dartboard.\n",
    "\n",
    "Then, the dartboard with a circle drawn inside it touching all its sides.\n",
    "\n",
    "And then, you throw darts at it. Randomly. That means some fall inside the circle, some outside. But assume that no dart falls outside the board.\n",
    "\n",
    "![boards](https://raw.githubusercontent.com/tirthajyoti/Stats-Maths-with-Python/master/images/boards.png)\n",
    "\n",
    "At the end of your dart throwing session, \n",
    "- You count the fraction of darts that fell inside the circle of the total number of darts thrown. \n",
    "- Multiply that number by 4.\n",
    "- The resulting number should be pi. Or, a close approximation if you had thrown a lot of darts.\n",
    "\n",
    "The idea is extremely simple. If you throw a large number of darts, then the **probability of a dart falling inside the circle is just the ratio of the area of the circle to that of the area of the square board**. With the help of basic mathematics, you can show that this ratio turns out to be $\\frac{\\pi}{4}$. So, to get $\\pi$, you just multiply that number by 4.\n",
    "\n",
    "The key here is to simulate the throwing of a lot of darts so as to make the fraction equal to the probability, an assertion valid only in the limit of a large number of trials of this random event. This comes from the [law of large number](https://en.wikipedia.org/wiki/Law_of_large_numbers) or the [frequentist definition of probability](https://en.wikipedia.org/wiki/Frequentist_probability).\n",
    "\n",
    "See also the concept of [Buffon's Needle](https://en.wikipedia.org/wiki/Buffon%27s_needle_problem)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 51,
   "metadata": {},
   "outputs": [],
   "source": [
    "from math import pi,sqrt\n",
    "import random\n",
    "import matplotlib.pyplot as plt\n",
    "import numpy as np"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Center point and the side of the square"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 52,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Center point\n",
    "x,y = 0,0\n",
    "# Side of the square\n",
    "a = 2"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Function to simulate a random throw of a dart aiming at the square"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 53,
   "metadata": {},
   "outputs": [],
   "source": [
    "def throw_dart():\n",
    "    \"\"\"\n",
    "    Simulates the randon throw of a dirt. It can land anywhere in the square (uniformly randomly)\n",
    "    \"\"\"\n",
    "    # Random final landing position of the dirt between -a/2 and +a/2 around the center point\n",
    "    position_x = x+a/2*(-1+2*random.random())\n",
    "    position_y = y+a/2*(-1+2*random.random())\n",
    "    \n",
    "    return (position_x,position_y)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 54,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(0.22132208501841943, 0.7844172623010937)"
      ]
     },
     "execution_count": 54,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "throw_dart()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Function to determine if the dart landed inside the circle"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 55,
   "metadata": {},
   "outputs": [],
   "source": [
    "def is_within_circle(x,y):\n",
    "    \"\"\"\n",
    "    Given the landing coordinate of a dirt, determines if it fell inside the circle\n",
    "    \"\"\"\n",
    "    # Side of the square\n",
    "    a = 2\n",
    "    \n",
    "    distance_from_center = sqrt(x**2+y**2)\n",
    "    \n",
    "    if distance_from_center < a/2:\n",
    "        return True\n",
    "    else:\n",
    "        return False"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 56,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "False"
      ]
     },
     "execution_count": 56,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "is_within_circle(1.9,1.9)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 57,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "True"
      ]
     },
     "execution_count": 57,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "is_within_circle(0.2,-0.6)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Now, throw a few darts"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 58,
   "metadata": {},
   "outputs": [],
   "source": [
    "n_throws = 10\n",
    "count_inside_circle=0"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 59,
   "metadata": {},
   "outputs": [],
   "source": [
    "for i in range(n_throws):\n",
    "    r1,r2=throw_dart()\n",
    "    if is_within_circle(r1,r2):\n",
    "        count_inside_circle+=1"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Compute the ratio of `count_inside_circle` and `n_throws`"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 60,
   "metadata": {},
   "outputs": [],
   "source": [
    "ratio = count_inside_circle/n_throws"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Is it approximately equal to $\\frac{\\pi}{4}$?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 61,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "2.8\n"
     ]
    }
   ],
   "source": [
    "print(4*ratio)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Not exactly. Let's try with a lot more darts!"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 62,
   "metadata": {},
   "outputs": [],
   "source": [
    "n_throws = 10000\n",
    "count_inside_circle=0\n",
    "\n",
    "for i in range(n_throws):\n",
    "    r1,r2=throw_dart()\n",
    "    if is_within_circle(r1,r2):\n",
    "        count_inside_circle+=1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 63,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "3.136\n"
     ]
    }
   ],
   "source": [
    "ratio = count_inside_circle/n_throws\n",
    "print(4*ratio)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Let's functionalize this process and run a number of times"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 64,
   "metadata": {},
   "outputs": [],
   "source": [
    "def compute_pi_throwing_dart(n_throws):\n",
    "    \"\"\"\n",
    "    Computes pi by throwing a bunch of darts at the square\n",
    "    \"\"\"\n",
    "    n_throws = n_throws\n",
    "    count_inside_circle=0\n",
    "    for i in range(n_throws):\n",
    "        r1,r2=throw_dart()\n",
    "        if is_within_circle(r1,r2):\n",
    "            count_inside_circle+=1\n",
    "            \n",
    "    result = 4*(count_inside_circle/n_throws)\n",
    "    \n",
    "    return result"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Now let us run this experiment a few times and see what happens."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 65,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Computed value of pi by throwing 3 darts is: 4.0\n",
      "Computed value of pi by throwing 10 darts is: 2.4\n",
      "Computed value of pi by throwing 31 darts is: 3.225806451612903\n",
      "Computed value of pi by throwing 100 darts is: 3.44\n",
      "Computed value of pi by throwing 316 darts is: 3.367088607594937\n",
      "Computed value of pi by throwing 1000 darts is: 3.152\n",
      "Computed value of pi by throwing 3162 darts is: 3.135989879822897\n",
      "Computed value of pi by throwing 10000 darts is: 3.1004\n",
      "Computed value of pi by throwing 31622 darts is: 3.141483777117197\n",
      "Computed value of pi by throwing 100000 darts is: 3.13996\n",
      "Computed value of pi by throwing 316227 darts is: 3.145550506439994\n",
      "Computed value of pi by throwing 1000000 darts is: 3.142276\n",
      "Computed value of pi by throwing 3162277 darts is: 3.140461129749228\n",
      "Computed value of pi by throwing 10000000 darts is: 3.1421584\n"
     ]
    }
   ],
   "source": [
    "n_exp=[]\n",
    "pi_exp=[]\n",
    "n = [int(10**(0.5*i)) for i in range(1,15)]\n",
    "for i in n:\n",
    "    p = compute_pi_throwing_dart(i)\n",
    "    pi_exp.append(p)\n",
    "    n_exp.append(i)\n",
    "    print(\"Computed value of pi by throwing {} darts is: {}\".format(i,p))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 66,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAf0AAAF1CAYAAADiLzM9AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOydaZhUxdWA3zMDMzCgbMOibKMiKhE/ECSAqCyKyBKjgBoJioiogMYlrrigBDERlaCiIigqGCC4IJsLmxoFCSgxiqIEYUBQQGVnYJbz/ajbTU9P99A909M9y3mf5z63b62n6t6+51bVqSpRVQzDMAzDKP8kJVoAwzAMwzDigyl9wzAMw6ggmNI3DMMwjAqCKX3DMAzDqCCY0jcMwzCMCoIpfcMwDMOoIJjSNwyjUERko4hsTLQcxUVElolI3OYoi8goEVER6RxFnAwvztSSk8yoyJjSNyoUInKqiDwlIl+KyG4ROSwiW0VkvohcKyJVEi1jvIm3MqzoeEp9WaLlMComlRItgGHECxF5AHgQ97G7AngZ2AfUBzoDk4EbgbYJErG00i3RAsSIq4C0OOb3NDADyIxjnoZRKKb0jQqBiNwLPARsBvqr6qchwvQGbo+3bKUdVf1fomWIBaoaV+WrqjuBnfHM0zCOhnXvG+UeEckARgHZQM9QCh9AVecBPULEv0xEPvSGAw6KyH9F5B4RSQ0RdqN3VBeRJ0VksxdnjYj83gtTSUTuFZHvRCRLRP4nIiNCpNXZ6woeJSIdRGSRJ8NeEXlXRAr0SIjIVC9ORmHp+erF69Y/z7vWgGNZcJmC0hrkhRskIl28IYK9IrLHGyo5LVQdi0hzEXldRH4Vkf0i8omI9ApML1S8EOn4x8tF5GoR+dyr5+0i8qKINAgRJ+JhDBFZLiKHRKRakPuHXr5TgtxbeO6vhJLRux4UkP95QfU9KoQMGSIyQ0R2es/JKu/D1DCKjLX0jYrANUBlYIaqfllYQFU9FHgtIo8A9+BabK/hhgMuAh4BLhSRC1Q1OyiZysD7QG1gDpAC/AF4XUS6A8OA3wILgUNAf+ApEdmhqjNDiPVbT4ZFwDNAM+BS4FwR6a6qH0VUCwXZhev9GAQ09X772BhhGr2Bi3FleQ5oAfQEzhKRFl5rF3D2FMDHuHqZD3wBnAi8CSwoYhluBboDM4F3gE64+91ZRH6rqjuKmO5ioD1wjpcuIpKGuxdQcMija0C8cKzB1fGDwCZgaoDfsqCwTYGVwAbgVVydXQ7MEZHzVXVp5EUxjABU1Q47yvWBexErMCTKeB28eJlAgwD3SsBcz+/eoDgbPfe5QGqA+zme+y/Av4GaAX4nAoeBz4PS6uzFUWBEkN/Fnvt3QFKA+1TPPSNEeXzpjQpyX+ZeBWHrYSOwMchtkJdWDtAtyG+s53dnmPtwY5D7RQHlHBThvRnlhT8MtA7ye9LzmxJNOYPCdvXSeCzA7ULP7T3vfFKA35ueW+MQMnYOSluBZWHyzQioiweD/Hz5L0jk/8mOsn1Y975RETjOO2+JMt5g7/wXVf3R56iqObix/zxgSJi4t2hAr4G61vj3QC3gLlXdFeC3AdcCbikiySHSWg9MDHRQ1TnAB7hW/zlRliuWzFDV4NbtJO/czucgIo1xinQ98HxgYFVdiOvFKAqvqurnQW6jgN3AlaGGYCLkEyCL/C36briPnAcDrhGRJNwH1XequrmI+QWzCfhLoIOqvov7AG0XMoZhRIApfaMiIN452mlpZ3rnJcEeqvot7iPiBBGpGeS9S0Mbv231zqtD+P0AJAMFxqKBj1Q1L4T7Mu/cOoRfvFgVws2n+GoFuLXyzsvDlOVfRcz/g2AHVd2N60qvAoS0LTgaqpqFU/ytRKSO59wV+LeqLgd+4sgHwZlATUI8J8VgjarmhnDfTP56NYyoMKVvVAR8yrZRlPFqeOdtYfy3BYXzsTtM+BzwK6WQfjh7gGB+CpOer/chOP94sivYwesJAfcR48MnY7iyhHM/GiVZN4txH4xdRKQW7uPK16uxBOgqIsIR5V/YeH60FKhXjxzsvW0UA3t4jIqArxUZ7Xxzn3IO1fqGI8MG4ZR8rKgfxt0nV2D+vlZ0KCPd4B6JeLLHO4crSzj3oxFN3USLr+V+PtAF974MVPrpwP/hnisFzLjOKPWY0jcqAi/hpuv1FZEWhQUMGgP2jRV3DhGuGa7n4PvA8fkSopM3bhyMT67AMe1fvXPjEOHDLTqUCxDGniBW+GTsEKYsnYqY7nnBDiJSAzeckAV8XcR0wRlc7sEp9a7AQWC55+dT/j2Bs4EvNGCmwlHII38viGHEDVP6RrlHVTfijLtSgPmh5rcDiEgP3NQzHy965/tEpG5AuGRgHO7/k2++dglxMm6anx8RuRin8NYDgVP2Vnrn64LCtwT+FCb9n71zk2JLGgbPwG0ZzvDw+iDZeuBa00VhoIgE2zSMwnXr/0ODpmBGgzem/iFO5v7Av3zpqer3uFkNf8Kt8hfNeP7PhP4oM4wSx+bpGxUCVX1ERCrhLK//LSKf4IzQfMvwnotTrqsC4nwiIn8D7gS+FJHZwH7cFLPTccMGj8VB/HeAx0XkIuA/HJmnnwVcG2QYNwc3je8PItII+BSnzC/2/C4Lkf5inFJ7Q0QW4Fq0m1T11RiXYzhulsJEEenJkXn6fT3ZLubI8ESkLAQ+FpFZOBuLTt6xEbg7BjIvxq1FUI+CY/aLgWsDfkeT5hUiMhdn1JkDfKiqHxZTVsM4KtbSNyoMqvowTlk/jWsJXgPcAfQC/oebftcpKM5duIV1vsOt3X4z7n9zH3CBqh6Og+if4rryU4ERuI+OJcC5wYrCszrvBszClXUETrFeCTwbJv3JuLn1NXAfOKM5osxihqquxa198CZumuEtuHnpl3DE7mJPyMjheRLXC9LKS+9U3FoFHVV1e7GFzq/Mg1vzPr8cXI9ApPwJ+Adu6t39uPruWmgMw4gRomqbaxlGacRbvnUp8JCqjkqsNCWLiEzHfZicqqrrIgg/Ctdr00VVl5WsdIZRfrCWvmEYcUFEksKsid8Nt8Ts2kgUvmEYRcfG9A3DiBcpwGYRWQp8g+sW/w1wAW453eEJlM0wKgSm9A3DiBfZuE15uuI2rknDbWT0T+DREMvpGoYRY2xM3zAMwzAqCDambxiGYRgVBFP6hmEYhlFBMKWfAEQkQ0RURKYmWpbSiogM8upoUKJlKeuIyEYR2ZhoOWKJiLQVkfdFZKf3nKxJtEzFRURGeWXpnGhZYkF5fO7KA2bIZxhGmUJEjgXm47bOfRVnDPhjoZGMmCMiy4DzVFWOFtYoPZjSTww/4Pb5Lund2coybwIrCL+trVFxaYdbFnekqj6SaGEMoyxhSj8BqGo2bp6yEQZvz3n7KDJCcbx33ppQKQyjDGJj+gkg3Ji+iEz13DNE5HoR+a+IZInITyIyydsyNFR6jURkgoh854X/RURWisj9QeE2esexIvKE9zvbW9LUF6aSiAwTkRUiskdEDojI5yIyItSWqN7Y++siskFEDnpxPhaRP4aR9USvLOu98L945XxOROoEpVtgTD+gDGki8piIZIrIIS+9u0SkQFejOP4kImu9+vlBRJ4WkRrRjjt6Mi0TkXSvHNu8/L8SkWvC1E9Y2wRfekFu/rFdEfmDiKz27sNW776leuG6erLsEZFfReTVwDoMkVcNr9w/ePWwVkRuDlVnXvjfishsEflRRA6LyGYReV5Ejg8Rdpknc4qIPCAi67x6mVpohR6J301E3vGehywR+VZEHg185r3/hQIve04veXke1fbDq0v16radiMz38lIRyfDCdPHu6VqvTg+KyJci8qCIVAmRZuB96ifuP3fAS3eGiDQMI0sbr6x7vXwWiUiH4tZPQFjfvajs3Yv/eXG+EZHrAsLdIO6/d1BEtojIQxJ62+Pg9H334TzvWgOOZSHCR/RflYD3oog0F5GZIrJdRPIkwM5BRE4WkVe85/iw9794RURODkrvei+94B0nB3vuByT/Vtp49zBLRKoGuP1ORBbLkf/6VhH5QETy7XxZVrCWfunkb8CFwFzgPaALbqvUZgRtzCFum9h3gdq4TT/ewC160gK3xejooLRTcBuH1PbS3gN876VV2cvzQmAd8BpuJ7cuwFO4BVUGBqX3LLDWy3sbUAe3x/irInKKqvo/PETkONwe5ccCC4DXceOyJ3jpPs2RbV4Lo7In+/G4XdZygN8Dj3rpPRQU/hngRlzLcBJu9bff4bqJK+MWjYmGmrjd4g4Ds708+wEvikieqr5cWOQouAm3uc5buG1puwO3ArVFZA4wAze2PQnoCPwRSPfiBJMCLPJkn+Fd9wX+DpxC0Gp44j5gXgAOAW8Dm3G7EA4B+ohIe1XNDJHP68BZuPvyFnDUTW9E5Hrcc7Qft1DPdtwGQ3d5eZ2tqruAXbh724ojuwb6DPgiNeTrANyD2+DnRVx9+TZNugu3Yc8nHLEZOBv3P+osIud72+0GMwz3PL0NfID7n1wO/J+ItArc3ldEOuLuQwruv7reK88ywmzPG0X9BDPDk2UB7hnvB0wSkWzgDOBqYB5u46DfAQ8AB4C/hpIjAN99GAQ0Jf//bWNQ2Gj/qwAn4TaZ+haYDlTF24hJRM7C1d8xuPpei7tnA4CLRaSbqvp2yvRtiNQN9yz78L1Dq+Keh2Ve2jWAM4GPVPWg5zYUeB5nMzIXZz9SD1d/1wATQ9ZQaUZV7YjzgdtZTIGpQe5TPfdMoEmAeyWcUlWgXYB7Ck5hK3BliHwaB11v9MIuAqqFCD/K838KSA5wT8btG6/AxUFxTgqRTgruD5cNNAxwv8lL408h4lQDqgZcD/LCDgpThgVB4evhXka7gMoB7ud44dcBNYNk9NXpxijunXrH5KA6aoF7oa0NCh+yHEHpLQtzH3YDpwW4pwJfAbm4j6PzAvySgPe9eK3C1Nm/gNQA99q43QUVt2Ofz705ThGuD7x/nl9XL/83g9yXeel8AaRHUZ9NcR8We3Cb7QT6TfTSnBRNnYbJp3PAvbs+TJgT8RYsC3If7cW7PMx92gO0DPJ7zfO7LMBNcMN6of5HfwqQr3Mx68d3L/4d9Myf6N3XX3HvjcD/Zk2cQtsBVIqwTpcBWoi/77mL9L+aEVAHj4RIT4CvPf8BQX6Xe+7fAEkB7ptwH0kS4LYV937KBUYHuF/spXF/gNtqr/7rhZAn4ue8NB3WvV86eVgDWlGqmgO85F22CwjXB/dHeVtVXwtORFU3h0n/dlXdH+jgdeuNwH3R3qoBLRrv9+14f7agPP4XIt/DuNZ1JdxXdjAHQ8TZr97XdYTcHBhe3Taqc3Dbw54SEO5q7zxGA1pDnoz3RJFfIAeA24LqaC2u9X+aiBxTxHSDmaCqXwfkcQiYiVPw81X1gwC/PGCad/l/YdK7RwNanar6C0d6ggKHJm7EtdD+pKo/BCagqktwLaw+Ycp5v6rujKRwHn/EfYA9rarBdi4jgb3AwOBu2GKwRlWfD+WhqhvUe5sHMd47XxgmzQmq+t8gN1/LMvD/2hH3bH6oqnOCwj+N+wALpjj1c3fQM78B9+FXE6fsfgjw24VryaYDIYclikGk/1UfPxG6B6AjrlW/XFWnB3qo6kxc2U4h//bYS4C6QEsAEWkBHIfrofuM/O8n3+/A7ZTBfcwX6A2M8jkvNVj3fulkVQg3nwKvFeDW3jsvjCLtLFxrLJjmuK7574D7JPQw70HcrAM/ItIE183YDWiC6zILJPAF8jbwCPCMiFyIG5b4GNc6jmY96N2quj6Ee6g6au2d/0VBVuD+0NHynaqG2vfdl39N3Mu4uIR6DnzGa6tD+Ple4o1C+OXguq2DWeadWwe4+caXz/O6U4Oph+v9aR5CjpUhwhfGmd65QNe2qv4qIp8D5+Je9v+JMu1QhJVPRKrhWtyX4Mp2DK516SOcMoz0/+or6wdBYVHVXBH5F65rO5Di1E9xnp9NIfyLQjT/VR//Cfw4DSBsXQS4d8I9yx8GuA3CvZ++4EjX/mJcg+k2ETlGVfd6fvvI/4xMBx4HvhKRmbh797Gq7ggjQ6nHlH7pJNT4nE85JQe41fTOPxA528MoWJ8B2Mm4fcrDUd33Q0ROxP1BagEf4cbuduO6zTJwrWx/C0RVN4lIO1y3aA/gUs9rs4iMU9UJEZYhVP1A6DryGTr9FBzYe9FGYkNQnPyLQ6jZCzkR+FUO4bdTQ49H++a3BxqE+Z6FO44iX/UQbtHOl/flG25qps+9Zhj/aAkpn2fPsgTXMv8S16OygyMtvAcJeJaDiPT/GvZZLES2ItePuhkw4eSK9vkpKkX5r4R7hopSF4Hj+k965y2q+q2ILAbuxH3c/hu34+MCr2cVAFV9QkR24uw2bgZuAVREPgDu0CP2A2UGU/plG98fKpruuHAtat9L4E1VvTRMmGBuwymIa1R1aqCHiPyBI13rRzJ33dWXi0glXDf0+bix/r+LyH5VnRJh3pHia5HXBzYEyZjsyR/NR1O05HnnAv81EYmVIouEdBFJDqH4ffvbByoB3+8aYXo0whJlj01gXg1w9grBHBcUrriEk+9inMJ/WVUHBXp4BqiFfQhHiq8M9cP4NwjhFu/6KQ0c7R0Vqp4gRF2o6lYRWYdT7Kk42w7f0Mq/cDYO5+OMiyF0j8orwCve/7UjridoMPCuiJzmDVeUGWxMv2yzwjuHstaOlm9wHxHtvVZPJDTzzq+H8DuvsIiqmqOqq1X1r8AfPOffR5hvNPi2a+0Uwq89Jf/h+6t3bhzCr20J5x1IJdwLK5jO3jlwW1vfc3VOSQoUlG/nYA/vJdsKNyT1dbB/jCnysxwFn4VLz/sADfWMlpb6CUWuJ0eseraORti6CHL/LMh9Ma5X6kZcL8BiAFU9gHvWu5G/2z8kqrpLVReo6nU4o+vaxOc/ElNM6Zdt5uIsZH/ntazzIWHmCYfC69J6Cve1PCFwnmpAesd5hjA+NnrnzkHhLsRN7QqO305EQrVyfG4HIpU3Cl7xziMl/5zvFJx9QUmzCtfav1JE0gLyr42bmhlPxgYafHky3OddvhQQ7mlct/aTItI8OBFxc/Fj9bKb5uV1k4g0C/IbjWuBTQszxhtLNnrnzoGO3hDW0aawRconuFkk54rIxUF+Iyg4ng+lp35C4RsaaxKn/D7G1V8nEekX6OFdn4ub5hdsv+Nrvd8TdO37fTpuyuLPBNlFiEgPr1cymHreuSTeWSWKde+XYVT1sIj0x42lv+bN512Bm/96Gu4LNpp7PBrX5X4Dzjp7Ca7rux5urP9snMXwWi/8RJzV9z9F5HUv7Om48fpZuGk0gVwJDPfGw9bjWsEn4WYhHOKIlXTMUNUPRGQSMBRnjPM67iXaB9cNuJUjXfAxR1W3ich03DoEa0RkPu5F3RNnbNS6sPgxZBtuTPpLEXkbN27bD/eRN1FVfYZPqOo3IjIYN4/9KxF5B/cyrYx7wZ+DG+8+tbhCqepGEbkFN9vjMxGZ5aV9Hs6g8BucoWhJMxf3TN4mIi1xrcomQG/cnP1iKzZVVRG5Fje18nUR8c3T9w1zvYP77wTGKS31E4rFQH/gDRFZgDP03aSqr5ZEZl79XY2rv5neWhXf4Cz2f48znr3Km8kSyFLcf7we8E3QjJTFOBujusDsEMNTM4Asz8hyI86w8xzcWhSrcdOfyxSm9Ms4qrpKRFoBd+O6+TviHv71RDkOqarZIvJ73DShQbgXXnXcS+Z74H6cNasv/Bci0gX4C06JVcJ9KV+KGyoIVvr/wCmejjhL3Kq4D4UZwOOq+mU08kbBjbiXw/W4D5qfcWv73wtsIfRUqVhyHc546w+4RXAygQnAY8BlJZy3D9/Y5SPAFbipWRtwi6Q8FRxYVaeJyH9wUzW74BYG2o/7SJqNM3SLCao6UUTWA3/GLRiUhrPufgw3XzucMVjMUNX9ItIVVx+dcS/2DbgP4Sco+CwXNZ+PvV6SMRwZlvvUy/NCgpS+Fyfh9ROGybh1BK7AGcRVwlm3l4jSB1DVT70ZJffhnuc+uPUF/oGbhrguRJxfxO3CeCYFx+w/xT3X1UL4gXuvXujF7YkbStmE+9B6Vt2S6mUKid7uxjDKB+KW7fwWmKGqBYZHDMMwyhs2pm+Ue0SkgQStKe6Nr/uGE96Mv1SGYRjxx7r3jYrALcAfxG0Gsg035acbbhGShbj1zA3DMMo9pvSNisD7OGOp7rhpNjm4bv0JwPgizC03DMMok9iYvmEYhmFUEMp1Sz89PV0zMjJimub+/fupVq1aTNMsjVg5yxdWzvKFlbN8Eetyrl69eqeq1g3lV66VfkZGBqtWxXZp5GXLltG5c+eYplkasXKWL6yc5QsrZ/ki1uUUkbAbJpn1vmEYhmFUEEzpG4ZhGEYFwZS+YRiGYVQQTOkbhmEYRgXBlL5hGIZhVBBM6RuGYRhGBcGUvmEYhmFUEBKq9EXkXhFREXn6KOFaisgHInJQRH4QkQdEROIlp2EYhmGUBxKm9EWkPW6f8S+OEu5Y3NrpPwFnATcDdwC3lbSMiWL69OlkZGSQlJRERkYG06dPP3okwzAMwzgKCVH6IlIDmA5cC/x6lOADgDTgalX9UlVfB/4K3FYeW/vTp09nyJAhbNq0CVVl06ZNDB061BS/YRiGUWwS1dKfBMxW1SURhO0AfKSqBwPc3gWOBzJKQLaEMnLkSLKysvK5HThwgJEjRyZIIsMwDKO8EPdd9kTkOuAGoIOqHvb2OP9SVUeECf8esEVVBwe4NQE2AR1VdXlQ+KHAUID69eu3mTFjRkzl37dvH9WrV49pmoF07dqVUPdERFiyJJJvpNhQ0uUsLVg5yxdWzvKFlbNodOnSZbWqtg3lF9cNd0TkFOAR4BxVPRxF1GAtKGHcUdVJuJ4E2rZtq7HerKGkN4Bo0qQJmzYV3CuhSZMmcd14wja6KF9YOcsXVs7yRTzLGe/u/Q5AOvCliOSISA5wHjDMu04NEedHoEGQWz3v/FPJiZoYxowZQ1paWj631NRUxowZkyCJDMMwjPJCvJX+W0BLoFXAsQqY4f0O1fpfDpwjIlUC3C4AtgIbS1LYRDBgwAAmTZqUT/FfffXVDBgwIIFSGYZhGOWBuCp9Vd3lWeD7D2A/8It3rSIyVkQWB0R7DTgATBWR00XkUuBu4AmNt0FCnBgwYAA33XST/7phw4YJlMYwDMMoL8R1TD9CjgNO8l2o6m4RuQB4Btcr8CvwOPBEYsSLDxkZGf7fGzduTJgchmEYRvkh4UpfVTsHXQ8KEea/wLlxEqlU0LRpU//vUIZ9hmEYhhEttvZ+KSVQ6VtL3zAMw4gFpvRLKYFKf/PmzeTl5SVQGsMwDKM8YEq/lFKtWjXS09MByM7OZtu2bQmWyDAMwyjrmNIvxZgxn2EYhhFLTOmXYsyYzzAMw4glpvRLMWbMZxiGYcQSU/qlmMDufWvpG4ZhGMXFlH4pxrr3DcMwjFhiSr8UY4Z8hmEYRiwxpV+KCW7pl9OtBgzDMIw4YUq/FFOjRg1q1qwJQFZWFtu3b0+wRIZhGEZZxpR+KcfG9Q3DMIxYYUq/lGNK3zAMw4gVpvRLOWbMZxiGYcQKU/qlHGvpG4ZhGLHClH4px1r6hmEYRqwwpV/KsZa+YRiGEStM6ZdybK6+YRiGEStM6Zdy6tSpQ7Vq1QDYu3cvv/76a4IlMgzDMMoqpvRLOSJiXfyGYRhGTIir0heR4SLyhYjs8Y7lItLrKHEu9MLtFZGdIjJHRJrHS+bSgBnzGYZhGLEg3i39LcBdwJlAW2AJ8JaInBEqsIicAMwBPgJaA+cDVYEFcZG2lGAtfcMwDCMWVIpnZqo6J8hppIjcCHQAvggRpQ1QGbhHVXMBRGQssERE0lV1Z4kKXEoIVPrW0jcMwzCKSsLG9EUkWUSuAKoDn4QJtgrIBoZ44Y8Brgb+XVEUPuTv3reWvmEYhlFUJN5TwESkJbAcqALsAwao6vxCwp8D/BNIx32kfA5cpKoht5wTkaHAUID69eu3mTFjRkzl37dvH9WrV49pmkdj7dq1DB8+HICTTz6ZSZMmlXieiShnIrByli+snOULK2fR6NKly2pVbRvSU1XjegApQDPcmP5YYCdwepiwDYBvgb/hxvTPBZZ5R9LR8mrTpo3GmqVLl8Y8zaOxbds2BRTQWrVqxSXPRJQzEVg5yxdWzvKFlbNoAKs0jF6M65g+gKoeBtZ7l6tE5CzgVuDaEMGHA/tV9U6fg4j8EdgMdAT+VcLilgrq1atHamoqhw4d4tdff2XPnj0ce+yxiRbLMAzDKGOUhnn6SUBqGL80IDfIzXddGmSPC0lJSTRp0sR/beP6hmEYRlGI9zz9R0XkHBHJEJGWniV+Z2C65z9WRBYHRJkPnCkiD4rIySJyJvASrqW/Op6yJxoz5jMMwzCKS7xbyw2AacA6YDFwFs4ob6Hnfxxwki+wqi4BrgQuxhnwvYuz5u+hqvvjKHfCsbn6hmEYRnGJ9zz9QdH6q+oMILYm+GUQW5XPMAzDKC4VZly8rGMtfcMwDKO4mNIvI9iqfIZhGEZxMaVfRjBDPsMwDKO4mNIvIxx//PFUquRMMLZv387BgwcTLJFhGIZR1jClX0ZITk6mcePG/mtr7RuGYRjRYkq/DGHGfIZhGEZxMKVfhrBpe4ZhGEZxCDtPX0Rm4fax/5/3uzBUVS+PrWhGMNbSNwzDMIpDYYvz1AUqe7/r4XZ5MxKIKX3DMAyjOIRV+qraJeB357hIYxSKde8bhmEYxcHG9MsQ1tI3DMMwikPESt/bFe81EVkvIvu982sickZJCmgcoVGjRiQluVu2detWDh8+nGCJDMMwjLJEREpfRH6P28q2NTAbuN87twZWef5GCZOSksLxxx8PgKqyefPmBEtkGIZhlCUiben/FZgDtFDVu1X1CVW9G2gBzAX+VlICGvmxLn7DMAyjqESq9BsDk1U1nwW/dz0JaBRrwYzQmDGfYRiGUVQiVfqrgN+E8Tsd+Cw24hhHw1r6hmEYRlEpbJ5+ILcBM0SkMvAWsB03d/8SYAhwhYik+QKr6oFYC2o4rKVvGIZhFJVIlf5K7zwWeCTAXbzzp0HhkyRmzH0AACAASURBVIsjlBEea+kbhmEYRSVSpT8YW5GvVGBK3zAMwygqESl9VZ1awnIYEdKkSRP/782bN5OTk0OlSpF+uxmGYRgVmbiuyCciw0XkCxHZ4x3LRaTXUeKIiNwiIt+IyCER2SYij8ZL5tJG1apVqV+/PgC5ubls3bo1wRIZhmEYZYV4L8O7BbgLOBNoCywB3jrKqn6PA8O8eKcBPYEPS1jOUo0Z8xmGYRhFIa79wqo6J8hppIjcCHQAvggOLyKnADcBZ6jq1wFen5eclKWfpk2b8umnznbSxvUNwzCMSEnYhjsikiwiVwDVgU/CBLsY2AD0EJENIrJRRF4WkXpxE7QUEmjMZy19wzAMI1LCtvRFZAkwTFW/EZGrgPmq+nNxMxSRlsByoAqwD7hEVf8bJviJQFPgCmAQbgbBOGCuiHRQ1bwQ6Q8FhgLUr1+fZcuWFVfkfOzbty/maUZL4EY7K1asKBF5SkM544GVs3xh5SxfWDlLAFUNeQDZQHvvdy7QLlzYaA4gBWiGG9MfC+wETg8TdhJO0TcPcGvuuf32aHm1adNGY83SpUtjnma0zJs3T7060PPPP79E8igN5YwHVs7yhZWzfGHlLBrAKg2jFwsb098M9BeRfbhFeE7wfof7eFgb4UfGYWC9d7lKRM4CbgWuDRF8G5Cjqt8GuH0H5ABNKLgoUIXADPkMwzCMolCY0h8LTARuwbUqXwsTTjz/oq7ClwSkhvH7GKgkIiep6v88txNxcldYC7bAMf3MzEzy8vJISkqYeYZhGIZRRgir9FX1BRF5GzgZN0VuOBBRaz4c3vz6+bhehGOAK4HOQC/PfyxuGKGbF2URbjOfF0XkFs9tPK6Fv6o4spRlqlevTu3atfnll184fPgwP/74I8cff3yixTIMwzBKOYVO2VPVn4CfROQhYI6qFnclmAbANO+8GzdN7yJVfdfzPw44KSD/PBHpDUzAfXgcBN4HbtMQRnwViYyMDH755RfATdszpW8YhmEcjUiX4X0IQERSgJZAbeAX4L/eGH1EqOqgaP1VdRvQP9I8KgpNmzbls8/cjsabNm2iQ4cOCZbIMAzDKO1EPBAsIncCP+F23HvXO/8kIneUkGxGIZgxn2EYhhEtEbX0vfH0scBzwEyc8q8PXA6MFZFDqjqhxKQ0CmC77RmGYRjREukyvMOBR1V1ZIDbOuBDEdkF3IwbdzfihK3KZxiGYURLpN37jYGlYfyWAY1iIo0RMYHd+9bSNwzDMCIhUqWfCXQP43eB52/EkeDufbcIk2EYhmGEJ9Lu/QnABBGpDczGjenXw1nVD8J17xtxpGbNmhx77LHs2bOHAwcOsHPnTurWrZtosQzDMIxSTEQtfVV9Grge6AEswC2Ms9C7vkFVJ5aYhEZIRMSM+QzDMIyoiHjKnqq+gBvbbwp08M6NVXVyCclmHAWbtmcYhmFEQ6Td+wB4u/ds9g4jwVhL3zAMw4gG26WlDGNK3zAMw4gGU/plGOveNwzDMKLBlH4Zxlr6hmEYRjSY0i/DWEvfMAzDiIaolL6I1BKRc0TkShGp5blVERH7eEgA6enpVK1aFYA9e/awa9euBEtkGIZhlGYiUtYikiwifwO2AB8ArwIneN6vAw+WjHhGYdhc/aIzffp0MjIy6Nq1KxkZGUyfPj3RIhmGYZQ4kbbQHwGuA0YAJwIS4DcH6BNjuYwIsS7+6Jk+fTpDhw71L1+8adMmhg4daorfMIxyT6RK/yrgblV9iYJz9P+H+xAwEoC19KPn3nvv5cCBA/ncDhw4wL333psgiQzDMOJDpEq/Jk65hyIFSI6NOEa0WEs/OlSVzMzQ+0NlZmby4IMPWj0ahlFuiVTpfwlcHMbvIuCz2IhjRIu19KPjiSeeKNT/4Ycf5sQTT6R79+7MnDmTQ4cOxUkywzCMkifSZXj/ArwuIlWBfwIKtBKRS3Ab8fyuhOQzjoIp/cj5xz/+wZ///OejhlNV3n//fd5//33q1KnDwIEDufbaazn99NPjIKVhGEbJEekue3OAK4HzcbvrCTAZt63uQFV9N5J0RGS4iHwhInu8Y7mI9Iow7skisldE9kUSvqJg3fuRsWTJEq6++mr/9SmnnELjxo39MyCmTp3KP//5T3r06IHIETvVn3/+mfHjx9OyZUvat2/PCy+8wN69exNRBMMwjGITzS57s1Q1AzgV6AS0AJqo6qwo8tsC3AWcCbQFlgBvicgZhUUSkRRgBvBhFHlVCBo0aEBKSgrgFNS+ffZNFMx//vMfLrnkErKzswFo0aIFy5cvJzMzkyVLlrBx40auvvpq+vXrx8KFC9m4cSMPPfRQvl4UgE8//ZShQ4fSoEEDBg8ezMcff4zbg8owDKNsEPWiOqr6rap+oqrfaJRvPFWdo6oLVXW9l85IYC9uq97C+CvwBW5owQggKSmJJk2a+K+tiz8/mzZt4qKLLmLPnj0ANGzYkIULF1KrVq2wcZo0acIDDzzAhg0beO+997jsssv8H1bgLP1feuklOnXqRIsWLRg3bhzbt28v8bIYhmEUF4lEb3sL8xSKqt4ZVcYiyUB/4BWgjar+N0y4XsB4XO9AX+BpVa1eSLpDgaEA9evXbzNjxoxoxDoq+/bto3r1sNknhNtvv53PPnO2lGPHjqV9+/bFTrM0ljNa9uzZw0033eS31q9WrRoTJkzgxBOPzDCNtJy7d+/m/fffZ/78+SGHUZKTk+nYsSO9evWibdu2JCeXrgkt5eF+RoKVs3xh5SwaXbp0Wa2qbUN6qupRD+D7EMcuIA/4FdgQSTpeWi2BfUCOl0avQsIeB2wF2nvXg4B9kebVpk0bjTVLly6NeZrFZfDgwYozrtRnnnkmJmmWxnJGw4EDB/Tss8/210tKSoouWbKkQLhoy5mXl6crVqzQ6667TqtXr+5PP/Bo1KiR3n///bphw4YYlab4lPX7GSlWzvKFlbNoAKs0jF6M1JDvhBBHTVy3fCYwINIvEGAd0ApoDzwLvCwi4cyipwHPquqKKNKvcAQa81n3PuTm5jJgwAA+/vhjv9srr7xCly5dip22iPDb3/6WSZMmsW3bNl588UU6duyYL8yWLVsYPXo0J554IhdccAEzZswgKyvLv/RvUlKSLf1rGEZCiHTKXkhU9VMReQx4GmgTYZzDwHrvcpWInAXcClwbInhX4DwR8a3tL0CSiOQAw1R1UnHkLy/YtL0jqCo333wzb775pt/tiSee4PLLL495XtWrV+eaa67hmmuu4euvv2bKlCm88sor7Nixwx9m0aJFLFq0iGrVqnHo0CFycnIA/Ev/AgwYEM03s2EYRtGJxe54PwOnFFOG1DB+LXG9Ar7jAeCg99uM+jxs2t4RHn30USZOnOi/vu2227j11ltLPN/TTjuNcePGsWXLFmbPnl1g6t/+/fv9Ct/HgQMHGDlyZInLZhiG4SOilr6IpIVwTgFOAx4GvoownUeB+bj1+4/Bzf3vDPTy/McC7VS1G4CqfhkUvy2QF+xe0bGWvuPll1/Ot37+FVdcwWOPPRZXGVJSUujbty99+/YlMzOTqVOn8uKLL4a9L+GWBDYMwygJIm3p78NNrQs8fgY+BhoAwyJMpwFunH4dsBg4C7hIVRd6/scBJ0WYluHRsGFDv7X4jz/+SFZWVoIlij/vvvsuQ4YM8V936dKFqVOnkpQUi86sohE49a9evXphwxiGYcSLSMf0B+MskwPJwi22s1JVsyNJRFUHFdN/KjA1krwqEpUqVaJRo0b+1mRmZibNmzdPsFTxY/Xq1fTt29fffd6yZUvefPNNUlPDjRrFl6SkJJ544gmGDh1aYHe/rl27JkgqwzAqIhEpfU/ZGqWYpk2b+pX+pk2bKozS37BhAz179mT//v0ANG7cmIULF1KjRo0ES5Yfn7Hevffem69L/+WXX6Zfv3707NkzUaIZhlGBSFzfpxFTKqIx344dO+jRo4d/NbxatWrxzjvv0LBhwwRLFpoBAwawadMm9u7dy5lnnglAXl4el112GatXr06wdIZhVATCKn0R2SEi2yM94im0UZCKZsy3f/9+evfuzXfffQdAamoqb7/9Ni1atEiwZEenevXqzJ8/33/PfGWpCPfNMIzEUlj3/jMUHMc3SimBSr+8t/RzcnK44oorWLlyJeAWzHnttdfo1KlTgiWLnAYNGrBgwQI6duzI7t27+fHHH+nZsycff/wxNWvWTLR4hmGUU8IqfVUdFUc5jGJSUVblU1WGDRvGvHnz/G4TJkzg0ksvTaBURaNFixa89dZbdO/enezsbNauXcsll1zCO++8U2qMEA3DKF/YmH45oaJ07z/88MO88MIL/uu7776bESNGJFCi4tG5c2emTp3qv162bBlDhgyxLXsNwygRIl6GV0Q64JbKbQ5UCfZX1XYxlMuIksaNGyMiqCo//PAD2dnZVK5cOdFixZTJkyczatQo//XAgQN55JFHEidQjLjyyivZuHGjf3W+adOmkZGRwejRoxMsmWEY5Y2IWvoicgHwIdAI6ATswC3Y839AHcBWyEswqampHHfccYCzCN+yZUuCJYot8+bN44YbbvBfd+/enSlTpuRb6rYsc88993Ddddf5r//yl78wefLkBEpkGEZ5JNLu/YeBv+Mtlwvcr6pdca3+bGBZ7EUzoqW8Ttv79NNPueyyy8jNzQWgdevWzJ49u1z1ZIgIEydOpEePHn63G264gXfffTeBUhmGUd6IVOm3ABYCeTiL/moAqroJGAXYriGlgPI4rv/tt9/Su3dvDh48CLgPmwULFnDMMcckWLLYU6lSJWbNmkWrVq0At0Vwv379WLNmTYIlMwyjvBCp0s8CktRZF20j//r4e3Dd/kaCKW9K/6effqJHjx7s3LkTgDp16vDOO+/QoEGDBEtWchxzzDHMnz+fxo0bA7Bv3z569erF5s2bEyyZYRjlgUiV/n84sn3uYuAeEblARM7Ddf3/tySEM6KjPHXv7927l549e/L9998DULVqVebNm8cppxRnF+eywfHHH8+CBQs49thjAdi6dSs9e/Zk9+7dCZbMMIyyTqRKfzxHFuq5F9gPvAssBeoBw2MvmhEt5aWln52dTf/+/fnss88At2HNzJkzad++fYIlix+nn346b7zxBpUquQk2X375Jf369ePw4cMJlswwjLJMREpfVReo6jPe7x+ANriWfyugmarawuGlgPLQ0ldVrrvuunwGbM8++yx9+vRJoFSJoVu3bkyZMsV/vWjRIoYOHWpz+A3DKDKRTtnrIgFzo9Txnap+oarW9CglBO7NvnnzZr+1e1nivvvu4+WXX/ZfP/DAAwwdOjSBEiWWq666ioceesh//fLLL/Pwww8nUCLDMMoykXbvLwa2isjfRaRjSQpkFJ20tDTq1q0LuPXpt23blmCJomPixIn5FtsZPHhwvsV4Kir3338/11xzjf961KhR+VbxMwzDiJRIlX5L4AXgQuBfIpIpIo+JSNuSE80oCmWti3/69OlkZGSQlJTE8OFHTEN69uzJc889V24W3ykOIsLzzz/PBRdc4He77rrrWLRoUQKlMgyjLBLpmP5XqvqAqp4KnAlMBy4BVorIehH5S0kKaUROWTLmmz59OkOHDmXTpk35xqlPOOEEZs2aVa4W3ykulStXZvbs2ZxxxhmA68np27cv//2vTZwxDCNyot5wR1XXqOo9qtoM+B1QFbgn5pIZRaIstfRHjhzJgQMHCrhnZ2dTrVq1BEhUujn22GOZP38+DRs2BGDPnj307NmTH374IcGSGYZRVoha6YtIbREZIiLvA28A1YHXYi6ZUSTKUks/MzMzpLspsfA0atQo34qEW7ZsoVevXuzZsyfBkhmGURaI1Hr/WBG5WkQW4Fbk+zvwK3AFUE9VB0aYznAR+UJE9njHchHpVUj4ziIyR0S2icgBL+7gSPKqqJQlpR842yASd8NxxhlnMHv2bP8c/v/85z/079+f7OzsBEtmGEZpJ9KW/nbgOeAwMAin6C9T1TdU9VAU+W0B7sLZBbQFlgBvicgZYcJ3xK321w84HXgWmCQiV0aRZ4WiLHXvjxkzpoChXlpaGmPGjEmQRGWH7t27M2nSJP/1e++9x4033mhz+A3DKJRKEYa7AXhDVYvVh6iqc4KcRorIjUAH4IsQ4YM3S39WRLoAfbEhhZAEtvQzMzNR1VJrAd+7d+98Sqpx48aMHTuWAQMGJFCqssM111zDxo0b/fP2p0yZQkZGBvfdd1+CJTMMo7QiiWoZiEgy0B94BWijqhGZIYvIO8AWVR0Sxn8oMBSgfv36bWbMmBEjiR379u2jevXqMU0z1vzud79j7969ALz++uvUrl076jTiUc5///vf3HnnnQCcfPLJ+Vqu8aIs3M/CUFUeffRR3nvvPb/bvffem296H5T9ckaKlbN8YeUsGl26dFmtqqGn1KtqXA/cnP99QA6wC+gVRdzeQDbQLpLwbdq00VizdOnSmKcZa1q1aqW4vRJ0xYoVRUojHuUcNWqUX85hw4aVeH6hKAv382gcOnRIu3bt6q/LypUr6+LFi/OFKQ/ljAQrZ/nCylk0gFUaRi9Gbb0fA9bh1uxvjxujf1lETj9aJBE5G9elf7OqrixZEcs2gV38pXlcf/ny5f7fHTp0SKAkZZuUlBRef/11fvOb3wBuyuOll17KV199lWDJDMMobcRd6avqYVVdr6qrVPUeYA1wa2FxRKQTsBB4QFWfjYecZZlAY77SasGfl5fHihUr/Nem9ItHzZo1WbBgAccddxwAu3fvpmfPnmVuKWbDMEqWRLT0g0kCUsN5isi5OIX/kKqOj5tUZZiyMG3v66+/9u8PX69ePU488cQES1T2adKkCfPnz/ePDWZmZtKrVy/27duXYMkMwygtxFXpi8ijInKOiGSISEsRGQt0xi3ri4iMFZHFAeE74xT+c8B0EWngHXXjKXdZoyxM2wvu2i+tMwzKGq1bt+af//wnycnJAHz++efUrVuXrl27kpGRwfTp0xMsoWEYiSSs0heR70VkQ6RHhPk1AKbhxvUXA2cBF6nqQs//OOCkgPCDgDTgz7hFgXzHv6MpZEWjLLT0bTy/5OjRowfPPntkFCwrKwtVZdOmTQwaNIjLL7+cF198kXnz5rFy5Uo2btzI/v37Eyhx2SRws6jy/EGViHL68oznx2pFKWfYKXsiMg5nDezjCpwCfh+3WE894AJgPzBDVe8sWVGjp23btrpq1aqYprls2TI6d+4c0zRjzc8//0x6ejoA1apVY+/evVG3pEu6nC1atODrr78G4IMPPuDcc88tsbwKoyzcz6JSo0aNqJbnTUtLo169egWOunXrFnBLT08nJSUlZDrTp09n5MiRZGZm0qRJE8aMGRO3tRfidT99m0UF7h1RpUoVHnnkES655BKSkpKiPpKTk0lKSir0vxpYt40bN2b06NFceumlZGVlkZWVxaFDhwr8DuUWqf/333/PN998Q15enl+GpKQkWrRoQePGjalcuTKVKlWK+lyY3/Lly3n++ec5dOjIum+pqan8+c9/pnv37v56iuV51qxZDBs2LN/9TEtL46mnnqJ///7k5eWRm5tLXl5ezH6/9957jBs3Ll8509LSmDRpUrH/LyISdspeRPP0ReReoDtuet3+APfqwDxgkaqWup32SkLpt2rVipo1a+Zzu+yyy/wPTM+ePQvEGTRoEIMGDWLnzp3069evgP+NN97I5ZdfzubNmxk4sOCKxrfffjt9+vRh3bp1XH/99QX877vvPs4//3zWrFnDLbfcgqry8ccfk5ubC8CCBQu46KKL+OSTT7j33nsLxB8/fjytWrVi0aJF/OUv7jbu2rXLX87nn3+eU045hblz5/L4448XiP/qq6/SuHFjZs6cma+F6WP27Nmkp6czdepUpk6dSnZ2Np988gngto3dvn076enpTJw4kVmzZhWIv2zZMgDGjRvHvHnz8vlVrVqVhQtdR9Ho0aNZvHhxPv86derw+uuvA3DPPffk62EAt3vd+++/D8Att9zCmjVr8vk3b97cv37A0KFD+fbbb/P5t2rVivHjnanJH//4R7Zs2ZLPv0OHDowdOxaAvn378vPPP+fz79atG/fffz8AF110EQcPHszn37t3b/785z8DhFRmhT17H3zwQYHwsaRSpUo0aNCAjIwMUlJS+Prrrzl8+DC//PJLvkWX0tLSuOOOO/z3MZBQz14g0T57gc8tFHz2glmwYAFpaWkRPXtvvPEGe/bsYe/evWzevLnEVz8UESpXrkxSUhI5OTnk5ubaiosVgKZNmxZ7WLYwpR/pinzDgaGBCh9AVfd5PQIvAKVO6VdURITU1FT/V+uPP/6YYInyE9j6rF69OmlpaQmUpvySmpqarxXhIzk5mfT0dA4fPkzNmjXJycnhp59+4vDhw1Gln5OTw5YtWwp86ARz4MAB/va3v9GwYUOqV69OtWrV/PsGlFby8vLYv38/Tz/9NCtWrOCdd94p8MFW0qhq1PfEKPuE24gsVkTa0t+Hmx//Ygi/IcCTqnpMCchXLCpq9z64FuL8+fMBtyrfpZdeGlX8kixnnz59/C32m266iQkTJpRIPpFQWDlzc3O5//77mT59Olu2bKF169bE4nnq1KkTVapUYdGiRcVOqzBCdUGH6z5UVfbu3cv27dv9x44dO/JdB7rv2LEjX5dvtDRr1ozWrVvTqlUr/9k33bA4FPW53bJlCytWrPAfq1evJisrK+p0kpOTadSoUYEu3WiOolCrVi2qVKlClSpVSE1NLfA7lFskYYcMGcL27dsL5FevXj2mTJlCTk4OOTk5ZGdnR3SOJMzcuXML9HiB+4g966yz8tVtrM7h7rWIUK1atQJDMLH4vXLlypAf5SXd0o90JbzXgJ9xG9+kem6puGV0fwZeiySdeB8VdUU+VdVhw4b5V2h74oknoo5/zjnnaKVKlXT79u1hw4wfP14Bffvtt6NK+7TTTvPL9tprr0UtWywp7H6+8MIL/tUCX3nlFV2wYEFM8jz77LO1W7duMUnraEybNk2bNm2qIqJNmzbVadOmxSTd3Nxc3bFjh3711Ve6dOlSnTlzpj799NNao0YN/72N9qhfv75eeOGFevfdd+vMmTN13bp1mpubG5Vckfw/Dxw4oB999JE+9thj2rdvX23YsGFE8qWkpGiHDh301ltv1REjRmjVqlXz+aelpRW7fvPy8jQ3N1ezs7P10KFDevDgQd2/f7/u3btXGzduHFKupk2bFivPwpg2bZqmpaXFvJyWZ8nmSSEr8kWq9GsAbwJ5QC5u+dxc7/otoEYk6cT7qMhK/29/+5v/Qbr55pujjv/www8roBMmTAgbpm3btpqenq6HDx+OON2cnBxNTU31y/b9999HLVssKex+Xn755Zqenh7zPA8dOhRVncWCeD23oV5kqampOmDAAL3qqqu0ZcuWmpycHPGHQLVq1bRjx446fPhwfeGFF3TVqlV68ODBiMuZl5en69ev12nTpumIESO0TZs2WqlSpYjyPuGEE/QPf/iD/v3vf9dPP/1Us7KyCpS1JD6owpEIxeTLN57ltDyLT7GVvj8w/AY3je4u4GqgRTTx431UZKU/c+ZM/4vh4osvjjr+e++9p7Vr19Z27dqF9F+3bp0COmLEiKjS/eKLL/xy1a1bV/Py8qKWLZYUdj+7dOmiJ510UvyEKUHi+dwe7UV28OBBXbVqlb7wwgs6fPhw7dixo1arVi3iD4FKlSppy5YtdeDAgfrEE0/o0qVL9fnnn/fnWa9ePe3Xr5/27t1b09PTI/646Ny5s9599906Z84c/fHHH+NWX9GQCMWUSMrK+7a4xHPt/YQr5pI8KrLSX7Fihf+F1qpVq6jjL126VG+44QYF9Ntvvy3gf9999ymgn376qd/tzTff1N69e2vDhg21cuXK2rBhQx02bJju3r3bH+a5557zy9WjRw+/e7gu75EjR2pycnIB93/84x/arl07rVq1qh577LHau3dv/fLLLyMq286dO/X666/XBg0aaOXKlfXUU0/Vxx9/3P8B8t1334VUDK+++mrYNAcMGKCpqam6adMm7d27t1avXl3r1Kmjw4YN03379uULG8/ufR+l/bnNzc3VdevW6cyZM/Xuu+/WHj16aP369Ys8VHC049RTT9VBgwbpc889p2vWrNHs7OxEV0FUlPb7GSusnEWjMKUfsQmtiNQDbgfaAo2AS1X1KxH5E7BSVZcXmoARV2KxKt/AgQN57rnnmD59OqNGjcrn99prr9G8eXPatWvnd5syZQopKSncdNNN1KpVi88//5zJkyfz5Zdf+qePBU6Za9OmTZHkevTRR7nnnnu49NJLueqqq9i3bx8TJ07k7LPPZtWqVTRr1ixs3KysLDp37szXX3/NjTfeSFJSEt999x233347mZmZjB8/ngYNGvDqq68yZswYfv31V8aNGwfA2WefXahceXl59OjRg9/85jf89a9/ZeXKlUycOJHMzEzmzp1bpLJWFJKSkmjevDnNmzfnsssu87tv27aNNWvW8Pnnn/vP69evjyrtmjVr0r59e//Rrl07atWqFesiGEbZINzXQOABtAN2A+uBKbjx/DM9v0eB2ZGkE++jIrf08/LytEqVKv6Wza5du6KK7yvnSSedpM2aNcvn969//UsBHT16dD73/fv3F0jnpZdeUjiyxW/z5s39Mr3xxhv+cJG29L///ntNTk7WBx54IF+4rVu3ao0aNfSqq64qtFxPPvmkAjp58mR/OfPy8vSSSy5REdFvvvnGH/a8886LuHt/wIABCujgwYPzud95550K6KJFi45a1pKkrDy3kbBnzx796KOP9KmnntLBgwcX2qL/5ptvojYGLAuUp/tZGFbOokEMttZ9ElgKNAeuBwKXjFrpfRQYpQgRiclyvH/84x9Zv359vh3xpk2bhogUmPblzQIM5gAAIABJREFUm2+vquzevZudO3fSqVMnAFavXs3PP/+cb3Gbli1bRi3P7Nmzyc3N5YorrmDnzp3+o3Llyvz2t79lyZIlhcafN28e6enpDBo0yO8mItxxxx2oqn+aY1G55ZZbQl4HLypkFJ1jjjmGTp06MWLECKZMmZLvOQ+kadOmnHLKKSQllYZ9xQyjdBDpv+FMYKKq5kG+pXnBTdmrF1OpjJgQK6UPTtEDHD58mFmzZtGpUydOOOGEfGHXrl1Lnz59qF69OjVr1qRu3bqcfPLJgFvhL/DDAdyypdHi+2ho0aIFdevWzXe89957IecUB7Jx40ZOPvlk/4Y0Plq0aAHA999/H7VMPkSE5s2b53M77rjjOOaYY0rtxkflgTFjxhRY4CktLY0xY8YkSCLDKL1EOqa/Gwi3s92JwE+xEceIJbFQ+s2aNaNDhw7MnDmT8ePHs2DBAn755ZcCywXv2rWL8847j2rVqjF69GiaNWtGWloahw8fplevXuTl5RVYAjeQcOuN+5YS9uFbvGThwoUhV3UrbqvOdvsre/h6nBK13r9hlCUiVfpzgIdEZDng0x4qIum4HfDeKAnhjOIRqy12Bw4cyLBhw3j33XeZNm0aqamp9O/fP1+YxYsXs3PnTt566618Bm9r1671//attx+KWrVqsXXr1gLuwXL7jPSaNGnib51HQ0ZGBl988QW5ubn5Wvu+zX8C6yxaVJVvv/0237DFtm3b2Lt3b7HSNY7OgAEDGDBgQJlZMdMwEkWkzaK7gT3AWuBDz+053Ba5B4EHYi+aUVxitcXu5ZdfTkpKCs888wzz5s2jT58+BTYd8rWwg5cRfeyxx/zuK1euDJtHs2bN+Oqrr9ixY4ffLTMzk7fffjtfuH79+pGcnMyDDz4YcsnSwPih6NOnDzt27OCVV17xu6kq48aNQ0To1atXofGPhm/zneDr4qZrGIYRCyJq6avqryLSHhgIdMNtp/sLMBl4RVULLiBsJJxYtfRr165Nz549eeuttwBC7gR4zjnnULt2bQYMGMBNN91ElSpVmDt3Ljt37gTgp59+8u/ZXrNmTXbt2pUv/pAhQ/j73/9O9+7dufbaa/n111959tlnOe200/LtfNesWTPGjh3LnXfeSfv27bnkkkuoU6cOmZmZLFy4kNatWzN58uSwZRk6dCiTJ09m6NChfP7554gIjz32GAsWLOBPf/oTp5xySpHrqXLlyixfvpz+/fvTpUsXVq5cycsvv0zPnj05//zzi5yuYRhGrIh4AFRVD6vqFFW9UlW7q+oVqvqCKfzSS6xa+nBE0depU4eLLrqogH96ejoLFiygadOmPPTQQzz00EM0aNCABQsWALB582Z/2JNOOqlA/BYtWjBt2jT27dvHbbfdxowZM5gwYQI9evQoEPaOO+5gzpw5VKtWjbFjx3LLLbcwc+ZMzjjjDK677rpCy1G1alWWLVvGtddey6xZs5g4cSLr169n3LhxPPnkk1HVSTBJSUm88847ZGVlcddddzF37lxuuOEGZs6cWax0DcMwYka4uXyBB25efrswfm2A3EjSifdRkefpq7p17gPXGQ81jz4csS7nwIED/XI8/vjjMU27OMSqnL4V+UorZem5LQ5WzvKFlbNoEIN5+oWZNFcGcor60WGUHMnJyTRp0sR/XdzWfnEINOLr2LFjwuQwDMOoyIQd0xeRJkBGgFNrEQmeWF0Ft/FO0Sc3GyVK06ZN2bBhA+CU/mmnnRZ3GbZv387//vc/AFJSUmjdunXcZTAMwzAKN+S7BniQI0taPhsm3EFgSCSZichw3Ip+GZ7TV8BfVDXsMmgi0hJ4Grfq3y/A88BorwvDOAqxMuYrDoGL8rRp04bU1NSEyGEYhlHRKax7fyLQEvg/XPf+AO868DgFqK2q/4gwvy24bXnPxG3cswR4S0TOCBVYRI4F3sct/nMWcDNwB3BbhPlVeGJpzFdUAhfl6dChQ0JkKGmmTZtGVlZWosUwDMMolLAtfVXdAewAEJETgG2qerg4manqnCCnkSJyI/D/7d17vBV1vf/x1xuRBEEUVBAvoKalB0sF0TAVyUteHpl2So38STcq83Q1yzyV/c6P9JTp8ZT+0spjdeqYmZVU5gUBNbwA3iUUU0MuoiQiyB0+54+ZvRsWa+89e+9Za+211vv5eMyDPTPfNfP5rLXZnzXf+c7MO4DHy7xkAtAPODci1gBPSjoA+IKkK3y03zEXfTMza5H3jnzbA29u7xalETG3zZVlSNoGeD/QH2jrVm3vAO5NC36L24F/IzlF4LEEHah19/6GDRu2uCmPi76ZWe3kLfpPsvWDdkpt08F6oPUc/f0kgwBXAadHxBNtNB9Kckoga2lm3VZFX9IkYBLAkCFDmD59ep6wclu1alXh26ykJUuWtP48f/783LEXlefTTz/NmjXJd7YhQ4Ywf/585s+f3+3tFqXePs+ucp6NxXk2lqrm2da1fNkJOKbMdDrJ4L6/Aqfm2U66rT7Am0nO6V8KLANGttH2DuDHJcuGk3wBOaKjfTX7dfoREevXr49evXq1XiO/du3aXK8rKs/vfe97rfs+88wzC9lmkert8+wq59lYnGdjqeZ1+nlvwzujjVW/kfT/gA8AuR4YHsm4gGfT2dmSDgM+D3y0TPOXSI7os1oe4+sn++Ww7bbbsvvuu7feEe/FF19sfWhNNfh8vplZz9G955AmpgGndTOGtq7huh84quT+AMcDi4EXurHPplLLwXwu+mZmPUcRRf8U4LUOWwGSLpN0lKQRkg6SdCkwDvh5uv5SSVMzL/kFsBq4QdJISWeQPPHPI/c7oVaD+V566SWefz4ZdrHddttx8MEHV23fZma2tVzd+5JuKrO4D/BWYD/gqzn3NxT47/TfFSSX6Z0UEben63cDWp/GEhErJB0PXA3MBpYD3wWuyLk/o3ZH+tmj/NGjR9OnT5+q7dvMzLaWd/T+LmWWrQXuBb4QEX/Ms5GImNjZ9ZGM7D86z/atvFod6btr38ysZ8k7kO/YSgdildMTjvRd9M3Maq/T5/SV2EXt3anHepRs0a/Wkf769euZPXt267yLvplZ7eUu+pJOljSTpFv/JWCtpJmSTqlYdFaI7ON1Fy1axMaNlX8S8mOPPdZ6L/oRI0YwdGjplZdmZlZtuYq+pE8AU0juoPdZktvnfjadvzVdbz3Udttt11p0N23axKJFiyq+z5kz/3FnZR/lm5n1DHmP9L8KXBcRJ0TEDyLilvTfE4AfAhdXLkQrQrUH82XP548dO7bi+zMzs47lLfqDgVvaWPdrYFAx4VilVHswnwfxmZn1PHmL/jSS++2XcwxwTzHhWKVU80h/8eLFLFiwAIC+ffvytre9raL7MzOzfPJep/+fwI8kDQZ+C7xMcg/804GTgI9JOrClcXTyMbtWedU80s8e5R922GFsu+22Fd2fmZnlk7fot9wx7xPpFED2kr0/pf8qXZfrMbtWPdUs+h7EZ2bWM+Ut+r45T52rZve+B/GZmfVM3X20rtWJ7JH+ggUL2Lx5M716FfG8pS2tW7eOOXPmtM4fccQRhe/DzMy6pit35OstqV/pVIngrDjbb789O++8MwAbNmxgyZIlFdnPI488wvr16wHYd9992XXXXSuyHzMz67y8N+cZKOkaSUtI7si3ssxkPVw1zuv7Uj0zs54r7zn9G0guzfsh8CywvlIBWeUMHz68tev9b3/7W0XOt7vom5n1XHmL/ruAT0TE/1QyGKusagzmy47c9yA+M7OeJe85/QXA6koGYpVX6e79F198sfW+/ttvvz0jR44sfB9mZtZ1eYv+hcC/Stqrw5bWY1X6SD/btT9mzBh6987bkWRmZtWQ95K9P0o6DnhW0gvAa2XajCk4NitYpY/0fT7fzKxny1X0JV0OfA6YhQfy1a3Soh8RSGrnFZ3jom9m1rPl7X/9GHBxRFxayWCssnbccUcGDhzIihUrWLNmDa+88kph19GvXbuWhx9+uHXeN+UxM+t58p7TXw3M6bBVByRdJGmWpNclvSJpiqQOR3tJOlHS/ZJWSlom6XeS9u9uPM2oUl38c+bMYcOGDQDsv//+rTcCMjOzniNv0b8KmKTu9wWPA64BxgLjgY3AXZIGtfUCSXsDvwPuBQ4BjgP6An/sZixNqVKD+dy1b2bW8+Xt3t8ZOBx4WtJ0th7IFxHx5Y42EhEnZuclnQOsAI4EprTxslHAtsBFEbEpfd2lwN2Sdo6IZTlzMCp3pO+ib2bW8+Ut+v9MclS+LXB8mfUBdFj0yxhA0tuwvJ02s4ENwMck/QjoB5wLzHLB77xKFP2I8ON0zczqgCKidjuXbgL2A0a3HMW30e4o4FckPQ69gEeAkyLi5TJtJwGTAIYMGTLqxhtvLDTmVatW0b9//0K3WU0zZszgkksuAZLi/K1vfatsu87k+dJLL3H22WcD0K9fP2699Va22WabQuKttHr/PPNyno3FeTaWovM89thj50TE6LIrI6ImE3AFsBjYp4N2Q4FngG+TnNM/GpieTr3ae+2oUaOiaNOmTSt8m9U0a9asIOmZiYMOOqjNdp3J8xe/+EXrNo877rgCoqyeev8883KejcV5Npai8wRmRxt1Mfct0yTtA3wJeCcwCHiVZHDd5RHxXGe+hUi6EjgLODbHaz8NvBERF2Ze/yHgRZIBgfd1Zt/NrnQgXxRwrb7P55uZ1Ye8j9YdBTwKvI/kBj0/Tf99H/CIpEPz7lDSVcAHgfERMS/HS/oBpV3/LfN5rz6w1ODBg+nXrx8AK1eu5LXXtrq5Yqe56JuZ1Ye8RfNykvPoIyLiIxFxUUR8BNg7XX55no1Iuhr4MHA2sFzS0HTqn2lzqaSpmZf9AThU0jck7Zd+wfgvkiP9bt87oNlIKnQw35o1a3j00Udb531THjOznitv0R8DfDsitnjSXjp/OcnlfHmcRzJifyqwJDNdkGmzG7BvZh93k/QMnEbyBeN2ktH8746IN3Lu1zKKvFZ/9uzZbNy4EYC3vvWt7LTTTt3anpmZVU7ec/prgMFtrBsErM2zkYjo8ORxREwss+xGoNhh+E2syCN9X6pnZlY/8h7p/wG4TNI7swvT+Utp+8Y61gMVeaSfPZ8/duzYbm3LzMwqK++R/hdIboU7Q9IrwFJg13SaCXyxMuFZJRR1pB8RHsRnZlZHchX9iPg78E5J7wYOIznvvgR4MCLuqGB8VgHZot+dI/3nn3+el19O7o80cOBADjjggO6GZmZmFZT7On2AiPgT8KcKxWJVku3e786RfvYo//DDD6dXL19BaWbWk7X5V1rSYEm/lnRiO21OTNsU81B2q4ohQ4bQp08fAF599VVWrlzZpe14EJ+ZWX1p79Dsc8A+QHvd93eQXKvvc/p1pFevXoWc1/cgPjOz+tJe0f8A8IP0Pr5lpeuuJbmG3upId4v+G2+8weOPPw4kN/w5/PC8t2owM7Naaa/oDwfm5tjGX4ARhURjVdPdy/ZmzZrFpk3J3ZAPPPBABg4cWFBkZmZWKe0V/TXADjm20T9ta3Wku0f6vlTPzKz+tFf0Hwbek2Mbp6VtrY50t+h7EJ+ZWf1pr+hfDXxU0rltNZD0f0geoPP9ogOzyupO935E8MADD7TOexCfmVl9aPM6/Yi4JX0M7n9JOp/k+vwFQAB7AScCo4ErI+I31QjWitOdI/1nn32WZcuWAbDTTjux//77FxqbmZlVRrs354mIL0qaTnL53gXAm9JV64A/A6dFxO8rGqFVxLBhw+jduzcbN25k6dKlrFmzhr59++Z6bfZ8/hFHHOGb8piZ1YkO/1pHxJSIeBfJI3F3S6cBEXGcC3796t27N3vssUfr/IIFC3K/1oP4zMzqU+5DtIjYGBFL02ljJYOy6uhqF7+LvplZfXK/bBPrymC+lStX8sQTTwDJTXnGjBlTgcjMzKwSXPSbWFeO9B966CE2b94MwEEHHcQOO+S5lYOZmfUELvpNrCtH+u7aNzOrXy76TawrR/ou+mZm9ctFv4l1tuhv3rx5i5vyuOibmdWXqhZ9SRdJmiXpdUmvSJoiaWSO10nS5yTNk7RO0hJJl1Uj5ka25557IgmARYsWsX79+nbbP/PMM7z66qsADB48mP3226/iMZqZWXGqfaQ/DrgGGAuMBzYCd0ka1MHrvgucB3wZOAA4GbincmE2hz59+jBs2DAgubXuwoUL221f2rXf8oXBzMzqQ7t35CtaRJyYnZd0DrACOBKYUu41kt4C/Avwtoj4S2bVI5WKs5mMGDGCRYsWAclgvn322afNtj6fb2ZW32p9Tn9AGsPydtqcBjwHvFvSc5JekPQTSbtWJcIG15nz+i76Zmb1rapH+mVcBTwK3N9Om32A4cBZwESSB/5cDkyR9I6I2JxtLGkSMAlgyJAhTJ8+vdCAV61aVfg2e4oZM2aw9957A1vnuWrVKp566ikAevXqxdq1axvifWjkzzPLeTYW59lYqppnRNRkAq4AFgP7dNDuOpJCv39m2f7pssPbe+2oUaOiaNOmTSt8m7V07bXXRvpexrnnntu6vDTP22+/vbXdwQcfXN0gK6jRPs+2OM/G4jwbS9F5ArOjjbpYk+59SVcCZwPjI+K5DpovATZGxDOZZfNJBgHuVaEQm0be7n137ZuZ1b+qF31JVwEfJCn483K85M9Ab0n7ZpbtQ3JqonMPgret5L0rX7bojx07toIRmZlZpVT7Ov2rgQ+THOUvlzQ0nfpn2lwqaWrmZXcBDwPXSzpE0iHA9cCDwOwqht+Q9trrH50lCxcuZNOmTVu18U15zMwaQ7WP9M8jGbE/laTbvmW6INNmN6D1qD6SgXqnAi+TXJt/O7AQOC1KBvFZ5/Xt25ddd00uhNi4cSOLFy/eqs28efNYsWIFALvssku7l/WZmVnPVe3r9Du8m0tETCyzbAnw/krEZEkX/8svvwwkXfx77rnnFut9Ux4zs8ZQ6+v0rQfoaDDfzJkzW392176ZWf1y0bcOB/N5EJ+ZWWNw0bd2j/SXL1/OX/6S3P24d+/ejB49uqqxmZlZcVz0rd0j/QcffLD157e//e3069evSlGZmVnRXPSt3SN935THzKxxuOjbFkV/wYIFbN78jyshPYjPzKxxuOgbAwYMYNCgQQCsW7eOpUuXArBp06Ytuvc9iM/MrL656BtQvot/7ty5rFy5EoChQ4du0cbMzOqPi74B5Qfz+aY8ZmaNxUXfgPJH+h7EZ2bWWFz0DShf9D2Iz8yssbjoG7B19/6KFSt45plngOSmPKNGjapRZGZmVhQXfQO2PtKfO3du6/yhhx5K3759axGWmZkVyEXfgK2P9LNF3137ZmaNwUXfANhxxx0ZMGAAAKtXr/YgPjOzBuSibwBI2qKL/69//Wvrzy76ZmaNwUXfWmW7+FsMGzaMPffcs/rBmJlZ4Vz0rVW5O+6NHTvWN+UxM2sQLvrWqtyRvrv2zcwah4u+tSp3pO+ib2bWOKpa9CVdJGmWpNclvSJpiqSRnXj9fpJWSlpVyTibVWnR79OnD4ceemiNojEzs6JV+0h/HHANMBYYD2wE7pI0qKMXSuoD3AjcU8kAm1m57v2bb765+oGYmVlF9K7mziLixOy8pHOAFcCRwJQOXv7vwOPADOCYigTY5O68884t5tevX8+kSZMAmDBhQi1CMjOzAtX6nP6ANIbl7TWSdApwKvCZagTVrC6++OKtlq1evbrscjMzqz+KiNrtXLoJ2A8YHRGb2mizGzAHOCMiHpA0Efh+RPRvo/0kYBLAkCFDRt14442Fxrxq1Sr69y+767o3fvx4yv0+SOLuu++uQUSV18ifZ5bzbCzOs7EUneexxx47JyJGl10ZETWZgCuAxcA+HbSbCnwtMz8RWJVnH6NGjYqiTZs2rfBt9hTDhw8PYKtp+PDhtQ6tYhr588xyno3FeTaWovMEZkcbdbEm3fuSrgTOBsZHxHMdNB8PfEPSRkkbgR8D26fzkyodazOZPHky/fr122JZv379mDx5co0iMjOzIlV1IB+ApKuAs4BxETEvx0sOKpk/DbgYGAMsKji8ptYyWO/iiy9mwYIF7LXXXkyePNmD+MzMGkRVi76kq4FzgPcCyyUNTVetiohVaZtLgTER8S6AiHiyZBujgc2ly60YEyZMYMKECUyfPp1x48bVOhwzMytQtbv3zyMZsT8VWJKZLsi02Q3Yt8pxmZmZNbxqX6ff4ZNbImJiB+tvAG4oJiIzM7PmUevr9M3MzKxKXPTNzMyahIu+mZlZk3DRNzMzaxIu+mZmZk2ipvferzRJrwB/S2cHkjzRL6t0WXa+rZ93BpZ1M7RysXSlXZ6cyi1zns6zK5xn59o5z/aX1VOeba3rqXkOj4hdyrZq6/68jTYB13W0LDvfzs9t3tO4O7F0pV2enJyn83SeztN5dq9dW+vqMc9m6t6fkmPZlBw/VyqWrrTLk1O5Zc6zWM6zc+2cZ/vLnGexisizrXV1l2dDd+9XgqTZ0dYjCxuI82wszrOxOM/GUs08m+lIvyjX1TqAKnGejcV5Nhbn2ViqlqeP9M3MzJqEj/TNzMyahIu+mZlZk3DRNzMzaxIu+mZmZk3CRb8gkk6W9LSk+ZLOq3U8lSLpVknLJd1c61gqRdKekqZLmivpMUln1DqmSpE0M83xSUlfr3U8lSSpl6RZDf67+4KkxyU9KmlareOpFEkjJN2d/h99StLOtY6paJIOSj/HlmmNpPd2e7sevd99knoDfwHGA38HZgPvioglNQ2sAiQdC/QHzo2If651PJUgaTdgSEQ8KmlXYA7wlohYXePQCidph4h4XdI2wH3ApyLi0VrHVQmSPg0cBfRu4N/dF4CREbGq1rFUkqQZwNci4h5JA4G1EbGu1nFViqQBwPMkt9d9ozvb8pF+McYAcyPixbQw/AY4tcYxVURETANW1jqOSoqIJS2FLyJeBpaT3Bu74UTE6+mPfdKpIaVf3s6gea77bliS/gnYEBH3AETEikYu+KnTgLu6W/DBRR8ASUen3daLJIWkiWXanCfpeUlrJc2RdFRm9TDgxcz8QmD3CofdaQXkWReKzFPSaGBbtvx8e4Si8pT0IPAyyR+VHneUX1Ce3wG+BmyuRsxdUVCeAcxIT2NMqErgnVRAnvsBKyX9TtIjkv5v1YLvhIL/3p4J/LKIuFz0E/2BJ4HPAmtKV0o6E7gK+BZwCDATuE3SXi1NymyzJ5436W6e9aKQPCUNBn4KfDR65nmwQvKMiMNJvqQeLGlkpYPugm7lKeloICJiZtUi7poiPs8jI2IU8B7gq5IOqnjUndfdPHsD44DPkPSyjpJ0euXD7rSi/g7tCBwB3FZIVN19sk+jTcAqYGLJsgeBH5Ysmw9cmv48FvhdZt1k4OO1zqXoPDPLxgE31zqHSuYJvAm4Bzin1jlU+vPMrPsycEGtcyk6T+ArwCLgBeAl4A3gx7XOpQqf53dKt9HTpi5+nkeQ9Eq1rPsU8M1a51KpzxP4CPDTomLxkX4HJPUBRgF3lKy6g6TYAzwE/JOSUd99gdOBP1Qvyu7LmWfdy5OnJAE3AHdHxM+qGmBBcua5o9JRz5K2A04A5lUzzu7Kk2dEXBYRu0fECOAs4LaI+GhVA+2mnJ/n9umALyT1JxlY/FQ14+yunH+HZgGDJQ1O/68eA8ytXpTd18m/t2cCNxa1bxf9ju0MbAMsLVm+FBgKEBEbgc8DU4EngP8fEYurGWQBOswTQNJdwK+AkyUtlPSO6oVYiDx5HknyH+29+sflMj2xm7Q9efIcBNwu6XGSK05mRMTvqxdiIXL93jaAPHkOAe6T9BjwAMnR4azqhViIPH9vN5H0Sk0DHgeWATdVMcYi5P17uzNwMHBnUTvuXdSGmkDpOV1ll0XEFIp/PnItdJTncdUNp2LazDMi7qNxvhC3l+dzJEcbjaDd39vWRhHTgelViKdSOvo83171iCqjo79DdwBvq2pEldFRnstIvswVplH+sFXSMmATWx817MrW39LqmfN0nvXIeTrPelSzPF30OxAR60luznJ8yarjSUZbNgTn6TzrkfN0nvWolnm6e5/WQS9vTmd7AXtJOhh4NSIWAFcAP5P0EPBn4JMk1+b/oBbxdpXzdJ44zx7LeTpPqpFnrS9l6AkTySVoUWa6IdPmPJJLftaRfEM7utZxO0/n6TydZ63jdp71lafvvW9mZtYkfE7fzMysSbjom5mZNQkXfTMzsybhom9mZtYkXPTNzMyahIu+mZlZk3DRNzMzaxIu+tZQJF0iKSTdXmbdzZKmVzGWcWksI6u1z86QdICkeyW9kcY5otYxZUk6X1JNbyQi6QOSJpZZPl3SzTUIyaxbXPStUZ0g6bBaB9HDfQfYEXgP8A5gSW3D6ZE+AEysdRBmRXHRt0b0Kslzti+udSCVJGm7bm7ircCdETE1Ih6IiHU599u3m/ttWn7vrNZc9K0RBfAt4D2SDmqrUXoqYFmZ5SHp/Mz8C5Iul/QVSUskrZD0XSVOlvSUpJWSfitppzK7Gibp92k3+gJJnyyzz3dKmiFptaS/S/qhpAGZ9RPTuMakXctrgC+1k9vBkqam21su6eeShqTrRqTd5vsCn0+3O72dbYWkL0j6D0mvAE+ky0+RdKeklyW9LukBSSeUe48lHZKuXy3pEUlHlbR7k6TvS3pN0quSrgS2LRPL3un7/Hr6nk+R9OaSNiHp8+ln9Pd0/xek686V9Fy6n+vb++Ik6QbgfcAx6TZD0iUlbT4o6dk0ntsk7ZFZNyJ9zQRJP5X0GjAlXbdN+t4skLQu/R36YOa149PXDsssu1/SJkk7ZpY9IWly+vOOkn4kabGktem2f9hWftakav1QAk+eipyAS0ieVd0LmAfcmFl3MzC9tG2ZbQRwfmb+BWAhcAvwbpIehACuJHlIxhnABGA58IPM68al7V4k+RJyInBsHyn3AAAGSElEQVRtuuzUTLsjSR648UvgZOAcYBFwc6bNxPR1fwUuAI4FDmnjPdgFeA24H3gv8KE0/seBPsCbgCNIuvN/nv58YDvvaaRtf5nmf3K6/HzgM2lex5M8NWwTcGTJe7w63feHgZOAB9LPqF+m3ZXAWuCLaZtb0pgj0+ZNwHPA08CZJAX5yfS9GlQS78L0vT4B+M902beBacCpJA86WQd8pZ289wXuBh5O36MjgD3SddPTz3UmcBpwFslz0P+Yef2IzHt3dfoejU/XTQY2AP+avn/XpW3PTtf3BdYDZ6bz/dL5NcAp6bJBwGbg3en89SS/82cCx6Sf+3W1/j/pqWdNNQ/Ak6ciJzKFnKRQbgL2T+e7U/SfBbbJLHsI2AjsnVn2bWBpZn5cuq3rSrZ/J/BAZv5eYFpJm/Hpa0dmcgngszneg8tIiv4OmWVjskUlk9flObYXwCMdtOlF8qju24HrS97jaCl26bKD02UtxWpwWsy+XLK9eWxZ9D+Zvuf7ZJbtkRbDi0rinVayrSUkX8qy78lNwIMd5LXF70xm+XRgBbBTZtnn0n33TedHpPO/KXntIOAN4Bsly/8IPJ2Zvx/4fub3YRlwI3BZuuw9JL/fO6TzTwL/Uqv/e57qY3L3vjWy/wYWABcVsK3pEbEpM/8s8EJEPF+ybBdJfUpe+5uS+VuAUWkXbz+SQXQ3SerdMgH3kRwJjip57R9yxDoGuCMiXm9ZEBEPkRT5d+Z4fTlb7VfSHpJ+ImkRSTHeQHJkvX9J0w0kRbLF3PTflq7wg4DtgN9l4t2cnU+NAR6OiOcy7RaSPIu8NK+pJdt6HpiTfU9IPq/dS/PqhFkRsTwz35JX6TZL37uRJEfuvypZ/ktgf0m7pvP3Ai2nQY5O52eULHssk9OjwJcknSep9DMwA3xO3xpYRGwkOfr+kKTh3dzcayXz69tYJpIu9KyXy8z3BnYGdgK2Aa4hKY4t0zqSc9p7lrx2aY5Yd2uj3VKSo8yu2GJ7knoBtwJjga+TnG44DLiNpIBnvZ4WXgAiYn36Y0u7oem/5d6nrM7klffz6s5gyHLbo8w2S2PerY3lLfMt40LuAUam5/CPIin69wKj07EILctanA/8luTzeFrSfEln5czFmoSLvjW660mKx5fLrFtLSYFuYyBed+1aZn4jSXftayRdwN8gKZql0/Ulr81z3fqSMvsEGEJyZUNXlO73zcAhJN3JP46IGRExm+RcdGe9lP5b7n3KqkRe1VD63rVcGlmay5D035Zc/kzyJXIcyXiCe4CngFXAu4BDyRT9iHgtIj4TEUOBtwMPAj+XdGAxaVgjcNG3hhbJZWiXAx/hH0dYLRYCAyRlu2NPoHinl5mfExGbIuINkoFtb4mI2WWmxV3Y34PAiSWj/w8jOcd8XxdzKNVS3Fsv80t7U47swraeIPkCdlpmW72y86kHSU6L7J1ptztJb0NReZXqbm9AOU+SDG58f8nyDwDPRMQrAOmpgyeBz5Ocu38kIoIk1wtJeovK5h0Rj5Nc3dGL5NJMMyD5pTFrdNcCXyUpDjMyy/9EMoDseknfBfYmGSxWtJPSy6pmkIz0P54tC9qFwFRJm0kGjq0E9gJOAS6OiGc6ub8rgE8Bt0v6d6A/yeC+J4BfdyeRjHkkX5q+K+lrwADgmyQj6TslIv4u6Trgm5I2khzNfjyNO+sGkh6b2yR9naQQXkLSY3Jt19Lo0DzgNEnvJcl3cRe/iLWKiFcl/Qfwr2m+s0l+L04Gzi5pfg/waeD2zJiSe0lurDQ/Ilp6SZB0H8n4kSdJehc+TjJg8KHuxGuNxUf61vAiYjXJJWGly5eRXPa1B8m50A8BHyxtV4CPkXTF/pbkcrFPR8StmTjuIxmUtQvwM5JruS8kuSQszzn8LaRHiseSHD3/D8nlYvcCx2fOp3dL2oNyBslpipuBfwMuZcsvVZ1xIcmpjK+TxLyY5MtL6T6PIynEPwZ+AvwNGBcRlerevwa4I41tFjCpoO1+neT9+hTwe5LP/0MRcWNJu5bu+3vKLCs9yr+f5CqPm0muTNgZOCkd7GgGgJLeIjMzM2t0PtI3MzNrEi76ZmZmTcJF38zMrEm46JuZmTUJF30zM7Mm4aJvZmbWJFz0zczMmoSLvpmZWZP4X1K7bv597wAbAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 576x360 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure(figsize=(8,5))\n",
    "plt.title(\"Computing pi with \\nincreasing number of random throws\",fontsize=20)\n",
    "plt.semilogx(n_exp, pi_exp,c='k',marker='o',lw=3)\n",
    "plt.xticks(fontsize=14)\n",
    "plt.yticks(fontsize=14)\n",
    "plt.xlabel(\"Number of random throws\",fontsize=15)\n",
    "plt.ylabel(\"Computed value of pi\",fontsize=15)\n",
    "plt.hlines(y=3.14159,xmin=1,xmax=1e7,linestyle='--')\n",
    "plt.text(x=10,y=3.05,s=\"Value of pi\",fontsize=17)\n",
    "plt.grid(True)\n",
    "plt.show()"
   ]
  }
 ],
 "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.7.0"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}