04_Matrices.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Matrices"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "from sympy import init_printing, Matrix, eye, zeros, diag, Symbol\n",
    "init_printing() # Pretty printing"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Make a Vector"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "\n",
    "\n",
    "# Column vector\n",
    "\n",
    "\n",
    "v = Matrix([\n",
    "    [1],\n",
    "    [2],\n",
    "    [3]\n",
    "])\n",
    "\n",
    "# Row vector\n",
    "u = Matrix([    [4,5,6]  ])\n",
    "\n",
    "\n",
    "display(v, u)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Make a Matrix"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "A = Matrix([\n",
    "    [1,2,3],\n",
    "    [4,5,6],\n",
    "    [7,8,9]\n",
    "])\n",
    "\n",
    "A"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Identity Matrices"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "eye(2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "eye(5)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Zero Matrices"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "zeros(3)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "zeros(6)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Diagonal Matrices\n",
    "\n",
    "You can pass in a series of values that will be positioned along the diagonal of a square matrix."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "diag(1,2,3)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "diag(4,5,6,7)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "diag(3.14, 2.71, 6.626)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Access Rows/Columns at a Time\n",
    "\n",
    "Once you have created a `Matrix`, you can access it's rows and columns with the appropriate method for that object.\n",
    "\n",
    "There are however two caveats to indexing matrices in python. \n",
    "1. Python is 'column-major', meaning if you access a two-dimensional array you specify **column first and then row**. This is the opposite way around as the way we do it with matrices.\n",
    "2. Numbering starts from `0` in python. So what we would usually refer to as row 1, we will from here on in refer to as row 0."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "A = Matrix([\n",
    "    [1,2,3],\n",
    "    [4,5,6],\n",
    "    [7,8,9]\n",
    "])\n",
    "\n",
    "# Row access\n",
    "display(A.row(0))\n",
    "\n",
    "# Column access\n",
    "display(A.col(0))\n",
    "\n",
    "# Element-specific indexing\n",
    "display(A[0,2]) # See, row first, then column!\n",
    "display(A[2,2])\n",
    "\n",
    "# You can also access the n-th element in the matrix if you fancy.\n",
    "display(A[4])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Want reduced-row-echelon-form (rref)?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "A = Matrix([\n",
    "    [1,2,3],\n",
    "    [4,5,6],\n",
    "    [7,8,9]\n",
    "])\n",
    "\n",
    "A.rref() # As simple as that!"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can see that two elements were returned. The first is the row-reduced form of $\\mathbf{A}$, the other is a `tuple` of the columns of $\\mathbf{A}$ that contain pivots."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Want the inverse of a matrix?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Resetting A (in case of unordered use of cells)\n",
    "A = Matrix([\n",
    "    [1,2,-1],\n",
    "    [2,-1,-1],\n",
    "    [-1,-1,3]\n",
    "])\n",
    "\n",
    "A.inv() # One-liner!"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Transpose of a matrix?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Resetting A (in case of unordered use of cells)\n",
    "A = Matrix([\n",
    "    [1,2,3],\n",
    "    [4,5,6],\n",
    "    [7,8,9]\n",
    "])\n",
    "\n",
    "display(A)\n",
    "display(A.T) # Note that this is an attribute and not a method (no need for parantheses to call a function)."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Now with row-reductions:\n",
    "\n",
    "Because `sympy` has automated a lot of the matrix processes, it seems that row operations are not particularly intuitive operations for whatever reason that may be.\n",
    "\n",
    "The best way to learn this is likely to closely follow the example below and ask questions."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Resetting A (in case of unordered use of cells)\n",
    "A = Matrix([\n",
    "    [1,2,-1],\n",
    "    [2,-1,-1],\n",
    "    [-1,-1,3]\n",
    "])\n",
    "\n",
    "display(A)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# R1 = R1 - 2*R0                          .  < That dot is above the row number you're subtracting.\n",
    "A.row_op(1, lambda one, row_el: one - 2*A[0,row_el])\n",
    "#                                            ^ This is a varaible to allow sympy to iterate over all row elements\n",
    "display(A)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# R2 = R2 + R0\n",
    "A.row_op(2, lambda two, row_el: two + A[0,row_el])\n",
    "display(A)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# R1 <--> R2\n",
    "A.row_swap(1,2)\n",
    "display(A)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# R2 = R2 + 5*R1\n",
    "A.row_op(2, lambda two, row_el: two + 5*A[1, row_el])\n",
    "display(A)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# R2 = (1/11)*R2\n",
    "A.row_op(2, lambda two, row_el: int((1/11)*A[2,row_el]))\n",
    "display(A)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# R1 = R1 - 2*R2\n",
    "A.row_op(1, lambda one, row_el: one - 2*A[2,row_el])\n",
    "display(A)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# R0 = R0 + R2\n",
    "A.row_op(0, lambda zero, row_el: zero+ A[2, row_el])\n",
    "display(A)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# R0 = R0 -2*R1\n",
    "A.row_op(0, lambda zero, row_el: zero -2*A[1, row_el])\n",
    "display(A)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Invert a Matrix to Solve \n",
    "$$\n",
    "\\mathbf{Ax}=\\mathbf{b} \\quad \\Rightarrow \\quad \n",
    "\\mathbf{x}=\\mathbf{A^{-1}b}\n",
    "$$"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Changing A\n",
    "A = Matrix([\n",
    "    [1,2,-1],\n",
    "    [2,-1,-1],\n",
    "    [-1,-1,3]\n",
    "])\n",
    "\n",
    "# Right-hand side of the equation\n",
    "b = Matrix([\n",
    "    5,\n",
    "    0,\n",
    "    8\n",
    "])\n",
    "\n",
    "# Solving the equation.\n",
    "x = A.inv()*b\n",
    "\n",
    "display(x)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Symbolic Matrices\n",
    "\n",
    "The useful thing about some calculations is specifying a value in your matrix so that some property is true.\n",
    "\n",
    "E.g.\n",
    "\n",
    "Determine the value of $c$ so that $\\mathbf{A}$ has a determinant of 0.\n",
    "\n",
    "$$\n",
    "\\mathbf{A} = \n",
    "\\begin{bmatrix}\n",
    "    c & 2 & 1 \\\\\n",
    "    5      & -1 & 2 \\\\\n",
    "    -3     & 1 & -1 \\\\\n",
    "\\end{bmatrix}\n",
    "$$"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "c = Symbol('c')\n",
    "\n",
    "A = Matrix([\n",
    "    [c,2,1],\n",
    "    [5,-1,2],\n",
    "    [-3,1,-1]\n",
    "])\n",
    "\n",
    "display('A is: ', A)\n",
    "\n",
    "\n",
    "display('Det(A) = ', A.det())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "So if:\n",
    "$$\\text{det}(\\mathbf{A}) = 0$$\n",
    "then $$-c=0 \\quad \\text{ and so } \\quad c=0$$ \n",
    "\n",
    "for the matrix to have a determinant of zero.\n",
    "\n",
    "A matrix with a zero determinant is called a 'singular' matrix."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Eigenvalues & Eigenvectors\n",
    "\n",
    "An eigenvalue and eigenvector are solutions to the following equation:\n",
    "\n",
    "$$\n",
    "\\mathbf{Ax}\n",
    "=\n",
    "\\lambda\n",
    "\\mathbf{x}\n",
    "$$\n",
    "\n",
    "Where $\\lambda$ is an _eigen**value**_ and $\\mathbf{x}$ is an _eigen**vector**_.\n",
    "\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "* The direction of an _eigen**vector**_ is unchanged, but its magnitude may not be.\n",
    "\n",
    "* So we can say that $\\mathbf{x}$ is a vector who's direction is unchanged by the linear transformation, $\\mathbf{A}$.\n",
    "\n",
    "* The _eigen**value**_ is how much the magnitude of the _eigen**vector**_ is changed."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The approach is to:\n",
    "\n",
    "1. Find the _eigen**values**_.\n",
    "2. Use each _eigen**value** to find one _eigen**vector**_."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Eigenvalues"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "A = Matrix([\n",
    "    [1,2],\n",
    "    [2,4]\n",
    "])\n",
    "\n",
    "e_vals = A.eigenvals()\n",
    "\n",
    "print('Type of object: \\t' + str(type(e_vals)))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "So the `eigenval` method has returned us a `dictionary` object. So, what are the `keys` and `values`?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "e_vals.keys()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "e_vals.values()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "display(e_vals)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The eigenvalues: $\\lambda_1 = 0$ and $\\lambda_2=5$ are the `keys` of the `dictionary`. The `values` in that `dictionary` are the 'multiplicity' of that eigenvalue, meaning how many times that eigenvalue occurred. For both values, their multiplicity is one."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Eigenvectors"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Resetting A (in case of unordered cells).\n",
    "A = Matrix([\n",
    "    [1,2],\n",
    "    [2,4]\n",
    "])\n",
    "\n",
    "A.eigenvects()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "From the `eigenvects` method we are returned a `list` containing `tuples`.\n",
    "\n",
    "We also have the information we previously gathered (the eigenvalues) as well as the requested eigenvectors.\n",
    "\n",
    "Interpretting the first element of the `list`:\n",
    "\n",
    "The eigenvalue $\\lambda_1 = 0$ with multiplicity $1$ is used to obtain the vector $\\begin{bmatrix} -2 \\\\ 1 \\end{bmatrix}$.\n",
    "\n",
    "And interpretting the second value in the `list`:\n",
    "\n",
    "The eigenvalue $\\lambda_2 = 5$ with multiplicity $1$ is used to obtain the vector $\\begin{bmatrix} 1/2 \\\\ 1 \\end{bmatrix}$."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# -------- Questions --------"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Question 0:\n",
    "\n",
    "Display the transpose of the following matrix:\n",
    "\n",
    "(Run the cell to obtain the matrix)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAMsAAAB9CAMAAAA7tmdpAAAAP1BMVEX///8AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAADFBd4eAAAAFHRSTlMAMquZdlQQQO0wRIlmIs273e98bAefp9AAAAAJcEhZcwAADsQAAA7EAZUrDhsAAAqXSURBVHgB7V3blqM4DHQgYWYDue3y/9+6koxdkq9M7+acTh/80BF2RXJhcyvktDutUgb3ueXpKTh3WsczlcvnUnELExhW5nL6YBro+mK4TOfzPF659X6e3GmsDtXpOYwzFwcLTrU1EW64hZrTGKz8E8h2cESEJd4sF4k0MIMXT8A5D7jVLH5+roODVQRPr8W52yuQeT6KKK5UyHZwRIQlXg2XC8WlyHf6c5+Hc+iAAO2f0bfR6QKWRWxbw4uNYRuO86vORSHbwRERlkQzXGaJND05urRW//gRm+k4g1UEr7xj3Ewzmcptude5KGQ7OCLCEveGy3V9TBSYR6ftTr7qrmf/6ZQVarbPaZXOn/3ZZXZ1Lhq5I3gptuHihnVdhIobluU8dk5v2MewEi7utY0Ln1CWW4OLRu4IjoiwLBeaDOudhsY5HqDbKqe0tHth25PmLVihLX762XVfaawnGsb6uGzzUJD94IgIi05D+vpyGa7P9RUZPPjIqZXJHwPUDKuAvRODy7zSuZFnd4MLXQYCUty0giMiLBp3zWXhiXpeXzIyZI5r41R2ltnDYWHxVlroknya6Xi58IxtceGLtyC9h1ZwRISVcPGXgRPtxEFN87Rz2zYuFbAqUDoOp0lOOk0u/G1C7giOiLAslzBew7IdiCM5rpWJrpO+wKph3X1wl4HL+hriya+IJmQ/OCLCIm9mjm2X58fN+WvbM06jPC7dkG6VsHKUW9jnKc7VxrUSyG5wRISVcrnKOWyhPXdjf+K90D+pukQusArYM7mc7nI/wa3+clPA8VEXkN3giAiLfJpxcSe6Y/T3lrdxHsb6DCO2r3DjCavUx3kcvUdqnB/rShO4UoDsBUdEWOTUcqlE+ZDqg8v3HKhjXI5xefceOObYu/fw1/wf4/K1/fbub/3scRE5LlHRsj2q2iHRZSiuuI3j4yHCwY2kwUe4iStgITQWGnUVgHDO7fm4yNNNoqJpT2KjXUl0GYoqJn7KmVk4kLvfCz/4V4rc6ovQSM/K8dm2AI5AOBdUxsXLcYmKljlEu5LoMhRVDHKv/SLF78zP/G6tSggQGqdh5MfLWgEQzgWbctnkuERFy9yiXUl0GYoqXtL3gZ73r6yJTB5eQiqhkYnXuQAI5+Iw5aLlOOhppdhe4dMSXQl1l04xFymNOaaExjYXABPnCRcjx0FFK/WSNDSu1mJeGUa1z203T08/oEUghMY2F6VIip/gPDn2jRynVLRibN+uxLwiiiuv/o3BbX42H1Wj0NjhAkVSOSfTjgvvtSD7BFWG8aUS2hOJrgR9BsnGzRASM6AWGlvHC0k6WpGEc8PFyHFaRcviUkVstxJdATp6YUVaXtWTrREaW1wM0CnnzOXX+kviWDlOq2iFDjrT3jqHnj2Vm2igDzkzl/xBaKTWFhcLVPvpb6XBGjnOqGiF2LadJbpKuXA0IrLK8V/lEqas12kaXAwwOJfYZo5JzSbHGRWt0M3QDomuAKIq/4aS3rF5v88o+6VwCI3U0uDiFDA6F2c5l02OMypaGpe2QzskugKIBuRJd2HjQFfMhe9eLmtVgoXQSLgWFwDhXGKnXKIcZ1S0QjdjOyS6Asrd/UtavvpfSfV7xPchORhC4zi81ntdFIxA5ZzdpVzyEJ9Tc3D5nmN1jMsxLu/eAz9tjk2Uk/QTyo2uwz8zf+yzR+enHS96js3zmTOp3H5hroOMwpySB3vDbyW8Hjq223GZ7nTvN9Jb/f3CXA8pz0qs4EEejMFTY1P4EgkvRcXtVBC0XB4c+UH6yn5hroOEMAd5MHbGGFD4EgnPoOIG4LHKcLnGB6X9wlwHCWEO8mAMnhrbU0si4aWouJ0+5BguD0mOjNi9whw/YtVkYghz4rYtH26dSyS82J/UaHJZ76f5HLMkdgtzroXUCt4mD6Z9itumc1HCi82pYeDUqMeFxF7evT45ZL8w10EqBa+VAMgd1Z3b9EGurhUNZ4zlIlLJGCSsvcIcvZSoS3hamAsSyp7OQcKroQ31lMv2QmGJaZY7hTnyU0UaYS7Kg5XuqR2tJLwK2A5j4BK0PufzxRY6kPcLcx2kEeaMPFjoIrhs+mABo6oA95Va63P+vRCPy35hro0Ms0pEFSsPqk4FM3bOSHihNfuM8K1FHy90w8w3/yOdmfcLcx2kEubIvRJMs55RReiclfBKSKkL8AAwXNxAF7SJz2P7hbkOEsIc5MEQO/vcOpdIeBksVLS5OBLuBnmFtV+Y6yCjMEc3eTEVMPTGfEaFL5HwDAgbER6r7LjE6o80Di7fc9iOcTnG5d174KfNsUPre/eM+XP/P22OaX3sz/fG9/mGHRcsnoVV7mts70l4nbS/1DvBFxEb0wa1rV0GbZKbLRcsnoWlnCgztnckPJP2l4pzyl8w73STPj0by9T4Rp5UCb8sOGqT8nXLBYtnYYUo9jO2dyQ8pP0VxDnrkrcWUbXmxnooShUUjCwLjtqkeLJckAkCKw/INbG9I+GZtL/0gSN3/RC/ly3dLG/nGriENinI/8xFvFQlPJv21+fiF9RdqwlAFE25TLTJhEtcudtbRmvb6xmAJu2vy2XyT9HXeqYJkYFLq00mxz4Wz8KSHZ/9Me2NDECT9tflcvMJgKfGLwdwhjD3hhb7Gm2Squy4SJexeBaWNGR/tvYgtmTtXKHT/nZwEXWjzSW6nLwGFLXJAhcsnoVV7GZY49uW8FTaX5fLrjlGA7At9vXJzlGbZC6/fv/l+4qVu7DKLEx7T8LzK3PJUZeL88f+pXXsbx3iTEJok1L5z2/1HtkfVbxyF1aZi27vSng00fwpvM+F32PxxO+/pGeX0Call2aOeSmO2cIqc9HtTQnPpP31ufhrJb9lrBe4hDYpaMMFi2dhlX3q9pDhV0SatL8+F/fke5jXlmVe9KgW+0KbzLnQK9ewchdW2aFqb0t4SPvLxbmC62mm399pUqGTMpYFR22SXZlxKfj+pKqDy/ccrWNcjnF59x445ti79/DX/B/j8rX99u5v/fBx2SO3qaW5HTiEOVi1AQrCXSdV0Kb9+Rfo7LIwLjvkNrU017XhEOZgVahE4a6TKpik/eFJMOeyR25TS3M7cAhzsCpconDXSRW0aX/qJwBzLnvkNrU0twOHMAerzAXCXSdV0KT96Z8AzLnskNtY0eGnWF6a24ZDmINVpuIS4a6eKojY5EmvOc647JRCpEOUeteDQ5iDVSZjhbtWqmCITZ9mzXHGZZfc5nvDqXc9eBTmINGVqRjhrpMqSB582p9Zc5yfx25eBu3IbdIfTr3rwqH1wSqxSYW7Rqogf92n/ek1x4Vzcm/SoB+SeteHB2GOYpnf44MjbyXCXT1VkOE+7c+sOS5w2S23bal3u9Q5LPGFlXKBcNdJFaQv+th2zXGJy065LaTe7YJjiS+slAuEu3aqIH1vi23WHFN1duxvr6bacptamttR5yDMwUpZ+G0Id51UQSwLli96NJs5lz1ym1qa24FD64Pl+5793ZtUmKT94ScAC1z2yG0q9a4DhzAHK6PhK6Jw104VVLHpWql+ArDApRLp+1cfXL7nGB3jcozLu/eAn2Py4wf+Pdy7A77Hf/wfELRskEvj5xve04H/0SvdtnJx/wLI4ZNKMG86FgAAAABJRU5ErkJggg==\n",
      "text/latex": [
       "$$\\left[\\begin{matrix}85 & 77 & 94 & 95 & 77\\\\47 & 79 & 23 & 81 & 82\\\\63 & 83 & 87 & 12 & 1\\\\55 & 57 & 91 & 0 & 96\\\\0 & 20 & 99 & 63 & 24\\end{matrix}\\right]$$"
      ],
      "text/plain": [
       "⎡85  77  94  95  77⎤\n",
       "⎢                  ⎥\n",
       "⎢47  79  23  81  82⎥\n",
       "⎢                  ⎥\n",
       "⎢63  83  87  12  1 ⎥\n",
       "⎢                  ⎥\n",
       "⎢55  57  91  0   96⎥\n",
       "⎢                  ⎥\n",
       "⎣0   20  99  63  24⎦"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "from sympy import init_printing, Matrix, eye, zeros, diag, Symbol, randMatrix\n",
    "init_printing()\n",
    "\n",
    "A = randMatrix(5,5)\n",
    "display(A)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Answer 0:\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Question 1:\n",
    "\n",
    "If you have a matrix:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAHgAAABLCAMAAACMVLPjAAAAPFBMVEX///8AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAo1xBWAAAAE3RSTlMAMquZdlQQQOkwRN2JZiK7ze9sVohKbgAAAAlwSFlzAAAOxAAADsQBlSsOGwAABE5JREFUaAXtWtui7BAMVdTeZ6q34///9UTcQivmYXe/nPHQ0rAWEURmxOQwSfFLaQl8QkxOaUjzL/GK1bNJ54mn3+IsPCslnhaprE9CHNqISTFKsFZb399NqX1nOl7kRtkVmyB7RbwG5TuY7dNnoQOdZI6XEOoQwnjLsA5K94nIDxiGWbZYryJW4avHOqzUqc4N5K7g474LIY2XnudNHfxU5CvWsdBZTBVxGCGqgzfxl0udOhcPI7tWUuQ7Qs4OeypERYxdeWmEwnzvsecRHgjUJy5yP4FCvFy0mysx6A+SXFetujbjjslqVSZ2SePAttcHyI3zswPEOK6bEdsV2+2gka1nM8YdvtaZjP7FmGFgswCGEzmlqu2IDSyvnHacwFzMGeNwhOqM87XwFiG8fAsj7hLrZHWeRWUbypwh47BHa9SIQiU2VUgR5SNVL3GKkb+7PhcUrw6nRQ94ozwY19wxLhPE4kRk1bMZiVYdRjx73i0tLzLQmE1yv+q9Td0vJzgvsH54hoFFAPqa4hzDtwmr9jebLA8biN/tMDXGNUfizcOtZ3ccEkzUeKveFtjalexYYSVf/JZ5piXaEG9piWzKShW1EvtYvaxS0mMcYXvvEhO5sXAYJt6bnavCf67QjPg5ohb5Q9xq5LHyR9WPqbYF/qi61chjZU7V4WBhqMFFt7J7kJSGd0At8ZzcGWgWnYIC0OTMCY7A1j/CcvUERMErYiOVTAe1EPoM3khu32aCOyDD2d0KSTkA1eCX00ln4m09BsQu+keE4y5bgAr41b0tMisGxMZhx3T3FhF7UYAKOEO8biPi6Jh1PcLIS4DeITbg8A9GDJdEj30EXzPSXF4U6B1i7/aPiOESLcRsk8N64cQPFOgN4hn9qYFxwURpPVl2jiuglvjr+0/ubZAZvOcMR4ytyALMKDlTA1Hiv991DCTIZumTO2W82mWgm8zBXZxqIEpcbSAAS2WjDQT97ql3vyJ9TEAUnCMO65QANFl9GGHwwtoI2mIC6hMrebpDhhuy3Z2L2RYolcGrJ7fz9LV9JyAKft1A2laPlVtVP0bUAn+IW408Vv6o+jHVtsAfVbcaeazcqjoF3geEwwB+bl+cfiZCD2dNDLxDO+p9Z5iUoQF8tiJx+pkIvciB98b7Tnz5nQP4o4rF6eci9CXwDgz08MyEKUO9DrZicfq5CH0JvP8UMXH6uQh9FXhnB0ID+GzFEI0Fp58LG9eBdxaPBvDZitnp5yL0deCdxcO5jgF8vmJy+tkIfRV45/E8cwzgDypGp59TtagD7/nCisOrHjJeUPF3mAGxb+idfi5CTwPv7HKqAvhvEHunn4vQk8A7v46rAD5LXJx+NkKfA++gHg6vCuBzFQVx+tkIfQq8i9r7ribYF0oAf1BRFKf/E6Gvr6kXpT71oXUEnuK54P6vxPirFT3cL5r5yQ/5DyjG/yFE6/ID+E+y3GDhH1C0Fv8AXGI6C8H4XfgAAAAASUVORK5CYII=\n",
      "text/latex": [
       "$$\\left[\\begin{matrix}75 & 62 & 20\\\\14 & 9 & 41\\\\61 & 51 & 90\\end{matrix}\\right]$$"
      ],
      "text/plain": [
       "⎡75  62  20⎤\n",
       "⎢          ⎥\n",
       "⎢14  9   41⎥\n",
       "⎢          ⎥\n",
       "⎣61  51  90⎦"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "A = randMatrix(3,3)\n",
    "display(A)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Find the `rref` and determine if it has an inverse."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Answer 1:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Question 2:\n",
    "\n",
    "\n",
    "$$\n",
    "    \\begin{bmatrix} \n",
    "        4 & 3 \\\\\n",
    "\t\t-10 & -2 \\\\\t\t\n",
    "    \\end{bmatrix}\n",
    "    \\begin{bmatrix} \n",
    "        x \\\\\n",
    "\t\ty \\\\\t\t\n",
    "    \\end{bmatrix}\n",
    "    =\n",
    "    \\begin{bmatrix} \n",
    "        -13 \\\\\n",
    "        5 \\\\\n",
    "    \\end{bmatrix}\n",
    "$$\n",
    "\n",
    "Solve for $x$ and $y$."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Answer 2:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Question 3:\n",
    "\n",
    "Determine the eigenvalues and eigenvectors of:\n",
    "\n",
    "$$\n",
    "\\mathbf{Q}\n",
    "=\n",
    "\\begin{bmatrix}\n",
    "0 & 1 \\\\\n",
    "-1 & 0 \\\\\n",
    "\\end{bmatrix}\n",
    "$$"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Answer 3:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Question 4:\n",
    "\n",
    "Show, using `if` statements and comparisons, that the following matrix is equal itself after being multiplied by an appropriately sized identity matrix.\n",
    "\n",
    "# Answer 4:\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAPQAAACWCAMAAADezVWdAAAANlBMVEX///8AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAABHL6OuAAAAEXRSTlMAMquZdlQQQN0iRO/NZom7fEVHPZwAAAAJcEhZcwAADsQAAA7EAZUrDhsAAA5HSURBVHgB7V3booMqDrW29fRqx///2Um4ZCUQhDNv47YPu5EsAkvQ2kXKni5beM3TH3g9ItdpumzXG73uf4Dz9GWm88akL3+BLzh+m6TX9bY2T8ZyXdf5mcNcrtnaeR8Ccf0hYAR9bgvBB6bocrut11fuXUn6/l6Ca/kQ5PrJsOJ9eX+n6fnOrB+/wu8dDoG44hAwgt58ha5ea6ks0wmjMuezY0gv83XeIukfw34tMvObg85pgG/vFi41zW9DoGFgivZZ51s+86qtZILOncaIRimPoSFNjlsk/drasQi1heorXRr0en4/fdJDoH8drfuBk+isoYPLI/R3mhqkf2EoE6Z6W7YQ5BZvgOs0QHoIxA0NATNolPRr+9EEXsN4UxsN0tvnst5w5Ze032mk+dbwfQ6QHgJxK0NAAc3f7+3avNtSuDTS07xtX+HcIL1sHz4r73zlc3f0K07sz0ag5Tb1SQ+BuIEhIEA8fs9N7sq6i9HOpKd12z7xbkUOf6SXLVza13Qrr4N9iOp93eik8L2zO72HQNzKELAA/fKlWndTRvo+vx7bO58dn/S0hUDf9kmkJ5vLStf0nSdXj/QQiLs8BCxB152bbhrpL1/8ty2PYYP0I1y0X56/7Rd9vC181nukh0AcZwioQLO6tXD9+pVIxyeKC89MfjVIxw/inZHmup95us/82t4zTffmawjEtYeAChTvp9f0aOG1H0kv8bN1mtMYNkhf0jXtBaKyL5+6i0yrkYeTaQjE7Q0BAyg+HMVZ6ffUjPT0Sw8fDdLTTPN2ad69b3QnXMINPrQVP7b9ZqV0CMToIWAAPZl1GABppTAS6Ve4c3/zdLSkr/N7+8RJsF6vc/sTkLzyKb7+ti1PnKJNHA6BGD4EzKDndZ2v8kmE1pIFOpf5uvOFo6p4xAI70kdk6HA6STsn5ZBF50gfclgdUudIOyflkEXnSB9yWB1S50jnk7Kr82dQEuWf1+vv135In+B/ruv6awlQrBRlPd6uJEh7YiBkBynuy4OevfkVYzgjva/zS9NRlF9Yk1jbCgv84SvRvS1LhO+JrMcXKwloMFkI2UHC/Y3LdltSTx3S+zq/9CHq7XP4jvNuKsbw34JuEXUoCQIDerxdSQAiWwjZQcJ9jV+ks2Jck+7o/LnppN6/g5g2N5cA4X+xLkcyaw5QvEOPNysJBYoPEbKDhDtOa1maq0nv6/zSiaS3f4LE0iZd+NvTW/R4u5Ig7cGQkB1k6X5lDaHSyHjJZlfnT42L3h6OHzsyFQOyf3mkW0mKot9Ej9crCRpgbQ7ZQRZurD1VI93T+WPT0Nv5+LW7cij+5/rYkTlEj1crCZanPgpNdpDWjQWOeqS7On9omgcMavcj3yF0t5QN/yqCu3JHE3o8VhIqkBTEkB2kdmdFlCNUI93X+amW1duvaclWelQY2v/OgnuBmbQen1cSSgyOc8gOUrlv6g7KpP/Z/kG4aUDnV3o7Vbx1OCf/M3xs/LLgrloMZqHH05JbW/CzTe4iaRU9BtJL/f+p0i8GdH6lt9OoM+fIqCQSjrM/Lo+1SOfZJ7IqrSQ0XzlkBOwhCRHdS34w4Tr19O7o/LEh/htF+Zj80c4HEH+EP2SFAIGClXI5SI8vVhIKHB1KyA5SuymDCoFq0h2dH1Wj3v6gJ9rr3Fw4fIr/y2sq900+LBEoWNDji5WEAkezSkJ2kNp975Ce9nX+1Imkt3/iU22TtPK/aIHgl1dLKyoT9Hi9klDjJhWyg1Tup16tcUbaaedgRSfpgw1ok8450s1TczDHOdIHG9AmHR7phXJp/9TrSU9IZ773Xxjy80Y2rszzfIhfdkhQ/4rM6EwTJN8Pi/1OlEZRp3GsCsDiSGakIaP3lHmuGr+Wf1ief6QELS4uXpJ83wsZvvql5PuckF/EksPs328cdGCFEIY0ZPSOMs9Vo9j/DTL/qrQY6Vk0JPm+ExJiPxLyi1DpEP5O46ADK8QwpCGjd5R5qprE/l8QOO5tbUcEkE5IiP0UPOW8JZb1W/J3GgcdWCGWIS0yemymrcyTP4n9UYV5tZQvkqhMl9shRexn/CDpTuOgAyv0xpCO/RtR5nMG/hIFiVdTEKEkVJV8PyT2D5PuN86EMh1l1aSTcr+vzGex/xll/ktb7VfJ9/shRewfJt1vnEJhIQJWTXpImc9i/zOO9A5pPtmTJN8Pif2j03uocdCBVZHOMjr3tanMi9g/NsMmJN83Q2qxf/CaHmkcdGCFz2kt9g8p80rsj/eSe/NGhuT7fyP2/6sbWbtxmjIi/MKapkLszzL6vjKvxP74I71v8yMLyff7Ia3YP0i617haiMjEeALbJzLI6B1lPlQNmPh80PxR5hTPNC8VdUJC7Kfgg6R7jcuqAIiFnptrGjJ6R5kPVWMG/iNkiTQzbZB83wkJsX+c9LTfOOjAqkkrGb2jzCMDf1npR9hNzvTkJsn3nZAQ+5GQH7pY/YF/v3HQgVWTrsIftMBM74NyrGidpKtTctCCc6QPOrAVrT870qfYX82FAxb82en9x5d1oMx3JzWy/2F5lSS7XqXuezgqy4F6vVAhzW5CVVgArctOb1HmCZQFdYtPR8j+h+UDseGPVvMdLALpXnhACYnUfQdW/ERA07GkRZmHoO6Gm5D9D8tFIrsear4LVCGlFz4QIWG5SLgLOpa0Fqn3vsgj+x+W267a8Meo+TVYBdK9qIEqJFL3HVix2ZCm87+RRvY/LLddlV1v1PwarALtk0ZIWHU4KrHuHdJKmdeoMiiy/2GVmHissusldd9FqkBmfaAGIySsGkUlxq3p2JFWyvyeToXsf1hus/zrJXaEDX+0ml+hdSDdiwqoQ+rgu0BytkmHmkmZ16giIrL/YRUQOUR2PVL3xQmjCiTrA8BkCyFhZZ95125Nx450qJKUeY0yoegAu/zAKjH5OGfXGzU/O/FeBsL6ADDZyiFJzE27CWVP8a7cmg6TFrEfyjxV1qgilsr+H/gdAFXm7Poydb+IiUCmFwUKh8jnhwWvspJb0zFif7zy07Y4GqViBDN+AvIuP7BKjD6m7Hqr5mtntBHI9KIGphLk88Nywcmt6ZjpDWWeqmtUGQ3Z/7BKTDxGdr1R82swAple1ECV+Y/gDozmPjeZNxvSdAxpKPMURKOqmNjlB1YF4gJk1xs138FKINMLB4iQsByYbps7olaeDGmlzENQd+Op7P/O7wCQXQ81vxMS6wNNYNpDCMFdJNyWjiXtVj1e4Un6eGPqMzpH2j8vxys9R/p4Y+oz4pGWZ28fcrxS8+x9PHo+o/Oa9s/L8UrtSKutdGRTHZczgOyOqUsuMBSm3wDsh5RAxZ47Xti8LGAT9isk3LAYZElDj4dVheIC69Y/vffg0W/reLj4Y4Fiz50aKMsCRcJ+iYQbVsAY0tDjYZWRwrF197bjj35bxw0bgcWeOzVS1heKhP0SCTesgDGkocfDKiOFY+Pubcef/KaOGzUDYxPN3DQsCxQJ+2VQuGEFjCENPR5WGSkcG3feQN9FUmHymzouVgfCnjsVFMsCRcJ+iYQbVsAY0mpr+31lXgFzjn/ZohzLhj+dkDYQ9tyRQNlQywKhCKn7GWHe4RbLklZ6vNnk3kQJB+LOOf41JJYov9RxsQpIs6O9m7peFuBASNh3w8INy5KGHg/LDQV3zvF3YVQIP+q4WAB5GcqFhMJyWQAJ+24duGExaUm0gR4Py40Et+T4uzi94Q/quFATSO+5U6HtsoBK2K+QVAA3rCn8qjb/gA56PCwvEoR7lePvApV/P6QCUqDdD34sCxBSJ+w7PYAbFo0E7VCSc07ypJq/sJxAmH3zV+X4u0j4OyEBpDhmz50qLpYFqPv8xNPeQwhuWFTBXNNxMHhre1hVm1xQuGPSvosMhcFf1HHRKZDZc6cCYlmgSNivkeEhkLcYQo4/gwxp6PGwqkhcULh7u+wHf1HHjZsCmT13aiCWBWRDnxpEJcjnhxWAhrTaSqejzGt33kDfbZgKs1/XcbEZSP+Jq729KNfM6wtFwn4ZFG5YAWNJl9UOenySPujAVrTOka5OyUELzpE+6MBWtM6Rrk7JQQvOkT7owFa07EhrDT8L6lUVLkDyvZXRXfBoYW6xs9+PWgxope6nFlUvaXdT/L8ASxp6vAjqfo8l+b6Q0T00cuphOThpsbffDxYDOpn9k/SyCGlIKz1eBHWne1QkyfeFjF6hkVMPqwKFAmmxs9/PhMWAKCfk/xNch5VeFiENaejxENTrUFwieeiFjO6hkbYGq8ahxc5+P0FqpK+XtBjQyexHL4uQhjT0eAjqdfe4REgXMrqHBlVYNa5osb3fT6zKiwE2db8OiV4GH0Ia0tDwS0G9jFgk34uMXuL4GFRh1Tjb4t5+P6FuWAwwqft1yOYWQpZ0FvtLQb0KaJPvIaNXQCoAVVgVzrTY2e+HKsfFgE5m/6R6aUJa0lmPLwX1qo+hQJLvIaN7QFCFVeHKFnf2+6G6WVrVqftVyFQgvURIQxp6vBXUG/Fy8r2S0T0kqMKqcWWLzf1+uKosBqjU/TpkLMm9VLsSGdJRpeX/Y2sE9TqcSb7XMnoNHZzeqsXOfj/chFkMaGb2o5dFSE06TxraN18L6g4RnXxvZHQHO3ZNqxb39/vhFuxiQDOzH70sQmrSSsNXgrrHIzy5xcGxMroHxqSGVePQYme/H6qaFwNM6n4dEr0sQhrSSo8XQb2ORSVIvi9kdA8NqrAcnLTY2e+HqubFgE5mP3pZhDSkldgvgrrTPSqS5PtCRq/RyKmHVaO4JEv4U2e/H2o9LwYgdd8NKb0sQlrSbtXjFZ6kjzemPqNzpP3zcrzSPzzS4X/FyHfk4w0tGD3iv8Whp7pbeLX/7Q3q/N9b9A2FX9N/AfPau9ytVTKSAAAAAElFTkSuQmCC\n",
      "text/latex": [
       "$$\\left[\\begin{matrix}16 & 94 & 4 & 4 & 45 & 18\\\\6 & 42 & 23 & 89 & 29 & 97\\\\23 & 40 & 53 & 81 & 10 & 22\\\\51 & 66 & 98 & 65 & 99 & 91\\\\88 & 84 & 48 & 47 & 62 & 22\\\\86 & 52 & 11 & 63 & 79 & 53\\end{matrix}\\right]$$"
      ],
      "text/plain": [
       "⎡16  94  4   4   45  18⎤\n",
       "⎢                      ⎥\n",
       "⎢6   42  23  89  29  97⎥\n",
       "⎢                      ⎥\n",
       "⎢23  40  53  81  10  22⎥\n",
       "⎢                      ⎥\n",
       "⎢51  66  98  65  99  91⎥\n",
       "⎢                      ⎥\n",
       "⎢88  84  48  47  62  22⎥\n",
       "⎢                      ⎥\n",
       "⎣86  52  11  63  79  53⎦"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "from sympy import randMatrix\n",
    "\n",
    "A = randMatrix(6,6)\n",
    "display(A)\n",
    "\n",
    "# Your code here"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Question 5:\n",
    "\n",
    "(i) Write a loop to create five square, random matrices. Determine if each of them have an inverse.\n",
    "\n",
    "(ii) Write another loop to create five rectangular, random matrices. Should these have an inverse? What happens when you attempt to invert them?\n",
    "\n",
    "What can we say about the `randMatrix` function?"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Bonus Question :\n",
    "\n",
    "What is meant by a 'pivot'?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "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.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}