Transformations_for_Car.ipynb

{
  "nbformat": 4,
  "nbformat_minor": 0,
  "metadata": {
    "colab": {
      "name": " Transformations for Car",
      "provenance": [],
      "include_colab_link": true
    },
    "kernelspec": {
      "name": "python3",
      "display_name": "Python 3"
    }
  },
  "cells": [
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "view-in-github",
        "colab_type": "text"
      },
      "source": [
        "<a href=\"https://colab.research.google.com/github/KapilKhanal/Advanced-Linear-Algebra-projects/blob/master/Transformations_for_Car.ipynb\" target=\"_parent\"><img src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open In Colab\"/></a>"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "ubqt1c0MDF-X",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "import numpy as np\n",
        "import matplotlib.pyplot as plt"
      ],
      "execution_count": 0,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "1nVTY4zaDNLG",
        "colab_type": "text"
      },
      "source": [
        "In this worksheet, we will generate matrix functions and use them to verify our solution to the car driving problem.\n",
        "The car will be represented by a vector $$\\vec{c} = \\left( \\begin{array}{c} x \\\\  y \\\\  1 \\\\ \\end{array} \\right)$$\n",
        "and the headlights will be represented by a vector\n",
        "$$ \\vec{l} = \\left( \\begin{array}{c} x \\\\ y \\\\ 1 \\\\ \\end{array} \\right)$$\n",
        "which will extend one unit in front of the car.\n",
        "\n",
        "##\n",
        "\n",
        "We will use the car and light definition to practice function declarations.  To place the car at (3,0) on the x-axis, we could specify:"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "_t_f3feYFzdY",
        "colab_type": "code",
        "outputId": "4d1aeb7a-672f-4ef9-a02e-22f9bbf0112a",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 153
        }
      },
      "source": [
        "import pprint\n",
        "car = np.array([3,0,1]) \n",
        "pprint.pprint(car)\n",
        "car = np.mat('3,0,1')\n",
        "pprint.pprint(car)\n",
        "car = np.mat('3,0,1').T\n",
        "pprint.pprint(car)\n",
        "car = np.matrix([[3],[0],[1]])\n",
        "pprint.pprint(car)\n"
      ],
      "execution_count": 16,
      "outputs": [
        {
          "output_type": "stream",
          "text": [
            "array([3, 0, 1])\n",
            "matrix([[3, 0, 1]])\n",
            "matrix([[3],\n",
            "        [0],\n",
            "        [1]])\n",
            "matrix([[3],\n",
            "        [0],\n",
            "        [1]])\n"
          ],
          "name": "stdout"
        }
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "OdtEsY3TETwL",
        "colab_type": "text"
      },
      "source": [
        "Instead, we will make it easy to specify where the car is by letting us define the car as a function of two inputs, the $x$ and $y$ location, like so:  $$location(x,y) = \\left( \\begin{array}{c} x \\\\ y \\\\ 1 \\\\  \\end{array} \\right)$$\n",
        "\n",
        "so that we can specify\n",
        "\n",
        "$$ car = location(3,0)   \\hspace{2cm} lights = location(3,1)$$\n",
        "\n",
        "by modifying the code block below:\n",
        "\n",
        "```\n",
        "def square(a):\n",
        "       \"\"\"Return the square of a.\"\"\"\n",
        "       return a ** 2\n",
        "```\n"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "hHdMlbyMGZG3",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "def location(a,b):\n",
        "       \"\"\"creating a vector at location (a,b,1)\"\"\"\n",
        "       return np.matrix([[a],[b],[1]])\n"
      ],
      "execution_count": 0,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "OrmBpkI990FS",
        "colab_type": "code",
        "outputId": "7f65a446-b63c-4c77-cae8-2f60b1fd18fd",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 119
        }
      },
      "source": [
        "#Verify that it works to define the car and the light direction:\n",
        "\n",
        "car = location(3,0)\n",
        "print(car)\n",
        "light = location(3,1)\n",
        "\n",
        "light"
      ],
      "execution_count": 18,
      "outputs": [
        {
          "output_type": "stream",
          "text": [
            "[[3]\n",
            " [0]\n",
            " [1]]\n"
          ],
          "name": "stdout"
        },
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "matrix([[3],\n",
              "        [1],\n",
              "        [1]])"
            ]
          },
          "metadata": {
            "tags": []
          },
          "execution_count": 18
        }
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "yZKKbf1UHaQw",
        "colab_type": "text"
      },
      "source": [
        "In your homework, you constructed several matrices by combining translation matrices and rotation matrices\n",
        "\n",
        "$$\n",
        "T_{c,d} = \\left( \\begin{array}{ccc} 1 & 0 & c \\\\ 0 & 1 & d \\\\ 0 & 0 & 1 \\\\ \\end{array} \\right) \\hspace{2cm} R_\\phi = \\left( \\begin{array}{ccc} \\cos \\phi & -\\sin \\phi & 0 \\\\ \\sin \\phi & \\cos \\phi & 0 \\\\ 0 & 0 & 1 \\\\ \\end{array} \\right)\n",
        "$$\n",
        "\n",
        "into complicated matrices which could rotate a car at a fixed position.  We're going to verify your results by combining matrices and checking the effect on the car's direction.\n",
        "\n",
        "First, let's define the two basic matrices $T(a,b)$ and $R(\\phi)$:"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "WJAryfZ0FARN",
        "colab_type": "code",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 68
        },
        "outputId": "bd7c03c7-098c-4b2f-e46c-b40b5d0a3e92"
      },
      "source": [
        "def T(c,d):\n",
        "  \"\"\"Return a 3x3 matrix which translates x by a and y by b\"\"\"\n",
        "  return np.matrix([[1,0,c], [0,1,d],[0,0,1]])\n",
        "\n",
        "def R(phi):\n",
        "  \"\"\"Return a 3x3 matrix which translates around the origin by angle phi\"\"\"\n",
        "  phi = np.deg2rad(phi)\n",
        "  return np.matrix([[np.cos(phi), -1*np.sin(phi), 0], [np.sin(phi), np.cos(phi), 0], [0, 0, 1] ])\n",
        "\n",
        "import math\n",
        "T(5,1)\n",
        "R(90)\n",
        "T(3,4) * T(3,4)"
      ],
      "execution_count": 19,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "matrix([[1, 0, 6],\n",
              "        [0, 1, 8],\n",
              "        [0, 0, 1]])"
            ]
          },
          "metadata": {
            "tags": []
          },
          "execution_count": 19
        }
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "8DNVwH4RNDDR",
        "colab_type": "text"
      },
      "source": [
        "Be Careful!!!  Just because something looks like a matrix does not mean it is a matrix!!!   In your group, develop some tests to make sure we have the results we expect.  Each group contributes one idea:"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "imccf5kBMx-Q",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        ""
      ],
      "execution_count": 0,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "EA9kTmUoRgHn",
        "colab_type": "text"
      },
      "source": [
        "Here's some documentation on the difference between arrays and matrices:\n",
        "\n",
        "https://docs.scipy.org/doc/numpy-1.13.0/reference/arrays.classes.html#index-0"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "y4qNC33EOR9F",
        "colab_type": "text"
      },
      "source": [
        "We'll test the result  of the translations on our original car and light position, which moved the car from $(3,0)$ to $(3,0)$ but from pointing 'north' to pointing 'east'  (e.g. pointing toward $(4,0)$ to $(3,1)$  )\n",
        "\n",
        "\n"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "JcgVR64fO8YV",
        "colab_type": "code",
        "outputId": "3e5d260e-b7a0-4be9-c0b2-4a76cac80946",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 68
        }
      },
      "source": [
        "car = location(3,0)\n",
        "light = location(3,1)\n",
        "\n",
        "new_Car = T(3,3) * R(270)  * car\n",
        "new_Light = T(3,3) * R(270)  * light\n",
        "\n",
        "\n",
        "new_Car"
      ],
      "execution_count": 40,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "matrix([[3.],\n",
              "        [0.],\n",
              "        [1.]])"
            ]
          },
          "metadata": {
            "tags": []
          },
          "execution_count": 40
        }
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "SxXmwY_YQIlc",
        "colab_type": "code",
        "outputId": "e15dcfe8-c649-4461-b9cb-083e3539c656",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 279
        }
      },
      "source": [
        "ax = plt.axes()\n",
        "\n",
        "ax.arrow(car[0,0], car[1,0],new_Car[0,0]- car[0,0],new_Car[1,0]- car[1,0], head_width=0.5, head_length=0.5, fc='lightblue', ec='black')\n",
        "ax.arrow(car[0,0], car[1,0],light[0,0]- car[0,0],light[1,0]- car[1,0], head_width=0.5, head_length=0.5, fc='lightblue', ec='red')\n",
        "\n",
        "plt.grid()\n",
        "\n",
        "plt.xlim(-5,5)\n",
        "plt.ylim(-5,5)\n",
        "\n",
        "plt.title('Initial Car Position (black) and final car position (red) ?',fontsize=10)\n",
        "\n",
        "plt.savefig('Initial Car Position (black) and final car position (red).png', bbox_inches='tight')\n",
        "plt.show()\n",
        "plt.close()"
      ],
      "execution_count": 42,
      "outputs": [
        {
          "output_type": "display_data",
          "data": {
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXIAAAEGCAYAAAB4lx7eAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0\ndHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAXfElEQVR4nO3de5hcdX3H8fcnCQhkEyiELpWlBrxA\nUQTMglZUdiFYQCJU0keIQqFipIKEmz4gPopaKliVolARsVKEukK4lquhsnIRSLLkgiFAQ4wCVS4K\nlA1Gcvn2j/ktMywz2ctM9swv+3k9zz6ZOWfO73zPd85+9jdnZjeKCMzMLF9jii7AzMzq4yA3M8uc\ng9zMLHMOcjOzzDnIzcwy5yA3M8vcqAxySb2DeMylknZNtz/fb90vhrsPSdtJ6pL0uKQeSbdIettg\na68yXrekRyUtknSvpJ2HMcbxko5Ot4+R9MaKda/2oV6SDpP0xXT7MknTqzymQ9JNwxx/haRJ/ZZt\nKukuSeOGV/WA+6z1PJ8kaamkKyV9WNIZjd5Hs5H0FUlT0+2TJW1Rse4WSVs1aD8n952vdYyxQtKk\n9Z0fknaR9AtJD0n6ef9zq6lExKj7Ano35ONrbQMIuA84vmLZ7sD7BzmmgDH9lnUD7en2TODGOnvz\n6ngboO+/ACal25cB06s8pgO4aZjjr+gbv9/yLwEfG8lzCXgEaNuQ+6hjvHEboheDeS4aUTuwuNox\nDOW4KuurdX4AuwA7pdtfA76wofs23K9ROSPvk2Z/3ZJmS3okzZ6U1nVLapd0LrC5pIWSrkzretO/\nLZL+W9KD6af2oQPsshNYHREX9y2IiEURcXetsSRNTjPuy4FfAjusZ/y7gLek7faXtCCN9e+S3pCW\nnyvpYUmLJX0jLTtb0ulphtwOXJmOd/O+PqTHHZnG+6Wk8yr62CvpnPSq4H5JrVV6/TbgTxHxXMXi\nqZLmS3pM0iFVttlb0n3pOH7R92pD0lhJ30h1LJb0mX7bbS7pVkmfTIuuBz5WrWGSrk+vjJZImjnQ\nMUnaMdX0kKR/qjHmxcBOwK2STkmvci5M6y6T9O10PMv7XpUM41xC0tHp+BdJ+lFaNk3SA6lnd1TU\nfbakH0m6F/hRv3E60qz05nSuXSxpTFr3uuc89f+ytOwhSadUHNt0SScBbwTulHRnWvfqqyVJp6Zt\nfynp5LRsskqvYL6fnoufStq8ymHvBzwYEWvSdt2S/lXSfGCWpG0lXSNpXvraJz1umzTmEkmXUpoU\n9al6fkTEIxGxPN19A7BqoOekMEX/JCniizTDoTT7exFoo3SZ6T7gfWldN+WZbm+N7ccBE9PtScAy\nQNW2SctOAs6vUVPVsYDJwDrgPTW2q6zzs8BPgM2AJ4C3peWXAycD2wCPVtS4Vfr3bOD0/uNV3qf0\njfkbYNtU68+Aw9JjApiWbn+dKjMX4FjgmxX3LwNuS31/K/BkqruDNCMHJpJmWcBU4Jp0+x+B2RXr\ntk7/rkj9ugM4umJfY4Fna/Svb9vNKf2g3GZ9xwTc2Dc2cEK157milr4Z3zHAhRXHfXU67l2BZcM8\nl94OPFaxj77j+LOK7Y7r63l6jnuAzauM1UEppHZKvZoDTK/1nANTgDkV229VcWzT+x9/5f207UPA\neKAFWALsmZ63NcAe6fFXAR+vUuuXgc/0Oz//reL+f1L+Hv5LYGm6/W3gi+n2h9Lz29e7mudHWv83\nwNK+42zGr1E9I0/mRsSTEbEOWEjphBosAf8saTGl8NgeeN1stAFj/Toi7l/PtldKWgjsA5wO7Az8\nKiIeS+v/A/gApR9aq4AfSPoI8PIQ6tsL6I6IZ6M0G7oyjQnwCtB3XbuH6j38C+DZfsuuioh1EfE/\nwHJKL2UrbQlcLemXwPmUwgtKof69VAcR8YeKbW4AfhgRl/ctiIi1wCuSJlSp6yRJi4D7Kb3aeesA\nx7QP8ON0+zUz2yG4Ph33w5Sf46GeS/sBV0d6hVPRgzbgdkkPUfrB/vaKbW6MiD/WGG9uRCxPvfox\n8D5qP+fLgZ0kfUfSgcD/DeHY3wdcFxErI6IXuBZ4f1r3q4hYmG4P5Tz6ScXtqcCF6fvhRmCipJZU\n9xUAEXEz8HzfBus7P9Irkx8AH46IF4ZwnCPKQQ5/qri9ltLMY7A+Rmm2MiUi9gCepjSrrGUJpRnJ\nUMdaOVAdEbFHRBwWEU/UelD6Ztyb0mz2EEoz4kZYHWnqQu0e/pHX96b/H/rpf/+rwJ0R8Q5gWpXt\nq7kXOFCS+i1/3UtjSR2UvvH/OiJ2BxZU7GN9x1TvHyiqPOf66hzquVTLdyjN/ncDPtVvjPWdRwM9\nF+UVEc9Tem+nGzgeuHQYdVYzmO/FaudR5XGNofTqdY/0tX36gTGQWpdO3gi8mCYbTctBPjirJW1S\nZfmWwDMRsVpSJ/CmAcb5GfCGftdi3ynp/cMYa30eBSZLeku6fxTw8zQz2TIibgFOofTN2N9LQLWZ\n61xgX5Xe6R8LHAn8fAg1LSVdv6/wd5LGSHozpZf1j/ZbvyXwVLp9TMXyOcCnlD5pIGnrinVfpDTb\nuqhvgaRtgOciYnWV8Z+PiJcl7QK8ZxDHcS9wRLpd9br7MA3nXPq7dGyVPajs2d8PYf97p+v/Y4CP\nAvdQ4zlP17rHRMQ1wBeAd1UZr9Z5dDdwmKQtJI0H/jYtG6xq51GlnwKvvmciaY908y5gRlp2EKVL\nUH2PqXV+QOlcOm0I9RXCQT44lwCLld7srHAl0J5exh5N6ZMKNaUZ3t9SepPvcUlLKL0b/ruhjjXA\nflZRuiZ9dRpvHXAxpW+sm9LL93uAU6tsfhlwsdKbnRVj/hY4A7gTWAT0RMQNQyjrLmDPfjPl31AK\ni1spfZKn/4zo68DXJC3gtbOzS9O2i9NlkRn9tptF6Q3qr6f7ncDNVWq6DRgnaSlwLqXLKwOZBZyQ\n+rr9IB4/WEM9l5YA51AK1kXAt9Kqsyk97z3AczU2r2YecCGloPwVpcsftZ7z7YHudPniCuDMKuNd\nAtzW92ZnRd0PUjrH5gIPAJdGxIIh1Hkr5Ut61ZxEqY+LJT1M6RUDlK6tfyB9z32E0vnTp9b5AaUf\njMcNob5CqPzq0WzDknQB8F8RcccI7/da4IyK9wysQrrEdHpEvO6TQ81I0nXA5xp1uWNjOD88I7eR\n9M/AFgM+qoEkbUrpzcVsv0ntdc6g9KZn3TaW88MzcjOzzHlGbmaWOQe5mVnmNsgfEhrIpEmTYvLk\nyUXs+lUrV65k/PjxhdbQLNyLMveizL0oa5Ze9PT0PBcR2/ZfXkiQT548mfnz5xex61d1d3fT0dFR\naA3Nwr0ocy/K3IuyZumFpF9XW+5LK2ZmmXOQm5llzkFuZpY5B7mZWeYc5GZmmXOQm5llzkFuZpY5\nB7mZWeYc5GZmmXOQm5llzkFuZpY5B7mZWeYc5GZmmXOQm5llzkFuZpY5B7mZWeYc5GZmmWtYkEsa\nK2mBpJsaNaaZmQ2skTPyWcDSBo5nZmaD0JAgl9QGfAi4tBHjmZnZ4Cki6h9Emg18DZgAnB4Rh1R5\nzExgJkBra+uUrq6uuvdbj97eXlpaWgqtoVm4F2XuRZl7UdYsvejs7OyJiPb+y8fVO7CkQ4BnIqJH\nUketx0XEJcAlAO3t7VH0/0jdLP8rdjNwL8rcizL3oqzZe9GISyv7AB+WtALoAvaTdEUDxjUzs0Go\nO8gj4syIaIuIycARwM8i4uN1V2ZmZoPiz5GbmWWu7mvklSKiG+hu5JhmZrZ+npGbmWXOQW5mljkH\nuZlZ5hzkZmaZc5CbmWXOQW5mljkHuZlZ5hzkZmaZc5CbmWXOQW5mljkHuZlZ5hzkZmaZc5CbmWXO\nQW5mljkHuZlZ5hzkZmaZc5CbmWXOQW5mljkHuZlZ5hzkZmaZc5CbmWXOQW5mljkHuZlZ5hzkZmaZ\nc5CbmWXOQW5mljkHuZlZ5hzkZmaZc5CbmWXOQW5mljkHuZlZ5hzkZmaZc5CbmWXOQW5mlrm6g1zS\nDpLulPSwpCWSZjWiMDMzG5xxDRhjDXBaRDwoaQLQI2lORDzcgLHNzGwAdc/II+K3EfFguv0SsBTY\nvt5xzcxscBp6jVzSZGBP4IFGjmtmZrUpIhozkNQC/Bw4JyKurbJ+JjAToLW1dUpXV1dD9jtcvb29\ntLS0FFpDs3AvytyLMveirFl60dnZ2RMR7f2XNyTIJW0C3ATcHhHfGujx7e3tMX/+/Lr3W4/u7m46\nOjoKraFZuBdl7kWZe1HWLL2QVDXIG/GpFQE/AJYOJsTNzKyxGnGNfB/gKGA/SQvT18ENGNfMzAah\n7o8fRsQ9gBpQi5mZDYN/s9PMLHMOcjOzzDnIzcwy5yA3M8ucg9zMLHMOcjNrTitXsrqjk3UXXFB0\nJU3PQW5mzSeCNZ+cye833Zy1/3QO3H130RU1tUb8GVszs4Za993vsmrhIh74zxuYNO9+3n3EEYzr\n6YHttiu6tKbkGbmZNZd581j3xS9x7/nfY+3mW/D0B/bj8cM+ypqPHgFr1hRdXVNykJtZ8/j971kz\nfTrzzz6PlZN3enXxkhNO5YV1sPasLxRYXPNykJtZc1i3jjUzPsaKAz7E/x5w0GvXjR3L/V+/kDVX\nXEHccEMx9TUxB7mZNYW1X/0qLz3/Ig+dcmbV9a9svQ33ffO7rD3uOFi+fISra24OcjMrXNx+O2sv\n/h73ffO7xLjan8H4w57tLJl5Ems+cjisWjWCFTY3B7mZFeuJJ1h79NHcf953WPXnrQM+/PGjPsEz\nbW9izQknjkBxeXCQm1lxXnmFNYdP59GjjuO5d793cNtIzPvKN3jlrruIH/5ww9aXCX+O3MwKs/a0\n01j5xz/x232nMmHZY69bP2b1K6zbZNOq2y4+8bPs9ZmTGPuud8Huu2/oUpuag9zMCrPu6WfYYvWf\n2Pezn666fpNHHwFg9c67VN++rQ2tWMEYB7mZWTE2ueon63+ABJ/+NJtcdNHIFJQpXyM3M8ucg9zM\nLHMOcjOzzDnIzcwy5yA3M8ucg9zMLHMOcjOzzDnIzcwy5yA3M8ucg9zMLHMOcjOzzDnIzcwy5yA3\nM8ucg9zMLHMOcjNrOsuWLePcc88F4MUXXyy4mubnIDezwkUES5Ys4ewvf5m37/ZO3vPefbhr8VIA\nXnjhhYKra37+jyXMrBARwYIFC7h69myunn0NvStXsvcBB3HkGV9h5z3bGTt2LPz48qLLzEJDglzS\ngcAFwFjg0og4txHjmtnGZd26dcydO/fV8A7EXgcczCfPOZ+37LYHkoouMUt1B7mkscBFwAHAk8A8\nSTdGxMP1jm1m+Vu7di333HMPV82ezbXXXstmLRPYa+pBnHzB93nTzrs6vBugETPyvYFlEbEcQFIX\ncCjgIDcbxXp7eznl1NO4/obr+bNtW2mfehBnfv/HtL35rUMaZ93YsRuowo2HIqK+AaTpwIERcVy6\nfxTw7og4sd/jZgIzAVpbW6d0dXXVtd969fb20tLSUmgNzcK9KHMvyurtxdq1a1m2bBkrX36Z8RMm\nMn7ilmwxYQJjxgz+MxYvPb6MrbZrZcLWWw+7jkZolvOis7OzJyLa+y8fsTc7I+IS4BKA9vb26Ojo\nGKldV9Xd3U3RNTQL96LMvShrRC/2339/nn76aa677jp+cMn3mD93Lrv/9ftoP+BgpnRMZfyEievd\n/sJzz+NLnz218Oek2c+LRnz88Clgh4r7bWmZmRmtra0cf/zx3HnHHH694ld84ojpPHLnLfxj516c\nd/zR/OyaLl56/g9Fl5m1RgT5POCtknaUtClwBHBjA8Y1s43M1ltvzbHHHsttN9/M/z71JLM+eSxP\nzLubEz/4Xs75xJHc3nU5Lzz3bNFlZqfuSysRsUbSicDtlD5++O8RsaTuysxsozZx4kRmzJjBjBkz\nWLlyJbfddhs/uXo2J3/ra+y4y18xZerB9Pq3OgelIdfII+IW4JZGjGVmo8/48eM5/PDDOfzww1m1\nahVz5szhqqtn8+JzT9PW1lZ0eU3Pv9lpZk1ls802Y9q0aUybNq3oUrLhv7ViZpY5B7mZWeYc5GZm\nmXOQm5llzkFuZpY5B7mZWeYc5GZmmXOQm5llzkFuZpY5B7mZWeYc5GZmmXOQm5llzkFuZpY5B7mZ\nWeYc5GZmmXOQm5llzkFuZpY5B7mZWeYc5GZmmXOQm5llzkFuZpY5B7mZWeYc5GZmmXOQm5llzkFu\nZpY5B7mZWeYc5GZmmXOQm5llzkFuZpY5B7mZWeYc5GZmmXOQm5llzkFuZpa5uoJc0r9IekTSYknX\nSdqqUYWZmdng1DsjnwO8IyLeCTwGnFl/SWZmNhR1BXlE/DQi1qS79wNt9ZdkZmZD0chr5P8A3NrA\n8czMbBAUEet/gHQHsF2VVWdFxA3pMWcB7cBHosaAkmYCMwFaW1undHV11VN33Xp7e2lpaSm0hmbh\nXpS5F2XuRVmz9KKzs7MnItr7Lx8wyAci6RjgU8D+EfHyYLZpb2+P+fPn17XfenV3d9PR0VFoDc3C\nvShzL8rci7Jm6YWkqkE+rs5BDwQ+B+w72BA3M7PGqvca+YXABGCOpIWSLm5ATWZmNgR1zcgj4i2N\nKsTMzIbHv9lpZpY5B7mZWeYc5GZmmXOQm5llzkFuZpY5B7mZWeYc5GZmmXOQm5llzkFuZpY5B7mZ\nWeYc5GZmmXOQm5llzkFuZpY5B7mZWeYc5GZmmXOQm5llzkFuZpY5B7mZWeYc5GZmmXOQm5llzkFu\nZpY5B7mZWeYc5GZmmXOQm5llzkFuZpY5B7mZWeYc5GZmmXOQm5llzkFuZpY5B7mZWeYc5GZmmXOQ\nm5llzkFuZpY5B7mZWeYc5GZmmWtIkEs6TVJImtSI8czMbPDqDnJJOwAfBH5TfzlmZjZUjZiRnw98\nDogGjGVmZkOkiOHnr6RDgf0iYpakFUB7RDxX47EzgZkAra2tU7q6uoa930bo7e2lpaWl0BqahXtR\n5l6UuRdlzdKLzs7Onoho7798wCCXdAewXZVVZwGfBz4YES8OFOSV2tvbY/78+YMqfEPp7u6mo6Oj\n0BqahXtR5l6UuRdlzdILSVWDfNxAG0bE1BoD7gbsCCySBNAGPChp74j4XZ31mpnZIA0Y5LVExEPA\nn/fdH8qM3MzMGsefIzczy9ywZ+T9RcTkRo1lZmaD5xm5mVnmHORmZplzkJuZZc5BbmaWOQe5mVnm\nHORmZplzkJuZZc5BbmaWOQe5mVnmHORmZplzkJuZZc5BbmaWOQe5mVnmHORmZplzkJuZZc5BbmaW\nuQH/8+UNslPpWeDXI77j15oE+L+lK3EvytyLMveirFl68aaI2Lb/wkKCvBlIml/tf6MejdyLMvei\nzL0oa/Ze+NKKmVnmHORmZpkbzUF+SdEFNBH3osy9KHMvypq6F6P2GrmZ2cZiNM/Izcw2Cg5yM7PM\nOcgBSadJCkmTiq6lKJL+RdIjkhZLuk7SVkXXNNIkHSjpUUnLJJ1RdD1FkbSDpDslPSxpiaRZRddU\nNEljJS2QdFPRtVQz6oNc0g7AB4HfFF1LweYA74iIdwKPAWcWXM+IkjQWuAg4CNgVOFLSrsVWVZg1\nwGkRsSvwHuCEUdyLPrOApUUXUcuoD3LgfOBzwKh+1zcifhoRa9Ld+4G2IuspwN7AsohYHhGvAF3A\noQXXVIiI+G1EPJhuv0QpwLYvtqriSGoDPgRcWnQttYzqIJd0KPBURCwqupYm8w/ArUUXMcK2B56o\nuP8kozi8+kiaDOwJPFBsJYX6V0qTvXVFF1LLuKIL2NAk3QFsV2XVWcDnKV1WGRXW14uIuCE95ixK\nL62vHMnarPlIagGuAU6OiP8rup4iSDoEeCYieiR1FF1PLRt9kEfE1GrLJe0G7AgskgSlSwkPSto7\nIn43giWOmFq96CPpGOAQYP8Yfb9g8BSwQ8X9trRsVJK0CaUQvzIiri26ngLtA3xY0sHAZsBESVdE\nxMcLrus1/AtBiaQVQHtENMNfOBtxkg4EvgXsGxHPFl3PSJM0jtKbvPtTCvB5wIyIWFJoYQVQaWbz\nH8AfIuLkoutpFmlGfnpEHFJ0Lf2N6mvk9hoXAhOAOZIWSrq46IJGUnqj90Tgdkpv7l01GkM82Qc4\nCtgvnQsL04zUmpRn5GZmmfOM3Mwscw5yM7PMOcjNzDLnIDczy5yD3Mwscw5yM7PMOcjNzDL3/xzo\nDI7AJm/rAAAAAElFTkSuQmCC\n",
            "text/plain": [
              "<Figure size 432x288 with 1 Axes>"
            ]
          },
          "metadata": {
            "tags": []
          }
        }
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "oi8gu2VRTfc0",
        "colab_type": "text"
      },
      "source": [
        "# Homework 1:  Construct the transformation that allows the car to turn right by 90 degrees from location (a,b)\n",
        "Construct the single matrix $R_{a, b, 90}$ as a product of the simple matrices $T_{c,d}$ and $R_{\\phi}$ constructed above.\n",
        "\n",
        "Verify the result by applying your new matrix to the car located at various locations.\n",
        "\n",
        "Produce a plot that shows the result of the transformation on a car located in at least two different quadrants.\n"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "wIeLNqrMQKwR",
        "colab_type": "code",
        "outputId": "1ec66565-9789-40c6-ab0f-01e81eed1b6d",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 68
        }
      },
      "source": [
        "def RabBy90(a,b, T, R):\n",
        "  \"\"\"Return a 3x3 matrix which rotates the car at location a,b by 90 degrees.\"\"\"\n",
        "  # We have to pass in the T and R functions if we are going to use them in this function.\n",
        "  return T(a,b) * R(-90) * T(-a, -b)\n",
        "\n",
        "RabBy90(3,0, T, R)"
      ],
      "execution_count": 43,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "matrix([[ 6.123234e-17,  1.000000e+00,  3.000000e+00],\n",
              "        [-1.000000e+00,  6.123234e-17,  3.000000e+00],\n",
              "        [ 0.000000e+00,  0.000000e+00,  1.000000e+00]])"
            ]
          },
          "metadata": {
            "tags": []
          },
          "execution_count": 43
        }
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "vSDzmp-KYR0u",
        "colab_type": "code",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 279
        },
        "outputId": "e51819a3-061a-47c1-d3e9-82170c27a537"
      },
      "source": [
        "#car 1 location at (3,2) heading toward (3,3)\n",
        "#car 2 location at (-2, -3) heading toward (-1, -3)\n",
        "\n",
        "ax = plt.axes()\n",
        "\n",
        "#car 1 old position:\n",
        "ax.arrow(2, 2, 1.0, 1.0, head_width=0.5, head_length=0.7, fc='lightblue', ec='black')\n",
        "#car 1 new position:\n",
        "ax.arrow(2, 2, 2.0, -1, head_width=0.5, head_length=0.7, fc='lightblue', ec='red')\n",
        "\n",
        "#car 2 old position:\n",
        "ax.arrow(2, 2, 2.0, -1, head_width=0.5, head_length=0.7, fc='lightblue', ec='red')\n",
        "#car 2 new position:\n",
        "ax.arrow(2, 2, 2.0, -1, head_width=0.5, head_length=0.7, fc='lightblue', ec='red')\n",
        "\n",
        "\n",
        "plt.grid()\n",
        "\n",
        "plt.xlim(-5,5)\n",
        "plt.ylim(-5,5)\n",
        "\n",
        "plt.title('Showing Car rotates 90 degrees ?',fontsize=10)\n",
        "\n",
        "plt.savefig('Showing Car rotates 90 degrees.png', bbox_inches='tight')\n",
        "plt.show()\n",
        "plt.close()"
      ],
      "execution_count": 23,
      "outputs": [
        {
          "output_type": "display_data",
          "data": {
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXIAAAEGCAYAAAB4lx7eAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0\ndHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAbWElEQVR4nO3de5hVZd3/8feH4SQMCogMJhZWWnh5\nSiY1zRoseNJQ0qw8pNlBKvOMeYgy86nsQS8rO9iPqDD0CTXlEJmC4CRqFCcFQfABPCRapIQwoMAM\n398fa8EepxkYmD2z9mI+r+va16y99lr3+q574LPvfe+19ygiMDOz/OqQdQFmZtYyDnIzs5xzkJuZ\n5ZyD3Mws5xzkZmY55yA3M8s5B/keQNIoSYslLZT0pKRj0/XPS+pTxONUSrqtiO2dL+lpSYskLZB0\nVbHa3sExB0g6p1jbNfOYJ0man57rHZI6pusl6TZJy9Pf3dHNaOuGtuin3SXpK+m/xWcl3ZB1Pe2F\ngzznJH0AGAYcHRFHAB8F/t4ax4qIuRFxaTHaknQycDkwNCIOB44DXt+F/TvuzmPAAKA5Ad3c7XZI\nUgfgDuCsiDgMeAH4XPrwycDB6W0EcHtLj7cLde2oj1piOfA+4HDgc5L6t9JxrB4Hef7tD7waEZsA\nIuLViHi53uOXpKPBRZLeCyCpt6RJ6ShwtqQj0vWLJPVMR4qvSTo/Xf9bSUMkVUmamq67QdKvJVVL\nWilpe8BL+pakZZIek/S7JkaQ1wFXbas1IjZFxC/T/S+UNEfSU5Luk9QtXT9O0i8k/RUYXb8xSRdI\nmiJpJjAjPYeb6434P5Nu+gPgxPSVyxXpyHtW2kfzJR3fxHZlaXtz0n77cnrc/SU9mm73tKQTG5zn\nvsDmiHg2vT8d+GS6PBz4bSRmAz0l7d+wo9JXXM9Kegx4T73175L0oKR56Tm8t9762el5f1dSTbq+\nKt1uCrAkXfdZSX9L6/9/ksrS9UMl/SXtk3sllafrfyBpSdoHtzSsNSIejojNgICOwOZGfvdWbBHh\nW45vQDnwJPAs8HPgw/Ueex64JF2+CBibLv8E+Ha6fBLwZLr8C+DjwGHAHOCX6fr/A7oDVcDUdN0N\nwBNAF6AP8BrQCXh/Wk9XoEe671WN1L0G2KeJc9q33vJ3653DOGAqUNbIPhcALwG90/ufJAnNMqAC\neJHkSW/7OaTbdQO6pssHA3PT5YbbjQC+mS53AeYCBwEjgVHp+jKgR4O6RDIKr0zv/xhYlC5PBT5Y\nb9sZ27art24QsCitc2+SEe9V9bY/OF0+FphZr92z0+WvADX1zmkDcFB6fyDwB6BTev/nwPnp7/NR\noHu6/hrgepInpWWA0vU9d/Dv8rfAzVn//2gvt9Z6eWVtJCJqJA0CTgQGA3dLujYixqWb3J/+nAec\nkS5/kHRUGBEzJe0raW9gFvAhkuC5HRgh6QDg3xGxQVLDw/8xklcCmyStJgnME4DJEfEm8KakP+zG\naR0m6btAT5InqofqPXZvRNQ1sd/0iFhT7xx/l277T0l/JnmSWddgn07ATyUdBdQBhzTR9lDgCEln\npvf3IQn+OcCvJXUCJkXEk/V3ioiQdBbwQ0ldgGnpcZrrRGBiRGwESEfTpCPk44F76/1euqQ/PwB8\nIl3+X6D+yPlvEfFcuvwRkieKOWkbewGrSaa5DgUeT9d3Bv5CMvX1JvCr9JXZ1MYKlnQayZPmBbtw\nntYCDvI9QBpW1UC1pEUkc7Dj0oc3pT/r2Pnv+1Hga8DbgVHA6cCZJAHfmE31lpvTfn2LSUJkZiOP\njQM+ERFPSbqAZCS5zYYdtLmjx5pyBfBP4EiSqcY3m9hOJK8MHvqPB6QPkbySGSfp1oj4bf3HI+Iv\nJIGMpKEUnixWAQfW27R/uq45OgBrI+KoZm6/Tf0+EnBHRFxXfwNJp5I8KZ7dcGdJx5A8AZwJXEzy\niq6hI4BpEbF1F2uz3eQ58pyT9B5JB9dbdRTJiHpHZgHnpvtXkcyxr4uIv5O8rD44IlYCjwFXkQR8\ncz0OnCqpazpqHNbEdjcBN0vql9bRWdKX0sd6AK+ko9xzd+HY9c0CPpPObe9H8krjb8D6tP1t9gFe\nSUPnPJLpERrZ7iHgq2lNSDpEUndJ7wD+Gcn8/ljgP648kdQ3/dmFZJriF+lDU4Dz0/n844DXI+KV\nBrs/CnxC0l6SegCnAkTEOuA5SZ9K25akI9N9ZlOYhz9rB300AzizXn290/OZDZwg6d3p+u7p+ZaT\nTIc9QPIEeGQT7U5Kz83aiEfk+VcO/ERST6CWZA51xE72uYFkOmAhsJHCVRQAf6UQZrNIAvex5hYT\nEXPSl/8LSUa6i2jkapSIeEBSBfCwktfvAfw6ffhbaR3/Sn/2aLh/M0wkmWJ4Km376oj4h6TXgDpJ\nT5GM/H8O3Kfkjd0HKYxYFzbY7sckV7LMT+v9F8n0RRXwdUlbgBqSOeaGvi5pGMnA6faI2PYq5AHg\nFJLf2Ubg8w13jIj5ku5Oz2M1yVTONucCt0v6JskU0YR0u8uBOyWNSs+p0auBImJJuu80JVfXbAG+\nFhGz01dCv0uffAC+SfLkNllSV5LR/JWNtUsyrbWRZD7d2sC2Ny3MikZSeTp3341kRDkiIuZnXVd7\nkfb7G/Xm58+OiOFZ12WtxyNyaw1jJB1KcuXKHQ7xNjeI5A1cAWuBL2Rcj7Uyj8jNzHLOb3aameWc\ng9zMLOcymSPv06dPDBgwIItDb7dhwwa6d++eaQ2lwn1R4L4ocF8UlEpfzJs379WI2K/h+kyCfMCA\nAcydOzeLQ29XXV1NVVVVpjWUCvdFgfuiwH1RUCp9IanRz4h4asXMLOcc5GZmOecgNzPLOQe5mVnO\nOcjNzHLOQW5mlnMOcjOznHOQm5nlnIPczCznHORmZjnnIDczyzkHuZlZzjnIzcxyzkFuZpZzDnIz\ns5xzkJuZ5ZyD3Mws54oW5JLKJC2QNLVYbZqZ2c4Vc0R+GfBMEdszM7NmKEqQS+oPfBwYW4z2zMys\n+RQRLW9E+j1wE9ADuCoihjWyzQhgBEBFRcWgCRMmtPi4LVFTU0N5eXmmNZQK90WB+6LAfVFQKn0x\nePDgeRFR2XB9x5Y2LGkYsDoi5kmqamq7iBgDjAGorKyMrP8idan8VexS4L4ocF8UuC8KSr0vijG1\ncgJwmqTngQnASZLuLEK7ZmbWDC0O8oi4LiL6R8QA4CxgZkR8tsWVmZlZs/g6cjOznGvxHHl9EVEN\nVBezTTMz2zGPyM3Mcs5BbmaWcw5yM7Occ5CbmeWcg9zM9kgRwcKFC1m5cmXWpbQ6B7mZ7TG2bNnC\nI488wqWXXcY7BhzEkUceyS233pp1Wa3OQW5mubZ+/XruvfdezjrnXPpW9OOiy6/kX3ThxDPO5oAD\nD+S7N96YdYmtrqjXkZuZtYVVq1YxefJk7p88hdlPPMHAoyt53+ChjP7ylezb72289s9XuO7Mk5l0\n/3307t0763JbnYPczEpeRPD0008zadIk7ps0medWrmDQhz7C0aecwQXfv41u5T22b1tXV8fPrrmU\nSy7+Gh/84AczrLrtOMjNrCTV1tYya9Ys7p80icmTp1C7dSuVJw3l9EuvZeCgY+jYqVOj+0365U8p\n71TGN0eNauOKs+MgN7OSsX79eh588EHumziJhx58kH4Hvp2jqoZw+W1jecchA5G0w/2Xzp/DtP/9\nDQvmzaOsrKyNqs6eg9zMMrdkyRIuveLKRue7m6vm9bX85OqLGTtmDP3792/FakuPg9zMMrd+/Xrm\nzZnDh049g/OuuZ4uXffapf0jgjHXf53Th5/G8OHDW6nK0uXLD80sc8ceeyzLlj5Dty0bufr0oSye\nM3uX9p9+z12sffnv3HrLLa1UYWnziNzMSkLfvn35/b33MHHiRC66+GLeVzWEc678xluuSGnMi/+3\njHtu+x+eeOwxunbt2kbVlhaPyM2spJx++uk8s3gxnTa+znmV7+Gv0//U5Lab3nyD20Z+ldE/+AHv\nfe9727DK0uIgN7OSc9NNNzFl0iQA7hp9Az+95lLW/fu1/9hu/OgbGXTkEXzxi19s6xJLioPczErG\nc889hyRGjx7NJZdcQkSw7JlnOHxAf64a/lEe/9MUIgKA2dMeYPETf2bsL8fs9LLEPZ3nyM0scxHB\n5z73OcaPHw8kH8F/29uSSw+7d+/Oj3/0Q84+6zNc8Pkv8Jc/TuK0L32Nsd+5lgem/oF99tkny9JL\ngkfkZpaphQsX0qFDB8aPH8/3v/99ImJ7iNd33HHH8dSTC/jo8cdw/XlncNXIKznuuOMyqLj0eERu\nZpmICIYMGcKMGTMAWLNmDb169drhPl26dOHG73yHa66+mm7durVFmbngEbmZtblZs2bRoUMHZsyY\nwdixY4mInYZ4fd27d2/38+L1eURuZm2mtraWww8/nKVLl1JeXs7q1avZa69d+xSn/SePyM2sTUyZ\nMoVOnTqxdOlSJk6cyPr16x3iReIRuZm1qjfeeIN+/fqxbt06DjnkEBYvXkzHjo6eYvKI3MxazW9+\n8xu6devGunXrqK6uZtmyZQ7xVuAeNbOiW7t27fY3LwcPHsyMGTP85mQr8ojczIpq9OjR20N8wYIF\nzJw50yHeyjwiN7OieOWVV7Z/kOfcc89l/PjxDvA24hG5mbXYFVdcsT3Ely9fzp133ukQb0MOcjPb\nbcuXL0cSP/rRj7jyyiuJCN71rndlXVa746kVM9tlEcHZZ5/N3XffDSTTKv369cu4qvbLI3Iz2yUL\nFiygQ4cO3H333dx8881EhEM8Yx6Rm1mzbN26laqqKmbNmgUklxj6K2RLQ4tH5JIOlPSIpCWSFku6\nrBiFmVnpqK6upqysjFmzZjFu3DgiwiFeQooxIq8FRkbEfEk9gHmSpkfEkiK0bWYZ2rJlCwMHDmTF\nihX06tWLl19+ud3+geNS1uIReUS8EhHz0+X1wDPAAS1t18yytXbtWjp37syKFSuYMmUKa9ascYiX\nqKLOkUsaALwP+Gsx2zWzNvb66/RcsYIAtg4cSIfBg7OuyHZA2/6QaYsbksqBPwPfi4j7G3l8BDAC\noKKiYtCECROKctzdVVNTQ3l5eaY1lAr3RYH7oqDm1Vcpf+GFt67s3BkGDoR29sVXpfLvYvDgwfMi\novI/HoiIFt+ATsBDwJXN2X7QoEGRtUceeSTrEkqG+6LAfVGwvS9qayM++cnYChH1b337RqxcmWmN\nbaVU/l0Ac6ORTC3GVSsCfgU8ExG3trQ9MysxZWXw+9+jbRF+4YUEwOrV8M53ggQ9e8LChVlX2m4V\n4wNBJwDnASdJejK9nVKEds2sFI0ZUwj1a69NQv311+HII5NQ79YNHn006yrblWJctfJYRCgijoiI\no9LbA8UozsxK3E03FUL9llsICd54Az784STUu3SBKVOyrnKP54/om1lxjByJtqYz6XfcQXTsCJs3\nw/DhICX3x43Luso9koPczIrv/PPRli1JqP/hD9C5M6qrg89/Pgn1Dh1g9Oisq9xjOMjNrHUNGwab\nNiWh/thj0K1bMh1zzTVJqEvJsu02B7mZtZ0TToANG5JQX7QIevZEkIzOt4X6hRdmXWXuOMjNLBuH\nHQb//ncS6i+8ABUVSaiPHVsI9dNPh7q6rCsteQ5yM8ve298O//hHEuqvvQYHHZSE+qRJyadIJaiq\ncqg3wUFuZqWld29YuTIJ9Y0bk5E7wJ//XAj1o49OrogxwEFuZqVsr72SufQIqK2F449P1i9YkFyj\nLsHBBycfSGrHHORmlg9lZfD444VQP+WU5FOly5cnXxEgQf/+yRRNO+MgN7P8KSuDP/6x8KnSc85J\nQn3VKth//yTU+/SBFSuyrrRNOMjNLP/uuqsQ6pdemoT6a6/Bu9+dhPree8O8eVlX2Woc5Ga2Z/nx\njwuhfuONyWWM69dDZWUS6nvtBdOn717bdXVs2v9tbN6nJ9x+e3HrbgEHuZntub71rcL3v/z858lX\nA7z5Jgwdmlyr3qkT3Htvs5ur+8AH2NylK0suHkntyJHUHnZ48iVhGXOQm1n78NWvJt/3EgH33EN0\n7Ihqa+HTn05Cvaxsx6Psb3+brYuX8MSYO1l5/peYMXkG61TGlop+cOedbXcejXCQm1n786lPFb7U\na9o06No1GblfdFHhS73++78L22/cSO3NtzDvf25jwzsOAmDD2wdQfc8feeaiy6m9cAR1739/Zte2\nO8jNrH0bMiSZHomAuXOhR49kjv366wtfFbB0Kc99+lxeHnLyW/ft0IHln/8KM+9/kJp1NWzZry9M\nnNjmp+AgNzPbZtAgWLcuCfVly2DffZOvCojg6au+2eRuNe88mJn3T+PZC75M7TnnsPXEE9v06wQc\n5GZmjTnkEOpOPZXNe+/DuvccmrwxugPRsSPLLrqc6runsmHVK2zp02f3r47ZRQ5yM7PGTJlCTLib\n2T/5FVs7dWz2buvecygPT57Bis+cT+1pw9k6dGirj84d5GZmDa1Zw5ZzP8vSr17Oq8cev8u7R+fO\nLLniWh4dfz9vLHkmmTt//PFWKDThIDcza2DLEUfwRt8KXh7yMTqt/TfJR0V33drDj2T6A4/y/Gmf\npHbIEOrOOKO4haaa/3rBzKydqOvZiy6rVjH408PosHkzs77/PU4483S29NiHTb16s2nfPry5X1/e\n3K8vm3v3YVOv3mzu1ZtNvfZNf/amrls3kNjapSuLrvsOq/7r4xx7+Vfo2Gc/Oj08HY46qmj1OsjN\nzBro+vSit66YOZNOU6fSadEiuj37LLz4Irz0HFufmsvWNzdRV1sLdXWothZt2ULZ5s0QW6kt78Hm\nvdPw77Mfq487gQOmPUDtB45HV15B2fe+V5R6HeRmZjvToUPyF4qqqt66Or01GqSrVtF53jw6L15M\n+fLl8NJL8MJytlb0pW7TZmKvbpQVqTwHuZlZazjggOR22mlvWb0t/IvJb3aameWcg9zMLOcc5GZm\nOecgNzPLOQe5mVnOOcjNzHLOQW5mlnMOcjOznHOQm5nlnIPczCznihLkkj4maZmk5ZKuLUabZmbW\nPC0OckllwM+Ak4FDgbMlHdrSds3MrHmKMSI/BlgeESsjYjMwARhehHbNzKwZivHthwcAf693/yXg\n2IYbSRoBjACoqKigurq6CIfefTU1NZnXUCrcFwXuiwL3RUGp90WbfY1tRIwBxgBUVlZGVYPv9W1r\n1dXVZF1DqXBfFLgvCtwXBaXeF8WYWlkFHFjvfv90nZmZtYFiBPkc4GBJB0nqDJwFTClCu2Zm1gwt\nnlqJiFpJFwMPAWXAryNicYsrMzOzZinKHHlEPAA8UIy2zMxs1/iTnWZmOecgNzPLOQe5mVnOOcjN\nzHLOQW5mlnMOcjOznHOQm5nlnIPczCznHORmZjnnIDczyzkHuZlZzjnIzcxyzkFuZpZzDnIzs5xz\nkJuZ5ZyD3Mws5xzkZmY55yA3M8s5B7mZWc45yM3Mcs5BbmaWcw5yM7Occ5CbmeWcg9zMLOcc5GZm\nOecgNzPLOQe5mVnOOcjNzHLOQW5mlnMOcjOznHOQm5nlnIPczCznHORmZjnnIDczy7kWBbmkmyUt\nlbRQ0kRJPYtVmJmZNU9LR+TTgcMi4gjgWeC6lpdkZma7okVBHhHTIqI2vTsb6N/ykszMbFcUc478\nC8CfitiemZk1gyJixxtIDwP9GnloVERMTrcZBVQCZ0QTDUoaAYwAqKioGDRhwoSW1N1iNTU1lJeX\nZ1pDqXBfFLgvCtwXBaXSF4MHD54XEZUN1+80yHdG0gXAl4GPRMTG5uxTWVkZc+fObdFxW6q6upqq\nqqpMaygV7osC90WB+6KgVPpCUqNB3rGFjX4MuBr4cHND3MzMiqulc+Q/BXoA0yU9KekXRajJzMx2\nQYtG5BHx7mIVYmZmu8ef7DQzyzkHuZlZzjnIzcxyzkFuZpZzDnIzs5xzkJuZ5ZyD3Mws5xzkZmY5\n5yA3M8s5B7mZWc45yM3Mcs5BbmaWcw5yM7Occ5CbmeWcg9zMLOcc5GZmOecgNzPLOQe5mVnOOcjN\nzHLOQW5mlnMOcjOznHOQm5nlnIPczCznHORmZjnnIDczyzkHuZlZzjnIzcxyzkFuZpZzDnIzs5xz\nkJuZ5ZyD3Mws5xzkZmY55yA3M8s5B7mZWc45yM3Mcq4oQS5ppKSQ1KcY7ZmZWfO1OMglHQgMBV5s\neTlmZrarijEi/yFwNRBFaMvMzHaRInY/fyUNB06KiMskPQ9URsSrTWw7AhgBUFFRMWjChAm7fdxi\nqKmpoby8PNMaSoX7osB9UeC+KCiVvhg8ePC8iKhsuH6nQS7pYaBfIw+NAr4BDI2I13cW5PVVVlbG\n3Llzm1V4a6murqaqqirTGkqF+6LAfVHgvigolb6Q1GiQd9zZjhHx0SYaPBw4CHhKEkB/YL6kYyLi\nHy2s18zMmmmnQd6UiFgE9N12f1dG5GZmVjy+jtzMLOd2e0TeUEQMKFZbZmbWfB6Rm5nlnIPczCzn\nHORmZjnnIDczyzkHuZlZzjnIzcxyzkFuZpZzDnIzs5xzkJuZ5ZyD3Mws5xzkZmY55yA3M8s5B7mZ\nWc45yM3Mcs5BbmaWcw5yM7Oc2+kfX26Vg0r/Al5o8wO/VR/Af5Yu4b4ocF8UuC8KSqUv3hER+zVc\nmUmQlwJJcxv7a9TtkfuiwH1R4L4oKPW+8NSKmVnOOcjNzHKuPQf5mKwLKCHuiwL3RYH7oqCk+6Ld\nzpGbme0p2vOI3Mxsj+AgNzPLOQc5IGmkpJDUJ+tasiLpZklLJS2UNFFSz6xramuSPiZpmaTlkq7N\nup6sSDpQ0iOSlkhaLOmyrGvKmqQySQskTc26lsa0+yCXdCAwFHgx61oyNh04LCKOAJ4Frsu4njYl\nqQz4GXAycChwtqRDs60qM7XAyIg4FDgO+Fo77ottLgOeybqIprT7IAd+CFwNtOt3fSNiWkTUpndn\nA/2zrCcDxwDLI2JlRGwGJgDDM64pExHxSkTMT5fXkwTYAdlWlR1J/YGPA2OzrqUp7TrIJQ0HVkXE\nU1nXUmK+APwp6yLa2AHA3+vdf4l2HF7bSBoAvA/4a7aVZOpHJIO9rVkX0pSOWRfQ2iQ9DPRr5KFR\nwDdIplXahR31RURMTrcZRfLS+q62rM1Kj6Ry4D7g8ohYl3U9WZA0DFgdEfMkVWVdT1P2+CCPiI82\ntl7S4cBBwFOSIJlKmC/pmIj4RxuW2Gaa6ottJF0ADAM+Eu3vAwargAPr3e+frmuXJHUiCfG7IuL+\nrOvJ0AnAaZJOAboCe0u6MyI+m3Fdb+EPBKUkPQ9URkQpfMNZm5P0MeBW4MMR8a+s62lrkjqSvMn7\nEZIAnwOcExGLMy0sA0pGNncAayLi8qzrKRXpiPyqiBiWdS0Ntes5cnuLnwI9gOmSnpT0i6wLakvp\nG70XAw+RvLl3T3sM8dQJwHnASem/hSfTEamVKI/IzcxyziNyM7Occ5CbmeWcg9zMLOcc5GZmOecg\nNzPLOQe5mVnOOcjNzHLu/wMQu7NjCmVRIAAAAABJRU5ErkJggg==\n",
            "text/plain": [
              "<Figure size 432x288 with 1 Axes>"
            ]
          },
          "metadata": {
            "tags": []
          }
        }
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "viVT6qimX6U0",
        "colab_type": "text"
      },
      "source": [
        "# Homework 2:  Construct the transformation that allows the car to turn right by phi degrees from location (a,b)\n",
        "\n",
        "Construct the single matrix  $R_{a,b,phi}$  as a product of the simple matrices  $T_{c,d}$  and  $R_\\phi$  constructed above.\n",
        "\n",
        "Verify the result by applying your new matrix to the car located at various locations.\n",
        "\n",
        "Produce a plot that shows the result of the transformation on a car located in at least two different quadrants.\n",
        "\n",
        "#Rotation at location..use the co-ordinate there"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "UI9gA8XgUgGl",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "def RabByPhi(a,b,phi, T, R):\n",
        "  \"\"\"Return a 3x3 matrix which rotates the car at location a,b by 90 degrees.\"\"\"\n",
        "  return T(a,b) * R(phi) * T(-a, -b)"
      ],
      "execution_count": 0,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "duhFQStqY2DR",
        "colab_type": "code",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 279
        },
        "outputId": "3b5dd6a1-5803-4599-e396-4b8459a27d13"
      },
      "source": [
        "#car 1 location at (3,2) heading toward (3,3)\n",
        "#car 2 location at (-2, -3) heading toward (-1, -3)\n",
        "#  Have one car rotate by 180? degrees, and the other rotate by some other amount of degrees.\n",
        "\n",
        "ax = plt.axes()\n",
        "\n",
        "#car 1 old position:\n",
        "ax.arrow(2, 2, 1.0, 1.0, head_width=0.5, head_length=0.7, fc='lightblue', ec='black')\n",
        "#car 1 new position:\n",
        "ax.arrow(2, 2, 2.0, -1, head_width=0.5, head_length=0.7, fc='lightblue', ec='red')\n",
        "\n",
        "#car 2 old position:\n",
        "ax.arrow(2, 2, 2.0, -1, head_width=0.5, head_length=0.7, fc='lightblue', ec='red')\n",
        "#car 2 new position:\n",
        "ax.arrow(2, 2, 2.0, -1, head_width=0.5, head_length=0.7, fc='lightblue', ec='red')\n",
        "\n",
        "\n",
        "plt.grid()\n",
        "\n",
        "plt.xlim(-5,5)\n",
        "plt.ylim(-5,5)\n",
        "\n",
        "plt.title('Showing Car rotates 90 degrees ?',fontsize=10)\n",
        "\n",
        "plt.savefig('Showing Car rotates 90 degrees.png', bbox_inches='tight')\n",
        "plt.show()\n",
        "plt.close()"
      ],
      "execution_count": 25,
      "outputs": [
        {
          "output_type": "display_data",
          "data": {
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXIAAAEGCAYAAAB4lx7eAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0\ndHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAbWElEQVR4nO3de5hVZd3/8feH4SQMCogMJhZWWnh5\nSiY1zRoseNJQ0qw8pNlBKvOMeYgy86nsQS8rO9iPqDD0CTXlEJmC4CRqFCcFQfABPCRapIQwoMAM\n398fa8EepxkYmD2z9mI+r+va16y99lr3+q574LPvfe+19ygiMDOz/OqQdQFmZtYyDnIzs5xzkJuZ\n5ZyD3Mws5xzkZmY55yA3M8s5B/keQNIoSYslLZT0pKRj0/XPS+pTxONUSrqtiO2dL+lpSYskLZB0\nVbHa3sExB0g6p1jbNfOYJ0man57rHZI6pusl6TZJy9Pf3dHNaOuGtuin3SXpK+m/xWcl3ZB1Pe2F\ngzznJH0AGAYcHRFHAB8F/t4ax4qIuRFxaTHaknQycDkwNCIOB44DXt+F/TvuzmPAAKA5Ad3c7XZI\nUgfgDuCsiDgMeAH4XPrwycDB6W0EcHtLj7cLde2oj1piOfA+4HDgc5L6t9JxrB4Hef7tD7waEZsA\nIuLViHi53uOXpKPBRZLeCyCpt6RJ6ShwtqQj0vWLJPVMR4qvSTo/Xf9bSUMkVUmamq67QdKvJVVL\nWilpe8BL+pakZZIek/S7JkaQ1wFXbas1IjZFxC/T/S+UNEfSU5Luk9QtXT9O0i8k/RUYXb8xSRdI\nmiJpJjAjPYeb6434P5Nu+gPgxPSVyxXpyHtW2kfzJR3fxHZlaXtz0n77cnrc/SU9mm73tKQTG5zn\nvsDmiHg2vT8d+GS6PBz4bSRmAz0l7d+wo9JXXM9Kegx4T73175L0oKR56Tm8t9762el5f1dSTbq+\nKt1uCrAkXfdZSX9L6/9/ksrS9UMl/SXtk3sllafrfyBpSdoHtzSsNSIejojNgICOwOZGfvdWbBHh\nW45vQDnwJPAs8HPgw/Ueex64JF2+CBibLv8E+Ha6fBLwZLr8C+DjwGHAHOCX6fr/A7oDVcDUdN0N\nwBNAF6AP8BrQCXh/Wk9XoEe671WN1L0G2KeJc9q33vJ3653DOGAqUNbIPhcALwG90/ufJAnNMqAC\neJHkSW/7OaTbdQO6pssHA3PT5YbbjQC+mS53AeYCBwEjgVHp+jKgR4O6RDIKr0zv/xhYlC5PBT5Y\nb9sZ27art24QsCitc2+SEe9V9bY/OF0+FphZr92z0+WvADX1zmkDcFB6fyDwB6BTev/nwPnp7/NR\noHu6/hrgepInpWWA0vU9d/Dv8rfAzVn//2gvt9Z6eWVtJCJqJA0CTgQGA3dLujYixqWb3J/+nAec\nkS5/kHRUGBEzJe0raW9gFvAhkuC5HRgh6QDg3xGxQVLDw/8xklcCmyStJgnME4DJEfEm8KakP+zG\naR0m6btAT5InqofqPXZvRNQ1sd/0iFhT7xx/l277T0l/JnmSWddgn07ATyUdBdQBhzTR9lDgCEln\npvf3IQn+OcCvJXUCJkXEk/V3ioiQdBbwQ0ldgGnpcZrrRGBiRGwESEfTpCPk44F76/1euqQ/PwB8\nIl3+X6D+yPlvEfFcuvwRkieKOWkbewGrSaa5DgUeT9d3Bv5CMvX1JvCr9JXZ1MYKlnQayZPmBbtw\nntYCDvI9QBpW1UC1pEUkc7Dj0oc3pT/r2Pnv+1Hga8DbgVHA6cCZJAHfmE31lpvTfn2LSUJkZiOP\njQM+ERFPSbqAZCS5zYYdtLmjx5pyBfBP4EiSqcY3m9hOJK8MHvqPB6QPkbySGSfp1oj4bf3HI+Iv\nJIGMpKEUnixWAQfW27R/uq45OgBrI+KoZm6/Tf0+EnBHRFxXfwNJp5I8KZ7dcGdJx5A8AZwJXEzy\niq6hI4BpEbF1F2uz3eQ58pyT9B5JB9dbdRTJiHpHZgHnpvtXkcyxr4uIv5O8rD44IlYCjwFXkQR8\ncz0OnCqpazpqHNbEdjcBN0vql9bRWdKX0sd6AK+ko9xzd+HY9c0CPpPObe9H8krjb8D6tP1t9gFe\nSUPnPJLpERrZ7iHgq2lNSDpEUndJ7wD+Gcn8/ljgP648kdQ3/dmFZJriF+lDU4Dz0/n844DXI+KV\nBrs/CnxC0l6SegCnAkTEOuA5SZ9K25akI9N9ZlOYhz9rB300AzizXn290/OZDZwg6d3p+u7p+ZaT\nTIc9QPIEeGQT7U5Kz83aiEfk+VcO/ERST6CWZA51xE72uYFkOmAhsJHCVRQAf6UQZrNIAvex5hYT\nEXPSl/8LSUa6i2jkapSIeEBSBfCwktfvAfw6ffhbaR3/Sn/2aLh/M0wkmWJ4Km376oj4h6TXgDpJ\nT5GM/H8O3Kfkjd0HKYxYFzbY7sckV7LMT+v9F8n0RRXwdUlbgBqSOeaGvi5pGMnA6faI2PYq5AHg\nFJLf2Ubg8w13jIj5ku5Oz2M1yVTONucCt0v6JskU0YR0u8uBOyWNSs+p0auBImJJuu80JVfXbAG+\nFhGz01dCv0uffAC+SfLkNllSV5LR/JWNtUsyrbWRZD7d2sC2Ny3MikZSeTp3341kRDkiIuZnXVd7\nkfb7G/Xm58+OiOFZ12WtxyNyaw1jJB1KcuXKHQ7xNjeI5A1cAWuBL2Rcj7Uyj8jNzHLOb3aameWc\ng9zMLOcymSPv06dPDBgwIItDb7dhwwa6d++eaQ2lwn1R4L4ocF8UlEpfzJs379WI2K/h+kyCfMCA\nAcydOzeLQ29XXV1NVVVVpjWUCvdFgfuiwH1RUCp9IanRz4h4asXMLOcc5GZmOecgNzPLOQe5mVnO\nOcjNzHLOQW5mlnMOcjOznHOQm5nlnIPczCznHORmZjnnIDczyzkHuZlZzjnIzcxyzkFuZpZzDnIz\ns5xzkJuZ5ZyD3Mws54oW5JLKJC2QNLVYbZqZ2c4Vc0R+GfBMEdszM7NmKEqQS+oPfBwYW4z2zMys\n+RQRLW9E+j1wE9ADuCoihjWyzQhgBEBFRcWgCRMmtPi4LVFTU0N5eXmmNZQK90WB+6LAfVFQKn0x\nePDgeRFR2XB9x5Y2LGkYsDoi5kmqamq7iBgDjAGorKyMrP8idan8VexS4L4ocF8UuC8KSr0vijG1\ncgJwmqTngQnASZLuLEK7ZmbWDC0O8oi4LiL6R8QA4CxgZkR8tsWVmZlZs/g6cjOznGvxHHl9EVEN\nVBezTTMz2zGPyM3Mcs5BbmaWcw5yM7Occ5CbmeWcg9zM9kgRwcKFC1m5cmXWpbQ6B7mZ7TG2bNnC\nI488wqWXXcY7BhzEkUceyS233pp1Wa3OQW5mubZ+/XruvfdezjrnXPpW9OOiy6/kX3ThxDPO5oAD\nD+S7N96YdYmtrqjXkZuZtYVVq1YxefJk7p88hdlPPMHAoyt53+ChjP7ylezb72289s9XuO7Mk5l0\n/3307t0763JbnYPczEpeRPD0008zadIk7ps0medWrmDQhz7C0aecwQXfv41u5T22b1tXV8fPrrmU\nSy7+Gh/84AczrLrtOMjNrCTV1tYya9Ys7p80icmTp1C7dSuVJw3l9EuvZeCgY+jYqVOj+0365U8p\n71TGN0eNauOKs+MgN7OSsX79eh588EHumziJhx58kH4Hvp2jqoZw+W1jecchA5G0w/2Xzp/DtP/9\nDQvmzaOsrKyNqs6eg9zMMrdkyRIuveLKRue7m6vm9bX85OqLGTtmDP3792/FakuPg9zMMrd+/Xrm\nzZnDh049g/OuuZ4uXffapf0jgjHXf53Th5/G8OHDW6nK0uXLD80sc8ceeyzLlj5Dty0bufr0oSye\nM3uX9p9+z12sffnv3HrLLa1UYWnziNzMSkLfvn35/b33MHHiRC66+GLeVzWEc678xluuSGnMi/+3\njHtu+x+eeOwxunbt2kbVlhaPyM2spJx++uk8s3gxnTa+znmV7+Gv0//U5Lab3nyD20Z+ldE/+AHv\nfe9727DK0uIgN7OSc9NNNzFl0iQA7hp9Az+95lLW/fu1/9hu/OgbGXTkEXzxi19s6xJLioPczErG\nc889hyRGjx7NJZdcQkSw7JlnOHxAf64a/lEe/9MUIgKA2dMeYPETf2bsL8fs9LLEPZ3nyM0scxHB\n5z73OcaPHw8kH8F/29uSSw+7d+/Oj3/0Q84+6zNc8Pkv8Jc/TuK0L32Nsd+5lgem/oF99tkny9JL\ngkfkZpaphQsX0qFDB8aPH8/3v/99ImJ7iNd33HHH8dSTC/jo8cdw/XlncNXIKznuuOMyqLj0eERu\nZpmICIYMGcKMGTMAWLNmDb169drhPl26dOHG73yHa66+mm7durVFmbngEbmZtblZs2bRoUMHZsyY\nwdixY4mInYZ4fd27d2/38+L1eURuZm2mtraWww8/nKVLl1JeXs7q1avZa69d+xSn/SePyM2sTUyZ\nMoVOnTqxdOlSJk6cyPr16x3iReIRuZm1qjfeeIN+/fqxbt06DjnkEBYvXkzHjo6eYvKI3MxazW9+\n8xu6devGunXrqK6uZtmyZQ7xVuAeNbOiW7t27fY3LwcPHsyMGTP85mQr8ojczIpq9OjR20N8wYIF\nzJw50yHeyjwiN7OieOWVV7Z/kOfcc89l/PjxDvA24hG5mbXYFVdcsT3Ely9fzp133ukQb0MOcjPb\nbcuXL0cSP/rRj7jyyiuJCN71rndlXVa746kVM9tlEcHZZ5/N3XffDSTTKv369cu4qvbLI3Iz2yUL\nFiygQ4cO3H333dx8881EhEM8Yx6Rm1mzbN26laqqKmbNmgUklxj6K2RLQ4tH5JIOlPSIpCWSFku6\nrBiFmVnpqK6upqysjFmzZjFu3DgiwiFeQooxIq8FRkbEfEk9gHmSpkfEkiK0bWYZ2rJlCwMHDmTF\nihX06tWLl19+ud3+geNS1uIReUS8EhHz0+X1wDPAAS1t18yytXbtWjp37syKFSuYMmUKa9ascYiX\nqKLOkUsaALwP+Gsx2zWzNvb66/RcsYIAtg4cSIfBg7OuyHZA2/6QaYsbksqBPwPfi4j7G3l8BDAC\noKKiYtCECROKctzdVVNTQ3l5eaY1lAr3RYH7oqDm1Vcpf+GFt67s3BkGDoR29sVXpfLvYvDgwfMi\novI/HoiIFt+ATsBDwJXN2X7QoEGRtUceeSTrEkqG+6LAfVGwvS9qayM++cnYChH1b337RqxcmWmN\nbaVU/l0Ac6ORTC3GVSsCfgU8ExG3trQ9MysxZWXw+9+jbRF+4YUEwOrV8M53ggQ9e8LChVlX2m4V\n4wNBJwDnASdJejK9nVKEds2sFI0ZUwj1a69NQv311+HII5NQ79YNHn006yrblWJctfJYRCgijoiI\no9LbA8UozsxK3E03FUL9llsICd54Az784STUu3SBKVOyrnKP54/om1lxjByJtqYz6XfcQXTsCJs3\nw/DhICX3x43Luso9koPczIrv/PPRli1JqP/hD9C5M6qrg89/Pgn1Dh1g9Oisq9xjOMjNrHUNGwab\nNiWh/thj0K1bMh1zzTVJqEvJsu02B7mZtZ0TToANG5JQX7QIevZEkIzOt4X6hRdmXWXuOMjNLBuH\nHQb//ncS6i+8ABUVSaiPHVsI9dNPh7q6rCsteQ5yM8ve298O//hHEuqvvQYHHZSE+qRJyadIJaiq\ncqg3wUFuZqWld29YuTIJ9Y0bk5E7wJ//XAj1o49OrogxwEFuZqVsr72SufQIqK2F449P1i9YkFyj\nLsHBBycfSGrHHORmlg9lZfD444VQP+WU5FOly5cnXxEgQf/+yRRNO+MgN7P8KSuDP/6x8KnSc85J\nQn3VKth//yTU+/SBFSuyrrRNOMjNLP/uuqsQ6pdemoT6a6/Bu9+dhPree8O8eVlX2Woc5Ga2Z/nx\njwuhfuONyWWM69dDZWUS6nvtBdOn717bdXVs2v9tbN6nJ9x+e3HrbgEHuZntub71rcL3v/z858lX\nA7z5Jgwdmlyr3qkT3Htvs5ur+8AH2NylK0suHkntyJHUHnZ48iVhGXOQm1n78NWvJt/3EgH33EN0\n7Ihqa+HTn05Cvaxsx6Psb3+brYuX8MSYO1l5/peYMXkG61TGlop+cOedbXcejXCQm1n786lPFb7U\na9o06No1GblfdFHhS73++78L22/cSO3NtzDvf25jwzsOAmDD2wdQfc8feeaiy6m9cAR1739/Zte2\nO8jNrH0bMiSZHomAuXOhR49kjv366wtfFbB0Kc99+lxeHnLyW/ft0IHln/8KM+9/kJp1NWzZry9M\nnNjmp+AgNzPbZtAgWLcuCfVly2DffZOvCojg6au+2eRuNe88mJn3T+PZC75M7TnnsPXEE9v06wQc\n5GZmjTnkEOpOPZXNe+/DuvccmrwxugPRsSPLLrqc6runsmHVK2zp02f3r47ZRQ5yM7PGTJlCTLib\n2T/5FVs7dWz2buvecygPT57Bis+cT+1pw9k6dGirj84d5GZmDa1Zw5ZzP8vSr17Oq8cev8u7R+fO\nLLniWh4dfz9vLHkmmTt//PFWKDThIDcza2DLEUfwRt8KXh7yMTqt/TfJR0V33drDj2T6A4/y/Gmf\npHbIEOrOOKO4haaa/3rBzKydqOvZiy6rVjH408PosHkzs77/PU4483S29NiHTb16s2nfPry5X1/e\n3K8vm3v3YVOv3mzu1ZtNvfZNf/amrls3kNjapSuLrvsOq/7r4xx7+Vfo2Gc/Oj08HY46qmj1OsjN\nzBro+vSit66YOZNOU6fSadEiuj37LLz4Irz0HFufmsvWNzdRV1sLdXWothZt2ULZ5s0QW6kt78Hm\nvdPw77Mfq487gQOmPUDtB45HV15B2fe+V5R6HeRmZjvToUPyF4qqqt66Or01GqSrVtF53jw6L15M\n+fLl8NJL8MJytlb0pW7TZmKvbpQVqTwHuZlZazjggOR22mlvWb0t/IvJb3aameWcg9zMLOcc5GZm\nOecgNzPLOQe5mVnOOcjNzHLOQW5mlnMOcjOznHOQm5nlnIPczCznihLkkj4maZmk5ZKuLUabZmbW\nPC0OckllwM+Ak4FDgbMlHdrSds3MrHmKMSI/BlgeESsjYjMwARhehHbNzKwZivHthwcAf693/yXg\n2IYbSRoBjACoqKigurq6CIfefTU1NZnXUCrcFwXuiwL3RUGp90WbfY1tRIwBxgBUVlZGVYPv9W1r\n1dXVZF1DqXBfFLgvCtwXBaXeF8WYWlkFHFjvfv90nZmZtYFiBPkc4GBJB0nqDJwFTClCu2Zm1gwt\nnlqJiFpJFwMPAWXAryNicYsrMzOzZinKHHlEPAA8UIy2zMxs1/iTnWZmOecgNzPLOQe5mVnOOcjN\nzHLOQW5mlnMOcjOznHOQm5nlnIPczCznHORmZjnnIDczyzkHuZlZzjnIzcxyzkFuZpZzDnIzs5xz\nkJuZ5ZyD3Mws5xzkZmY55yA3M8s5B7mZWc45yM3Mcs5BbmaWcw5yM7Occ5CbmeWcg9zMLOcc5GZm\nOecgNzPLOQe5mVnOOcjNzHLOQW5mlnMOcjOznHOQm5nlnIPczCznHORmZjnnIDczy7kWBbmkmyUt\nlbRQ0kRJPYtVmJmZNU9LR+TTgcMi4gjgWeC6lpdkZma7okVBHhHTIqI2vTsb6N/ykszMbFcUc478\nC8CfitiemZk1gyJixxtIDwP9GnloVERMTrcZBVQCZ0QTDUoaAYwAqKioGDRhwoSW1N1iNTU1lJeX\nZ1pDqXBfFLgvCtwXBaXSF4MHD54XEZUN1+80yHdG0gXAl4GPRMTG5uxTWVkZc+fObdFxW6q6upqq\nqqpMaygV7osC90WB+6KgVPpCUqNB3rGFjX4MuBr4cHND3MzMiqulc+Q/BXoA0yU9KekXRajJzMx2\nQYtG5BHx7mIVYmZmu8ef7DQzyzkHuZlZzjnIzcxyzkFuZpZzDnIzs5xzkJuZ5ZyD3Mws5xzkZmY5\n5yA3M8s5B7mZWc45yM3Mcs5BbmaWcw5yM7Occ5CbmeWcg9zMLOcc5GZmOecgNzPLOQe5mVnOOcjN\nzHLOQW5mlnMOcjOznHOQm5nlnIPczCznHORmZjnnIDczyzkHuZlZzjnIzcxyzkFuZpZzDnIzs5xz\nkJuZ5ZyD3Mws5xzkZmY55yA3M8s5B7mZWc45yM3Mcq4oQS5ppKSQ1KcY7ZmZWfO1OMglHQgMBV5s\neTlmZrarijEi/yFwNRBFaMvMzHaRInY/fyUNB06KiMskPQ9URsSrTWw7AhgBUFFRMWjChAm7fdxi\nqKmpoby8PNMaSoX7osB9UeC+KCiVvhg8ePC8iKhsuH6nQS7pYaBfIw+NAr4BDI2I13cW5PVVVlbG\n3Llzm1V4a6murqaqqirTGkqF+6LAfVHgvigolb6Q1GiQd9zZjhHx0SYaPBw4CHhKEkB/YL6kYyLi\nHy2s18zMmmmnQd6UiFgE9N12f1dG5GZmVjy+jtzMLOd2e0TeUEQMKFZbZmbWfB6Rm5nlnIPczCzn\nHORmZjnnIDczyzkHuZlZzjnIzcxyzkFuZpZzDnIzs5xzkJuZ5ZyD3Mws5xzkZmY55yA3M8s5B7mZ\nWc45yM3Mcs5BbmaWcw5yM7Oc2+kfX26Vg0r/Al5o8wO/VR/Af5Yu4b4ocF8UuC8KSqUv3hER+zVc\nmUmQlwJJcxv7a9TtkfuiwH1R4L4oKPW+8NSKmVnOOcjNzHKuPQf5mKwLKCHuiwL3RYH7oqCk+6Ld\nzpGbme0p2vOI3Mxsj+AgNzPLOQc5IGmkpJDUJ+tasiLpZklLJS2UNFFSz6xramuSPiZpmaTlkq7N\nup6sSDpQ0iOSlkhaLOmyrGvKmqQySQskTc26lsa0+yCXdCAwFHgx61oyNh04LCKOAJ4Frsu4njYl\nqQz4GXAycChwtqRDs60qM7XAyIg4FDgO+Fo77ottLgOeybqIprT7IAd+CFwNtOt3fSNiWkTUpndn\nA/2zrCcDxwDLI2JlRGwGJgDDM64pExHxSkTMT5fXkwTYAdlWlR1J/YGPA2OzrqUp7TrIJQ0HVkXE\nU1nXUmK+APwp6yLa2AHA3+vdf4l2HF7bSBoAvA/4a7aVZOpHJIO9rVkX0pSOWRfQ2iQ9DPRr5KFR\nwDdIplXahR31RURMTrcZRfLS+q62rM1Kj6Ry4D7g8ohYl3U9WZA0DFgdEfMkVWVdT1P2+CCPiI82\ntl7S4cBBwFOSIJlKmC/pmIj4RxuW2Gaa6ottJF0ADAM+Eu3vAwargAPr3e+frmuXJHUiCfG7IuL+\nrOvJ0AnAaZJOAboCe0u6MyI+m3Fdb+EPBKUkPQ9URkQpfMNZm5P0MeBW4MMR8a+s62lrkjqSvMn7\nEZIAnwOcExGLMy0sA0pGNncAayLi8qzrKRXpiPyqiBiWdS0Ntes5cnuLnwI9gOmSnpT0i6wLakvp\nG70XAw+RvLl3T3sM8dQJwHnASem/hSfTEamVKI/IzcxyziNyM7Occ5CbmeWcg9zMLOcc5GZmOecg\nNzPLOQe5mVnOOcjNzHLu/wMQu7NjCmVRIAAAAABJRU5ErkJggg==\n",
            "text/plain": [
              "<Figure size 432x288 with 1 Axes>"
            ]
          },
          "metadata": {
            "tags": []
          }
        }
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "EwXluSqZZIDX",
        "colab_type": "text"
      },
      "source": [
        ""
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "9uiX7BiQkc4b",
        "colab_type": "text"
      },
      "source": [
        "# Putting the car in the garage.\n",
        "\n",
        "With the car starting at $(3,0)$ and pointing toward $(3,1)$, we want to make the car go to $(-2,-2)$ pointing toward $(-2,-3)$.\n",
        "\n",
        "Draw the car's route out on the whiteboard, making sure to note it's location at each corner.\n",
        "\n",
        "Then, iteratively apply the matrices constructed above to verify your sequence of matrices puts the car into the garage appropriately.\n",
        "\n",
        "You may want to generate a plot of the car's position and path using the plotting code given above.\n",
        "\n"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "6rf8U2Bhk_wa",
        "colab_type": "code",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 68
        },
        "outputId": "3360735b-6603-4504-c11e-c616d4831ef0"
      },
      "source": [
        "#car 1 location at (3,0) heading toward (3,1)\n",
        "#car 1 final location at (-2, -2) heading toward (-2, -3)\n",
        " np.set_printoptions(precision=0)\n",
        "car = location(3,0)\n",
        "light = location(3,1)\n",
        "\n",
        "#apply matrices in order to get the car and the light to the correct locations\n",
        "#Car\n",
        "T(0,1)*RabByPhi(-2,-3,90,T,R)*T(-2,0)*RabByPhi(0,-3,-90,T,R)*T(0,-3)*RabByPhi(0,0,90,T,R)*T(-3,0)*RabByPhi(3,0,90,T,R) * car\n",
        "\n",
        "#light\n",
        "T(0,1)*RabByPhi(-2,-3,90,T,R)*T(-2,0)*RabByPhi(0,-3,-90,T,R)*T(0,-3)*RabByPhi(0,0,90,T,R)*T(-3,0)*RabByPhi(3,0,90,T,R) * light"
      ],
      "execution_count": 53,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "matrix([[-2.],\n",
              "        [-1.],\n",
              "        [ 1.]])"
            ]
          },
          "metadata": {
            "tags": []
          },
          "execution_count": 53
        }
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "sGRmNznBlAFL",
        "colab_type": "text"
      },
      "source": [
        "Once you have verified you've completed this task, create the single matrix from all of your transformations:\n"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "aJCtEa3elH_w",
        "colab_type": "code",
        "outputId": "43bd4628-022b-4fe2-8d79-4af60c38e52d",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 68
        }
      },
      "source": [
        "#create the full matrix as either a function or as a multiplication of matrices\n",
        "\n",
        "#fullmatrix\n",
        "\n",
        "def fullMatrix(a,b, T, R):\n",
        "  \"\"\"Describe this matrix\"\"\"\n",
        "  # We have to pass in the T and R functions if we are going to use them in this function.\n",
        "  return T(a,b) * R(-90) * T(-a, -b)\n",
        "\n",
        "fullMatrix(3,0, T, R)\n"
      ],
      "execution_count": 27,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "matrix([[ 6.123234e-17,  1.000000e+00,  3.000000e+00],\n",
              "        [-1.000000e+00,  6.123234e-17,  3.000000e+00],\n",
              "        [ 0.000000e+00,  0.000000e+00,  1.000000e+00]])"
            ]
          },
          "metadata": {
            "tags": []
          },
          "execution_count": 27
        }
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "oKnhjUqXCgDE",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        ""
      ],
      "execution_count": 0,
      "outputs": []
    }
  ]
}