Eigenvalues.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Eigenvalues"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In this chapter we shift focus away from solving linear systems, and look closer at the effect of matrix multiplication.  We restrict our attention now to square matrices, which define linear transformations from $\\mathbb{R}^n$ to $\\mathbb{R}^n$.  In this context we will study special values called eigenvalues, and corresponding vectors called eigenvectors, that can be used to analyze the effect of a corresponding matrix.\n",
    "\n",
    "Given a square $n\\times n$ matrix $A$, a scalar $\\lambda$ is called an **eigenvalue** of $A$ if there exists some nonzero vector $V$ in $\\mathbb{R}^n$ such that $AV=\\lambda V$.  The vector $V$ is the **eigenvector** associated with $\\lambda$.  The equation states that when an eigenvector of $A$ is multiplied by $A$, the result is simply a multiple of the eigenvector.  In general, there may be multiple eigenvalues associated with a given matrix, and we will label them as $\\lambda_1$, $\\lambda_2$, etc., to keep an orderly notation.  We will label eigenvectors in a similar way in order to track which eigenvectors are associated with which eigenvalues.  \n",
    "\n",
    "We will visualize examples in $\\mathbb{R}^2$.  "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Example 1:  Matrix representing horizontal stretch\n",
    "\n",
    "Let's consider first the following matrix.\n",
    "\n",
    "$$\n",
    "\\begin{equation}\n",
    "A = \\left[ \\begin{array}{cc} 2 & 0 \\\\ 0 & 1 \\end{array}\\right]\n",
    "\\end{equation}\n",
    "$$\n",
    "\n",
    "We saw in [Chapter 3](Planar_Transformations.ipynb) that multiplication by this matrix has the effect of a horizontal stretch."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "tags": [
     "hide-input"
    ]
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQAAAAD8CAYAAACYVXqwAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/d3fzzAAAACXBIWXMAAAsTAAALEwEAmpwYAAAer0lEQVR4nO2dbWxcZ5XH/yd+ieMWE0UppqlJGpEG2tDEaSLobkWXRrAqy4oVSIsoS6nYaPOBdmER2hVsukILWq2Q2MIXhLZqQ0DsgngpUGC1tEBChIBmHdfGiTPN2LE7+CUehjhMJ5lMMp6zH8bj2L6eO3fmnnufM/b5SZZfcuc+vzuT+9xzz33O8xAzwzCMtck61wKGYbjDOgDDWMNYB2AYaxjrAAxjDWMdgGGsYawDMIw1TKvETohoHMArAOYAFJl5v8R+DcOIFpEOYJ4HmDkjuD/DMCLGbgEMYw1DEiMBiWgMwCwABvCfzPzkCtscAnAIADo7O/fdcccdqLS9bt06FItFtLW14dq1a+jo6MCVK1fQ2dnp+X716lW0t7fj+vXraG1tRalUquwfpVIJra2tuHbtGtavX4+rV69iw4YNnn3k83l0dHTg2rVraG1txcDAAPbs2QMAKJVKaGlpQbFYRHt7e819FAoFtLe3o1gsYt26dWqOaW5uDuvWrbNjWqPHNDAwkGHmW2qeu0IdwBZmniKi1wB4HsDfM/OJatvffffdPDQ0FLpdKYgImoZEj46O4vWvf71rjSVoczIff4joVJBcnMgtADNPzX9PA/gegDf7bd/aKpl6WH1s2rTJtYIHbU7mI0PoDoCIbiKiV1V+BvDnAE77vaYSOhkrc+XKFdcKHrQ5mY8MEpfibgDfI6LK/v6bmf9XYL9rlso9oSa0OZmPDKE7AGY+D2BPPa+Z7yyMKrS1tblW8KDNyXxkcNJt2S2AP7lczrWCB21O5iODkw7AkoD+bN682bWCB21O5iODkw7g2rVrLpptGiYmJlwreNDmZD4yOOkA1q9f76LZpmHHjh2uFTxoczIfGZx0AFevXnXRbNNw5swZ1woetDmZjwwiIwHrZf/+/dzX1xd7u9XQNhLQMMIS60jAemnWQRNxcerUKdcKHrQ5mY8MFgHAIgBj9WERQBOj8Wqizcl8ZLAIABYBGKsP1RFAPp930WzToKlUuoI2J/ORwUkH0NHR4aLZpmHnzp2uFTxoczIfGWwkoEJSqZRrBQ/anMxHBqsFUEh3d7drBQ/anMxHBicdwNzcnItmm4ZLly65VvCgzcl8ZHDSATTr5AlxoTFHos3JfGSwM9Ew1jA2IYhCNBZLaXMyHxmcdAAtLS0umm0aNm7c6FrBgzYn85HBSQdQLBZdNNs0zMzMuFbwoM3JfGRw0gG0t7e7aLZp2Lp1q2sFD9qczEcGmxBEIefOnXOt4EGbk/nIYMVAsGIgY/WhuhjIyoH90Vhaqs3JfGSwCAAWARirD4sAmhiNVxNtTuYjg0UAsAjAWH3EHgEQUQsRvUhEP6q1rU0I4s/g4KBrBQ/anMxHBslbgI8BOBtkw2YtnIiLXbt2uVbwoM3JfGQQ6QCIqAfAuwA8FWT7QqEg0eyqZWRkxLWCB21O5iODVATwRQD/BCBQlY+NBPSnp6fHtYIHbU7mI0PoDoCI/hJAmpl906BEdIiI+oiob3p6GplMBtPT05icnMTs7CxGR0eRz+cxPDyMUqmE/v5+ADeyq/39/SiVShgeHkY+n8fo6ChmZ2cxOTmJyv7Gx8eRy+WQSCRQLBYX7ssq+6h8HxoaQqFQQDKZRDabBQCk02mk02mkUilks1kkk0kUCoWFyR6X72NwcBDFYhGJRAK5XA7j4+Nix5RMJkMfUyqVEj2mZDLp/HNafEznzp1z/jktPqaXXnpJxedUOaaghH4KQET/DuBhAEUAHQC6ADzDzB+s9pre3l4eGBgI1a4k2p4CZDIZdctNa3MyH39iewrAzJ9i5h5mvh3A+wH83O/kn39N2GZXNdevX3et4EGbk/nIYDMCKUTjhCnanMxHBtHpeZn5OIDjtbazOQH96ezsdK3gQZuT+chgE4Io5OLFi64VPGhzMh8ZnHQAbW1tLpptGrZs2eJawYM2J/ORwVYGUsjY2JhrBQ/anMxHBisGgr7HgKVSSV2eRJuT+fhj5cBNjKYxEhW0OZmPDBYBQF8EYBhhsQigidE4uYQ2J/ORwSIAWARgrD4sAmhiKsUomtDmZD4yWAQAfRGAtowyoM/JfPxRHQHYwiD+JBIJ1woetDmZjwy2NJhCtm/f7lrBgzYn85HBSQfQrKWTcTE1NeVawYM2J/ORwUkH0NoqWoS46ti0aZNrBQ/anMxHBicdQLPWTseFxqck2pzMRwY9aUtjAU3Z5AranMxHBifWROSi2aZBY7m0NifzkcFuARSSy+VcK3jQ5mQ+MlgSUCGaZpetoM3JfGSwCUEUMjEx4VrBgzYn85HBSQewfv16F802DTt27HCt4EGbk/nIYEOBFXLmzBnXCh60OZmPDFYMBH3FQIYRFtXFQM06aCIuNE4uoc3JfGSwCAAWARirD4sAmhiNVxNtTuYjg0UAsAjAWH2ojgDy+byLZpuGyrrwmtDmZD4yhO4AiKiDiE4S0SARnSGif631mo6OjrDNrmp27tzpWsGDNifzkUEiAigAOMDMewD0AniQiO71e4GNBPQnlUq5VvCgzcl8ZAjdAXCZSiVE2/yX7w211QL4093d7VphgatXgW9+E7jlFj1OgKL36I9/BL7wBXQ36SxXIjkAImohogEAaQDPM/MLK2xziIj6iKjvwoULyGQymJ6exuTkJGZnZzE6Oop8Po/h4WGUSqWFaZYr2dX+/n6USiUMDw8jn89jdHQUs7OzmJycxPT0NDKZDMbHx5HL5ZBIJFAsFjE4OLhkH5XvQ0NDKBQKSCaTyGazAIB0Oo10Oo1UKoVsNotkMolCobBwb7d8H4ODgygWi0gkEsjlchgfHxc7prGxsdDHlEqlQh/T2bPj+Ju/yeBLX5rG+fNjzj+nxcc0Ojrq9nPq68Pk4cOYfuc7kXn6aYyMjTn7nFY6psAws9gXgI0AjgF4k992e/bsYU1gIZDRwczMjGsFnppi3rGDmYj5Ax/Q4bQYJz7XrjF/+9vM99zDvGEDc0sLc1sb8+c+p+79AdDHAc5Z0VicmS8R0XEADwI4LblvIz5GRoC3vhXIZICbbgL+7u9cGyng6lVg61Ygnwcqtf8bNgDvehfwj/8I/P73bv0aROIpwC1EtHH+5w0A3g7Ad5J0mxDEH5fFUv39wP79wMwMUCwC7e3A/ffrK+CK3aelBXjNa4DKILaWFmD7duBrXwOI1L0/QZGIAG4F8FUiakG5Q/kWM//I7wUtLS0Cza5eNm7c6KTdiQngT/8UKBTKv7e2Ah/6ELBunTunasTu8+tfA4sr/rq6gOeeK0cBLnyEkHgK8Ftm3svMu5n5Tcz8mVqvKRaLYZtd1czMzDhp99WvBu5d9AC3vR348IfdOlUjVp+DB4E/+zNgyxbg+98Hbr4Z+PGPgdtuc+MjiJPncbYykD9bt2510u6rXlW+2gPlC1t3N7B7t1unasTic/ly+WQHgCeeAD7+8fLPr7zixicCbEIQhZw7d85JuydPAseOAUePAqOjwK9+5d6pGpH7nDhx4+RPJm+c/K58IsKKgWDFQADAXL76EwFrPkd78CBw5Eg55E+lygm/JkN1MZCVA/vjorT0wIHy9xWiWwD6yl0j8bl8udwDHjlSDvknJwOf/Nren6BYBACLAE6eBN7ylnLo/8gjrm0cceJEOdEHlEP+Jp3ks4JFAE1MnFcT5vLJT+R/8mu7won6LM7yF4sNnfza3p+gWASAtR0BPPAAcPx4eXDbTTe5tomZaln+VYDqCMAmBPGnUkgSNSdPlk/+o0drn/xxOQUltE+dWf7IfRzhJALYt28fawqZtEUAxWIx8pLperP+cTjVQyifCLL82t4f1RFAoTLW1FiRkZGRyNuolfVfThxO9dCQT4gsfyQ+CrCRgArp6emJdP/1hP4Vonaql7p9Is7ya3t/guIkArBaAH8ymUxk+w6a9Y/TqRHq8hHI8ov6KMJJBLBunZN+p2m4uZKcioB6Q/8KUTo1QiCfGLP82t6foDjpADQl3DRyPaL55RoJ/StE5dQoNX1iHtij7f0Jil2KFRLFhCmNhv5ROoXB1yeGkL8uH8XYLYBCOjs7xffZaOhfIQqnMKzo43Bgj7b3JyiWBHTMl7/8ZXzkIx9Z+P3xxx/HhyuzcAgRJvSvcPHiRUml0Hh8hAf2hPV57LHHsG3btiV/W+mzfvjhh2Pxq0qQmUOlv/bu3Ss1+akIcDgr8OXLl7mnp4dnZ2f5hz/8Iff29vIf/vAHsf2XSsxAeXbfMFy5ckVGSIglPn/7t+WD3LKFuVh07nP+/HnetWsXv/a1r+VsNrvw95U+66jeVwScFdhJBGArA92gs7MTDz30EA4fPoyPfvSj+M53voMLFy6I7T9s6F9hbGwsvIwgY2NjkQ7sachnnk9/+tN4/PHHcdddd+HMonkEV/qsN8zPKeiMIL2E9Ne+ffsi6fUaBY7XBTh79iwD4B/84AfMzDw3Nyey3xdeKF8Yjx4Nvy8pJynmjh8vHxzAnEy61ll4f06fPs179+7lUqnEjz76KD/11FNLtlv+WUcFNEcAVg68lM985jO45ZZbFnIjAwMDofcZNuu/HAknMQ4exMAzz8Sa5a9F5f05fPgwPvvZz4KIcOedd+L06aXLYyz/rJ0TpJeQ/rII4Aaf//zn+T3veQ8fO3aM77//frH9vu1t5YtjLie2S/fkcjeu+k884drGw29+8xvu6Ojgbdu28bZt27i7u5vf8Y53LPx7VJ/1SiBgBOCkA7jzzjujPPa6cdUB/OxnP+Ndu3YtJIp6e3v5xRdf5L6+vlD7lQz9K4R1Cs0vfrEk5Hfus4y+vj4+cOAA//SnP13424ULF/jWW29l5uqfdVSo7gAsAmB++eWXefv27Tw8PLzwt6985Sv8yCOPhNqvVNZfFQqy/LV47rnn+L777vP8vaurK7LP2g/VHYBFAP6cOnWq4ddGFfqHcWoYn5DfiY8P2nyCdgBORgI266ipuOjt7W3odRIDfqrRqFPD1BjLH7tPDbT5BMUWBlFIIuG7tuqKsHDWX8KpYQKM5Y/VJwDafIIisTrw64joGBGdJaIzRPSxWq+xCUH82b59e92vkRrwU41GnOqmjoE9sfjUgTafoEhEAEUAn2DmOwHcC+BRIrrL7wXNWjoZF1NTU3VtH2XoX6Fep7qpcyx/5D51os0nKBKrA08zc//8z68AOAvgNr/XaJo8USObNm0KvG3UoX8jTnXTQPlupD4NoM0nKKI5ACK6HcBeAC+s8G+HiKiPiPpmZmaQyWQwPT2NyclJzM7OYnR0FPl8HsPDwyiVSujv7wdwY8GF/v5+lEolDA8PI5/PY3R0FLOzs5icnMT09DQymQzGx8eRy+WQSCRQLBYXpmqu7KPyfWhoCIVCAclkEtlsFgCQTqeRTqeRSqWQzWaRTCZRKBQwNDS04j4GBwdRLBaRSCSQy+UwPj4udkwTExOBj+lznyt/7+tbekypVEr0mCYmJuQ/p6EhFDdswGBbG/DEEzj17LNAS0vVz0n6mCT/77388ssN/d+L6pgCE+RRQZAvADcDOAXgvbW23b17d4QPQOoHyh4DTk1NBdouigE/1QjqFJhlA3uc+4REmw/irAUgojYA3wXwX8z8TIDtJZpdtbS1tdXcJq7Qvx6nwAjM2CPqI4A2n6BIPAUgAE8DOMvMTwR5TbNOnxQXuVyu5jZRZ/2XE8SpJoLluyI+gmjzCYpENu4+AA8DGCKigfm//TMz/0/VRi0J6MvmzZt9/z2OrP9yajnVRHiSztA+wmjzCYrEU4BfMjMx825m7p3/qnryAzYhSC0mJiaq/lvcoX8Qp5pEMElnKJ8I0OYTFCeX4vXr17totmnY4XOCxB36V/BzqkqEk3Q25BMh2nyCYkOBFbJ4GqnFuAj9K1RzqkrEk3TW7RMx2nyC4mR14P3793NfX1/s7VZD2+rAK8F1rubrlAhW3zXqQ/XqwDYlmD8rLZ3uKvSvEGg59xgn6dS0vDygzycoFgFAfwRw8mQ58Xf0aLyJv7qIeSkuwx+LAJqYxVcTV1l/PycPDpbi0nbF1eYTFIsAoDsCeOCBcuIvl4s/8VcTh0txGf6ojgDy+byLZpuGShGIy6z/cipOCzheisvj4xhtPkFxEgHs27ePNYVM2iKAQqGA9vb1qrL+hULhxvgNBVn+JT4K0OajOgKwkYD+pFIp51n/5aRSKVVLcaVSKSftVkObT1CcjAS0WgB/0uluNaF/he6JCWDnzvIvCrL83d3dTttfjjafoDg5E+fm5lw02xQwA//yL5dA1KXnkd/Bg7g0NoYuRQN7Ll26hK6uLtcaC2jzCYqTW4B165w02xQcOADMznboCP0Xhfwd73uf05B/OR0dHa4VlqDNJygWiyuikvX/+tcVhP7LB/Y04dXNqI2TS7FNCOJl8YCft77VcbHUCgN7tBVwmY8MTiKAFiVhpCYWZ/3n5ja6kfAZ2LNxoyOnKpiPDE4iADVroyth+YCfmZmZ+CVqDOxx4uSD+chgA4HgdiDQSmW+sQ8qCTCwR9tAF/PxR/VAoGa9X4qClQb8nDt3Lp7G6xjYE5tTQMxHBic5gA0bNrhoVh3Vxvrffffd0TdeZ/luLE51YD4yWDmwI/zKfCO/PWqgfFfTLRtgPlJYOTDc5ACclPla+e6aQXUOYK1HALXKfCO5moQs39V2hTMfGSwCQLwRgJPJPRWU7xrxojoCWMsTggQp862sLBsawfJdMSchzEcGJ08BmrVwIixBZ/jZtWtX+MaEJ+kUcRLEfGRwEgEUCgUXzTqlnsk9R0ZGwjUWwSSdoZ2EMR8ZpJYHP0JEaSI6HWT79vZ2iWabinpm+Onp6WmskQhn7GnYKSLMRwapCOAogAeDbrzWagHqndwzk8nU30jEk3Q25BQh5iODSAfAzCcAXAzc6BqaEKSRef1vrpzIQYlhXv66nSLGfGSI7UwkokNE1EdEfel0GplMBtPT05icnMTs7CxGR0eRz+cxPDyMUqmE/v5+ADeer/b396NUKmF4eBj5fB6jo6OYnZ3F5OQkpqenkclkMD4+jlwuh0QigWKxuJCZreyj8n1oaAiFQgHJZBLZbBYAkE6nkU6nkUqlkM1mkUwmUSgUFqZ7Xr6PwcFBFItFJBIJ5HI5jI+Pr3hM7353Hh/84DD++Mfgx3ThwoVgx3T5Mk59/OPAkSMY+upXUTh/Hsnz55HNZpFKpUSP6cKFC84/p8XHNDU1Jfo5hT2mycnJ0Mck+X8vMMws8gXgdgCng2y7e/du1kT5bZDnhReYAeajR+t73cTERO2NfvGL8s4B5mSyMUFppxgxH38A9HGAc9HmBIyIMEt6dXZ2+m/gYCmumk4xYz4y2IQgERFmXv+LF6ukUxzOy1/VyRHmI4PUY8BvAPg1gDcQ0QQRHfTbvq2tTaJZtYRd0mvLli3ePzpeimtFJ4eYjwxSTwEeYuZbmbmNmXuY+Wm/7VfzykASq/mOjY0t/YODkL+mk2PMRwYbCiyMxJJeb3zjG8s/KCrfXXBSgvnIYOXAgkit5jswMOA85F/RSRHmI4OVA0OmHFi0zNfKd42QqC4HXo0RgMhqvvNZ/lNdXc5X312OtgkvzEcGiwAQPgI4ebKc+Dt6tPHEn3T5rrG2sQggJiSy/suz/P3zw5M1URkeqwXzkcHJU4BmHTW1EqFC/ypZ/t7eXhE3SbQ5mY8MtjBICEJl/X2y/IlEQkpRDG1O5iODkxzAPffcw5pCpkZyAKGy/jWy/Pl8Xt3iKdqczMcf1TmA69evu2hWlIZC/4Bj+aempmQkBdHmZD4yOMkBtLY6aVaMhkL/OrL8mzZtCqsojjYn85HBSQRQim1CfHkayvrXOZZf41MSbU7mI0NzX4odUFfo3+BYfo3zJWhzMh8ZnHQAROSi2dDUFfqHGNijsVxam5P5yGC3AAGpK/QPWb6by+UaF40IbU7mI4MlAQMSKPQXKt/dvHlzQ6+LEm1O5iODkwig2SYECRT6C5bvTkxMNPzaqNDmZD4yOBkItG/fPtZUPeU3ECjQgB/h8t1isaguStLmZD7+qB4I1ExDgX1D/4gm6Txz5kzofUijzcl8ZLByYFSPAHzLfK1811CM6gigGQZN+Gb9I56kU9PtUQVtTuYjg0UAWDkCeOCBcuIvl1uU+FM0Sadh+GERQAhWzPrHOEmnxquJNifzkcEiACyNAFbM+tsknUaToToCyOfzLpoNxJKsv6OluCqrwmpCm5P5yGALgyxiSeh/yl2Wf+fOnbG1FRRtTuYjg40EnGdJ1v+E26W4UqlUrO0FQZuT+cggtTjog0T0EhGNENEna22vacRUhQMHgE5cRondrL67mO7u7tjbrIU2J/ORIXQHQEQtAL4E4J0A7gLwEBHd5feaubm5sM2KM3f8BC5Dx1Jcly5dctZ2NbQ5mY8MEhHAmwGMMPN5Zr4G4JsA/sq3UUWTJ1QegpyA29V3F6MxR6LNyXxkkDgTbwPwu0W/T8z/rSkYfnYEAFD49/9QtRSXYcSBxM34StP7eAYXENEhAIcW/S7QtBwdn/oE8KlPuNYwjFiR6AAmALxu0e89ADxzJDPzkwCeBIC9e/fyiy++KNC0DBKrA0uSzWbR1dXlWmMJ2pzMx5+gF1iJW4D/A3AHEW0nonYA7wfwrN8LisWiQLOrl5mZGdcKHrQ5mY8MoSMAZi4S0WMAfgKgBcARZvYtjm5vbw/b7Kpm69atrhU8aHMyHxlE0vHM/D/MvJOZX8/M/1Zr+2aaEMQF586dc63gQZuT+chgxUDQlwMwjLCoLgbSXg7sGo2lpdqczEcGiwBgEYCx+rAIoInReDXR5mQ+MlgEAIsAjNWH6ghA84QgGhgcHHSt4EGbk/nIYAuDQF8EoG2RCUCfk/n4ozoCKBQKLpptGkZGRlwreNDmZD4yOOkAbCSgPz09Pa4VPGhzMh8ZnHQAVgvgTyaTca3gQZuT+cjgpAPQNCGIRm6urD+gCG1O5iODkzNRU8JNI9evX3et4EGbk/nIYJdihZSqrkPuDm1O5iOD3QIopLOz07WCB21O5iODJQEVcvHiRdcKHrQ5mY8MTjqAtrY2F802DVu2bHGt4EGbk/nIYCsDKWRsbMy1ggdtTuYjgxUDQd9Q4FKppC5Pos3JfPxRPRTYyoH9GRgYcK3gQZuT+chgEQD0RQCGERaLAJoYTZWSFbQ5mY8MFgHAIgBj9WERQBPT39/vWsGDNifzkcEiAOiLALRllAF9Tubjj+oIwBYG8SeRSLhW8KDNyXxksAlBFLJ9+3bXCh60OZmPDE46gGYtnYyLqSnP4srO0eZkPjI46QA0TZ6okU2bNrlW8KDNyXxkCNUBENFfE9EZIioRUc2EQ4VmrZ2OC41PSbQ5mY8MYSOA0wDeC+CEgIsxj6ZscgVtTuYjQ6hYnJnPAuXHaPVQ7/ZrDY3l0tqczEeG2G7GiegQgEPzvxaI6HRcbQdgMxFpmtZ1MwBNPoA+J/Px5w1BNqrZARDRTwG8doV/OszMPwhqw8xPAnhyfp99QQYpxIX51Eabk/n4Q0SBRtrV7ACY+e3hdQzD0EhzZi4MwxAh7GPA9xDRBIA/AfBjIvpJwJc+GabdCDCf2mhzMh9/Avk4KQYyDEMHdgtgGGsY6wAMYw3jrANodBhxBB4PEtFLRDRCRJ905THvcoSI0lrGSBDR64joGBGdnf+sPubYp4OIThLR4LzPv7r0qUBELUT0IhH9yLULABDROBENEdFArceBLiMA58OIiagFwJcAvBPAXQAeIqK7XPkAOArgQYftL6cI4BPMfCeAewE86vj9KQA4wMx7APQCeJCI7nXoU+FjAM66lljGA8zcW2tsgrMOgJnPMvNLrtqf580ARpj5PDNfA/BNAH/lSoaZTwBQs8YUM08zc//8z6+g/J/8Noc+zMy5+V/b5r+cZrGJqAfAuwA85dKjUdZ6DuA2AL9b9PsEHP4H1wwR3Q5gL4AXHHu0ENEAgDSA55nZqQ+ALwL4JwCaSlwZwHNEdGp+CH5VIq0FkBpGHCErVSXZc9FlENHNAL4L4B+YOevShZnnAPQS0UYA3yOiNzGzk5wJEf0lgDQznyKit7lwqMJ9zDxFRK8B8DwRJeajSw+RdgBNMIx4AsDrFv3eA6A5p3aJCCJqQ/nk/y9mfsa1TwVmvkREx1HOmbhKmt4H4N1E9BcAOgB0EdHXmfmDjnwAAMw8Nf89TUTfQ/lWd8UOYK3fAvwfgDuIaDsRtQN4P4BnHTupgcp1208DOMvMTyjwuWX+yg8i2gDg7QCczcbJzJ9i5h5mvh3l/zs/d33yE9FNRPSqys8A/hw+HaTLx4CNDiMWg5mLAB4D8BOUE1zfYuYzcXtUIKJvAPg1gDcQ0QQRHXTlMs99AB4GcGD+kdLA/NXOFbcCOEZEv0W5836emVU8elNEN4BfEtEggJMAfszM/1ttYxsKbBhrmLV+C2AYaxrrAAxjDWMdgGGsYawDMIw1jHUAhrGGsQ7AMNYw1gEYxhrm/wEUMtjZkkamPgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "fig, ax = plt.subplots()\n",
    "\n",
    "options = {\"head_width\":0.1, \"head_length\":0.2, \"length_includes_head\":True}\n",
    "\n",
    "ax.arrow(0,0,2,3,fc='b',ec='b',**options)\n",
    "ax.arrow(0,0,4,3,fc='r',ec='r',**options)\n",
    "\n",
    "ax.set_xlim(-1,5)\n",
    "ax.set_ylim(-1,5)\n",
    "ax.set_aspect('equal')\n",
    "ax.set_xticks(np.arange(-1,6,step = 1))\n",
    "ax.set_yticks(np.arange(-1,6,step = 1))\n",
    "\n",
    "ax.text(0.8,2.2,'$X$')\n",
    "ax.text(4,2.2,'$AX$')\n",
    "\n",
    "ax.axvline(color='k',linewidth = 1)\n",
    "ax.axhline(color='k',linewidth = 1)\n",
    "\n",
    "ax.grid(True,ls=':')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "$$\n",
    "\\begin{equation}\n",
    "X = \\left[ \\begin{array}{c} 2  \\\\ 3 \\end{array}\\right] \\hspace{1cm}\n",
    "AX = \\left[ \\begin{array}{c} 4  \\\\ 3 \\end{array}\\right]\n",
    "\\end{equation}\n",
    "$$"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "For this example it is possible to deduce the eigenvalues and eigenvectors since the effect of the matrix is rather simple.  Any vector that does not have its direction changed by this transformation is an eigenvector.  In particular, we notice that any vector along the $x_1$-axis will be scaled by the stretch, but the direction will be unchanged.  Let's choose the unit vector to label as $V_1$\n",
    "\n",
    "$$\n",
    "\\begin{equation}\n",
    "V_1 = \\left[ \\begin{array}{c} 1  \\\\ 0 \\end{array}\\right] \\hspace{1cm}\n",
    "\\end{equation}\n",
    "$$\n",
    "\n",
    "We can calculate and observe that $AV_1 = 2V_1$  This tells us that $\\lambda_1 = 2$ is the eigenvalue associated with $V_1$.  It is important to note that given $\\lambda_1 = 2$, there is a whole collection of eigenvectors that satisfy the equation $AV = 2V$.  Indeed if $V=kV_1$, where $k$ is any scalar, we have \n",
    "\n",
    "$AV = A(kV_1) = k(AV_1) = k(2V_1) = 2(kV_1) = 2V$.\n",
    "\n",
    "We now consider if any other vectors whose directions are unchanged by the transformation. \n",
    "After a bit of thought, we realize that any vector along the $x_2$-axis is entirely unaffected by the matrix multiplication.  Let's label $V_2$ as the unit vector along that axis.\n",
    "\n",
    "$$\n",
    "\\begin{equation}\n",
    "V_2 = \\left[ \\begin{array}{c} 0  \\\\ 1 \\end{array}\\right] \\hspace{1cm}\n",
    "\\end{equation}\n",
    "$$\n",
    "\n",
    "Since $AV_2 = V_2$, it must be that the associated eigenvalue is $\\lambda_2 = 1.$  Again we note that there is a whole collection of eigenvectors associated with $\\lambda_2$, since any multiple of $V_2$ is also unchanged by the transformation."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Example 2:  Projection matrix\n",
    "\n",
    "For the next example, let's consider a matrix which projects vectors orthogonally onto its one-dimensional column space.\n",
    "\n",
    "$$\n",
    "\\begin{equation}\n",
    "B = \\left[ \\begin{array}{cc} 0.2 & -0.4 \\\\ -0.4 & 0.8 \\end{array}\\right]\n",
    "\\end{equation}\n",
    "$$\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[[-0.6]\n",
      " [ 1.2]] \n",
      "\n",
      "[[0.]\n",
      " [0.]]\n"
     ]
    }
   ],
   "source": [
    "B = np.array([[0.2, -0.4],[-0.4, 0.8]])\n",
    "X = np.array([[1],[2]])\n",
    "Y = np.array([[2],[1]])\n",
    "\n",
    "print(B@X,'\\n')\n",
    "print(B@Y)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "tags": [
     "hide-input"
    ]
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQAAAAD8CAYAAACYVXqwAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/d3fzzAAAACXBIWXMAAAsTAAALEwEAmpwYAAApuElEQVR4nO2de3yU1Z3wvychIQm3GLkYCCiC4VKRQPD21lUsVN2uVm33Um29rmJ1t9q17XatLm3lbf3USl/crrXaKqK1Vnpx8dJtrfVe5RpzkTCQRMIkIRICCWFgmMlkzvvHMGkGyOSZzDNzfpOc7+eTDxnyXL7P78xzznnOcy5Ka43FYhmeZJkWsFgs5rAZgMUyjLEZgMUyjLEZgMUyjLEZgMUyjLEZgMUyjHEtA1BKZSulPlBKvezWMS0WS2pxswZwF7DNxeNZLJYU40oGoJQqAf4O+Lkbx7NYLOlhhEvHWQX8OzCmvw2UUsuAZQAFBQXlZ5xxBtFeiFlZWYRCIXJycggGg+Tl5XH48GEKCgqO+/fIkSPk5ubS3d3NiBEjCIfD0eMTDocZMWIEwWCQkSNHcuTIEfLz8487ht/vJy8vj2AwSPvhHg74Q8yckE/W0WNkZ2cTCoXIzc0d8BiBQIDc3FxCoRBZWVkirmnEiBH09PSQlRXJ3924purqasrKyobUNQ3FdIpeU2VlZbvWesJAN27SGYBS6nKgTWu9RSm1uL/ttNaPA48DzJs3T1dWViZ7alfYfyjIrsZGFnyi1LRKLw0NDcyYMcO0RgxKKT744APTGr1Ii5E0H6XULifbuVED+CTwWaXUZ4A8YKxS6hda6y/1e9IRblU8kqdoVC5qciSj7O4Jk5Nt/sVIUVGRaQXxSIuRNB+nJP1t11rfo7Uu0VqfBnwBeD3ezQ/0Vp2kcOjQIW5Zs5nl67aaVgHg8OHDphXEIy1G0nycIqcoNkh2djZzi8cwriDXtApA7zOhpX+kxUiaj1NczQC01m8Cbw60nVLKzdMmTU5ODndfMsu0Ri85OTmmFcQjLUbSfJxiJNuS9gjg8/kA0Frz9o69tHT6RfhY+kdajKT5OMVIBiCpERBg/PjxAOz1BbhlzWaefq9RhI+lf6TFSJqPU4xkAMFg0MRp+6W5uRmAiWPyePbWc/ma4ceBqI+lf6TFSJqPU5SJKcHKy8v1li1b0n7e/giFQsfVSrTWxtoqTuRjGqUUkqaPkxYjaT5KqS1a60UDbWekBnDkyBETp+2XrVtjX/9VN3dy+Y/fpbnDzKudY30sxyMtRtJ8nGIkA8jPzzdx2n6ZP39+zOfxo0cS1rD3YECEj+V4pMVImo9TjGQA0jpNHPs4Mrkwn9/feQELpp0kwsdyPNJiJM3HKUbaABYtWqQ3b96c9vMmSndPmMqmTs4+LTO7ebqJtDYAS3xEtwFIrwFEWfXaDq55fD2tB9LbLyBTS5N0Ii1G0nycYmsAcWg94GdbaxcXz5oorvdiurE1gMxCdA3A7zfb0+5YampqTvj/xePy+dTsSWm/+fvzsfwVaTGS5uMUIxlAXl6eidP2S2lp/LkA1rzXyA//6EmTzcA+FnkxkubjFNsTEPB6vXH/Xt/m48OWLsLh9FSBB/KxyIuRNB+nGOm6JKnHFMCkSZPi/n35FXPTOlHIQD4WeTGS5uMUIzWAnp4eE6ftl87Ozrh/j978B/zd7POlvnPQQD4WeTGS5uMUIxmAtMkTnLRJHOnuYcnKN3no1R0ifIY70mIkzccpsurigsnLyeYbl87izCnjTKtYLK5hJAOQNiGI08FJ/3T2tBSbRJA2WEoi0mIkzccpRuri2dnZJk7bL4WFhY633X8oyIqXa1M6UjARn+GKtBhJ83GKkQwgFAqZOG2/7Nmzx/G2R7p7+OUGL+s/2i/CZ7giLUbSfJxiJwQBAoEAI0eOdLx95+EghSmcQThRn3QgrSuwtBhJ8xHdFVja89KOHYm17Edv/gOHu1Ohk7DPcERajKT5OMVOCALMmzcv4X3+54MWznvgzylpCxiMz3BDWoyk+TjFDgdmcEM5z5lexBfOmUpejvsNmpIej6QiLUbSfJxihwNbHCGtDcASH9FtAEOhBhDF83EXz210dyBIppYm6URajKT5OMVIBlBQUGDitP1SXl4+6H1/ucHLg3/wcDjo3qvNZHyGC9JiJM3HKUlnAEqpPKXURqVUlVJqq1LquwPtI21CkKqqqkHv+9Wlpbzx9cUU5LrXqTIZn+GCtBhJ83FK0m0AKjJdziittU8plQO8C9yltV7f3z7S+gG4tahDd0/YlWHD0haZAHltANJiJM0nbW0AOkJ0ZcScoz9xvymBgJn59vujvr4+qf211tyyZhPL130owmc4IC1G0nyc4kobgFIqWylVCbQBf9Jab4i3fW5u6nrRDYaSkpKk9ldKMbd4LDMnjhHhMxyQFiNpPk5xJQPQWvdorcuAEuAcpdSZx26jlFqmlNqslNrc2tpKe3s7ra2ttLS00NHRQUNDA36/n9raWsLhMBUVFcBfW1crKioIh8PU1tbi9/tpaGigo6ODlpYWosdrbGzE5/Ph8XgIhUK9z2XRY0T/rampIRAIUFdXR1dXFzt27KCtrY22tja8Xi9dXV3U1dURCAR6J3s89hhVVVWEQiE8Hg8+n4/PlY7kytljXLmmurq6pK/J6/UmfU2NjY296QQYT6e+17Rjx46kr8nN79727dtFpFP0mpziej8ApdS3gUNa64f626asrExXVla6et5kaG9vd2V5Z601b9e1M3PiaKYUDr63o1s+biKtDUBajKT5pK0NQCk1QSlVePT3fGApEHcKXUlfJIDubnf69O/1Bbh1zWaefq9RhM9QRlqMpPk4xY1my2JgjVIqm0iGslZr/bILx00bbk1QMnFMHs/eei7zSwpF+AxlpMVImo9Tks4AtNbVwIJE9pE2J6CbHZOi6whqrQe9oIi0jlISkRYjaT5OsROCAPv3uzu5R3VzJ5f/+N1BjxR022coIi1G0nycYiQDyMnJMXHafpk8ebKrx5swZiRaw96Dg+vv4LbPUERajKT5OMWuDATs3LnT1eMVj8vnlTsvYMG0k0T4DEWkxUiaj1Ps2oDA7NmzXT+mUorunjCbGhOvGqbCZ6ghLUbSfJxihwMDqeqTsOq1HVzz+HpaDyQ2+ElSHwmpSIuRNB+n2AlBUkjrAT+ejw+yuHRC2pcYdxtpHYEs8bETgiRAqkYmFo/L5+JZExO++SWNlJSKtBhJ83GKrQGkgTXvNdJ28AjfuDQznxPB1gAyDVsDSIDo4I9UUd/mY+vuLsJhZzdQqn2GAtJiJM3HKbYGQKQbZyp7JyY6UUiqfQaDtBqAtBhJ8xFdA5C2MIjHE3fsUtJEb/4D/m72+QbuHJRqn6GAtBhJ83GKkQxA2oQg06dPT/k5jnT3sGTlmzz06sAryKTDJ9ORFiNpPk4xkgFIGzq5e/fulJ8jLyebb1w6i+vOO1WET6YjLUbSfJxiZBZDSZMnAhQVFaXlPP909jRH26XLJ5ORFiNpPk4xUgOQNnY6nW8l9vkC3P9SLU37+z+ntLckEpEWI2k+TpHTbGmQdLbeBkJhntvoZcPO/scISGpNloq0GEnzcYoRa2ndYtM5PHlyYT7v3/Mp/r68/1lkpQ2XNs2jjz7KHXfc0fv5vvvu46677jJodDyZmmb2EQDw+XwDb+QihQWRtyAHDp+4MTTdPtK54YYbeOmll+js7OTll1/mlVde4f777zetFUOmppmRDEBaI6CJ2Vz/54MWzn3gtRPOGiRpdlkJFBQUcM0113Dvvfdy55138pvf/EbcPPyZmmZ2QhCgubk57ec89/QirjlnGnk52SJ8pHPzzTfzk5/8hFWrVjFjxgxxMZLm4xQjRfHIkSNNnLZfZs6cmfZzFo/L59tXfOKEfzPhI53777+fCRMm9M4nKS1G0nycYrsCA1u3bjV2bs/HXTy30RvzfyZ9JLJy5UqOHDnC2rVrefjhhwF5MZLm4xQjNYD8/MGvmpMK5s+fb+zcz6738nL1bq4sm9y7xLhJH2m8/vrrrF69mvfff58xY8bQ1dVFZWUlZWVlptViyNQ0s8OBMTuZw799upQ3vr649+Y37SMJr9fLLbfcwq9//WvGjIksvHrXXXexatUqcTGS5uMUOxxYEIkOG04n0oYDW+IjejiwrQHEorXmljWbWL5uqwifTEBajKT5OMVIG4C0ZZTKy8uNnl8pxdzisb0dhEz7ZALSYiTNxylGagB+f2LTZKea6DrsJrn7klncfEFkTLkEH+lIi5E0H6e4sTz4VKXUG0qpbUqprUqpATtpS1sYpLS01LQCEHkUeHN7G6MnOhs2PJyRkmZRpPk4xY0aQAj4mtZ6DnAe8C9KqbnxdpDWE9Dr9Q68URrY6wuw7OktPPpaZr5TTidS0iyKNB+nuLE8eCvQevT3g0qpbcAUoLbfkwobCzBp0iTTCgBMHJPHL289l+njZL4JkML//i/MmycjzaJI+Q4liqvfNKXUacACYMMJ/rZMKbVZKbX5448/pr29ndbWVlpaWujo6KChoQG/309tbS3hcLh3muVo62pFRQXhcJja2lr8fj8NDQ10dHTQ0tJCa2sr7e3tNDY24vP58Hg8hEIhqqqqYo4R/bempoZAIEBdXR1dXV00NDTQ1tZGW1sbXq+Xrq4u6urqCAQCvc92xx6jqqqKUCiEx+PB5/PR2NjoyjXNGKdo9u5i9+7dSV2T1+t19ZoA4+nk9Xr5zW/a+Na32qisbDCaTsdeU319vYh0il6TU1zrB6CUGg28BXxPa/27eNuWlZVpSWuptbW1MXHiRNMavbxZvZMH32jm8evLKTlJxhsTCf0ANmyAxYvhyBFoaGjj9NPlpJm071Ba+wEopXKA3wLPDnTzWwbm5FE5KAXtPlltJSZpbIRLL43c/CUlMHq0aaOhQdIP4yoyvc8TwDat9Y+c7CNtQhBpg5PG5YR5+SsXiJs5yRSdnZGS/+DByOe/+Rt5aSbNxylu1AA+CVwHfEopVXn05zPxdsjOPn4MvEkKCwtNK8RQWFiIUopgKMzGOHMHDge6u+Ezn4HWVgiHYdQoWLJEZpplIklnAFrrd7XWSmt9lta67OjP7+PtEx3TLYU9e/aYVogh6vPwn3dw7c/W03pAVsepdPKXv8D69dD3zfH558tNs0zDrgwETJsmq+NN1Oe6807jZzcs4pSxsjpOpZOLLoKf/CTye35+pBYwe7bcNMs07IQgwI4dAy/XlU6iPqeMy+PiWROHdVuAUnD33ZHf166Fn/0MsrLkplmmYYcDZwBP/WUne30BvnHpbGMOpl4DNjfD1Knw4otwxRVpP33GYocDJ4C0oZzH+nzUfoja3V2Ew8NvPH60i/2xN7/0NMsUbA0gA5AwUYiJGoAt/QePrQEkgLTc+1if6M1/4HA3+3wBE0pG6K/0B/lplikYyQDshCDxOZHPke4elvzoTR56NTMbmxKluRn8/kjpfyIyIc0yATshCPQO3JDCiXzycrL590tnc/35pxowSj/xSn/IjDTLBIyMy5U2IcgnPnHiBTpM0Z/PP549Nc0mZhio9IfMSTPpGKkBBAKynmPr6+tNK8QQz2efL8D9L9XStF9WO4qbDFT6Q2almWRsT0AQt9BkPJ9AKMxzG71sGKJjBJyU/pBZaSYZIxmAtLEA7e3tphViiOczuTCf9fcs4e/LM/MLNxBOSn/IrDSTjJEMICtL1pRXo4UNLh/IZ1xBDhB5LTiUcFr6Q+almVSM3ImmZ5Y5lu5uWTeSE58XPmjm3Adeo7lj6LQFOC39ITPTTCKyimJDSJugxInPeaefzLXnnEp+jqy5FQZLIqU/ZGaaScTIa0BpjwDSOiY58Skel8/yK+LOvp5RJFL6Q2ammURsIyCwf7+sFvWoz2OPPUZxcTFlZWXMnDmTq666imAwyKOPPsodd9wBwLbWLq6++Stcd911JpWTItHSH+SmWcahtU77z4IFC7QkDh8+bFohhqjPHXfcoZ988kmttdY9PT16xowZetOmTfrQoUO6pKREd3R06H/61n/r/FNm6PbOrpQ6Rb4qqSE/X+tEDy81zaQAbNYO7kUjNQBpKwPt3LnTtEIMUZ+amhoWLFgARDqaaK0pLS2loKCAa665hnvvvZf1v3yI9/70MiePG2NSedAMpvQHuWmWadjhwEQacCS1S0R9Tj75ZKZMmUJ3dzfNzc288sorXHjhhQB4PB7mzJnDunXr+OxnPwukdthwqoYDFxREMoBEDy01zaRghwMngKRFSiDi09TUxMSJE6murmbbtm088sgjrFixoneb+++/nwkTJhAKhdBa889PbWL5usxaU3CwpT/ITLNMxA4HBhYuXGhaIYaFCxdSXV3N3Ll/beWfP38+bW1tAKxcuZIjR46wdu1aHn74YZRSfGLyWM6YmFmdURJt+e+LxDTLRGwNAHmTOWzZsoWamhrmzJkDRBpq16xZw9KlS3n99ddZvXo1a9asYfHixXR1dVFZWcndl8zi5qlZ8M47hu2dkUzpDzLTLCNx0lLo9k95eXlqmj6HENdee62eMmWKLisr0wsWLNC333673rVrl54+fbqura3t3W716tX6hhtu0FprHX5opX7jjHN189sbXPfB5bcAg2n5tzgHh28BjGQAc+bMSeW1J8yWLVtMK8QwWJ89b/xFn/G1F/T3L7td6127XHVyMwNoaop88158cfDHGCppliqcZgD2LQDyWnAH7RMIsGnWOczfvZ3cqVPggw9g7FhXnNx8CzDYlv++DJk0SxGi3wJIWxjE4/GYVohh0D4jR3L2qBC53QF0SwtcfjkI63WZ7LN/lCGTZoZxa3nwJ5VSbUqpD51sL21CkOnTp5tWiCEpn4svprr4DP7umgdp3rELli1Lrqh1mWRa/vsypNLMIG7VAJ4CLnO6sbShk7t37zatEENSPhddxAQdRGlNe1YerF4Ne/e6J5cEbpX+MMTSzCCujAbUWr+tlDrN8UlHGBmE2C9FRUWmFWJIyuf88ynet5uX13wVBfDeezBxoltqSeFW6Q9DLM0MkrY2AKXUMqXUZqXU5j179tDe3k5raystLS10dHTQ0NCA3++ntraWcDhMRUUF8Nf3qxUVFYTDYWpra/H7/TQ0NNDR0UFLSwutra20t7fT2NiIz+fD4/EQCoV6p2qOHiP6b01NDYFAgLq6Orq6uvB6vbS1tdHW1obX66Wrq4u6ujoCgQA1NTUnPEZVVRWhUAiPx4PP56OxsdG1a2pubh78Nfn9dP3Xf9H0xhu0LCjnj799zZVrApK6purqFs48s5V164ZOOvX97u3atUvUNTnFtbcAR2sAL2utzxxo2/nz52tJ86i3trZSXFxsWqMXt3wevPNHPD5yBu/cs5TiolFJHSvZtwButPz3ZaimmVs4fQtgpC4ubbnrnJwc0woxuOVz3Xdv45zyCzllz2/h6addOeZgcPPZP8pQTbN0Y+Q1oLTpk3w+n2mFGNzyKT5pFIuXLEQ98wz09LhyzMHg5rN/lKGaZunGrdeAzwHvA7OUUs1KqX+Ot720RsDx48ebVojBVZ/HHuOphZfz4J0/cu+YCZCK0h+GeJqlEVcyAK31NVrrYq11jta6RGv9RLztpU0I0tzcbFohBld9srJo+JtL2NZ2iHB3+jsFpaL0hyGeZmnESFE8cuRIE6ftl5kzZ5pWiMFtn+U/vJ0RuTmo/I/S2haQqtIfhn6apQvbFRjYulXWRBpu++TkjEDdcgsH1v6O9gPpG4qdqtIfhn6apQs7GGiYcCTQzSe/vpZL2McDP74z4f0TfQ3Y3AxTp0ZK/1RkAJb4iB4MZCcEiU8qfPJG5vBN/RHXP///0vJGIJWlPwyPNEsHtgYwnAiHITsbrrsu4baARGoAtvQ3j60BJIC03DtlPllZ7Lv1Dr67O4/m9tS9t0516Q/DKM1SjJ0UFCgvLzetEEMqfYI/eJDnz7qE9fc+mJLjp7Llvy/DKc1SiZEMwO/3mzhtv0QHXUghlT7FJ43ifd/r/P3jK1LSFpCO0h+GV5qlEiMZQF5enonT9ktp9FsrhFT7jPvpfwPQefNtrh43XaU/DL80SxV2aTDA6/WaVogh5T5ZWbxwx3c4b/zfutoWkK7SH4ZhmqUIIxmAtLEAkyZNMq0QQzp8zvu/X+fayj+Q//V/c+V46Sz9YXimWSowkgH0GByZdiI6OztNK8SQDp/ik0ax/HQ4ec3PXWkLSGfpD8MzzVKBkQxA0vTJIK9NIm0+jz3Gtgmn8dy/rBh42ziku/SHYZxmLiPrTrSkl6wsnrv2bn6YW8phf2DQh0l36W9xDyMP49ImBJE2OCmdPl9dcStfmzCegs4/DGqkoInSH4Z3mrmJkRpAdna2idP2S2FhoWmFGNLpUzQmj3HXXQPPPEN3MPHp2k2V/sM5zdzESAYQErZazZ49e0wrxJBuH/3Tn/LPn1/O8q89mtB+pkp/sGnmFkYyAGkrA02bNs20Qgzp9lHZ2Xxi8lhmvfenhN4ImHz2H+5p5hZ2QhBgx44dphViMOFz96p/48aKl+Gmmxxtb7L0B5tmbmGHA1t60bfeypuvf0DpxjeZcvLomL8dOxzY7Xn+Le5ihwMngLShnKZ89v7wYW67+j6eXv7TuNuZLv3Bpplb2BqAJYbNt3+Tsx7/EbnBI5HJQ47StwZgS3/52BpAAkjLvU36LHrkAXLDIXQ/bQESSn+waeYWtgZgOY6qL3+DewJTefyHN1MyPtIWEK0B2NI/MxBdA5A2IYikhUrBvM/EB75Dlta0f/O+mP+XUvqD+RgdizQfpxipAZSXl2tJVaZQKCRqiLIIn1tvhZ//HEIhyM5GKUV+vhZT+ouIUR+k+YiuAQQCgx94kgrq6+tNK8QgwuexxwhmjWDDbf/e+19SSn8QEqM+SPNxiluLg16mlNqulKpXSv3HQNtL6wlYUlJiWiEGET5ZWaz6lx/wxZMupHX/od7/ljLiT0SM+iDNxylJ11mUUtnAI8CngWZgk1LqRa11bX/7SBsL0N7ezujRowfeME1I8bnuu7dxTvmFjPnot4Cc0h/kxCiKNB+nuFEDOAeo11p/pLUOAr8Crox7UmETgkhLOCk+xSeNYvGShYz+3TOAnNIf5MQoijQfp7jRajEFaOrzuRk499iNlFLLgGV9Prtwaks6sWk29HCjKD7Rt+K4dmKt9eNa60Va60VnnXUWWmsxP83NzcYdJPvc90JkzvuenrBxF6kxkubjFDcygGZgap/PJcDuuCcV9gggbaUiaT7Lr5gLQFaWnBqAtBhJ83GKG3fiJuAMpdR0pVQu8AUgbnORtEbA/fv3m1aIQZpPTnbka9J5OEi7T8YrXGkxkubjlKQzAK11CPhX4I/ANmCt1nprvH1ycnKSPa2rTJ482bRCDNJ8oixZ+RYrX91uWgOQFyNpPk5xpS6utf691rpUaz1Da/29gbaXtjLQzp07TSvEIM0nyjf/djY3/J/TTGsA8mIkzccpdjAQkVmKJbVLSPOB4ycEMY20GEnzEd0VWNpw4MrKStMKMUjz6cs+X4DvvrSV5g6zaSgtRtJ8nGIkA5DWYrpw4ULTCjFI8+lLIBTm+U1NbPjIbKOXtBhJ83GKrQEgbzIHaT59mVyYz/v3LOHz5Wb7vkuLkTQfp9g2AIsjTtQG0Hk4SGGBrIFdlgi2DSABKioqTCvEIM3nRLzwQTPnPfBnY20B0mIkzccptg0AKCsrM60QgzSfE3He6SfzxXNPJT/HzDJv0mIkzccpdmEQwOPxmFaIQZrPiSgel89/Xj6Xk0ePNHJ+aTGS5uMUuzQYMH36dNMKMUjziUft7i5+ucGb9vNKi5E0H6cYyQC6uxNfhTaV7N4dd+xS2pHmE49fbfLy0KvbORxM7/gOaTGS5uMUI7MYSpo8EaCoqMi0QgzSfOLx1aWlfO3TsyjITW+aSouRNB+nGKkBhMNhE6ftF2lvJaT5xKNoVC7jCiKDu7p70peu0mIkzccpcjovG0RSH26Q5zMQWmtuWr2R5eviDgJ1FWkxkubjFCN1cWlTS0kbnizNZyCUUswrKaSoIH3e0mIkzccp9hEA8Pl8phVikObjhLs/XcqNn0xfS7i0GEnzcYqRDEBaI+D48eNNK8QgzccpWmve8LTR0pn6pd+kxUiaj1OMZADSJgRpbm42rRCDNB+n7D0Y4LZntvD0+40pP5e0GEnzcYpdGxB567pJ8wHnE4JsbtzPWSWF5I5IbdkiLUbSfEQPBpLWFXjr1vS1XjtBmk8iLDqtiNwRWSmfPUhajKT5OMVIBpCfn2/itP0yf/580woxSPNJlKqmTj7zX++mdKSgtBhJ83GKHQ6MvMkcpPkkysSxI8nOgn2+1LX1SIuRNB+n2AlBLI6QNimoJT6i2wBsDSA+0nwGSzAUZsNH+1JybGkxkubjFFsDsDhiMDWAB//g4fG3P+Ldb36KU8blpcjMciJE1wD8/tR3FEmEmpoa0woxSPMZLNedfypP3Hg2k8a6P2lITU0Nfr+fiy66iJ6eHgAee+wxiouLKSsrY+bMmVx11VUEg0GCwSAXXnhhSpeky9Q0M/LiMi9PVmlQWlpqWiEGaT6DpXhcPsXj3Hvj09TUxBNPPEFLSwunnnoqJ510Ep/73OfIzo5MS1ZdXc33v/99brrpJsLhMKWlpVRXV7No0SKWLFnC888/zxe/+EXXfPqSqWlmewICXm/6Z7SJhzSfZFn9l508+IfkpszasGEDX/nKV7j99tv5/Oc/z7p163j22We58sore7epqalhwYIFANTX16O17r0xr7rqKp599tmkHOKRsWmW5Brk/wBsBcLAIqf7lZWVaUkcOHDAtEIM0ny01jryVRkc971Qo29avVH39IQHtX9XV5eeMWOGbm1t1Vpr7ff7dV1dnZ40aVLMdkVFRXrevHl69uzZevTo0fqtt97q/VsoFNLjx48f9DUMhLQ0AzZrB/disjWAD4HPAW8nslP0mU0KnZ2dphVikOaTLN++Yi5P3ng2WVmDGwb+61//moULF3LKKacAkUfIjo4OCgsLe7dpampi4sSJVFdXs23bNh555BFWrFjR+/fs7Gxyc3M5ePBgUtfSH5maZkllAFrrbVrrhNeLljZ5grQ2CWk+yTIiO5LenYeDtPsCCe+/fft2ZsyYEfN/48aNi+lSXl1dzdy5c3s/z58/n7a2tph9AoFAymKbqWkm6060DFn8wR4+tfItVr6acHnBjBkzqK6u7v384YcfEgwG6enp6c0EampqmDNnDhB5rF2zZg1Lly7t3Wffvn1MmDAhYyfuSBUDZgBKqdeUUh+e4OfKgfY95jjLlFKblVKb9+zZQ3t7O62trbS0tNDR0UFDQwN+v5/a2lrC4XDvSivRDhYVFRWEw2Fqa2vx+/00NDTQ0dFBS0sLra2ttLe309jYiM/nw+PxEAqFqKqqijlG9N+amhoCgQB1dXV0dXXR0tJCW1sbbW1teL1eurq6qKurIxAI9L7eOfYYVVVVhEIhPB4PPp+PxsZG166ptbU16Wvyer2uXhOQ1DXt3/sxXz7/FK6cU5jwNV1wwQUUFhZy0UUXceONN/Luu+/S3d3N4sWL+dWvfkUgEOCdd97hqaeeYtasWZSXl9PU1MSKFSt6r+kXv/gFS5cudTWd+n73mpubRaRT9Joc46ShYKAf4E1sI6BrSPPROrlGwFRw4MABXVFRob/0pS852v7qq6/WHo8npT6SIE2NgIMilR0yBsOePXtMK8QgzcdN9vkCfOfFrTTtT647+J49e1iwYAEXX3zxgI3KwWCQq666ilmzZiV1zoF8MpGkMgCl1NVKqWbgfOAVpdQfnewnbWWgadOmmVaIQZqPmwRCYdZubmLjzv1JHScao5tvvrm3I1B/5Obmcv311yd1Pqc+mUaybwFe0FqXaK1Haq0naa0vdbKftAlBduzYYVohBmk+bjK5MJ/371nC58tLkjqOtBhJ83GKHQxkcUQqhgN3Hg5SWCCrNjhUED0YyA4Hjo80n1TwwgfNnPfAnwc9a5C0GEnzcYqtAVgc4XYNoPWAn5+/s5M7Fs8wtsT4UMbWABJAWu4tzScVFI/L5z8vnzvom19ajKT5OMXWACyOSNWUYLW7u6hs6uTaczOzFV0qomsA0iYEifbakoI0n1Ty3EYvD726ncPBxPqGSIuRNB+n2IVBkLeogzQfSF0NYP+hINlZinH5ifXRlxYjaT6iawCBQOIjwlJJfX29aYUYpPmkkqJRub03f3eP80VjpcVImo9TjGQA0noClpQk1ynFbaT5pBqtNTet3sjydc5X15EWI2k+TrFjAYD29nbTCjFI80k1SinmlRQya9Jox/tIi5E0H6cYeWiRNiHI6NHOv3jpQJpPOrj704lNqiktRtJ8nGLkTjTR8BiP7u5u0woxSPNJF1pr3vC00dI58FsiaTGS5uMUWUWxIcJh541P6UCaT7rYezDAbc9s4en3GwfcVlqMpPk4xT4CAAUFBaYVYpDmky4mjs3juWXnclZJ4YDbSouRNB+n2EZAYP/+5Mamu400n3RSfmoROdlZAz4mSouRNB+nGMkApE3MOHnyZNMKMUjzSTdVTZ387cPvxB0pKC1G0nycYlcGAnbu3GlaIQZpPulm4tiR5GRnsc/X//dEWoyk+TjFDgYi0oAjqV1Cmg+krivwYJEWI2k+orsCSxsOXFlZaVohBmk+pgiGwqz/aN8J/yYtRtJ8nGIkA5DWYrpw4ULTCjFI8zHFqtd28KWfb+DjA8fPISktRtJ8nGJrAMibzEGajymuP/80nrjxbCaNPX7SEGkxkubjFNsGYHGEtDYAS3xsG0ACRJeCkoI0H9M8+e5OHvyDJ+b/pMVImo9TjPQElNYGUFZWZlohBmk+pmncd4iWDj/hsO5dYlxajKT5OMVIDUDawiAej2fgjdKINB/TLL98Lk/ceHbvzQ/yYiTNxyl2QhBg+vTpphVikOZjmhHZka9p5+Eg7b7IbFLSYiTNxylGMgBpQyd3795tWiEGaT4S8Ad7+NTKt1j5amQJLmkxkubjFCNtAJImTwQoKioyrRCDNB8J5Odm8x+XzWb+1EJAXoyk+Tgl2dWBf6iU8iilqpVSLyilCp3sJ23stLS3EtJ8pPCPZ09l1iljAHkxkubjlGQfAf4EnKm1PgvYAdyTvFL6kdSHG+T5SGKfL8B3XtzKHp+sx8hMTbNklwd/VWsdHdy/HnA0NapSauCN0oi04cnSfCQR7AmzdnMTNR/LWlwmU9PMtZ6ASqmXgOe11r/o5+/LgGVHP54JfOjKid1hPCBpWldpPiDPyfrEZ5bWesxAGw2YASilXgNOOcGf7tVarzu6zb3AIuBz2kGOopTa7KSbYrqwPgMjzcn6xMepz4DN8VrrpQOc6AbgcmCJk5vfYrHIIan3cUqpy4BvAhdprTOzGdRiGcYk23T538AY4E9KqUql1E8d7vd4kud1G+szMNKcrE98HPkYGQ5ssVhkkJkvLy0WiyvYDMBiGcYYywCUUiuOdiGuVEq9qpQyOrH6YLs1p9DnH5RSW5VSYaWUsddLSqnLlFLblVL1Sqn/MOXRx+dJpVSbUkpEPxKl1FSl1BtKqW1H0+suwz55SqmNSqmqoz7fjbuD1trIDzC2z+93Aj815XLU4RJgxNHffwD8wLDPHGAW8CawyJBDNtAAnA7kAlXAXMNxuRBYCHxo0qOPTzGw8OjvY4h0iTcWI0ABo4/+ngNsAM7rb3tjNQCtdVefj6MAo62RepDdmlPos01rvd2kA3AOUK+1/khrHQR+BVxpUkhr/TYgZh0urXWr1rri6O8HgW3AFIM+WmvtO/ox5+hPv/eW0TYApdT3lFJNwBeB5SZdjuFm4H9NSwhgCtDU53MzBr/c0lFKnQYsIFLqmvTIVkpVAm3An7TW/fqkNANQSr2mlPrwBD9XAmit79VaTwWeBf41lS5OfI5ucy8QOupk3McwJxq1Zd8bnwCl1Gjgt8BXj6ndph2tdY/WuoxILfYcpdSZ/W2b0pk59ADdiPvwS+AV4Nsp1BnQJ93dmhOIjymagal9PpcAmTn1TQpRSuUQufmf1Vr/zrRPFK11p1LqTeAy+hl8Z/ItwBl9Pn4WMDqrYp9uzZ/VtltzlE3AGUqp6UqpXOALwIuGnUShImPbnwC2aa1/JMBnQvQNllIqH1hKnHvLWE9ApdRvibRyh4FdwJe11i1GZCI+9cBIILoY3Xqt9ZcN+lwN/BiYAHQClVrrSw14fAZYReSNwJNa6++l2+EYn+eAxUSG3+4Bvq21fsKgzwXAO0ANke8ywLe01r835HMWsIZIemUBa7XW9/e7vakMwGKxmMf2BLRYhjE2A7BYhjE2A7BYhjE2A7BYhjE2A7BYhjE2A7BYhjE2A7BYhjH/H+fvcya1+Od7AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, ax = plt.subplots()\n",
    "x=np.linspace(-6,6,100)\n",
    "\n",
    "options = {\"head_width\":0.1, \"head_length\":0.2, \"length_includes_head\":True}\n",
    "\n",
    "ax.arrow(0,0,1,2,fc='b',ec='b',**options)\n",
    "ax.arrow(0,0,-0.6,1.2,fc='r',ec='r',**options)\n",
    "ax.plot(x,-2*x,ls=':')\n",
    "\n",
    "ax.set_xlim(-3,3)\n",
    "ax.set_ylim(-2,4)\n",
    "ax.set_aspect('equal')\n",
    "ax.set_xticks(np.arange(-3,4,step = 1))\n",
    "ax.set_yticks(np.arange(-2,5,step = 1))\n",
    "\n",
    "ax.text(0.8,2.2,'$X$')\n",
    "ax.text(-1,1.2,'$BX$')\n",
    "ax.text(0.85,-1,'$\\mathcal{C}(B)$')\n",
    "\n",
    "ax.axvline(color='k',linewidth = 1)\n",
    "ax.axhline(color='k',linewidth = 1)\n",
    "\n",
    "ax.grid(True,ls=':')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In this example we can again determine the eigenvectors and corresponding eigenvalues by considering the transformation represented by the matrix.  Any vectors in $\\mathcal{C}(B)$ will be unchanged by the projection.  Again we will label one such vector as $V_1$.\n",
    "\n",
    "$$\n",
    "\\begin{equation}\n",
    "V_1 = \\left[ \\begin{array}{r} -1  \\\\ 2 \\end{array}\\right] \\hspace{1cm}\n",
    "\\end{equation}\n",
    "$$\n",
    "\n",
    "Since $BV_1 = V_1$ the eigenvalue corresponding to $V_1$ is $\\lambda_1 = 1$.\n",
    "\n",
    "The other eigenvalue for this matrix may not be as easy to discover by examining the directions of various choices of $X$ and the corresponding directions of $BX$.  If we closely consider the definition of eigenvalues however, we realize that if there exists any nonzero vector $V$ such that $BV=0$, then that vector must be an eigenvector with zero as the corresponding eigenvalue.  To make the statement in terms of the fundamental subspaces, zero is an eigenvalue of the matrix $B$ if there are any nonzero vectors in $\\mathcal{N}(B)$.  \n",
    "\n",
    "In this example, since $B$ represents an orthogonal projection, $\\mathcal{N}(B)$ is the subspace orthogonal to $\\mathcal{C}(B)$.  We can label any such vector as $V_2$ \n",
    "\n",
    "$$\n",
    "\\begin{equation}\n",
    "V_2 = \\left[ \\begin{array}{r} 2  \\\\ 1 \\end{array}\\right] \\hspace{1cm}\n",
    "\\end{equation}\n",
    "$$\n",
    "\n",
    "The eigenvalue corresponding to $V_2$ is $\\lambda_2 = 0$."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Exercises\n",
    "\n",
    "**Exercise 1:** Determine the eigenvalues and corresponding eigenvectors of the following matrix by considering the transformation that it represents.  Check your answers with a computations.\n",
    "\n",
    "$(a)$\n",
    "\n",
    "$$\n",
    "\\begin{equation}\n",
    "A = \\left[ \\begin{array}{cc} 1 & 0 \\\\ 0 & -3 \\end{array}\\right]\n",
    "\\end{equation}\n",
    "$$\n",
    "\n",
    "$(b)$\n",
    "\n",
    "$$\n",
    "\\begin{equation}\n",
    "B = \\left[ \\begin{array}{cc} 1 & 1 \\\\ 0 & 1 \\end{array}\\right]\n",
    "\\end{equation}\n",
    "$$\n",
    "\n",
    "$(c)$\n",
    "\n",
    "$$\n",
    "\\begin{equation}\n",
    "C = \\left[ \\begin{array}{cc} \\cos(\\frac{\\pi}{2}) & -\\sin(\\frac{\\pi}{2}) \\\\ \\sin(\\frac{\\pi}{2}) & \\cos(\\frac{\\pi}{2}) \\end{array}\\right]\n",
    "\\end{equation}\n",
    "$$\n",
    "\n",
    "$(d)$\n",
    "\n",
    "$$\n",
    "\\begin{equation}\n",
    "D = \\left[ \\begin{array}{cc} -0.6 & -0.8 \\\\ -0.8 & 0.6 \\end{array}\\right]\n",
    "\\end{equation}\n",
    "$$"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "**Exercise 2:** Find the eigenvalues and the corresponding eigenvectors of the following matrix $R$ that represents the reflection transformation about the line $ x_1 = x_2 $.\n",
    "\n",
    "\n",
    "$$\n",
    "\\begin{equation}\n",
    "R = \\left[ \\begin{array}{cc} 0 & 1 \\\\ 1 & 0 \\end{array}\\right]\n",
    "\\end{equation}\n",
    "$$"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "## Code solution here"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "**Exercise 3:** Find a matrix that represents a vertical and horizontal stretch by a factor of $2$. Then, find the eigenvalues and the eigenvectors associated with those eigenvalues. (You may have to take a look at the Planar Transformations section  [Planar Transformations](Planar_Transformations.ipynb))."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "## Code solution here"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "**Exercise 4:** Find a matrix that represents reflection about the $x_1-$axis and find its eigenvalues and eigenvectors."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "## Code solution here"
   ]
  }
 ],
 "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.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}