0. Python_Intro (optional).ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "#Python Introduction\n",
    "##### Created by Professor Eric Larson, contact at eclarson@smu.edu for errors or questions\n",
    "\n",
    "In this notebook, we will introduce the basics of programming in python as well as the installation of different modules and libraries. Debugging in python can be a pain, but there are many different ways of debugging depending on what you are most comfortable with. \n",
    "\n",
    "We will go over:\n",
    "+ Analysing and debugging errors \n",
    "+ Setting breakpoints with pdb \n",
    "+ Using an IDE with python that can look at memory, set breakpoints, etc.\n",
    "\n",
    "\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "slideshow": {
     "slide_type": "subslide"
    }
   },
   "source": [
    "![Guido Von Rossum](https://raw.githubusercontent.com/eclarson/DataMiningNotebooks/master/data/guido.png)\n",
    "\n",
    "Python was founded by a guy named Guido Von Rossum in 1989. I will let him explain his motivations.\n",
    "From wikipedia:\n",
    ">Over six years ago, in December 1989, I was looking for a \"hobby\" programming project that would keep me occupied during the week around Christmas. My office ... would be closed, but I had a home computer, and not much else on my hands. I decided to write an interpreter for the new scripting language I had been thinking about lately: a descendant of ABC that would appeal to Unix/C hackers. I chose Python as a working title for the project, being in a slightly irreverent mood (and a big fan of Monty Python's Flying Circus).\n",
    "\n",
    "        -Guido van Rossum in 1996\n",
    "\n",
    "So let that be a lesson! Python is not named after the snake, it is after Monty Python. However, even though the language is a bit quirky, do not let that fool you. It is absolutely made for people to use as a serious programming language and tool. Whether you have novice programming skills or are a sesoned veteran at OOP. \n",
    "\n",
    "For example, take a look at the most popular programming languages as identified by the number of projects in GitHub, and the number of related questions on stack overflow. Notice the other langauges next to python, like java, php, c++, ruby, and all the heavy hitters. Its also on every top ten list of programming languages needed to get a job in programming. Check it out. Graphic credit given to RedMonk.\n",
    "![Python](https://raw.githubusercontent.com/eclarson/DataMiningNotebooks/master/data/python_ranking.png)\n",
    "\n",
    "Moreover, python is one of the fastest growing languages used in scientific computing and data science. According to a recent survey, the three most helpful languages or packages for doing data science are R, SAS, and python. You will likely make use of all three for a given project or set of projects. See http://www.r-bloggers.com/data-science-toolbox-survey-results-surprise-r-and-python-win/ for a high level description of the survey.\n",
    "\n",
    "# Running Python\n",
    "Python can be run in various ways. We will cover how to run it from the interpreter prompt, command line scripts, and from iPython notebook."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "# this is a comment\n",
    "\n",
    "# place this at the top of python file to enable running as >>./filename.py\n",
    "#! /usr/bin/python\n",
    "\n",
    "# otherwise you can run with >>python filename.py {or} >>python27 filename.py\n",
    "\n",
    "# those commands can be run from a terminal window"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "hello world\n"
     ]
    }
   ],
   "source": [
    "# hello world from iPython Notebook\n",
    "print (\"hello world\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Syntax\n",
    "Okay, so lets first dive into some of the syntax in python programming."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Variable type examples:\n",
      "<class 'int'>\n",
      "<class 'int'>\n",
      "<class 'float'>\n",
      "<class 'bool'>\n",
      "True\n",
      "False\n"
     ]
    }
   ],
   "source": [
    "# ============================================\n",
    "# Variables Example\n",
    "int_val = 8\n",
    "long_val = 23423423235 # in python 3 longs are integers!!\n",
    "float_val = 2.546923764237654239\n",
    "bool_val = False\n",
    "\n",
    "print (\"Variable type examples:\")\n",
    "print (type(int_val))\n",
    "print (type(long_val))\n",
    "print (type(float_val))\n",
    "print (type(bool_val))\n",
    "\n",
    "# testing for the type of a variable\n",
    "temp = isinstance(float_val,float)\n",
    "print (temp)\n",
    "print (isinstance(float_val,int))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "Arithmetic examples:\n",
      "2.6666666666666665\n",
      "2.6666666666666665\n",
      "2.6666666666666665\n",
      "False\n",
      "False\n",
      "True\n",
      "8.0\n",
      "5\n",
      "81\n"
     ]
    }
   ],
   "source": [
    "#============================================\n",
    "# Arithmetic and casting\n",
    "print (\"\\nArithmetic examples:\")\n",
    "print (8 / 3)\n",
    "print (float(8) / 3)\n",
    "print (float(8) / float(3))\n",
    "\n",
    "print (True and False)  # logicals\n",
    "print (8 == 3        )  # logical equallty\n",
    "print (5 <= 6        )  # logical comparison\n",
    "\n",
    "print (2.0*4.0)  # multiplication\n",
    "print (65%6   ) # remainder, modulus\n",
    "print (3**4   ) # 3 to the fourth power"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "String example:\n",
      "A string is double or single quotes\n",
      "Three quote means that the string goes over\n",
      "multiple lines\n",
      "This also spans multiple lines but has no newline\n"
     ]
    }
   ],
   "source": [
    "#============================================\n",
    "# strings example\n",
    "print (\"\\nString example:\")\n",
    "str_val = \"A string is double or single quotes\"\n",
    "str_val_long = '''Three quote means that the string goes over\n",
    "multiple lines'''\n",
    "str_val_no_newline = '''This also spans multiple lines \\\n",
    "but has no newline'''\n",
    "\n",
    "print (str_val)\n",
    "print (str_val_long)\n",
    "print (str_val_no_newline)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "A\n",
      "tr\n",
      "s\n",
      "uotes\n",
      "A string is double or single quotes\n"
     ]
    },
    {
     "ename": "TypeError",
     "evalue": "'str' object does not support item assignment",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mTypeError\u001b[0m                                 Traceback (most recent call last)",
      "Cell \u001b[0;32mIn[6], line 9\u001b[0m\n\u001b[1;32m      6\u001b[0m \u001b[38;5;28mprint\u001b[39m (str_val[\u001b[38;5;241m-\u001b[39m\u001b[38;5;241m5\u001b[39m:]) \u001b[38;5;66;03m# the last five elements\u001b[39;00m\n\u001b[1;32m      7\u001b[0m \u001b[38;5;28mprint\u001b[39m (str_val[\u001b[38;5;241m0\u001b[39m:\u001b[38;5;241m5\u001b[39m] \u001b[38;5;241m+\u001b[39m str_val[\u001b[38;5;241m5\u001b[39m:]) \u001b[38;5;66;03m# print the first five elements, then from the fifth and on\u001b[39;00m\n\u001b[0;32m----> 9\u001b[0m \u001b[43mstr_val\u001b[49m\u001b[43m[\u001b[49m\u001b[38;5;241;43m5\u001b[39;49m\u001b[43m]\u001b[49m \u001b[38;5;241m=\u001b[39m \u001b[38;5;124m'\u001b[39m\u001b[38;5;124mG\u001b[39m\u001b[38;5;124m'\u001b[39m \u001b[38;5;66;03m# this is an error, strings are immutable once they are set\u001b[39;00m\n",
      "\u001b[0;31mTypeError\u001b[0m: 'str' object does not support item assignment"
     ]
    }
   ],
   "source": [
    "# string can be accessed in a variety of different ways\n",
    "\n",
    "print (str_val[0]  ) # initial element \"0th\" element\n",
    "print (str_val[3:5]) # elements 3 and 4, but not 5\n",
    "print (str_val[-1] ) # the last element in the string\n",
    "print (str_val[-5:]) # the last five elements\n",
    "print (str_val[0:5] + str_val[5:]) # print the first five elements, then from the fifth and on\n",
    "\n",
    "str_val[5] = 'G' # this is an error, strings are immutable once they are set"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "A string is double or single quotesA string is double or single quotes\n",
      "True\n",
      "Eric\n",
      "a string is double or single quotes\n",
      "A STRING IS DOUBLE OR SINGLE QUOTES\n",
      "['this', ' is', ' separated', ' by', ' commas']\n"
     ]
    }
   ],
   "source": [
    "# some common operations for strings\n",
    "print (str_val*2) # mutliply is like adding the number of times, here it repeats the string twice\n",
    "print ('Python' > 'Java') # compare the strings alphabetically \n",
    "\n",
    "print (\"eric\".capitalize()) # the dot operator works like it does in most other OOP languages\n",
    "print (str_val.lower())\n",
    "print (str_val.upper())\n",
    "print (\"this, is, separated, by, commas\".split(',')) # this results is returned as a list, which we need to talk about!"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Tuples, Lists, Sets, and Dictionaries "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[45, 67, 'not a number', 'A string, appended as a new element in the list']\n",
      "[45, 67, 'not a number', 'A string, appended as a new element in the list', ['a list', 'within another list', 442]]\n",
      "['a list', 'within another list', 442]\n",
      "['A string, appended as a new element in the list', ['a list', 'within another list', 442]]\n"
     ]
    }
   ],
   "source": [
    "# A List is one of the most powerful tools in python from which\n",
    "# most abstract types get created and implemented\n",
    "\n",
    "# a list is very much like the mutable version of a tuple\n",
    "# it can hold any type of information\n",
    "a_list = [45, 67, \"not a number\"]\n",
    "\n",
    "# we can add to a list through the append function\n",
    "a_list.append(\"A string, appended as a new element in the list\")\n",
    "\n",
    "print (a_list)\n",
    "\n",
    "# Lists can have other lists in them\n",
    "tmp_list = [\"a list\", \"within another list\", 442]\n",
    "a_list.append(tmp_list)\n",
    "\n",
    "print (a_list)\n",
    "\n",
    "# all of the indexing we learned from before still works with lists\n",
    "print (a_list[-1])\n",
    "print (a_list[-2:])\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(45, 67, 'not a number')\n",
      "not a number\n"
     ]
    },
    {
     "ename": "TypeError",
     "evalue": "'tuple' object does not support item assignment",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mTypeError\u001b[0m                                 Traceback (most recent call last)",
      "Cell \u001b[0;32mIn[9], line 10\u001b[0m\n\u001b[1;32m      7\u001b[0m \u001b[38;5;28mprint\u001b[39m (a_tuple[\u001b[38;5;241m2\u001b[39m])\n\u001b[1;32m      9\u001b[0m \u001b[38;5;66;03m# but you cannot change a tuple, its immutable!!\u001b[39;00m\n\u001b[0;32m---> 10\u001b[0m \u001b[43ma_tuple\u001b[49m\u001b[43m[\u001b[49m\u001b[38;5;241;43m2\u001b[39;49m\u001b[43m]\u001b[49m \u001b[38;5;241m=\u001b[39m \u001b[38;5;124m'\u001b[39m\u001b[38;5;124mhey\u001b[39m\u001b[38;5;124m'\u001b[39m \u001b[38;5;66;03m# this will give you an error!!\u001b[39;00m\n",
      "\u001b[0;31mTypeError\u001b[0m: 'tuple' object does not support item assignment"
     ]
    }
   ],
   "source": [
    "# tuples are immutable lists and are designated by commas\n",
    "# you can store ANYTHING in side a tuple, its basically a complex object container\n",
    "a_tuple = 45, 67, \"not a number\" \n",
    "print (a_tuple)\n",
    "\n",
    "# you can access a tuple with square brackets\n",
    "print (a_tuple[2])\n",
    "\n",
    "# but you cannot change a tuple, its immutable!!\n",
    "a_tuple[2] = 'hey' # this will give you an error!!"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "{'orange', 'banana', 'apple', 'pear'}\n",
      "True\n",
      "False\n",
      "{'c', 'a', '!', 'd', 'r', 'b'}\n",
      "{'d', 'r', 'b', '!'}\n",
      "{'c', 'm', 'z', 'a', 'l', '!', 'd', 'r', 'b'}\n",
      "{'a', 'c'}\n",
      "{'m', 'z', 'l', '!', 'd', 'r', 'b'}\n"
     ]
    },
    {
     "ename": "AttributeError",
     "evalue": "'frozenset' object has no attribute 'add'",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mAttributeError\u001b[0m                            Traceback (most recent call last)",
      "Cell \u001b[0;32mIn[10], line 22\u001b[0m\n\u001b[1;32m     19\u001b[0m \u001b[38;5;28mprint\u001b[39m (a \u001b[38;5;241m^\u001b[39m b )  \u001b[38;5;66;03m# letters in a or b but not both\u001b[39;00m\n\u001b[1;32m     21\u001b[0m a_immutable \u001b[38;5;241m=\u001b[39m \u001b[38;5;28mfrozenset\u001b[39m(a)\n\u001b[0;32m---> 22\u001b[0m \u001b[43ma_immutable\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43madd\u001b[49m(\u001b[38;5;124m'\u001b[39m\u001b[38;5;124me\u001b[39m\u001b[38;5;124m'\u001b[39m)   \u001b[38;5;66;03m# the set is immutable, so we cannot add to it, this will give an error!\u001b[39;00m\n",
      "\u001b[0;31mAttributeError\u001b[0m: 'frozenset' object has no attribute 'add'"
     ]
    }
   ],
   "source": [
    "# Sets, taken from the Python sets tutorial\n",
    "# https://docs.python.org/2/tutorial/datastructures.html#sets\n",
    "basket = ['apple', 'orange', 'apple', 'pear', 'orange', 'banana']\n",
    "fruit = set(basket) # create a set without duplicates\n",
    "print (fruit)\n",
    "print ('orange' in fruit )# fast membership testing\n",
    "print ('crabgrass' in fruit)\n",
    "\n",
    "# Demonstrate set operations on unique letters from two words\n",
    "a = set('abracadabra')\n",
    "b = set('alacazam')\n",
    "a.add('!') # also add the some punctuation\n",
    "\n",
    "# set operations\n",
    "print (a     )   # unique letters in a\n",
    "print (a - b )  # letters in a but not in b\n",
    "print (a | b )  # letters in either a or b\n",
    "print (a & b )  # letters in both a and b\n",
    "print (a ^ b )  # letters in a or b but not both\n",
    "\n",
    "a_immutable = frozenset(a)\n",
    "a_immutable.add('e')   # the set is immutable, so we cannot add to it, this will give an error!"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "{'dog': 4, 'cat': 4, 'human': 2}\n",
      "4\n",
      "2\n",
      "===============\n",
      "{'dog': 4, 'human': 'two', 'bird': 'two and some wings', 45: 'a key that is not a string'}\n"
     ]
    }
   ],
   "source": [
    "# Dictionaries map keys to values.\n",
    "\n",
    "# Here we setup a key as a string and the value as a number\n",
    "num_legs = { 'dog': 4, 'cat': 4, 'human': 2 }\n",
    "\n",
    "# You access Subscripts via the \"key\"\n",
    "print (num_legs)\n",
    "print (num_legs['dog'])\n",
    "print (num_legs['human'])\n",
    "\n",
    "print ('===============')\n",
    "# Entries can be added, updated, or deleted.\n",
    "# Again, these are just containers for any memory type\n",
    "num_legs['human'] = 'two'\n",
    "num_legs['bird'] = 'two and some wings'\n",
    "num_legs[45] = 'a key that is not a string' # notice that the key does not need to be a string\n",
    "# the key just needs to be some immutable memory\n",
    "\n",
    "del num_legs['cat']\n",
    "print (num_legs)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# ADTs, Loops, and Conditionals"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "for loop output:\n",
      "DataMining <class 'str'>\n",
      "This statement is in the loop\n",
      "8 <class 'int'>\n",
      "This statement is in the loop\n",
      "23423423235 <class 'int'>\n",
      "This statement is in the loop\n",
      "2.546923764237654 <class 'float'>\n",
      "This statement is in the loop\n",
      "False <class 'bool'>\n",
      "This statement is in the loop\n",
      "This statement is outside the loop\n"
     ]
    }
   ],
   "source": [
    "\n",
    "#============================================\n",
    "# for loop example with list\n",
    "print (\"\\nfor loop output:\")\n",
    "\n",
    "list_example = [int_val, long_val, float_val, bool_val]\n",
    "list_example.insert(0, \"DataMining\")\n",
    "\n",
    "# notice that the loop ends with a colon and \n",
    "# is designated by the tab alignment\n",
    "for val in list_example:\n",
    "    print (str(val) + ' ' + str(type(val)))\n",
    "    print (\"This statement is in the loop\")\n",
    "    \n",
    "print (\"This statement is outside the loop\")\n",
    "    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0\n",
      "1\n",
      "3\n",
      "B. exiting for loop without break\n"
     ]
    }
   ],
   "source": [
    "import random\n",
    "val = 0\n",
    "for i in range(0, random.randint(1, 10) ):\n",
    "    val += i\n",
    "    print (val)\n",
    "    if val>20:\n",
    "        print (' A. leaving the loop on break')\n",
    "        break # break out of loop\n",
    "else: # this else belongs to the for loop\n",
    "    print ('B. exiting for loop without break')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "DataMining \t is at index \t 0\n",
      "8 \t is at index \t 1\n",
      "23423423235 \t is at index \t 2\n",
      "2.546923764237654 \t is at index \t 3\n",
      "False \t is at index \t 4\n"
     ]
    }
   ],
   "source": [
    "# more for loop examples\n",
    "\n",
    "# you can also get the index using the enumerate example\n",
    "for index, val in enumerate(list_example):\n",
    "    print (str(val), '\\t is at index \\t', index)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "What is your name?  It is lancelot.\n",
      "What is your quest?  It is the holy grail.\n",
      "What is your favorite color?  It is blue.\n"
     ]
    }
   ],
   "source": [
    "\n",
    "# this is a classic example for zipping, provided by the official python tutorial\n",
    "# notice the references to Monty Python\n",
    "\n",
    "# say you have two lists of equal size, that you would like to \n",
    "# loop through without indexing, you can use the \"zip\" function\n",
    "questions = ['name', 'quest', 'favorite color']\n",
    "answers = ['lancelot', 'the holy grail', 'blue']\n",
    "for q, a in zip(questions, answers):\n",
    "    print ('What is your %s?  It is %s.' % (q, a))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "dict_keys(['dog', 'human', 'bird', 45])\n",
      "dog => 4\n",
      "human => two\n",
      "bird => two and some wings\n",
      "45 => a key that is not a string\n",
      "===============\n",
      "dog => 4\n",
      "human => two\n",
      "bird => two and some wings\n",
      "45 => a key that is not a string\n",
      "===============\n",
      "human => two\n",
      "beast is not present.\n",
      "cat is not present.\n",
      "dog => 4\n",
      "45 => a key that is not a string\n"
     ]
    }
   ],
   "source": [
    "# Looping through dictionaries\n",
    "# Get all the keys.\n",
    "print (num_legs.keys())\n",
    "for k in num_legs.keys():\n",
    "    print (k, \"=>\", num_legs[k])\n",
    "\n",
    "print ('===============')\n",
    "# you can also use the iter_items function\n",
    "for k, v in num_legs.items():\n",
    "    print (k, \"=>\", v)\n",
    "\n",
    "print ('===============')\n",
    "# Test for presence of a key.\n",
    "for t in [ 'human', 'beast', 'cat', 'dog', 45 ]:\n",
    "    print (t,end=' ')\n",
    "    if t in num_legs:\n",
    "        print ('=>', num_legs[t])\n",
    "    else:\n",
    "        print ('is not present.')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[0, 1, 16, 81, 256, 625, 1296, 2401, 4096, 6561]\n",
      "['NAME', 'QUEST', 'FAVORITE COLOR']\n",
      "{0: 0, 1: 1, 2: 16, 3: 81, 4: 256, 5: 625, 6: 1296, 7: 2401, 8: 4096, 9: 6561}\n",
      "260\n"
     ]
    }
   ],
   "source": [
    "# list comprehensions and dictionary comprehensions\n",
    "\n",
    "# comprehensions are like shorthand for loops where each\n",
    "# pass through the loop saves the result in a list or dictionary\n",
    "\n",
    "# for example, imagine we want to take every element in a range to the fourth power\n",
    "times_four = [x**4 for x in range(10)]\n",
    "print (times_four)\n",
    "# this save the fourth power of the numbers 0 through 9\n",
    "# the result of each pass through the loop saves a value\n",
    "# in other languages this is also known as \"lambda functions\"\n",
    "\n",
    "# you can also call functions inside a comprehension, and take \n",
    "# advantage of the data type. For example, if each element is a string\n",
    "# we can use the string member methods\n",
    "questions = ['name', 'quest', 'favorite color']\n",
    "quest_upper = [x.upper() for x in questions]\n",
    "print (quest_upper)\n",
    "\n",
    "# you can also do comprehensions with dictionaries\n",
    "times_four = {x:x**4 for x in range(10)}\n",
    "# notice that the comprehension is wrapped in curly braces and the key for the value \n",
    "# is given inside the code itself, followed by a colon\n",
    "print (times_four)\n",
    "\n",
    "# you can also nest comprehensions inside each other, but this is not always good coding practice\n",
    "# especially if it makes the code less readable\n",
    "\n",
    "# Finally, all of the enumerate, zipping, and slicing we performed also applies to \n",
    "# list comprehensions\n",
    "x_array = [10, 20, 30]\n",
    "y_array = [7, 5, 3]\n",
    "# this prints the sum of the multiplication of the arrays\n",
    "print (sum(x*y for x,y in zip(x_array, y_array)) ) \n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "Stack Example:\n",
      "['LIFO', 0, 1, 2, 3, 4]\n",
      "=============\n",
      "4\n",
      "=============\n",
      "['LIFO', 0, 1, 2, 3]\n",
      "=============\n"
     ]
    }
   ],
   "source": [
    "#============================================\n",
    "# array as a stack\n",
    "print (\"\\nStack Example:\")\n",
    "list_example = []\n",
    "list_example.append('LIFO')\n",
    "\n",
    "for i in range(0, 5):\n",
    "    list_example.append(i)\n",
    "\n",
    "print (list_example)\n",
    "print (\"=============\")\n",
    "val = list_example.pop()\n",
    "print (val)\n",
    "print (\"=============\")\n",
    "print (list_example)\n",
    "print (\"=============\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "Queue Example:\n",
      "deque([9, 8, 7, 6, 5, 'FIFO'])\n",
      "=============\n",
      "FIFO\n",
      "=============\n",
      "deque([9, 8, 7, 6, 5])\n",
      "=============\n",
      "5\n",
      "6\n",
      "7\n",
      "8\n",
      "9\n"
     ]
    }
   ],
   "source": [
    "#============================================\n",
    "# array as a queue\n",
    "print (\"\\nQueue Example:\")\n",
    "from collections import deque # this is an import, we will get back to that later\n",
    "# essentially this is a set of utility functions shipped with python\n",
    "\n",
    "q_example = deque()\n",
    "q_example.appendleft(\"FIFO\")\n",
    "for i in range(5, 10):\n",
    "    q_example.appendleft(i)\n",
    "\n",
    "print (q_example)\n",
    "print (\"=============\")\n",
    "val = q_example.pop()\n",
    "print (val)\n",
    "print (\"=============\")\n",
    "print (q_example)\n",
    "print (\"=============\")\n",
    "\n",
    "# pop and print each element\n",
    "while len(q_example) > 0:\n",
    "    print (q_example.pop())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "Conditional Example:\n",
      "a is true\n",
      "b is false\n"
     ]
    }
   ],
   "source": [
    "#============================================\n",
    "# conditional example\n",
    "print (\"\\nConditional Example:\")\n",
    "\n",
    "a, b = True, False\n",
    "\n",
    "if a:\n",
    "    print (\"a is true\")\n",
    "elif a or b:\n",
    "    print (\"b is true\")\n",
    "else:\n",
    "    print (\"neither a or b are true\")\n",
    "\n",
    "# conditional assignment\n",
    "val = \"b is true\" if b else \"b is false\"\n",
    "print (val)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "I. Everybody is a five!\n",
      "II. These are the same object!\n",
      "III. Wish we had the same objects...\n",
      "IV. Everybody is a five!\n",
      "V. Wish we had fives...\n",
      "VI. Got through nested conditional\n"
     ]
    }
   ],
   "source": [
    "# more on conditionals, \"is\" versus ==\n",
    "\n",
    "# I. the traditional == works as expected\n",
    "a=5\n",
    "b=5\n",
    "if a==b:\n",
    "    print (\"I. Everybody is a five!\")\n",
    "else:\n",
    "    print (\"I. Wish we had fives…\" )\n",
    "    \n",
    "# II. the \"is\" function is for object comparison, much like comparing pointers\n",
    "a=327676\n",
    "b=a\n",
    "if a is b:\n",
    "    print (\"II. These are the same object!\")\n",
    "else:\n",
    "    print (\"II. Wish we had the same objects...\" )\n",
    "    \n",
    "# III. while these have the same value, they are not the same memory\n",
    "a=327676\n",
    "b=327675+1\n",
    "if a is b:\n",
    "    print (\"III. These are the same object!\")\n",
    "else:\n",
    "    print (\"III. Wish we had the same objects...\" )\n",
    "\n",
    "# IV. you would expect this to say wish we had fives,\n",
    "# but small integers like this are cached so right now they do point to the same memory\n",
    "a=5\n",
    "b=4+1\n",
    "if a is b:\n",
    "    print (\"IV. Everybody is a five!\")\n",
    "else:\n",
    "    print (\"IV. Wish we had fives...\" )\n",
    "    \n",
    "# V. but if we change the memory, that caching gets released\n",
    "b = b*2.0\n",
    "b = b/2.0\n",
    "if a is b:\n",
    "    print (\"V. Everybody is a five!\")\n",
    "else:\n",
    "    print (\"V. Wish we had fives...\")\n",
    "    \n",
    "# you can also perform nested conditionals, like bounding\n",
    "if 5 < 8 < 6: # not true because 8 is not less than 6\n",
    "    print ('VI. How did we get here')\n",
    "elif 4 < 18 < 22:\n",
    "    print (\"VI. Got through nested conditional\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Pythonic\n",
    "Now that we have gone through a consierable number of examples and data types, we need to talk about what it means to write code in python that is \"pythonic\". In a nutshell, there are many ways to implement something in python. The \"best\" method is one that is readable, simple, and brief. \n",
    "\n",
    "There many different coding “styles”\n",
    "+ “best” styles get the distinction of “pythonic”\n",
    "    + ill formed definition\n",
    "    + changes as the language matures\n",
    "+ pythonic code is:\n",
    "    + simple and readable\n",
    "    + uses dynamic typing when possible\n",
    "    + uses underscores to separate names\n",
    "    + names of variables are distinctive\n",
    "+ …or to quote Tim Peters…, run the next cell\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "The Zen of Python, by Tim Peters\n",
      "\n",
      "Beautiful is better than ugly.\n",
      "Explicit is better than implicit.\n",
      "Simple is better than complex.\n",
      "Complex is better than complicated.\n",
      "Flat is better than nested.\n",
      "Sparse is better than dense.\n",
      "Readability counts.\n",
      "Special cases aren't special enough to break the rules.\n",
      "Although practicality beats purity.\n",
      "Errors should never pass silently.\n",
      "Unless explicitly silenced.\n",
      "In the face of ambiguity, refuse the temptation to guess.\n",
      "There should be one-- and preferably only one --obvious way to do it.\n",
      "Although that way may not be obvious at first unless you're Dutch.\n",
      "Now is better than never.\n",
      "Although never is often better than *right* now.\n",
      "If the implementation is hard to explain, it's a bad idea.\n",
      "If the implementation is easy to explain, it may be a good idea.\n",
      "Namespaces are one honking great idea -- let's do more of those!\n"
     ]
    }
   ],
   "source": [
    "import this # the zen of python"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Pop Quiz!!!\n",
    "Which of these is more pythonic: You want to add the numbers between zero and N. \n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "8128\n",
      "8128\n",
      "8128.0\n"
     ]
    }
   ],
   "source": [
    "N = 128\n",
    "\n",
    "sum_value = 0\n",
    "for i in range(N):\n",
    "    sum_value += i\n",
    "print (sum_value)\n",
    "\n",
    "# OR\n",
    "print (sum(range(N)))\n",
    "\n",
    "# OR\n",
    "print (N*(N-1)/2 )"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Functions and Classes"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "Function Example:\n",
      "upper case\n",
      "data mining\n"
     ]
    }
   ],
   "source": [
    "#============================================\n",
    "print (\"\\nFunction Example:\")\n",
    "\n",
    "# create and call a function\n",
    "# the function can be defined almost anywhere in file, as long as it is defined before it gets used\n",
    "def make_strings_lowercase(str_input):\n",
    "    assert isinstance(str_input, str)  # test the type of input\n",
    "    return str_input.lower()\n",
    "\n",
    "# now we are back on the main execution\n",
    "print (make_strings_lowercase(\"UPPER CASE\"))\n",
    "print (make_strings_lowercase(\"Data Mining\"))\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[1, 2, 3, 4, 5]\n",
      "====================\n",
      "[1, 2, 3, 4, 5]\n",
      "[1, 2, 3, 4, 5]\n",
      "a cool X value\n",
      "[1, 2, 3, 4, 5]\n",
      "a cool X value\n",
      "a cool Y value\n",
      "====================\n",
      "(4, 1024)\n"
     ]
    }
   ],
   "source": [
    "# more functions examples\n",
    "def show_data(data):\n",
    "    # print the data\n",
    "    print (data)\n",
    "\n",
    "some_data = [1,2,3,4,5]\n",
    "show_data(some_data)\n",
    "\n",
    "# you can also define default values for the functions\n",
    "print ('====================')\n",
    "def show_data(data,x=None,y=None):\n",
    "    # print the data\n",
    "    print (data)\n",
    "    if x is not None:\n",
    "        print (x)\n",
    "    if y is not None:\n",
    "        print (y)\n",
    "\n",
    "some_data = [1,2,3,4,5]\n",
    "show_data(some_data);\n",
    "show_data(some_data,x='a cool X value')\n",
    "show_data(some_data,y='a cool Y value',x='a cool X value') \n",
    "\n",
    "# as well as have multiple return types in the function\n",
    "print ('====================')\n",
    "def get_square_and_tenth_power(x):\n",
    "    return x**2,x**10\n",
    "\n",
    "print (get_square_and_tenth_power(2) )\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "name: Heart has 1 beats per second\n",
      "This is a long string meant to be so long that the memory for it is not cached by python\n",
      "True\n",
      "False\n"
     ]
    }
   ],
   "source": [
    "# This is a class that inherits from a generic object\n",
    "class BodyPart(object):\n",
    "    kind = \"This is a long string meant to be so long that the memory for it is not cached by python\"\n",
    "    # this is a class variable, shared across all instances\n",
    "    def __init__(self,name):\n",
    "        self.name = name # the name attribute is unique to each instance of the BodyPart class\n",
    "\n",
    "# now define a class that sub classes from the defined BodyPart CLass\n",
    "class Heart(BodyPart):\n",
    "    def __init__(self,rate=60,units=\"minute\"):\n",
    "        self.rate = rate\n",
    "        self.units= units\n",
    "        super(Heart,self).__init__(\"Heart\")\n",
    "\n",
    "    def print_rate(self):\n",
    "        print (\"name: \" + str(self.name) + \" has \" + str(self.rate) + \" beats per \" + self.units) \n",
    "        \n",
    "\n",
    "my_heart = Heart(1,\"second\")\n",
    "my_heart.print_rate()\n",
    "\n",
    "generic_part = BodyPart(\"Foot\")\n",
    "print (my_heart.kind)\n",
    "print (my_heart.kind is generic_part.kind) # true, these are the same memory location\n",
    "\n",
    "# take the following for example, these are not the same object\n",
    "a = \"This is a long string meant to be so long that the memory for it is not cached by python\"\n",
    "b = \"This is a long string meant to be so long that the memory for it is not cached by python\"\n",
    "print (a is b) # not the same memory location"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "4 -1\n",
      "A: 4.5\n",
      "B: 0.6666666666666666\n",
      "C: 3\n",
      "D: 8\n",
      "\n",
      "That's odd, nothing went wrong.\n"
     ]
    }
   ],
   "source": [
    "# error handling and exceptions\n",
    "# Python exception handling.\n",
    "\n",
    "# example reissued from: http://sandbox.mc.edu/~bennet/python/code/exc_py.html\n",
    "\n",
    "# Choose two random integers.\n",
    "import random\n",
    "i = random.randrange(0, 8)\n",
    "j = random.randrange(-1, 6)\n",
    "print (i, j)\n",
    "\n",
    "# Get a nice little array, then try a bunch of dangerous stuff.\n",
    "some = [3, 10, 0, 8, 18];\n",
    "try:\n",
    "    # We try to execute this block.\n",
    "    den = some[j] / i\n",
    "    print (\"A:\", den)\n",
    "    frac = (i + j) / den\n",
    "    print (\"B:\", frac)\n",
    "    if frac < 2:\n",
    "        k = 3\n",
    "    else:\n",
    "        k = 'mike'\n",
    "    print (\"C:\", k)\n",
    "    print (\"D:\", some[k])\n",
    "    # This is the catch block.\n",
    "except ZeroDivisionError:\n",
    "    print( \"\\nDivision by zero.\")\n",
    "except TypeError as err:\n",
    "    # The detail provides extra information about the exception.\n",
    "    print( \"\\nSome type mismatch:\", err)\n",
    "except IndexError as err:\n",
    "    print (\"\\nSome value is out of range:\", err)\n",
    "except:\n",
    "    # Except without an exception name catches any exception.\n",
    "    print( \"\\nSomething else went wrong.\")\n",
    "\n",
    "# An else attached to an except block is run if no exception occurrs.\n",
    "else:\n",
    "    print (\"\\nThat's odd, nothing went wrong.\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "Could not read file: [Errno 2] No such file or directory: 'some_file.txt'\n"
     ]
    },
    {
     "ename": "NameError",
     "evalue": "name 'file' is not defined",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mNameError\u001b[0m                                 Traceback (most recent call last)",
      "Cell \u001b[0;32mIn[28], line 14\u001b[0m\n\u001b[1;32m     11\u001b[0m     \u001b[38;5;28mprint\u001b[39m (data)\n\u001b[1;32m     12\u001b[0m \u001b[38;5;28;01mfinally\u001b[39;00m:\n\u001b[1;32m     13\u001b[0m     \u001b[38;5;66;03m# this always gets called, close the file if its open\u001b[39;00m\n\u001b[0;32m---> 14\u001b[0m     \u001b[38;5;28;01mif\u001b[39;00m \u001b[38;5;129;01mnot\u001b[39;00m \u001b[43mfile\u001b[49m\u001b[38;5;241m.\u001b[39mclosed:\n\u001b[1;32m     15\u001b[0m         file\u001b[38;5;241m.\u001b[39mclose()\n",
      "\u001b[0;31mNameError\u001b[0m: name 'file' is not defined"
     ]
    }
   ],
   "source": [
    "# the with command and opening a file\n",
    "\n",
    "# the regular way of opening a file, lots of error checking\n",
    "try:\n",
    "    file = open(\"some_file.txt\")\n",
    "    data = file.read()\n",
    "except IOError as err:\n",
    "    print (\"\\nCould not read file:\", err)\n",
    "else:\n",
    "    print (\"Read successfully, file contents:\")\n",
    "    print (data)\n",
    "finally:\n",
    "    # this always gets called, close the file if its open\n",
    "    if not file.closed:\n",
    "        file.close()\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [
    {
     "ename": "FileNotFoundError",
     "evalue": "[Errno 2] No such file or directory: 'some_file.txt'",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mFileNotFoundError\u001b[0m                         Traceback (most recent call last)",
      "Cell \u001b[0;32mIn[29], line 7\u001b[0m\n\u001b[1;32m      1\u001b[0m \u001b[38;5;66;03m# the \"file\" class actually has a built in __init__ and __exit__\u001b[39;00m\n\u001b[1;32m      2\u001b[0m \u001b[38;5;66;03m# so we can really on their class implemntation to close the file propoerly\u001b[39;00m\n\u001b[1;32m      3\u001b[0m \u001b[38;5;66;03m# and handle any exceptions gracefully\u001b[39;00m\n\u001b[1;32m      4\u001b[0m \n\u001b[1;32m      5\u001b[0m \u001b[38;5;66;03m# in that respect, the \"with\" statement is a protocol that the \"file class\" adopts\u001b[39;00m\n\u001b[0;32m----> 7\u001b[0m \u001b[38;5;28;01mwith\u001b[39;00m \u001b[38;5;28;43mopen\u001b[39;49m\u001b[43m(\u001b[49m\u001b[38;5;124;43m\"\u001b[39;49m\u001b[38;5;124;43msome_file.txt\u001b[39;49m\u001b[38;5;124;43m\"\u001b[39;49m\u001b[43m)\u001b[49m \u001b[38;5;28;01mas\u001b[39;00m file:\n\u001b[1;32m      8\u001b[0m     data \u001b[38;5;241m=\u001b[39m file\u001b[38;5;241m.\u001b[39mread()\n\u001b[1;32m      9\u001b[0m     \u001b[38;5;28mprint\u001b[39m (\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mRead successfully, file contents:\u001b[39m\u001b[38;5;124m\"\u001b[39m)\n",
      "File \u001b[0;32m~/miniconda3/envs/mlenv2023/lib/python3.8/site-packages/IPython/core/interactiveshell.py:284\u001b[0m, in \u001b[0;36m_modified_open\u001b[0;34m(file, *args, **kwargs)\u001b[0m\n\u001b[1;32m    277\u001b[0m \u001b[38;5;28;01mif\u001b[39;00m file \u001b[38;5;129;01min\u001b[39;00m {\u001b[38;5;241m0\u001b[39m, \u001b[38;5;241m1\u001b[39m, \u001b[38;5;241m2\u001b[39m}:\n\u001b[1;32m    278\u001b[0m     \u001b[38;5;28;01mraise\u001b[39;00m \u001b[38;5;167;01mValueError\u001b[39;00m(\n\u001b[1;32m    279\u001b[0m         \u001b[38;5;124mf\u001b[39m\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mIPython won\u001b[39m\u001b[38;5;124m'\u001b[39m\u001b[38;5;124mt let you open fd=\u001b[39m\u001b[38;5;132;01m{\u001b[39;00mfile\u001b[38;5;132;01m}\u001b[39;00m\u001b[38;5;124m by default \u001b[39m\u001b[38;5;124m\"\u001b[39m\n\u001b[1;32m    280\u001b[0m         \u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mas it is likely to crash IPython. If you know what you are doing, \u001b[39m\u001b[38;5;124m\"\u001b[39m\n\u001b[1;32m    281\u001b[0m         \u001b[38;5;124m\"\u001b[39m\u001b[38;5;124myou can use builtins\u001b[39m\u001b[38;5;124m'\u001b[39m\u001b[38;5;124m open.\u001b[39m\u001b[38;5;124m\"\u001b[39m\n\u001b[1;32m    282\u001b[0m     )\n\u001b[0;32m--> 284\u001b[0m \u001b[38;5;28;01mreturn\u001b[39;00m \u001b[43mio_open\u001b[49m\u001b[43m(\u001b[49m\u001b[43mfile\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[43margs\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[43mkwargs\u001b[49m\u001b[43m)\u001b[49m\n",
      "\u001b[0;31mFileNotFoundError\u001b[0m: [Errno 2] No such file or directory: 'some_file.txt'"
     ]
    }
   ],
   "source": [
    "# the \"file\" class actually has a built in __init__ and __exit__\n",
    "# so we can really on their class implemntation to close the file propoerly\n",
    "# and handle any exceptions gracefully\n",
    "\n",
    "# in that respect, the \"with\" statement is a protocol that the \"file class\" adopts\n",
    "\n",
    "with open(\"some_file.txt\") as file:\n",
    "    data = file.read()\n",
    "    print (\"Read successfully, file contents:\")\n",
    "    print (data)\n",
    "\n",
    "# is the file closed? Let's check\n",
    "print (file.closed)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1. Just initialized body part with name Lungs\n",
      "2. Building up from \"with\" command\n",
      "3. Hi, my name is: Lungs\n",
      "4. Exit was called, no errors\n"
     ]
    }
   ],
   "source": [
    "# you can also define your own classes that adopt the \"with\" protocol,\n",
    "# if you are interested check out the example \n",
    "class BodyPart(object):\n",
    "    def __init__(self,name):\n",
    "        self.name = name\n",
    "        print ('1. Just initialized body part with name', name)\n",
    "        \n",
    "    def __enter__(self):\n",
    "        print ('2. Building up from \"with\" command' )\n",
    "        return self\n",
    "        \n",
    "    def __exit__(self, type, value, traceback):\n",
    "        if value is not None:\n",
    "            print ('4. An error occurred,', value)\n",
    "        else:\n",
    "            print ('4. Exit was called, no errors')\n",
    "            \n",
    "    def print_self(self):\n",
    "        # 5/0 # uncomment to raise an error\n",
    "        print ('3. Hi, my name is:', self.name)\n",
    "\n",
    "with BodyPart(\"Lungs\") as bp:\n",
    "    bp.print_self()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Importing Libraries\n",
    "Python shipd with many libraries already packaged. You simply import them using the \"import\" command. However, you can also install many third party libraries, of which there are many for data wrangling. In the next examples we will use some built in libraries, as well as ones that are third party.\n",
    "\n",
    "FOr installing the third party libraries, you will need to install a package mananger. There are three main ones that you may want to use on your system:\n",
    "+ Pip, https://pypi.python.org/pypi/pip\n",
    "+ Homebrew, http://brew.sh\n",
    "+ Macports, https://www.macports.org\n",
    "    \n",
    "I have had the best luck using pip (even on mac, after installing brew). But you will need to explore the options and decide which one is right for your system. If you are using windows, pip is the way to go. If you do not want the setup to affect other libraries installed on your system, consider setting up a virtual environment. http://docs.python-guide.org/en/latest/dev/virtualenvs/\n",
    "\n",
    "Now go our and install the Numpy and Matplotlib systems on your system! It might be as simple as:\n",
    "     > pip install numpy\n",
    "     > pip install matplotlib"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0.9998741276738751\n",
      "-1.2246467991473532e-16\n"
     ]
    }
   ],
   "source": [
    "# with the package installed, we can import the libraries and use them.\n",
    "\n",
    "# The numpy package is great for working with groups of numbers as matrices\n",
    "# and has many functionalities that are comparable to Matlab and Octave\n",
    "import numpy as np\n",
    "t = np.linspace(0, 1, 100)  # 100 equally spaced array of numbers from 0.0 to 1\n",
    "y = np.sin(2 * np.pi * t)  # take sine of the value\n",
    "\n",
    "print (np.max(y))\n",
    "print (np.median(y))\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjgAAAGdCAYAAAAfTAk2AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABaxklEQVR4nO3deVhU9f4H8PfMwMyAwgCyDAiIiIoLiqIiLqXJFZdfanlLyz3TFm2zTW6llZUt3m638ma5W5YtV8vUSEPNDUVRDBVXQNYBAZlhX2bO74/BKa6ioAxnlvfrec7T4+Gcw/sctPPhO99FIgiCACIiIiIbIhU7ABEREVFLY4FDRERENocFDhEREdkcFjhERERkc1jgEBERkc1hgUNEREQ2hwUOERER2RwWOERERGRzHMQOIAaDwYDc3Fy4uLhAIpGIHYeIiIiaQBAElJaWws/PD1Lpzdto7LLAyc3NRUBAgNgxiIiI6DZkZWXB39//psfYZYHj4uICwPiAXF1dRU5DRERETaHT6RAQEGB6j9+MXRY41z6WcnV1ZYFDRERkZZrSvYSdjImIiMjmsMAhIiIim8MCh4iIiGwOCxwiIiKyOSxwiIiIyOawwCEiIiKbwwKHiIiIbA4LHCIiIrI5LHCIiIjI5pi1wNm3bx/uvfde+Pn5QSKR4Mcff7zlOXv37kXfvn2hUCgQEhKCdevWXXfM8uXLERQUBKVSicjISCQmJrZ8eCIiIrJaZi1wysvL0bt3byxfvrxJx6enp2Ps2LEYPnw4kpOT8eyzz+LRRx/Fr7/+ajrm22+/xYIFC7B48WIcP34cvXv3RkxMDAoKCsx1G0RERGRlJIIgCK3yjSQSbNmyBRMmTGj0mJdffhnbt2/HqVOnTPsmT56MkpISxMXFAQAiIyPRv39/fPrppwAAg8GAgIAAPPXUU1i4cGGTsuh0OqhUKmi1Wq5FRUREZCWa8/62qMU2ExISEB0d3WBfTEwMnn32WQBATU0NkpKSEBsba/q6VCpFdHQ0EhISGr1udXU1qqurTX/W6XQtG5ysksEgIKOoHGc1pSgsq0ZljR5VtQZU1upRpzegXVsF/NyUULsq4efmBLVKCUcZu60REVkDiypwNBoNfHx8Guzz8fGBTqdDZWUlrl69Cr1ef8Njzp492+h1ly5dijfeeMMsmcl6VNXqse/8Fey7cAVncnU4qylFRY2+yec7OcrQL8gdg0M8MbiTJ7r7uUImvfWKtkRE1PosqsAxl9jYWCxYsMD0Z51Oh4CAABETUWupqtXj9/NXsCMlD/GpBSirrmvwdYWDFKFqF7R3d4LSQQalXAalgwwOMgmulFYjT1uJPG0V8rRVqKzVY/+FQuy/UAgAUDk5YkyYGlMiO6Bne5UYt0dERI2wqAJHrVYjPz+/wb78/Hy4urrCyckJMpkMMpnshseo1epGr6tQKKBQKMySmSzT1fIarD2UgfWHMqCtrDXt91UpEdNDjT6Bbujh54qgdm3g0ISPnQwGARcKynDwYiEOXSrCkbQiaCtr8U1iFr5JzELvADdMiQzEvb384CSXmfPWiIioCSyqwImKisKOHTsa7Nu1axeioqIAAHK5HBEREYiPjzd1VjYYDIiPj8f8+fNbOy5ZoAJdFVbuT8PGI5mmj5/8VEqMDvPFmDBf9Alwg/Q2PlaSSiXoqnZBV7ULHhnSEXV6A45mXMU3iZn45VQeTmaV4GRWCd7ZkYr5w0MwLaoDFA4sdIiIxGLWAqesrAwXL140/Tk9PR3Jycnw8PBAYGAgYmNjkZOTgw0bNgAAHn/8cXz66ad46aWX8Mgjj2D37t347rvvsH37dtM1FixYgBkzZqBfv34YMGAAPvroI5SXl2PWrFnmvBWycOXVdfh3/AWsO5SBmjoDAKCHnyvmDQ9BTA91i/eVcZBJEdWpHaI6tUNhWXf8kJSNr49kIrO4Am9tT8X6hAy8PCoUY8N8IZGwnw4RUWsz6zDxvXv3Yvjw4dftnzFjBtatW4eZM2ciIyMDe/fubXDOc889hzNnzsDf3x+vvfYaZs6c2eD8Tz/9FB988AE0Gg3Cw8Px8ccfIzIyssm5OEzctvx2Jh+Lt55GTkklAKBfB3fMuycEw7p4tWpxUac34IekbHy46zwKSo2j9sID3PD6uB4ID3BrtRxERLaqOe/vVpsHx5KwwLENedpKvLH1DOJOawAA/u5OeHN8Dwzv6i1qq0lFTR1W7kvH5/suoaJGD5lUgnnDQ/DUPSEcZk5EdAdY4NwCCxzr9/PJXMRuTkFZdR1kUgnmDA3GMyM6W1QH34LSKry9PRU/JecCAHr5q/CvSeHo5NVW5GRERNaJBc4tsMCxXjV1BryzIxXrDmUAAPoEuuGd+8LQzddyf44/n8zFK1tSoKuqg9JRin+M6YZpAzuwbw4RUTOxwLkFFjjWKaekEvM2HkdyVgkA4MlhnbDgb12aNMxbbBptFV784aRpDp0J4X54d2IvKB0tp8WJiMjSNef9bflvBiIA+85fwdiP9yM5qwSuSgesntEPL40KtYriBgDUKiXWzxqARf/XHTKpBD8m5+KhlYdxpbT61icTEVGzWcfbgezaD0nZmLXuKEoqahHWXoXtTw/FiG4+tz7RwkilEjwypCM2PDIAKidHnMgswfhPD+BMLtdGIyJqaSxwyGIJgoAVv1/CC9+fhN4g4L4+7fH941EI8HAWO9odGRziiS1PDkKwZxvkaqvw9xWHsOtM/q1PJCKiJmOBQxbJYBCwZFsq3v3FuIjqY3cF458P9LaZPivBXm2x5cnBGBLiiYoaPR7/KglbT+aKHYuIyGawwCGLU1NnwLPfJmPNwXQAwKtjuyF2TLfbWmLBkqmcHbF2Vn/c37c99AYBz246gR+SssWORURkEyxqLSqiOr0Bz2w6gV9OaeAglWDZA70xoU97sWOZjaNMimV/7w2FgxTfJGbhxR9OolZvwEMDAsWORkRk1VjgkMXQGwQ8//1J/HJKA7lMis+nR2B4V2+xY5mdVCrB2xPC4CiTYkPCZcRuTkGt3oDpUUFiRyMislr8iIosgsEg4JUtKfgpORcOUgn+M6WvXRQ310ilErwxrgceHdIRALDop9P46vBlkVMREVkvFjgkOkEQ8Oa2M9h0NAtSCfDvyX0Q3d36hoHfKYlEglfGdsMTwzoBAF776RR+SckTORURkXVigUOie//Xc1h3KAMSCbDsgd4Y28tX7EiikUgkeCmmKx4aEAhBAJ7ZlIzDaUVixyIisjoscEhUXx2+jM/2XgIAvD0hDPf39Rc5kfgkEgnemtATI7v7oEZvwJz1xzgZIBFRM7HAIdH8fv4KFm89DQB4YWQXPBzJkUPXyKQSfPxQHwwI8kBpdR1mrE1EVnGF2LGIiKwGCxwSxVmNDvM2HofeIGBiX3/MGx4idiSLo3SUYeWMfujq44IrpdWYsSYR2spasWMREVkFFjjU6gpKqzB73TGUVddhYLAHlt4fBonEtibxaykqJ0esf2QA/FRKpBWW45lNJ6A3CGLHIiKyeCxwqFVV1ugxZ/0x5JRUItizDVZMjYDcgX8Nb0atUuKL6f2gcJBi77krWLbznNiRiIgsHt8s1GoEQcDCzX/gZLYW7s6OWDOzP9yc5WLHsgo926vw/t97AQA+23sJP3PdKiKim2KBQ63mq8OX8VNyLmRSCVZMjUCQZxuxI1mV8eHt8djdwQCAF384iVM5WpETERFZLhY41CqSs0rw5rYzAIDY0aGIDG4nciLr9FJMKO7u4oWqWgMe+zIJRWXVYkciIrJILHDI7IrLa/DkV0mo1QsY3VON2fXLEVDzyaQSfDy5D4LaOSOnpBLPfpsMAzsdExFdhwUOmZXeIOCZTSeQq61CR882eP/vvThi6g6pnB3xxfR+UDpKsf9CIT7flyZ2JCIii8MCh8zq4/gL2H+hEEpHKT6b2hcuSkexI9mELj4ueP3eHgCAZTvPIenyVZETERFZFhY4ZDaH04rw8e4LAIB37gtDqNpV5ES2ZVL/ANzb2w96g4CnvzkBbQUnASQiuoYFDpmFtrIWC75NhiAAD/bz5xpTZiCRSPDOfT0R6GHsj7Nw8x8QBPbHISICWOCQmSz66RRytVXo0M4Zi+s/SqGW56J0xKcP94GjTIJfTmnw1ZFMsSMREVkEFjjU4n5KzjHNd/PRpHC0UTiIHcmm9fJ3w8ujQgEAS7adwYX8UpETERGJjwUOtajsqxV4dcspAMDT93RGn0B3kRPZh9lDOuLuLl6oqTPghe9Pok5vEDsSEZGoWOBQi9EbBCz47iRKq+vQN9AN84Z3EjuS3ZBIJHhvYi+4Kh1wMluLFb9fEjsSEZGoWOBQi1m5Pw2J6cVoI5fho0l94CDjX6/WpFYp8cZ4Y3+nf8dfwJlcnciJiIjE0ypvoOXLlyMoKAhKpRKRkZFITExs9Nhhw4ZBIpFct40dO9Z0zMyZM6/7+qhRo1rjVqgRl66U4cNd5wEAi8f1QGA7Z5ET2acJ4e0xsrsPavUCFnyXjJo6flRFRPbJ7AXOt99+iwULFmDx4sU4fvw4evfujZiYGBQUFNzw+M2bNyMvL8+0nTp1CjKZDA888ECD40aNGtXguG+++cbct0KNMBgEvPzDH6ipM+DuLl54IIJDwsUikUjw9n1h8Ggjx1lNKT6OvyB2JCIiUZi9wPnwww8xZ84czJo1C927d8eKFSvg7OyMNWvW3PB4Dw8PqNVq07Zr1y44OztfV+AoFIoGx7m7szOrWDYkZODY5atoI5fhnfvDuBSDyLxcFHhrQk8AwGe/X0JyVom4gYiIRGDWAqempgZJSUmIjo7+8xtKpYiOjkZCQkKTrrF69WpMnjwZbdq0abB/79698Pb2RteuXfHEE0+gqKio0WtUV1dDp9M12KhlZBVX4P1fzwEAFo7phvZuTiInIgAYE+aLcfWzHL/8wx+o5agqIrIzZi1wCgsLodfr4ePj02C/j48PNBrNLc9PTEzEqVOn8OijjzbYP2rUKGzYsAHx8fF477338Pvvv2P06NHQ6/U3vM7SpUuhUqlMW0BAwO3fFJkIgoDYzSmoqNFjQEcPTBkQKHYk+os3xvWARxs5zuWX4gsuyElEdsaih7msXr0aYWFhGDBgQIP9kydPxrhx4xAWFoYJEyZg27ZtOHr0KPbu3XvD68TGxkKr1Zq2rKysVkhv+74/lo0DFwuhcJDivYm9IJXyoylL4t5GjlfHdgNgXPT0clG5yImIiFqPWQscT09PyGQy5OfnN9ifn58PtVp903PLy8uxadMmzJ49+5bfJzg4GJ6enrh48eINv65QKODq6tpgoztTUFqFJdvPAAAW/K0LOnq2ucUZJIb7+rTH4JB2qK4z4JUtp7hWFRHZDbMWOHK5HBEREYiPjzftMxgMiI+PR1RU1E3P/f7771FdXY2pU6fe8vtkZ2ejqKgIvr6+d5yZmuad7akorapDL38VZg/pKHYcaoREIsHbE8KgcJDiwMVC/JicI3YkIqJWYfaPqBYsWICVK1di/fr1SE1NxRNPPIHy8nLMmjULADB9+nTExsZed97q1asxYcIEtGvXrsH+srIyvPjiizh8+DAyMjIQHx+P8ePHIyQkBDExMea+HQJw6FIhfkzOhUQCvD0hjBP6WbggzzZ4ekRnAMCSbakoLq8RORERkfmZfRXESZMm4cqVK1i0aBE0Gg3Cw8MRFxdn6nicmZkJqbThC/LcuXM4cOAAdu7ced31ZDIZ/vjjD6xfvx4lJSXw8/PDyJEjsWTJEigUCnPfjt2rqTPgtR+Na01NG9gBYf4qkRNRU8wZGoyfknNwPr8M7+xIxbIHeosdiYjIrCSCHX4or9PpoFKpoNVq2R+nmf6z9yLejzsHz7ZyxD8/DConR7EjURMlXS7GxM+M0zN8O3cgIoPb3eIMIiLL0pz3Nz9boCbLvlphmhn3lbHdWNxYmYgOHniofij/4q2nueI4Edk0FjjUZG/8fAZVtQZEdvTAhPD2Yseh2/BiTFeonBxxVlOKrxMzxY5DRGQ2LHCoSeJT87HrTD4cpBK8NaEnl2OwUh5t5HhhZBcAwD93nmeHYyKyWSxw6Jaq6/R4c5txzpvZQzuis4+LyInoTjwc2QHdfF2hrazFB/XLbBAR2RoWOHRLaw9m4HJRBbxdFHj6ns5ix6E7JJNK8Ma4HgCATUczkZKtFTkREVHLY4FDN3WltBqf7jbOEP3yqFC0UZh9ZgFqBQM6emB8uB8EAVi89RQMBrsbTElENo4FDt3UP3eeQ1l1HXr7q3BfH3YstiX/GNMNbeQyHM8swZYTnOGYiGwLCxxq1KkcLb49ZlyYdNG93bmYpo3xcVXiqfoZjt+LO4uKmjqRExERtRwWOHRDgiDgzW1nIAjAuN5+iOjgIXYkMoNZg4MQ6OGMgtJqfLEvTew4REQthgUO3VDcKQ0S04uhdJRi4ehQseOQmSgcZHh5lPHn+/nvacjXVYmciIioZbDAoetU1erx9o5UAMDcuzrBz81J5ERkTmPC1Ijo4I7KWj3+uZPDxonINrDAoeusP5SB7KuVULsq8fjdwWLHITOTSCR4ZWw3AMD3SdlIzdOJnIiI6M6xwKEGSipqsHyPcVj4CzFd4SznsHB70DfQHWN7+UIQgHd2pMIO1+AlIhvDAoca+M/eS9BV1SFU7cJh4XZm4ahQyGVS7L9QiL3nr4gdh4jojrDAIZPsqxVYdygDAPDy6FDIOCzcrgR4OGPm4CAAwDvbU7naOBFZNRY4ZPLhzvOoqTMgKrgdhnXxEjsOiWDe8BC4OTviQkEZfkjKFjsOEdFtY4FDAIAzuTpsSTbOZhs7JpSrhdsplZMjnqpfb+zf8RdQVasXORER0e1hgUMAgHfjzkIQgP/r5Yte/m5ixyERTYkMRHs3J+Rpq/BlwmWx4xAR3RYWOISDFwux7/wVOMokeDGmq9hxSGRKRxmeiTa24izfexG6qlqRExERNR8LHDsnCALe/eUsAGBKZAd0aNdG5ERkCe7v0x4h3m1RUlGLlVzCgYisEAscO/fraQ1ScrRoI5fhqXtCxI5DFsJBJsULI42teav2p+NKabXIiYiImocFjh3TGwT8c+d5AMDsIR3Rrq1C5ERkSWJ6+KB3gBsqa/X4dPcFseMQETULCxw7tvVkDi4UlMFV6YDZQ7kkAzUkkUjwcn2frK8TM5FZVCFyIiKipmOBY6dq9QZ89Jvxt/LH7u4ElZOjyInIEg0K8cTQzp6o1Qv412/nxY5DRNRkLHDs1A9J2bhcVAHPtnLMqp+9luhGro2s+zE5BxcLSkVOQ0TUNCxw7FBVrR4fxxtbb54cFsIFNemmevm7YWR3HwgCTK1+RESWjgWOHfomMRN52ir4qpR4ODJQ7DhkBZ6N7gIA2J6Sh7ManchpiIhujQWOnamoqcPyPRcBAE/d0xlKR5nIicgadPdzxZgwNQQB+DdbcYjICrDAsTMbEi6jsKwGHdo544F+/mLHISvyzIgukEiAX05pcDpXK3YcIqKbYoFjR8qr6/BF/ay0T9/TGY4y/vip6bqqXfB/vfwAsC8OEVk+vuHsyFeHL6O4vAZB7ZwxPtxP7DhkhZ4Z0RlSCbDrTD5SstmKQ0SWq1UKnOXLlyMoKAhKpRKRkZFITExs9Nh169ZBIpE02JRKZYNjBEHAokWL4OvrCycnJ0RHR+PCBf5GeTMVNX+23sy/pzMc2HpDtyHEuy3Gh7cHAM6LQ0QWzexvuW+//RYLFizA4sWLcfz4cfTu3RsxMTEoKCho9BxXV1fk5eWZtsuXLzf4+vvvv4+PP/4YK1aswJEjR9CmTRvExMSgqqrK3Ldjtb46fBlF5ca+NxPYekN34OkRnSGTSrD7bAFOZF4VOw4R0Q2ZvcD58MMPMWfOHMyaNQvdu3fHihUr4OzsjDVr1jR6jkQigVqtNm0+Pj6mrwmCgI8++givvvoqxo8fj169emHDhg3Izc3Fjz/+aO7bsUoVNXX4/Pf61pvhIWy9oTvS0bMN7utjbMW5Np8SEZGlMeubrqamBklJSYiOjv7zG0qliI6ORkJCQqPnlZWVoUOHDggICMD48eNx+vRp09fS09Oh0WgaXFOlUiEyMrLRa1ZXV0On0zXY7MnGw5koKq9BoIez6cVEdCfmDQ+BVALsOXeFfXGIyCKZtcApLCyEXq9v0AIDAD4+PtBoNDc8p2vXrlizZg1++uknfPXVVzAYDBg0aBCys7MBwHRec665dOlSqFQq0xYQEHCnt2Y1Kmv0+HzfJQBsvaGW09GzDcb1Nn7U+QlXGiciC2Rxb7uoqChMnz4d4eHhuPvuu7F582Z4eXnh888/v+1rxsbGQqvVmrasrKwWTGzZNh4xznsT4OGE+/qy9YZazvx7QiCRADvP5CM1z75aRYnI8pm1wPH09IRMJkN+fn6D/fn5+VCr1U26hqOjI/r06YOLF42z7147rznXVCgUcHV1bbDZg6paPVb8pe8N572hlhTi7YIxPX0BAJ/Wz45NRGQpzPrGk8vliIiIQHx8vGmfwWBAfHw8oqKimnQNvV6PlJQU+Poa/0fasWNHqNXqBtfU6XQ4cuRIk69pL749moXCsmq0d3PC/X05azG1vPn3hAAAdqTkcaVxIrIoZv+VfsGCBVi5ciXWr1+P1NRUPPHEEygvL8esWbMAANOnT0dsbKzp+DfffBM7d+5EWloajh8/jqlTp+Ly5ct49NFHARhHWD377LN46623sHXrVqSkpGD69Onw8/PDhAkTzH07VqOmzoDPfzf2vXn87mC23pBZdPN1Na00vnzPJbHjEBGZOJj7G0yaNAlXrlzBokWLoNFoEB4ejri4OFMn4czMTEilf758r169ijlz5kCj0cDd3R0RERE4dOgQunfvbjrmpZdeQnl5OebOnYuSkhIMGTIEcXFx100IaM9+PJGDXG0VvFwUeKCf/XSqptb31D2dsfNMPn5KzsEzIzojyLON2JGIiCARBEEQO0Rr0+l0UKlU0Gq1Ntkfp05vQPSHvyOjqAKvjOmGOXcFix2JbNystYnYc+4KHuznj/f/3lvsOERko5rz/ubnFjZoe0oeMooq4ObsiIcjA8WOQ3bgqRGdAQCbj+cgp6RS5DRERCxwbI7BIOA/9X0hHhncEW0UZv8Ukgh9A90xqFM71BkErKxf84yISEwscGzMb6n5OJdfirYKB8yIChI7DtmRJ4cZR1RtOpqJorJqkdMQkb1jgWNDBEHA8vr5SKZFdYDK2VHkRGRPBoe0Qy9/FapqDVh7MEPsOERk51jg2JADFwtxMlsLpaMUs4d0FDsO2RmJRIInh3UCAKxPyEBpVa3IiYjInrHAsSGf7TX2vZncPxCebRUipyF7NLK7Gp282qC0qg4bj2SKHYeI7BgLHBtxMqsEhy4VwUEq4bBwEo1UKsHjdxtbcVbtT0dVrV7kRERkr1jg2IgV9bMWjwv3Q3s3J5HTkD0bH94efiolCsuq8X1StthxiMhOscCxAWlXyhB3WgMApt+eicQid5CaWhG/2HcJdXqDyImIyB6xwLEBX+xLgyAA0d280cXHRew4RJjcPxAebeTIKq7Etj/yxI5DRHaIBY6Vy9dVYfPxHABsvSHL4SSXYdagIADGj0/tcEUYIhIZCxwrt+ZAOmr0BvTr4I5+QR5ixyEymRbVAc5yGc5qSrHvQqHYcYjIzrDAsWLaylrTUFy23pClcXOWY3J/41pon9d3giciai0scKzYxiOXUVZdhy4+bXFPqLfYcYiu88iQIMikEhy6VISUbK3YcYjIjrDAsVJVtXqsOZABAHjsrk6QSiXiBiK6AX93Z9zbyxcA8Pk+tuIQUethgWOlfjyRg8KyaviqlBgX7id2HKJGzb3L+PHpjpQ8ZBZViJyGiOwFCxwrZDAI+GJ/GgBg9pCOcJTxx0iWq7ufK+7q4gWDAKw6kCZ2HCKyE3wzWqHdZwuQdqUcLgoHTOofIHYcolt6vH7iv++OZaGorFrkNERkD1jgWKFrrTcPRwbCRekochqiW4vq1A5h7VWoqjVgQ8JlseMQkR1ggWNlkrNKkJheDAepBLMGdxQ7DlGTSCQSzK1vxdmQkIHKGi7CSUTmxQLHyqzcZ2y9GRfuB7VKKXIaoqYb3VONAA8nXK2oxQ/HuQgnEZkXCxwrkllUgV9OGdf1ufbbMJG1cJBJ8Uh9q+OaA+kwGLh8AxGZDwscK7L6QBoMAnBXFy+Eql3FjkPUbA/2C4Cr0gHpheX4LTVf7DhEZMNY4FiJq+U1+O6YsVl/7lC23pB1aqNwwMORHQAAq/ani5yGiGwZCxwrsfHIZVTW6tHd1xWDQ9qJHYfots0cFAQHqQSJGcU4mVUidhwislEscKxAdZ0e6w4Zh9bOuasjJBIuy0DWS61SYlxv4+zbK/dz4j8iMg8WOFZga3IuCsuqoXZVYmwYl2Ug6/do/cesO1LykFXM5RuIqOWxwLFwgiBg9QFjX4UZg4Igd+CPjKxfdz/jR60GAVh7MEPsOERkg/i2tHCHLhXhrKYUTo4yPDwgUOw4RC3mWivOt0czoa2sFTkNEdkaFjgWblV9H4UH+/lD5cxlGch2DOvihc7ebVFeo8emxEyx4xCRjWGBY8EuFpRiz7krkEjAZRnI5kgkEjw61Pj3ev2hDNTpDSInIiJb0ioFzvLlyxEUFASlUonIyEgkJiY2euzKlSsxdOhQuLu7w93dHdHR0dcdP3PmTEgkkgbbqFGjzH0brW71gQwAwMjuPgjybCNuGCIzGB/eHp5t5cjVVuGXUxqx4xCRDTF7gfPtt99iwYIFWLx4MY4fP47evXsjJiYGBQUFNzx+7969eOihh7Bnzx4kJCQgICAAI0eORE5OToPjRo0ahby8PNP2zTffmPtWWlVRWTU216/X8ygn9iMbpXSUYUr9xH/XOtMTEbUEsxc4H374IebMmYNZs2ahe/fuWLFiBZydnbFmzZobHr9x40Y8+eSTCA8PR2hoKFatWgWDwYD4+PgGxykUCqjVatPm7u5u7ltpVRuPZKK6zoDe/ir062Bb90b0V1MHdoBcJkVyVgmSLl8VOw4R2QizFjg1NTVISkpCdHT0n99QKkV0dDQSEhKadI2KigrU1tbCw8Ojwf69e/fC29sbXbt2xRNPPIGioqJGr1FdXQ2dTtdgs2TVdXpsSDBO7Dd7aDAn9iOb5uWiwPhw4/xOa9iKQ0QtxKwFTmFhIfR6PXx8fBrs9/HxgUbTtM/bX375Zfj5+TUokkaNGoUNGzYgPj4e7733Hn7//XeMHj0aer3+htdYunQpVCqVaQsICLj9m2oFP5/MQ2FZNXxVSozuqRY7DpHZza7vbPzLqTxkX+XEf0R05yx6FNW7776LTZs2YcuWLVAqlab9kydPxrhx4xAWFoYJEyZg27ZtOHr0KPbu3XvD68TGxkKr1Zq2rKysVrqD5vvfif0cZRb9IyJqEaHqPyf+W38oQ+w4RGQDzPr29PT0hEwmQ35+foP9+fn5UKtv3jKxbNkyvPvuu9i5cyd69ep102ODg4Ph6emJixcv3vDrCoUCrq6uDTZLdTitGKl5Ojg5yjC5v2W3NBG1pNlDjK04mxKzUFZdJ3IaIrJ2Zi1w5HI5IiIiGnQQvtZhOCoqqtHz3n//fSxZsgRxcXHo16/fLb9PdnY2ioqK4Ovr2yK5xbTmoLH1ZmJEe7g5y0VOQ9R6hnXxRrBXG5RW1+H7Y5bbykpE1sHsn38sWLAAK1euxPr165GamoonnngC5eXlmDVrFgBg+vTpiI2NNR3/3nvv4bXXXsOaNWsQFBQEjUYDjUaDsrIyAEBZWRlefPFFHD58GBkZGYiPj8f48eMREhKCmJgYc9+OWV0uKsdvqcbWLk7sR/ZGKpWY/t6vPZgBvUEQORERWTOzFziTJk3CsmXLsGjRIoSHhyM5ORlxcXGmjseZmZnIy8szHf/ZZ5+hpqYGf//73+Hr62vali1bBgCQyWT4448/MG7cOHTp0gWzZ89GREQE9u/fD4VCYe7bMat1hzIgCMDwrl7o5NVW7DhErW5i3/ZQOTkis7gC8an5tz6BiKgREkEQ7O7XJJ1OB5VKBa1WazH9cUqrahG1dDfKquvw5ewBGNrZS+xIRKJ495ezWPH7JUQFt8M3cweKHYeILEhz3t8comMhvjuWjbLqOnT2boshIZ5ixyESzfSoDpBJJUhIK0JqnmXPWUVElosFjgXQGwSsO2TsXPzIkI6c2I/smp+bE0bVz/+09iAn/iOi28MCxwLsOpOPrOJKuDs74r4+7cWOQyS6RwYHAQB+TM5FUVm1uGGIyCqxwLEA14aGPxwZCKWjTOQ0ROLrG+iO3v4q1NQZ8PWRTLHjEJEVYoEjstO5WiSmF8NBKsG0gUFixyGyCBLJn0PGvzx8GTV1BpETEZG1YYEjsnUHMwAAo8N8oVYpb34wkR0ZE+YLbxcFCkqrsSMl79YnEBH9BQscERWVVeOnk7kAgFn1fQ6IyEjuIMW0gR0AGD/GtcMZLYjoDrDAEdHXRzJRU2dAb38V+gS4iR2HyOI8HBkIuYMUf2RrcTzzqthxiMiKsMARSa3egC8PXwZgXJaBQ8OJrteurQLje/sBANbUf5xLRNQULHBEsiMlDwWl1fB2UWBMmPUvEkpkLtc6G8ed0iBPWylyGiKyFixwRLK2/rfRqQM7QO7AHwNRY7r7uSKyowf0BgFf1bd6EhHdCt+sIjiReRXJWSWQy6R4aECg2HGILN61TvhfH8lEVa1e3DBEZBVY4Ihg3aEMAMC9vf3g5WLdK6ATtYbobj5o7+aEqxW12JqcK3YcIrICLHBaWb6uCtv/MM7pwaHhRE3jIJNiepRxyPjaQxkcMk5Et8QCp5VtPHwZdQYB/YPc0bO9Suw4RFZjUv8AKB2lSM3TITG9WOw4RGThWOC0ouo6PTbWr6tzbWQIETWNm7Mc9/XxB/Dnx7xERI1hgdOKtp3MQ1F5DXxVSozs7iN2HCKrM3NQEADg19MaZF+tEDcMEVk0FjitRBAE02+d06I6wEHGR0/UXF3VLhgc0g4GAaaJMomIboRv2VZyPPMqUnK0kDtIMbk/h4YT3a6Zg4wf725KzEJlDYeME9GNscBpJdcm9psQ7gePNnJxwxBZsXtCvRHo4QxtZS22nMgROw4RWSgWOK0gT1uJX05pAAAz6vsQENHtkUklpiHj6zlknIgawQKnFWw8nAm9QcCAjh7o4ceh4UR36oF+AXBylOFcfikOp3HIOBFdjwWOmVXV6vFNYv3QcLbeELUIlZMj7u/bHgCw7lC6yGmIyBKxwDGzbX8Yh4b7qZT4G4eGE7WYa0PGd53J55BxIroOCxwzMg4NN/52OZVDw4laVGcfFwwJ8eSQcSK6Ib5xzSjp8lWcytFBwaHhRGZxrdP+t0c5ZJyIGmKBY0bXJvYbz6HhRGZxT6g3AjycUFJRi5+SOWSciP7EAsdMNNoqDg0nMjOZVILpA4MAGH+h4JBxIrqGBY6ZbDxy2Tg0PIhDw4nM6cH6IeNnNaU4wlXGiageCxwzqK7T4+v6VcPZekNkXirnvwwZr58xnIiIBY4ZNFg1vAeHhhOZ27VfJHae0SCnpFLcMERkEVqlwFm+fDmCgoKgVCoRGRmJxMTEmx7//fffIzQ0FEqlEmFhYdixY0eDrwuCgEWLFsHX1xdOTk6Ijo7GhQsXzHkLTSYIAtYnZAAApg7sAEcODScyuy4+f1llPIFDxomoFQqcb7/9FgsWLMDixYtx/Phx9O7dGzExMSgoKLjh8YcOHcJDDz2E2bNn48SJE5gwYQImTJiAU6dOmY55//338fHHH2PFihU4cuQI2rRpg5iYGFRVVZn7dm7pRFYJ/si+tmp4gNhxiOzGjKggAMCmo5moquWQcSJ7JxHMPOwgMjIS/fv3x6effgoAMBgMCAgIwFNPPYWFCxded/ykSZNQXl6Obdu2mfYNHDgQ4eHhWLFiBQRBgJ+fH55//nm88MILAACtVgsfHx+sW7cOkydPvmUmnU4HlUoFrVYLV1fXFrpTo6e/OYGtJ3Px9wh/LHugd4tem4gapzcIuPuDPci+Won3J/bCg/wFg8jmNOf9bdYWnJqaGiQlJSE6OvrPbyiVIjo6GgkJCTc8JyEhocHxABATE2M6Pj09HRqNpsExKpUKkZGRjV6zuroaOp2uwWYOBboq7EjJA/DnNPJE1Dr+usr4Wg4ZJxLN+fxSPPblMRxOKxI1h1kLnMLCQuj1evj4NOxo6+PjA41Gc8NzNBrNTY+/9t/mXHPp0qVQqVSmLSDAPL/ZbTySiTqDgH4d3NGzPYeGE7W2B/sFQOkoRWqeDkczroodh8gurTuUgV9P52N9/WS3YrGLHrCxsbHQarWmLSsryyzfZ2Jff8wZ2hFz7wo2y/WJ6ObcnOW4r49xyLjY/3MlskfailpsOW6cVVzsaVLMWuB4enpCJpMhPz+/wf78/Hyo1eobnqNWq296/LX/NueaCoUCrq6uDTZzCGznjFfGdsfIHjfOQUTmd+1/qnGnNcjlkHGiVvXdsSxU1uoRqnZBZEcPUbOYtcCRy+WIiIhAfHy8aZ/BYEB8fDyioqJueE5UVFSD4wFg165dpuM7duwItVrd4BidTocjR440ek0ish+halcMDPaA3iBg4xEOGSdqLXqDgA2HMwAY+6FKJBJR85j9I6oFCxZg5cqVWL9+PVJTU/HEE0+gvLwcs2bNAgBMnz4dsbGxpuOfeeYZxMXF4Z///CfOnj2L119/HceOHcP8+fMBABKJBM8++yzeeustbN26FSkpKZg+fTr8/PwwYcIEc98OEVmBa538v0nM4pBxolay+2wBsooroXJyxPjw9mLHgYO5v8GkSZNw5coVLFq0CBqNBuHh4YiLizN1Es7MzIRU+medNWjQIHz99dd49dVX8Y9//AOdO3fGjz/+iJ49e5qOeemll1BeXo65c+eipKQEQ4YMQVxcHJRKpblvh4isQHQ3H7R3c0JOSSV+PpmLB/pxyDiRuV3r9zZ5QACc5DJxw6AV5sGxROacB4eILMNney/hvbiz6OHnim1PDRG9uZzIll3IL8Xf/rUPUgmw76Xh8Hd3Nsv3sZh5cIiIxDK5fwAUDlKcztUh6TKHjBOZ07Uliv7W3cdsxU1zscAhIpvk3kaOCfX9ANZyyDiR2Wgra/HfJMsYGv5XLHCIyGaZhoyf0iBPyyHjRObwff3Q8K4+LogKbid2HBMWOERks7r7uWJAx/oh44czxY5DZHP0BgFfHjZOxzDDAoaG/xULHCKyabNMQ8a5yjhRS9t7rgCXiyrgqnTAhD5+YsdpgAUOEdm0v3X3gZ9KiaLyGmz7I0/sOEQ2ZV19/7aHBgTCWW72mWeahQUOEdk0B5kUU+tXGV/PVcaJWsyF/FLsv1AIqQSYOrCD2HGuwwKHiGze5P6BUDhIkZKjxfFMDhknagnXhoZHd/NBgIdlDA3/KxY4RGTzPNrIMT7c2D9g7cEMccMQ2QBtZS02168aPnNwkLhhGsECh4jswl+HjGu0VeKGIbJy3x/LQkWN5Q0N/ysWOERkF3r4qTCgowfqDAK+OsxVxolul94gYEOC8d/QzMGWNTT8r1jgEJHduDZk/GsOGSe6bXvOFiCzuAIqJ0fTbOGWiAUOEdmNv3U3rjJeXF6Dn0/mih2HyCpdGxo+ub9lrBreGBY4RGQ3HGRSTKsfMr72IIeMEzXXhfxSHLhouUPD/4oFDhHZlcn9A6B0lOJMng5HMzhknKg5rrXe/K27ZQ4N/ysWOERkV9yc5bivj7HfwLpD6SKnIbIe2oo/h4Zb0qrhjWGBQ0R259r/nH89nY+cEq4yTtQU3x7LRGWtHqFqyx0a/lcscIjI7oSqXTGoUzvjSsgJHDJOdCt1egPWHzL+W5llwUPD/4oFDhHZpZn1rTibjmaisoZDxolu5rdUY2unu7Mjxlvw0PC/YoFDRHZpRDcf+Ls7oaSiFj8m54gdh8iiralf4uThyEAoHS13aPhfscAhIrskk0owIyoIALD2YDqHjBM14nSuFonpxZBJJRY/NPyvWOAQkd16sH8AnOUynM8vw6FLRWLHIbJI6+pbb0b3VMNX5SRumGZggUNEdkvl5Ii/R/gDMLbiEFFDRWXV+Kl+1u9ZgzuKnKZ5WOAQkV27NmQ8/mwBMgrLxQ1DZGG+ScxETZ0Bvf1V6BvoJnacZmGBQ0R2rZNXWwzv6gVB+HOWViICavUGfHnY8lcNbwwLHCKye9ea3n9IykZpVa3IaYgsw46UPOTrquHlosDYMD+x4zQbCxwisntDO3sixLstyqrr8N2xbLHjEIlOEASsOWDslzZtYAfIHayvXLC+xERELUwikWDW4CAAwPpDGdAbOGSc7NvxzKs4ma2F3EGKhyMDxY5zW1jgEBEBuL+PP1ROjsgsrsDuswVixyES1ZoDGQCACeF+8GyrEDfMbWKBQ0QEwEkuw0MDjL+pcsg42bPsqxX45VQeAOCRIdY1NPyvWOAQEdWbHtUBMqkEhy4V4UyuTuw4RKL4MuEyDAIwOKQdQtWuYse5bWYtcIqLizFlyhS4urrCzc0Ns2fPRllZ2U2Pf+qpp9C1a1c4OTkhMDAQTz/9NLRabYPjJBLJddumTZvMeStEZAf83JwwqqcaALCGrThkh8qr6/BNYiYA4BErm9jvf5m1wJkyZQpOnz6NXbt2Ydu2bdi3bx/mzp3b6PG5ubnIzc3FsmXLcOrUKaxbtw5xcXGYPXv2dceuXbsWeXl5pm3ChAlmvBMishez65vktybnoqC0SuQ0RK3rv8ezoauqQ1A7Zwzv6i12nDviYK4Lp6amIi4uDkePHkW/fv0AAJ988gnGjBmDZcuWwc/v+jH1PXv2xH//+1/Tnzt16oS3334bU6dORV1dHRwc/ozr5uYGtVptrvhEZKf6BrqjT6AbTmSW4KvDmVjwty5iRyJqFQaDgLX1607NGtwRUql1Tez3v8zWgpOQkAA3NzdTcQMA0dHRkEqlOHLkSJOvo9Vq4erq2qC4AYB58+bB09MTAwYMwJo1a266EnB1dTV0Ol2DjYioMddacTYevoyqWr3IaYhax97zBUgvLIeL0sG0Rps1M1uBo9Fo4O3dsHnLwcEBHh4e0Gg0TbpGYWEhlixZct3HWm+++Sa+++477Nq1CxMnTsSTTz6JTz75pNHrLF26FCqVyrQFBAQ0/4aIyG6M6qFGezcnFJXXYGtyrthxiFrF6vqJ/Sb3D0Abhdk+4Gk1zS5wFi5ceMNOvn/dzp49e8fBdDodxo4di+7du+P1119v8LXXXnsNgwcPRp8+ffDyyy/jpZdewgcffNDotWJjY6HVak1bVlbWHecjItvlIJNixqAOAIydjW/WQkxkC1LzdDh4sQhSyZ8L0Fq7Zpdozz//PGbOnHnTY4KDg6FWq1FQ0HCyrLq6OhQXF9+y70xpaSlGjRoFFxcXbNmyBY6Ojjc9PjIyEkuWLEF1dTUUiusnJFIoFDfcT0TUmEn9A/HRbxdwVlOKgxeLMKSzp9iRiMzmWuvN6DBf+Ls7i5ymZTS7wPHy8oKXl9ctj4uKikJJSQmSkpIQEREBANi9ezcMBgMiIyMbPU+n0yEmJgYKhQJbt26FUqm85fdKTk6Gu7s7ixgiajEqJ0c82C8A6w5lYPWBNBY4ZLMKdFX4KTkHAPCoFU/s97/M1genW7duGDVqFObMmYPExEQcPHgQ8+fPx+TJk00jqHJychAaGorExEQAxuJm5MiRKC8vx+rVq6HT6aDRaKDRaKDXGzv6/fzzz1i1ahVOnTqFixcv4rPPPsM777yDp556yly3QkR2atbgIEgkwJ5zV3CxoPE5vIis2ZeHL6NWLyCigzv6BLqLHafFmLUX0caNGzF//nyMGDECUqkUEydOxMcff2z6em1tLc6dO4eKigoAwPHjx00jrEJCQhpcKz09HUFBQXB0dMTy5cvx3HPPQRAEhISE4MMPP8ScOXPMeStEZIc6tGuD6G4+2HUmH2sOpuOd+8LEjkTUoipr9Pjq8GUAttV6AwASwQ57z+l0OqhUKtMQdCKixhxOK8LkLw5D4SDFoYX3oJ2VLjxIdCMbj1zGK1tOIcDDCXtfGA6Zhc9905z3N9eiIiK6iciOHghrr0J1nQFfHc4UOw5RizEYBFPn4lmDOlp8cdNcLHCIiG5CIpHg0aHGpvsvD2dw4j+yGXvPFyDtSjlcFA54sL/tzQ/HAoeI6BbGhPmivZsTCstq8OOJHLHjELWIVfuNrTcPRQairQ1M7Pe/WOAQEd2Co0yKWYODAAAr96fBYLC7rotkY07nanHoUhFkUonNTOz3v1jgEBE1waT+AXBROODSlXLsPV9w6xOILNi11pvRPY3LktgiFjhERE3gonTE5AHGfgor96WLnIbo9uWWVOLnk8Y11h67q5PIacyHBQ4RURPNHGwcaZKQVoRTOVqx4xDdlrUH01FnEBAV3A5h/iqx45gNCxwioiZq7+aEsWG+AIBV+9NETkPUfLqqWnyTaFxweu5dwSKnMS8WOEREzTBnqPGl8PMfecgtqRQ5DVHzfHMkE2XVdejs3RZ3d7n1upLWjAUOEVEzhPmrMDDYA3qDgDUH2BeHrEdNnQFrD2YAAObcFQypjU3s979Y4BARNdNjdxs7Zn6TmAltZa3IaYia5ueTudDoquDtosD4cD+x45gdCxwiomYa1sULXX1cUF6jx8Yjl8WOQ3RLgiBgZX2/sZmDg6BwkImcyPxY4BARNZNEIsFjdxv74qw9yOUbyPLtu1CIs5pSOMtlmDKgg9hxWgULHCKi23Bvbz/4qpS4UlrN5RvI4n2x7xIA44SVKmdHkdO0DhY4RES3wVEmxewhxkU4v+DyDWTBUrK1OHjRuCzDI4M7ih2n1bDAISK6TZMHBMJF6YC0K+X4LTVf7DhEN7Tid2Przb29fBHg4SxymtbDAoeI6Da1VThg2kBjf4bP93HiP7I86YXl+OVUHgDg8WG2uyzDjbDAISK6AzMHB0EukyLp8lUcyygWOw5RA1/sS4NBAIZ39UKo2lXsOK2KBQ4R0R3wdlHi/r7tAQArfmcrDlmOAl0V/puUDQB4YliIyGlaHwscIqI7NOeuYEgkwG+p+TinKRU7DhEAYM3BDNToDegb6Ib+Qe5ix2l1LHCIiO5QJ6+2GN1TDeDPDp1EYtJV1WLjYeMklE8MC4FEYtvLMtwICxwiohbwZP1HAFtP5iKruELkNGTvvj6SidL6RTVHhHqLHUcULHCIiFpAz/YqDO3sCb1BwBccUUUiqqrVY3X9QrCP3d3J5hfVbAwLHCKiFnKtFefbY1koKK0SOQ3Zq83Hc3CltBq+KiXG9bb9RTUbwwKHiKiFDAz2QJ9AN9TUGbDmQIbYccgO1ekNpn5gc4YGQ+5gv695+71zIqIWJpFITK04Xx2+DG1lrciJyN78/EcuMosr0K6NHA8NCBQ7jqhY4BARtaARod7o4tMWZdV1+Kp+FAtRazAYBPxnj7H15pEhHeEkl4mcSFwscIiIWpBUKsET9VPirzmQjsoavciJyF7sPKPBhYIyuCgdMC2qg9hxRMcCh4iohd3byw/+7k4oKq/BpqOZYschOyAIAj7dcxEAMHNQEFyVjiInEh8LHCKiFuYgk+Lxu42tOJ//nobqOrbikHn9fv4KTuXo4OQow6zBHcWOYxFY4BARmcED/fyhdlVCo6vC98eyxY5DNm55fevNlMhAeLSRi5zGMpi1wCkuLsaUKVPg6uoKNzc3zJ49G2VlZTc9Z9iwYZBIJA22xx9/vMExmZmZGDt2LJydneHt7Y0XX3wRdXV15rwVIqJmUTjI8PjdwQCAz/ZeQk2dQeREZKuOpBXhaMZVyGVSzLkrWOw4FsOsBc6UKVNw+vRp7Nq1C9u2bcO+ffswd+7cW543Z84c5OXlmbb333/f9DW9Xo+xY8eipqYGhw4dwvr167Fu3TosWrTInLdCRNRskwcEwrOtAjklldhygq04ZB7X+t480M8fPq5KkdNYDrMVOKmpqYiLi8OqVasQGRmJIUOG4JNPPsGmTZuQm5t703OdnZ2hVqtNm6urq+lrO3fuxJkzZ/DVV18hPDwco0ePxpIlS7B8+XLU1NSY63aIiJpN6SjDY/W/US/fcwl1erbiUMs6nnkV+y8UQiaVmPp9kZHZCpyEhAS4ubmhX79+pn3R0dGQSqU4cuTITc/duHEjPD090bNnT8TGxqKi4s+F6xISEhAWFgYfHx/TvpiYGOh0Opw+ffqG16uuroZOp2uwERG1hikDjX0iMosrsPXkzX+5I2quf/92AQAwsW97BHg4i5zGspitwNFoNPD2briCqYODAzw8PKDRaBo97+GHH8ZXX32FPXv2IDY2Fl9++SWmTp3a4Lp/LW4AmP7c2HWXLl0KlUpl2gICAm73toiImsVZ7oBHhxpHtXy65yL0BkHkRGQrTmRexe/nr0AmlWD+8M5ix7E4zS5wFi5ceF0n4P/dzp49e9uB5s6di5iYGISFhWHKlCnYsGEDtmzZgkuXLt32NWNjY6HVak1bVlbWbV+LiKi5pkcFQeXkiLQr5diRkid2HLIR/443tt7c36c9Atux9eZ/OTT3hOeffx4zZ8686THBwcFQq9UoKChosL+urg7FxcVQq9VN/n6RkZEAgIsXL6JTp05Qq9VITExscEx+fj4ANHpdhUIBhULR5O9JRNSS2ioc8MjgjvjXb+fxye4LGBvmC6lUInYssmLJWSXYe66+9eaeELHjWKRmFzheXl7w8vK65XFRUVEoKSlBUlISIiIiAAC7d++GwWAwFS1NkZycDADw9fU1Xfftt99GQUGB6SOwXbt2wdXVFd27d2/m3RARtY6Zg4Ow6kAazueXYXtKHu7t7Sd2JLJi//7tPADgvj7t0aFdG5HTWCaz9cHp1q0bRo0ahTlz5iAxMREHDx7E/PnzMXnyZPj5Gf9h5+TkIDQ01NQic+nSJSxZsgRJSUnIyMjA1q1bMX36dNx1113o1asXAGDkyJHo3r07pk2bhpMnT+LXX3/Fq6++innz5rGVhogslsrJEXOGGkdUffTbefbFodt2MqsEe6613gxn601jzDoPzsaNGxEaGooRI0ZgzJgxGDJkCL744gvT12tra3Hu3DnTKCm5XI7ffvsNI0eORGhoKJ5//nlMnDgRP//8s+kcmUyGbdu2QSaTISoqClOnTsX06dPx5ptvmvNWiIju2KzBxr44l66U42eOqKLbdK3vzfhwPwR5svWmMRJBEOzu1widTgeVSgWtVttgjh0iInNbvuciPvj1HDp6tsGu5+6Cg4wr5lDT/ZFdgnGfHoRUAsQ/Pwwd7azAac77m/+yiIha0YxBQfBoI0d6YTl+TGYrDjXPh7uMfW8mhLe3u+KmuVjgEBG1orYKB9Psxh/HX0AtZzemJjqWUWwaOfVMNOe9uRUWOERErWxaVAd4tjXObrz5ONeoolsTBAEf/HoOAPBgP3+OnGoCFjhERK3MWe5gWjfo4/iLXGmcbungxSIcSS+GXCbFU/ew9aYpWOAQEYlg6sAO8HIxrjT+7THOrk6NEwQBH+w0tt5MGRgIPzcnkRNZBxY4REQiUDrKTHOYfBx/AZU1epETkaX6LbUAJ7NK4OQow5PDOO9NU7HAISISyUMDAuHv7oQrpdVYeyhd7DhkgQwGAf+sb72ZOTgIXi6c0LapWOAQEYlE7iDF8yO7AABW7L0EbUWtyInI0mxLycNZTSlc/jL6jpqGBQ4RkYjG9W6PULULdFV1WLHvkthxyILU6Q34qH7emzl3BcPNWS5yIuvCAoeISEQyqQQvjOwKAFh7MB35uiqRE5Gl+PZYFtIKy+HRRo5HhnQUO47VYYFDRCSyEd28EdHBHVW1Bnxcv84Q2beKmjp89Jvx78JT94SgrcJB5ETWhwUOEZHIJBIJXh4VCgDYdDQL6YXlIicisa3an44rpdUI9HDGlMgOYsexSixwiIgswICOHhje1Qv6v4yaIftUWFaNz3839sd6IaYr5A58Vd8OPjUiIgvxQoyxL862P/KQnFUibhgSzSfxF1Beo0cvfxX+L8xX7DhWiwUOEZGF6OGnwv192gMA3tmeCkEQRE5ErS2jsBwbj2QCABaODoVUKhE5kfVigUNEZEFeiOkKhYMUiRnF+PV0vthxqJV9sPMc6gwChnX1wqBOnmLHsWoscIiILIifmxPmDDVO6PbuL6lciNOOnMwqwfY/8iCRwNTpnG4fCxwiIgvz+LBO8GwrR0ZRBTYeuSx2HGoFgiDgre1nAAD39/FHN19XkRNZPxY4REQWpq3CAc/9zbiEw7/jL3AJBzuwPSUPRzOuwslRhhdiuogdxyawwCEiskCT+gWgs3dblFTU4tM9nPzPllXV6rF0x1kAwON3d4KvyknkRLaBBQ4RkQVykEnxj7HdAADrD11GZlGFyInIXFbuS0NOSSX8VErM5YKaLYYFDhGRhRrWxQtDO3uiRm/A2zvOiB2HzECjrcJ/9hon9Xt5dCic5DKRE9kOFjhERBZKIpHg1bHdIZNK8OvpfOy/cEXsSNTC3v/1LCpr9Yjo4I5xvf3EjmNTWOAQEVmwrmoXTI8yrkX0+tbTHDZuQ5KzSrD5eA4AYNH/dYdEwkn9WhILHCIiC/dsdBe0ayPHpSvlWH8oQ+w41AIEQcCbP58GANzftz16B7iJG8gGscAhIrJwKidH08Rv/46/gAJdlciJ6E7993gOjmeWwFku46R+ZsICh4jICvw9wh+9A9xQVl2Hd+POih2H7kBJRQ2W7kgFADw9ojN8XJUiJ7JNLHCIiKyAVCrBG+N6AAA2H89B0uVikRPR7frg13MoKq9BZ++2eGRwR7Hj2CwWOEREViI8wA0P9vMHACzeehp6A1cbtzbJWSX4OtG4WviSCT0hd+Br2Fz4ZImIrMhLo0LhonTAqRwdOxxbGb1BwKs/pkAQgPv7tMfA4HZiR7JpLHCIiKyIZ1sFFo42dkr9585zyC2pFDkRNdXXRy7jVI4OLkoHxI7pJnYcm2fWAqe4uBhTpkyBq6sr3NzcMHv2bJSVlTV6fEZGBiQSyQ2377//3nTcjb6+adMmc94KEZHFeKh/IPp1cEd5jR6LfjoNQeBHVZbuSmk13v/1HADgpZiu8HJRiJzI9pm1wJkyZQpOnz6NXbt2Ydu2bdi3bx/mzp3b6PEBAQHIy8trsL3xxhto27YtRo8e3eDYtWvXNjhuwoQJ5rwVIiKLIZVK8M79YXCUSfBbaj5+Pa0ROxLdwlvbz6C0qg5h7VV4OLKD2HHsgoO5Lpyamoq4uDgcPXoU/fr1AwB88sknGDNmDJYtWwY/v+unpJbJZFCr1Q32bdmyBQ8++CDatm3bYL+bm9t1xxIR2YsuPi547K5O+HTPRSzeehqDQjzhqnQUOxbdwO6z+fgpORdSCfDWhJ6QSTljcWswWwtOQkIC3NzcTMUNAERHR0MqleLIkSNNukZSUhKSk5Mxe/bs6742b948eHp6YsCAAVizZs1Nm2irq6uh0+kabERE1m7+PSEIaueMfF01Pog7J3YcugFdVS3+sfkUAODRocGcsbgVma3A0Wg08Pb2brDPwcEBHh4e0Gia1py6evVqdOvWDYMGDWqw/80338R3332HXbt2YeLEiXjyySfxySefNHqdpUuXQqVSmbaAgIDm3xARkYVROsrw9n1hAICvjlxG0uWrIiei/7V0x1lodFUIaueM56K7iB3HrjS7wFm4cGGjHYGvbWfP3vksm5WVlfj6669v2Hrz2muvYfDgwejTpw9efvllvPTSS/jggw8avVZsbCy0Wq1py8rKuuN8RESWYHCIJ+7v2x6CALz83z9QVasXOxLVO3SpEN/Uz3nz3sRecJLLRE5kX5rdB+f555/HzJkzb3pMcHAw1Go1CgoKGuyvq6tDcXFxk/rO/PDDD6ioqMD06dNveWxkZCSWLFmC6upqKBTX90xXKBQ33E9EZAteG9sd+y8U4mJBGf658xxeGdtd7Eh2r6KmDgv/mwIAmDowEJGc86bVNbvA8fLygpeX1y2Pi4qKQklJCZKSkhAREQEA2L17NwwGAyIjI295/urVqzFu3Lgmfa/k5GS4u7uziCEiu+TeRo537w/D7PXHsOpAOkb2UKN/kIfYsezaP3eeR2ZxBfxUSi6mKRKz9cHp1q0bRo0ahTlz5iAxMREHDx7E/PnzMXnyZNMIqpycHISGhiIxMbHBuRcvXsS+ffvw6KOPXnfdn3/+GatWrcKpU6dw8eJFfPbZZ3jnnXfw1FNPmetWiIgs3ohuPnggwh+CADz/3UmUV9eJHcluHcsoxpqD6QCAd+4PgwtHt4nCrPPgbNy4EaGhoRgxYgTGjBmDIUOG4IsvvjB9vba2FufOnUNFRUWD89asWQN/f3+MHDnyums6Ojpi+fLliIqKQnh4OD7//HN8+OGHWLx4sTlvhYjI4r12b3f4qZTILK7Au79wxXExlFbV4tlvkyEIwMS+/hjW1fvWJ5FZSAQ7nAJTp9NBpVJBq9XC1dVV7DhERC3mwIVCTF1tnIrjq9mRGNLZU+RE9mXBd8nYfDwH/u5O+OWZoWy9aWHNeX9zLSoiIhsypLMnpg4MBAC89MNJaCtrRU5kP7b9kYvNx3MglQAfTQpncSMyFjhERDYmdnQ3BHo4I1dbhdjNf3CtqlaQp63EK1uME/rNGx6CfuzkLToWOERENqaNwgEfP9QHDlIJdqRo8HX9XCxkHgaDgOe/M7aW9fZX4ekRncWORGCBQ0Rkk8ID3EzDk9/8+QzOarhEjbmsPpCOQ5eK4OQow78mhcNRxlerJeBPgYjIRs0e0hHDunqhus6A+V+fQEUNh463tGMZxXgvzjhibdG93RHs1fYWZ1BrYYFDRGSjpFIJ/vlAb/i4KnCxoAyvbz0tdiSbUlhWjXlfH0edQcC9vf0wuT/XObQkLHCIiGxYu7YKfDSpDyQS4Ltj2fjxRI7YkWyC3iDg6W9OIF9XjRDvtnj3/jBIJBKxY9FfsMAhIrJxUZ3a4al7jB1fYzen4HSuVuRE1u/DXedw6FIRnOUyrJjaF20UzV75iMyMBQ4RkR14ZkRn3NXFC5W1eszdkISismqxI1mt387kY/meSwCAdyf2Qoi3i8iJ6EZY4BAR2QGZVIJPJvdBUDtn5JRUYt7Xx1GrN4gdy+pcLirHgu+SAQAzBwVhXG8/cQNRo1jgEBHZCZWzI1ZO74c2chkOpxXj7e2pYkeyKtqKWsxadxS6qjr0CXTDP8Z0EzsS3QQLHCIiO9LZxwX/mhQOAFh3KAPfHc0SN5CVqKkz4LGvjiHtSjn8VEp8PjUCcge+Qi0ZfzpERHZmZA81novuAgB49cdTOJJWJHIiyyYIAv6xJQWH04rRVuGA1TP7w9tVKXYsugUWOEREduipe0IwuqcaNXoDHt1wDKl5nOm4Mf/Zewk/JGVDKgE+ebgPuvnefBVrsgwscIiI7JBUKsG/JoWjf5A7SqvqMGNNIrKKK8SOZXG2/ZGLD349BwB4fVwPDO/qLXIiaioWOEREdkrpKMOq6f3R1ccFBaXVmLEmEcXlNWLHshi/n7+CBd+eBADMGhyE6VFB4gaiZmGBQ0Rkx1TOjlj/yAC0d3NCWmE5Zq07ivJqrlmVcKkIczccQ43egDFharw6trvYkaiZWOAQEdk5tUqJ9Y8MgLuzI05mleCxL5NQWaMXO5Zoki5fxez1R1FdZ8CIUG98NKkPZFIuw2BtWOAQERFCvNtizcz+cJbLcOBiIR6x05aclGwtZq5JREWNHkNCPLF8Sl8OB7dS/KkREREAoE+gOzY8MgBtFQ5ISCvCjDWJKK2qFTtWqzmdq8W0NUdQWl2HAUEe+GJ6BJSOMrFj0W1igUNERCb9gjzw1aORcFU64Njlq5i2OhHaStsvchIuFWHy54dRUlGL8AA3rJnVH85yLqBpzVjgEBFRA+EBbvh6zkC4OTsiOasEU1YdtunFOXek5Blbq6rrENnRAxtmG1uxyLqxwCEiouv0bK/CprkD0a6NHKdydBi//CDOaUrFjtXivkzIwLyvj6NGb8Donmqsf2QAXJWOYseiFsACh4iIbihU7YrvHo9Ch3bOyL5aiYmfHcLus/lix2oRBoOAZb+ew2s/nYYgAFMHBuLTh/uyz40NYYFDRESN6uTVFj8+ORgDgz1QVl2H2euPYdX+NAiCIHa021ZcXoNZ647i0z0XAQAL/tYFS8b35FBwG8MCh4iIbsq9jRwbHonEQwMCIAjAW9tT8cL3f6DMCoeRn8i8iv/7eD9+P38FSkcp/vlAbzw9ojMkEhY3toYFDhER3ZLcQYp37gvDa//XHVIJ8N/j2Rjz7/1IulwsdrQmEQQB6w6m48HPE5CrrUJHzzbY8uRgTIzwFzsamYlEsOZ2xtuk0+mgUqmg1Wrh6spVYYmImiPhUhFe+P4kckoqIZUATwzrhGdGdLHYCfGyiivw2k+nsPfcFQDAmDA13pvYCy7sTGx1mvP+ZoHDAoeIqNl0VbV4fetpbD6eAwDo4eeKpfeHoZe/m7jB/qJWb8Cq/en4d/x5VNUaIJdJ8fLoUDwyOIgfSVkpFji3wAKHiKhl7EjJwz+2pKCkwjgZ4PhwP7wwsisCPJxFzXUsoxivbDmFc/nGoe0Dgz3w9n1h6OTVVtRcdGdY4NwCCxwiopZToKvCu3FnseVEDgQBkMukmDk4CPOGhUDl3LofAx1JK8Kney5i/4VCAIC7syNeHdsd9/dtz1YbG9Cc97fZPjB9++23MWjQIDg7O8PNza1J5wiCgEWLFsHX1xdOTk6Ijo7GhQsXGhxTXFyMKVOmwNXVFW5ubpg9ezbKysrMcAdERNQU3q5KfPhgOH6ePwSDQ9qhRm/AF/vSMPi93Xj1xxScydWZ9fsLgoA95wrwwIpDmPTFYey/UAiZVIJJ/QIQ//wwTIzwZ3Fjh8zWgrN48WK4ubkhOzsbq1evRklJyS3Pee+997B06VKsX78eHTt2xGuvvYaUlBScOXMGSqUSADB69Gjk5eXh888/R21tLWbNmoX+/fvj66+/bnI2tuAQEZmHIAjYe/4K3t1x1vTxEAD0CXTDlMgO+Ft3H6ic7rxVRxAEpORosf2PPGxPyUP21UoAxtajB/r547G7OiGwnbgfk1HLs6iPqNatW4dnn332lgWOIAjw8/PD888/jxdeeAEAoNVq4ePjg3Xr1mHy5MlITU1F9+7dcfToUfTr1w8AEBcXhzFjxiA7Oxt+fn5NysQCh4jIvARBQEJaETYeycSvpzSoMxhfNVIJENZehahOnhgc0g59A93RpgnrPtXqDUi7Uo4zeVqkZOuw84zGVNQAgLNchocHBGLOXcHwcVWa7b5IXM15f1vMamLp6enQaDSIjo427VOpVIiMjERCQgImT56MhIQEuLm5mYobAIiOjoZUKsWRI0dw33333fDa1dXVqK7+c6E4nc68zaVERPZOIpFgUCdPDOrkiYLSKnx/LBubj2fj0pVynMzW4mS2Fit+vwQAcHN2hNpVCT83J/i4KgBIUFWrR1WtHpW1elwprcaF/DLU6A0NvoeTowz3dPPG2DBfDO/qDSc5l1mgP1lMgaPRaAAAPj4+Dfb7+PiYvqbRaODt7d3g6w4ODvDw8DAdcyNLly7FG2+80cKJiYioKbxdlJg3PATzhocgT1uJhEtFOHixCIcuFSJPW4WSilqUVNTi7C0W82yrcEA3Xxd083XFwOB2LGropppV4CxcuBDvvffeTY9JTU1FaGjoHYVqabGxsViwYIHpzzqdDgEBASImIiKyT74qJ9zf1x/39zXOIKyrqoVGW4XckkpotFXQ6KoglUigdJTCyVEGpaMMKidHhKpd4e/uBCnXi6ImalaB8/zzz2PmzJk3PSY4OPi2gqjVagBAfn4+fH19Tfvz8/MRHh5uOqagoKDBeXV1dSguLjadfyMKhQIKheK2chERkfm4Kh3hqnREFx8XsaOQjWlWgePl5QUvLy+zBOnYsSPUajXi4+NNBY1Op8ORI0fwxBNPAACioqJQUlKCpKQkREREAAB2794Ng8GAyMhIs+QiIiIi62O2eXAyMzORnJyMzMxM6PV6JCcnIzk5ucGcNaGhodiyZQsAY4e0Z599Fm+99Ra2bt2KlJQUTJ8+HX5+fpgwYQIAoFu3bhg1ahTmzJmDxMREHDx4EPPnz8fkyZObPIKKiIiIbJ/ZOhkvWrQI69evN/25T58+AIA9e/Zg2LBhAIBz585Bq9WajnnppZdQXl6OuXPnoqSkBEOGDEFcXJxpDhwA2LhxI+bPn48RI0ZAKpVi4sSJ+Pjjj811G0RERGSFuFQD58EhIiKyChaxVAMRERGRWFjgEBERkc1hgUNEREQ2hwUOERER2RwWOERERGRzWOAQERGRzWGBQ0RERDaHBQ4RERHZHBY4REREZHPMtlSDJbs2ebNOpxM5CRERETXVtfd2UxZhsMsCp7S0FAAQEBAgchIiIiJqrtLSUqhUqpseY5drURkMBuTm5sLFxQUSiaRFr63T6RAQEICsrCyuc2VGfM6tg8+5dfA5tw4+59ZjrmctCAJKS0vh5+cHqfTmvWzssgVHKpXC39/frN/D1dWV/4BaAZ9z6+Bzbh18zq2Dz7n1mONZ36rl5hp2MiYiIiKbwwKHiIiIbA4LnBamUCiwePFiKBQKsaPYND7n1sHn3Dr4nFsHn3PrsYRnbZedjImIiMi2sQWHiIiIbA4LHCIiIrI5LHCIiIjI5rDAISIiIpvDAuc2LF++HEFBQVAqlYiMjERiYuJNj//+++8RGhoKpVKJsLAw7Nixo5WSWrfmPOeVK1di6NChcHd3h7u7O6Kjo2/5cyGj5v59vmbTpk2QSCSYMGGCeQPaiOY+55KSEsybNw++vr5QKBTo0qUL/9/RBM19zh999BG6du0KJycnBAQE4LnnnkNVVVUrpbVO+/btw7333gs/Pz9IJBL8+OOPtzxn79696Nu3LxQKBUJCQrBu3Tqz54RAzbJp0yZBLpcLa9asEU6fPi3MmTNHcHNzE/Lz8294/MGDBwWZTCa8//77wpkzZ4RXX31VcHR0FFJSUlo5uXVp7nN++OGHheXLlwsnTpwQUlNThZkzZwoqlUrIzs5u5eTWpbnP+Zr09HShffv2wtChQ4Xx48e3Tlgr1tznXF1dLfTr108YM2aMcODAASE9PV3Yu3evkJyc3MrJrUtzn/PGjRsFhUIhbNy4UUhPTxd+/fVXwdfXV3juuedaObl12bFjh/DKK68ImzdvFgAIW7ZsuenxaWlpgrOzs7BgwQLhzJkzwieffCLIZDIhLi7OrDlZ4DTTgAEDhHnz5pn+rNfrBT8/P2Hp0qU3PP7BBx8Uxo4d22BfZGSk8Nhjj5k1p7Vr7nP+X3V1dYKLi4uwfv16c0W0CbfznOvq6oRBgwYJq1atEmbMmMECpwma+5w/++wzITg4WKipqWmtiDahuc953rx5wj333NNg34IFC4TBgwebNactaUqB89JLLwk9evRosG/SpElCTEyMGZMJAj+iaoaamhokJSUhOjratE8qlSI6OhoJCQk3PCchIaHB8QAQExPT6PF0e8/5f1VUVKC2thYeHh7mimn1bvc5v/nmm/D29sbs2bNbI6bVu53nvHXrVkRFRWHevHnw8fFBz5498c4770Cv17dWbKtzO8950KBBSEpKMn2MlZaWhh07dmDMmDGtktleiPUetMvFNm9XYWEh9Ho9fHx8Guz38fHB2bNnb3iORqO54fEajcZsOa3d7Tzn//Xyyy/Dz8/vun9U9Kfbec4HDhzA6tWrkZyc3AoJbcPtPOe0tDTs3r0bU6ZMwY4dO3Dx4kU8+eSTqK2txeLFi1sjttW5nef88MMPo7CwEEOGDIEgCKirq8Pjjz+Of/zjH60R2W409h7U6XSorKyEk5OTWb4vW3DI5rz77rvYtGkTtmzZAqVSKXYcm1FaWopp06Zh5cqV8PT0FDuOTTMYDPD29sYXX3yBiIgITJo0Ca+88gpWrFghdjSbsnfvXrzzzjv4z3/+g+PHj2Pz5s3Yvn07lixZInY0agFswWkGT09PyGQy5OfnN9ifn58PtVp9w3PUanWzjqfbe87XLFu2DO+++y5+++039OrVy5wxrV5zn/OlS5eQkZGBe++917TPYDAAABwcHHDu3Dl06tTJvKGt0O38ffb19YWjoyNkMplpX7du3aDRaFBTUwO5XG7WzNbodp7za6+9hmnTpuHRRx8FAISFhaG8vBxz587FK6+8AqmUbQAtobH3oKurq9labwC24DSLXC5HREQE4uPjTfsMBgPi4+MRFRV1w3OioqIaHA8Au3btavR4ur3nDADvv/8+lixZgri4OPTr1681olq15j7n0NBQpKSkIDk52bSNGzcOw4cPR3JyMgICAlozvtW4nb/PgwcPxsWLF00FJACcP38evr6+LG4acTvPuaKi4roi5lpRKXCZxhYj2nvQrF2YbdCmTZsEhUIhrFu3Tjhz5owwd+5cwc3NTdBoNIIgCMK0adOEhQsXmo4/ePCg4ODgICxbtkxITU0VFi9ezGHiTdDc5/zuu+8Kcrlc+OGHH4S8vDzTVlpaKtYtWIXmPuf/xVFUTdPc55yZmSm4uLgI8+fPF86dOyds27ZN8Pb2Ft566y2xbsEqNPc5L168WHBxcRG++eYbIS0tTdi5c6fQqVMn4cEHHxTrFqxCaWmpcOLECeHEiRMCAOHDDz8UTpw4IVy+fFkQBEFYuHChMG3aNNPx14aJv/jii0JqaqqwfPlyDhO3VJ988okQGBgoyOVyYcCAAcLhw4dNX7v77ruFGTNmNDj+u+++E7p06SLI5XKhR48ewvbt21s5sXVqznPu0KGDAOC6bfHixa0f3Mo09+/zX7HAabrmPudDhw4JkZGRgkKhEIKDg4W3335bqKura+XU1qc5z7m2tlZ4/fXXhU6dOglKpVIICAgQnnzySeHq1autH9yK7Nmz54b/v732bGfMmCHcfffd150THh4uyOVyITg4WFi7dq3Zc0oEge1wREREZFvYB4eIiIhsDgscIiIisjkscIiIiMjmsMAhIiIim8MCh4iIiGwOCxwiIiKyOSxwiIiIyOawwCEiIiKbwwKHiIiIbA4LHCIiIrI5LHCIiIjI5rDAISIiIpvz/7TueDTNFdUPAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Matplotlib has a number of plotting tools that are extremely powerful for\n",
    "# exploring data. The plotting functionality is simliar to R, but many times\n",
    "# it is easier of more customizable. \n",
    "# Note: everthying you know in \"R\" will come back full circle. Python and R play very \n",
    "#   nicely with one another and python can call R code, share memory with it, and use its plotting functions\n",
    "#   ...with the correct libraries installed! Check out Rpy2, http://rpy.sourceforge.net \n",
    "\n",
    "%matplotlib inline\n",
    "\n",
    "#============================================\n",
    "# plotting example\n",
    "from matplotlib import pyplot as plt\n",
    "\n",
    "plt.plot(t, y)  # plot them\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "2.6.0\n"
     ]
    }
   ],
   "source": [
    "# this library ships with python and is for accessing SQL\n",
    "# You will learn how to use the module more extensively in the future!\n",
    "# here, we will just import it and get its version\n",
    "import sqlite3\n",
    "\n",
    "print (sqlite3.version)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {},
   "outputs": [],
   "source": [
    "# this is \n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Debugging\n",
    "\n",
    "Debugging in python is much like any other programming language. If you use the command line to run programs, you need to be able to catch and interpret runtime errors. There is no compilation process (that we have talked about), so all errors will occur at runtime--that includes things like spelling errors. Because of this, it is critical to run your code in varying contexts. That also means exploring ways that it could fail by changing the values of different variables.\n",
    "\n",
    "One easy way to do this is by the using the pdb package. You can run it from the iPython notebook, but it gets quite long on the output trace. You will likely want to copy this test over to another document (for example, debugging.py) and run it from the command line (for example, >python debugging.py). The pdb library allows you to set break points in the files and also step through the program, changing the variable values. \n",
    "\n",
    "+ import pdb\n",
    "    + pdb.set_trace()\n",
    "+ command line arguments\n",
    "    + s(tep), c(ontinue), n(ext), w(here), l(ist), r(eturn), j(ump)\n",
    "+ and much more… like print, p, pp\n",
    "+ can set numbered break points by running from python window\n",
    "    + python -m pdb your_function.py "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "> \u001b[0;32m/var/folders/8f/szs_zmy96hd71pg9d355wt9h0000gp/T/ipykernel_54612/3434866926.py\u001b[0m(4)\u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n",
      "\u001b[0;32m      2 \u001b[0;31m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0m\u001b[0;32m      3 \u001b[0;31m\u001b[0mj\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;36m0\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0m\u001b[0;32m----> 4 \u001b[0;31m\u001b[0;32mfor\u001b[0m \u001b[0mi\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mrange\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;36m5\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0m\u001b[0;32m      5 \u001b[0;31m    \u001b[0mj\u001b[0m \u001b[0;34m+=\u001b[0m \u001b[0;36m5\u001b[0m\u001b[0;34m*\u001b[0m\u001b[0mi\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0m\u001b[0;32m      6 \u001b[0;31m    \u001b[0mpdb\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mset_trace\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0m\n",
      "ipdb> c\n",
      "> \u001b[0;32m/var/folders/8f/szs_zmy96hd71pg9d355wt9h0000gp/T/ipykernel_54612/3434866926.py\u001b[0m(4)\u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n",
      "\u001b[0;32m      2 \u001b[0;31m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0m\u001b[0;32m      3 \u001b[0;31m\u001b[0mj\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;36m0\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0m\u001b[0;32m----> 4 \u001b[0;31m\u001b[0;32mfor\u001b[0m \u001b[0mi\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mrange\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;36m5\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0m\u001b[0;32m      5 \u001b[0;31m    \u001b[0mj\u001b[0m \u001b[0;34m+=\u001b[0m \u001b[0;36m5\u001b[0m\u001b[0;34m*\u001b[0m\u001b[0mi\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0m\u001b[0;32m      6 \u001b[0;31m    \u001b[0mpdb\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mset_trace\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0m\n"
     ]
    }
   ],
   "source": [
    "import pdb\n",
    "\n",
    "j = 0\n",
    "for i in range(0,5):\n",
    "    j += 5*i\n",
    "    pdb.set_trace()\n",
    "    \n",
    "print (j)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "There is still plenty more to learn about python, but now you have the basics. One topic that may be intersting to you is the concept of decorators, discussed here:\n",
    "    http://simeonfranklin.com/blog/2012/jul/1/python-decorators-in-12-steps/\n",
    "    \n",
    "Essentially, decorators are ways of manipulating functions with other functions. Mostly, this means recycling the functionality of one function to create a new function. \n",
    "\n",
    "There is also much to be learned about iterators, yield, etc. Python is an active language and it is still evolving!"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.8.16"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}