{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# First Step in SageMath\n", "\n", "The objective of this class is to present tools and classical data structures to interact with SageMath." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " ## Find help\n", " \n", " To find help, just type the name of the mathematical concept/object that you wish to use (with first letter in upper and lower case) followed by the question mark and press shift+Enter." ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": true }, "outputs": [], "source": [ "Permutations?" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "If you don't know the name of the mathematical function, you can guess it by writing the first letters and hit the TAB key." ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": false }, "outputs": [], "source": [ "Permu" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "You can also see the source code of a program by putting two question marks." ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [ "Permutations??" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " - You can find numerous tutorials in the page Help of the interface jupyter. There are and complete API of the available objects.\n", "\n", " - Try to find a programming tutorial of Python, a tutorial on linear algebra, a tutorial to use symmetric functions and the API on binary trees.\n", " \n", "### Reference\n", " \n", " http://doc.sagemath.org/html/en/tutorial/tour.html" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Polynomials" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " - Roots of a polynomial\n", " \n", " Look in the documentation how to write a polynomial and define the following one.\n", " $$\n", " P = x^3- 3\\cdot x^2 + 1.\n", " $$\n" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false }, "outputs": [], "source": [ "P = x^3 - 3*x^2 + 1" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Draw this polynomial on the interval $[1,3]$." ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAw8AAAJICAYAAAAn9ti2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzs3Xuc1nP+//HHVZMONEOpkGPOhy0Km5yPlTLlEDO7VLtO\nyzqzbMuy2LVfh12LRQ4bUSaRajogJVFYFDmsQznESgk1k06quX5/vH/ZtaKZ5rqu93V43G+3biTN\n56m9dmae1/v9er8TyWQyiSRJkiStQ4PYASRJkiTlBsuDJEmSpFqxPEiSJEmqFcuDJEmSpFqxPEiS\nJEmqFcuDJEmSpFqxPEiSJEmqFcuDJEmSpFrJeHlIJpNUV1fj3XSSJElSbsl4eVi8eDElJSUsXrw4\n04+WJEmSVA9uW5IkSZJUK5YHSZIkSbVieZAkSZJUK5YHSZIkSbVS7/Lw3HPPUVpaStu2bWnQoAGV\nlZWpyCVJkiQpy9S7PCxZsoQ999yT22+/nUQikYpMkiRJkrJQUX0/QLdu3ejWrRuAdzdIkiRJecyZ\nB0mSJEm1YnmQJEmS8kS6NwJZHiRJkqQ80aUL3H57+j5+vWce1teOO+5IIpGgbdu2tG3bFoDy8nLK\ny8tjRZIkSZJy1kcfwYsvwiWXpO8Z0crDrFmzKC4ujvV4SZIkKa88/jgUFcERR6TvGfUuD0uWLGH2\n7NnfnrT0wQcfMHPmTFq0aMFWW21V74CSJEmS1m38eDjgACgpSd8z6j3z8Morr7DXXnvRqVMnEokE\nF198MR07duSqq65KRT5JkiRJ67B8OUyaBEcfnd7n1Hvl4eCDD6ampiYVWSRJkiSthylTYNmy9JcH\nT1uSJEmSctz48bD11rDbbul9juVBkiRJynHjx0OPHpBIpPc5lgdJkiQph82aBbNnp3/LElgeJEmS\npJw2fjw0bgyHHpr+Z1keJEmSpBw2bhwcfDBsuGH6n2V5kCRJknLU4sXwzDPQs2dmnmd5kCRJknLU\nxImwcmUYls4Ey4MkSZKUo8aOhV13hXbtMvM8y4MkSZKUg2pq/nNEa6ZYHiRJkqQcNGMGzJuXuXkH\nsDxIkiRJOWncOCgpgS5dMvdMy4MkSZKUg8aOhW7doFGjzD3T8iBJkiTlmHnz4JVXMjvvAJYHSZIk\nKeeMHw+JBHTvntnnWh4kSZKkHDNuHHTuDJtumtnnWh4kSZKkHLJiBUyYkNlTltawPEiSJEk55Lnn\n4OuvMz/vAJYHSZIkKaeMHQtbbgnt22f+2ZYHSZIkKUckk6E89OgRBqYzzfIgSZIk5Yj33oP3348z\n7wCWB0mSJClnjBsHTZrAYYfFeb7lQZIkScoRY8fCoYdCs2Zxnh+tPJSVlVFaWkpFRUWsCJIkSVLO\nqKoKJy3F2rIEUBTrwcOGDaO4uDjW4yVJkqScMmECrFoV54jWNdy2JEmSJOWAceNgjz1gm23iZbA8\nSJIkSVmupgbGj4+76gCWB0mSJCnrvfwyLFgQd94BLA+SJElS1hs7FjbZBDp3jpvD8iBJkiRluTFj\noHt3KIp23FFgeZAkSZKy2Jw5MHMm9OoVO4nlQZIkScpqY8ZAo0bQtWvsJJYHSZIkKatVVsLBB0NJ\nSewklgdJkiQpa1VXwzPPQGlp7CSB5UGSJEnKUk8+CStXwjHHxE4SWB4kSZKkLFVZCe3bw7bbxk4S\nWB4kSZKkLLRqFYwblz1blsDyIEmSJGWladNg4ULLgyRJkqR1qKyEzTeHTp1iJ/kPy4MkSZKUZZLJ\nUB569oQGWfQdexZFkSRJkgTw7rswe3Z2bVkCy4MkSZKUdSoroWlTOPzw2Em+y/IgSZIkZZnKSjjq\nqFAgsonlQZIkScoiCxbA889n35YlsDxIkiRJWWXcuPDXHj3i5lgby4MkSZKURSoroXNnaNMmdpLv\nszxIkiRJWWL5cnjySTjmmNhJ1s7yIEmSJGWJyZNh6dLsnHcAy4MkSZKUNSoroV072G232EnWzvIg\nSZIkZYE1t0qXlkIiETvN2kUrD2VlZZSWllJRURErgiRJkpQ1ZsyAuXOzd8sSQFGsBw8bNozi4uJY\nj5ckSZKySmUlbLwxHHBA7CQ/zG1LkiRJUhaorISjj4ZGjWIn+WGWB0mSJCmyjz+G117L3iNa17A8\nSJIkSZGNHQtFRdCtW+wkP87yIEmSJEVWWQkHHxxmHrKZ5UGSJEmKqLoann46u09ZWsPyIEmSJEU0\nYQKsXJn98w5geZAkSZKiqqyEPfaA7baLnWTdLA+SJElSJCtXhmHp3r1jJ6kdy4MkSZIUybPPwsKF\nlgdJkiRJ6zBqFGy1FXTsGDtJ7VgeJEmSpAiSyVAeeveGRCJ2mtqxPEiSJEkRTJ8O//537mxZAsuD\nJEmSFMWoUbDJJnDQQbGT1J7lQZIkSYpg5Mhwt0NRUewktWd5kCRJkjLsvffgX/+CY4+NnaRuLA+S\nJElSho0eDU2bwlFHxU5SN5YHSZIkKcNGjgzFoVmz2EnqxvIgSZIkZdBnn8GLL+beliWwPEiSJEkZ\nVVkJDRpAz56xk9Sd5UGSJEnKoFGjwvGsLVvGTlJ3lgdJkiQpQ6qrYdKk3LoY7r9ZHiRJkqQMGT8e\nVq6EXr1iJ1k/lgdJkiQpQ0aNgo4dYZttYidZP5YHSZIkKQNWrAgrD7m6ZQksD5IkSVJGPP00LF6c\nm0e0rhGtPJSVlVFaWkpFRUWsCJIkSVLGjBoF228Pu+8eO8n6SySTyWQmH1hdXU1JSQlVVVUUFxdn\n8tGSJElSFKtXQ9u2cMopcOONsdOsP7ctSZIkSWn2z3/C/Pm5Pe8AlgdJkiQp7UaOhNatoXPn2Enq\nx/IgSZIkpVEyGcpDaSk0bBg7Tf1YHiRJkqQ0euMNeP99OP742Enqz/IgSZIkpdGIEVBSAocdFjtJ\n/VkeJEmSpDQaMQKOOQY22CB2kvqzPEiSJElp8u678NZb+bFlCSwPkiRJUtqMGAEbbghdu8ZOkhqW\nB0mSJClNRoyAo4+Gpk1jJ0kNy4MkSZKUBh9+CDNm5M+WJbA8SJIkSWnx2GPQuHFYecgXlgdJkiQp\nDUaMCLMOzZvHTpI6lgdJkiQpxT79FF54Ib+2LIHlQZIkSUq5kSOhqCjc75BPLA+SJElSio0YAYcf\nDptsEjtJalkeJEmSpBRasACefTb/tiyB5UGSJElKqVGjwl97946bIx0sD5IkSVIKjRgBBx0ErVrF\nTpJ6lgdJkiQpRRYuhEmT8nPLElgeJEmSpJQZMwZWrYJjj42dJD0sD5IkSVKKjBgB++0HbdvGTpIe\nlgdJkiQpBRYvhiefzN8tS2B5kCRJklJi/HhYsQKOOy52kvSxPEiSJEkpMGIEdOwI220XO0n6WB4k\nSZKkelq2LKw85POWJYhYHsrKyigtLaWioiJWBEmSJCklJkyAJUvyvzwkkslkMpMPrK6upqSkhKqq\nKoqLizP5aEmSJCkt+vaFGTPgzTdjJ0kvty1JkiRJ9fDNN1BZmf+rDmB5kCRJkupl4kSoqoITToid\nJP0sD5IkSVI9DB8Ou+wCe+wRO0n6WR4kSZKk9bRiBYwaBSeeCIlE7DTpZ3mQJEmS1tNTT4UtSyee\nGDtJZlgeJEmSpPU0fDjsthvsvnvsJJlheZAkSZLWw4oVMHp04aw6gOVBkiRJWi8TJkB1NfTpEztJ\n5lgeJEmSpPUwfHjYrrTbbrGTZI7lQZIkSaqj5csLb8sSWB4kSZKkOnvySVi8uLC2LIHlQZIkSaqz\n4cPhJz+BXXeNnSSzLA+SJElSHSxbBpWVhbfqAJYHSZIkqU6efBK+/tryIEmSJGkdhg+H9u1hl11i\nJ8k8y4MkSZJUS2u2LBXaKUtrWB4kSZKkWnr8cViypDC3LIHlQZIkSaq14cOhQwfYaafYSeKwPEiS\nJEm1sHQpjBlTuFuWIIXl4fbbb2e77bajadOmdO7cmZdffjlVH7reKioqYkeQ1puvX+UyX7/KVb52\ntTbjx4cCke1bltL5+k1JeXj44Ye5+OKLufrqq3n11Vfp0KEDXbt25YsvvkjFh683PwEol/n6VS7z\n9atc5WtXa/PII7DXXrDjjrGT/LisLw8333wzZ555Jn379mWXXXZh4MCBNGvWjEGDBqXiw0uSJElR\nLVkCY8cW9pYlSEF5WLlyJdOnT+fwww//9p8lEgmOOOIIXnjhhfp++LyTbe9kZFMes2S/bPpzMcva\nZVMWgE8//TR2hG9l05+NWdYum7L42l27bMoCmc2zri1L2fRnk87Xb73LwxdffMHq1atp06bNd/55\nmzZtmDdvXn0/fN7JphcWZFces2S/bPpzMcvaZVMW8BuwH2KWtcumLL521y6bskBm8wwfDh07wvbb\nx8+yLul8/Ral6wMnk0kSiURKP97ixYvX6/euWrWK6urqlGWpj2zKAtmVxyxrZ5a1M8vaZVMWCJ+7\nsyVPNv3ZmGXtsimLr921y6YskLk8X38dtiz99rfwQ4/Lpj+b+rx+mzdv/qPfwyeSyWRyfYNB2LbU\nrFkzRowYQWlp6bf/vH///lRVVTFy5Mjv/PvV1dWUlJTQunVrEokEbdu2pW3btgCUl5dTXl6+1ues\n+X2SJEmS0qOqqori4uIf/PV6rzw0atSITp06MWnSpG/LQzKZZNKkSZx33nk/+PtmzZr1o8H+V/Pm\nzamqqqpvXEmSJKlOysth3jyYPDl2kvRr3rz5j/56SrYtXXTRRfTr149OnTqx7777cvPNN7N06VL6\n9++fig8PhCHsupQNSZIkqb4WLYKJE+H//g/8VjRF5eHEE0/kiy++4Morr2T+/PnsueeePPnkk7Rq\n1SoVH16SJEmK4rHHYOVKOOmk2EmyQ71nHupqzezCuvZTSZIkSbEdeSSsWlUYW5ZqIyWXxEmSJEn5\nZt48ePrpMPOgwPIgSZIkrcUjj0CDBnD88bGTZI+8LQ/PPfccpaWltG3blgYNGlBZWRk7kvQddX2N\nTpkyhQYNGnznR8OGDfn8888zlFhatz//+c/su+++FBcX06ZNG4499ljee++92LGkb63Pa3Tw4MHf\nfs5d8/m3WbNmGUqsmCoqoGtXaNkydpIfN3DgQDp06EBJSQklJSV06dKFJ554Ii3PytvysGTJEvbc\nc09uv/32lF5WJ6XK+rxGE4kEs2bNYt68ecybN4/PPvuM1q1bpzmpVHvPPfcc5557Lv/85z+ZOHEi\nK1eu5KijjmLZsmWxo0nA+r9GS0pKvv3cO2/ePObMmZOhxIrlo4/ghRdyY8vSVlttxfXXX8/06dOZ\nPn06hx12GL169eLtt99O+bPSdsN0bN26daNbt25AuHdCyjbr+xpt1aqVhw0oa40fP/47P7///vtp\n3bo106dP54ADDoiUSvqP9X2NJhIJT5EsMA8/DE2bQq9esZOsW48ePb7z8z/+8Y/ceeedvPjii+y6\n664pfVberjxI+SiZTLLnnnuyxRZbcNRRR/H888/HjiT9qEWLFpFIJGjRokXsKNJa1fY1+vXXX7Pt\nttuy9dZb07t3b/71r39lKKFiqaiAY46BjTaKnaRuampqGDZsGEuXLmW//fZL+ce3PEg5YvPNN+eu\nu+5ixIgRPPbYY2y11VYccsghvPbaa7GjSWuVTCa54IILOOCAA9htt91ix5G+p7av0Z133plBgwZR\nWVnJ0KFDqampoUuXLnz66acZTKtMevttmDkzN7YsrfHmm2/SvHlzGjduzNlnn83IkSPZZZddUv6c\nvN22JOWbnXbaiZ122unbn3fu3Jn333+fm2++mcGDB0dMJq3d2Wefzb/+9S+mTZsWO4q0VrV9jXbu\n3JnOnTt/+/P99tuPXXfdlbvvvpurr7463TEVQUUFlJRA9+6xk9TeLrvswsyZM1m0aBEjRoygb9++\nPPvssykvEK48SDls3333Zfbs2bFjSN9zzjnnMH78eJ555hk233zz2HGk76nPa7SoqIi99trLz795\nKpkM5eG446Bx49hpaq+oqIh27drRsWNH/vSnP9GhQwduueWWlD/H8iDlsNdee81vzJR1zjnnHEaP\nHs3kyZPZeuutY8eRvqe+r9GamhrefPNNP//mqenTYfZsKCuLnaR+ampqWLFiRco/bt5uW1qyZAmz\nZ8/+9hSbDz74gJkzZ9KiRQu22mqryOmkdb9GBwwYwNy5c7/dknTLLbew3Xbbsfvuu7N8+XLuuece\nJk+ezFNPPRXzP0P6jrPPPpuKigoqKyvZcMMNmT9/PhCOuWzSpEnkdFLtXqP9+vWjbdu2XHfddQBc\ne+21dO7cmR122IFFixZxww03MGfOHE477bRo/x1Kn4oKaN0aDjssdpLau/zyy+nevTtbbbUVixcv\nZujQoUyZMoUJEyak/Fl5Wx5eeeUVDj30UBKJBIlEgosvvhgInxAGDRoUOZ207tfovHnz+OSTT779\n97/55hsuvvhi5s6dS7NmzWjfvj2TJk3ioIMOivWfIH3PwIEDSSQSHHLIId/55/fddx99+/aNE0r6\nL7V5jX7yySc0bNjw219buHAhZ5xxBvPmzWOTTTahU6dOvPDCC2kZRlVcNTXhiNY+faAoh75Lnj9/\nPn379uWzzz6jpKSE9u3bM2HCBA5LQwNKJDN8CUJ1dTUlJSVUVVV5Vr0kSZKyxpQpcMghMHUq7L9/\n7DTZyZkHSZIkibBlaeutIQ3XI+QNy4MkSZIK3sqV8OijYVC6gd8h/yD/aCRJklTwnnoKvvwyty6G\ni8HyIEmSpII3bBjssgt06BA7SXazPEiSJKmgLVsGI0eGVYdEInaa7GZ5kCRJUkEbNw6+/totS7Vh\neZAkSVJBq6iATp1gxx1jJ8l+lgdJkiQVrEWLYOxYVx1qy/IgSZKkgvXoo7BqleWhtiwPkiRJKlhD\nhsDhh8MWW8ROkhuilYeysjJKS0upqKiIFUGSJEkF7OOPYcoUOPnk2ElyR1GsBw8bNozi4uJYj5ck\nSVKBe+ghaNoUjj02dpLc4bYlSZIkFZxkMmxZ6t0bmjePnSZ3WB4kSZJUcF5/Hd56yy1LdWV5kCRJ\nUsEZMgQ23RSOPDJ2ktxieZAkSVJBWb06zDuUlUGjRrHT5BbLgyRJkgrKM8/A3LluWVoflgdJkiQV\nlCFDYIcdYN99YyfJPZYHSZIkFYylS2HEiLDqkEjETpN7LA+SJEkqGGPGwOLF8POfx06SmywPkiRJ\nKhhDhkDnzmHbkurO8iBJkqSCsGABPPGEg9L1YXmQJElSQRg+PPz1xBPj5shllgdJkiQVhCFDoFs3\naNUqdpLcVRQ7gCRJkpRus2fDiy/CsGGxk+Q2Vx4kSZKU94YOhebN4ZhjYifJbZYHSZIk5bVkEh58\nEI4/Hpo1i50mt1keJEmSlNemTYP334d+/WInyX2WB0mSJOW1wYNhm23goINiJ8l9lgdJkiTlrWXL\nwhGtp5wCDfzOt978I5QkSVLeGj0aqquhb9/YSfKD5UGSJEl5a/Bg6NIFdtwxdpL8YHmQJElSXpo7\nFyZMcFA6lSwPkiRJyktDh0KjRnDiibGT5A/LgyRJkvJOMhm2LPXuDRtvHDtN/rA8SJIkKe/MmAFv\nveWWpVSzPEiSJCnvDB4Mm20GRx4ZO0l+iVYeysrKKC0tpaKiIlYESZIk5aFvvoGHHoKTT4aiothp\n8ksimUwmM/nA6upqSkpKqKqqori4OJOPliRJUgEYNQqOPRbeeAP22CN2mvzitiVJkiTllcGDoWNH\ni0M6WB4kSZKUN774AsaNc1A6XSwPkiRJyhsVFeGY1vLy2Enyk+VBkiRJeWPwYOjRA1q1ip0kP1ke\nJEmSlBfeegumT3fLUjpZHiRJkpQXHngAWrYMKw9KD8uDJEmSct7q1TBkSJh12GCD2Gnyl+VBkiRJ\nOW/iRJg71y1L6WZ5kCRJUs677z7YfXfo1Cl2kvxmeZAkSVJO+/JLGDkSfvlLSCRip8lvlgdJkiTl\ntIcegpoaOPnk2Enyn+VBkiRJOW3QICgthdatYyfJf5YHSZIk5awZM+C118KWJaWf5UGSJEk5a9Ag\n2Hxz6No1dpLCYHmQJElSTlq+HIYOhf79oagodprCYHmQJElSTho5EhYtgl/8InaSwmF5kCRJUk4a\nNAgOPBB23DF2ksLhAo8kSZJyzkcfhVul778/dpLC4sqDJEmScs7998NGG8EJJ8ROUlgsD5IkScop\nNTVw331QVgYbbhg7TWGxPEiSJCmnTJoEH38Mp54aO0nhsTxIkiQppwwaBLvuCj/9aewkhcfyIEmS\npJzx1VfhiNZf/hISidhpCo/lQZIkSTnjoYdg9Wo45ZTYSQqT5UGSJEk5Y9Ag6NkT2rSJnaQwRSsP\nZWVllJaWUlFRESuCJEmScsirr4Yfv/xl7CSFK5FMJpOZfGB1dTUlJSVUVVVRXFycyUdLkiQph517\nLjz6KHzyCRR51XEUbluSJElS1lu+HIYOhX79LA4xWR4kSZKU9UaNgoUL4Re/iJ2ksFkeJEmSlPXu\nvhsOOgh23jl2ksLmoo8kSZKy2nvvweTJMGRI7CRy5UGSJElZ7d57YZNN4PjjYyeR5UGSJElZa8UK\nuO++MCjdpEnsNLI8SJIkKWuNGgVffAFnnBE7icDyIEmSpCx2991w4IGw666xkwgcmJYkSVKWmjUL\nnn4aHnwwdhKt4cqDJEmSspKD0tnH8iBJkqSs8803YVC6b19o2jR2Gq1heZAkSVLWGT0aFixwUDrb\nWB4kSZKUde6+Gw44AHbbLXYS/TcHpiVJkpRV3n8fJk6EBx6InUT/y5UHSZIkZZV77oGNN4YTToid\nRP/L8iBJkqSs4aB0drM8SJIkKWtUVsLnn8Ppp8dOorWxPEiSJClr3H03dOkCe+wRO4nWxoFpSZIk\nZYUPPoCnnoL774+dRD/ElQdJkiRlhXvvDYPSJ54YO4l+iOVBkiRJ0a1cCYMGwSmnOCidzSwPkiRJ\nim7MGJg/30HpbGd5kCRJUnR33QWdO8NPfhI7iX6MA9OSJEmKatYsmDABBg+OnUTr4sqDJEmSoho4\nEFq0cFA6F0QrD2VlZZSWllJRURErgiRJkiJbujTcKH3qqdCkSew0WpdEMplMZvKB1dXVlJSUUFVV\nRXFxcSYfLUmSpCyzpjjMmgXbbx87jdbFbUuSJEmK5s47oVs3i0OucGBakiRJUbz8cvgxZkzsJKot\nVx4kSZIUxZ13wjbbQPfusZOotiwPkiRJyrivvoKKCjjzTGjYMHYa1ZblQZIkSRl3//2wenUYllbu\nsDxIkiQpo2pqwpalPn2gdevYaVQXDkxLkiQpoyZOhNmzw+qDcosrD5IkScqoO+6A9u2hS5fYSVRX\nlgdJkiRlzMcfh6NZzz4bEonYaVRXlgdJkiRlzN13w4Ybws9/HjuJ1oflQZIkSRnxzTdw773Qty9s\ntFHsNFoflgdJkiRlxMiRMH8+nHVW7CRaX5YHSZIkZcTf/w4HHwy77x47idaXR7VKkiQp7V59FaZO\nhUcfjZ1E9eHKgyRJktLutttgq62gV6/YSVQflgdJkiSl1YIF8NBD8OtfQ5H7XnKa5UGSJElpdc89\n4U6H006LnUT1ZXmQJElS2qxcCXfeGe51aNkydhrVl+VBkiRJaTNqFPz733DuubGTKBXqXR5GjhxJ\nt27daNWqFQ0aNOD1119PRS5JkiTlgVtvhYMOgg4dYidRKtS7PCxZsoQDDjiA66+/nkQikYpMkiRJ\nygNrjmc977zYSZQq9Z53P/nkkwGYM2cOyWSy3oEkSZKUHzyeNf848yBJkqSU++KLcDzr2Wd7PGs+\nsTxIkiQp5TyeNT/VqTw89NBDNG/enObNm1NcXMy0adPSlStlnngCTjwRVq2KnUSSJKkwrFoFd9wB\nP/sZbLpp7DRKpTotIvXq1YvOnTt/+/O2bduu94N33HFHEokEbdu2/fbjlJeXU15evt4fc22KimDk\nyHCj4cCBoQFLkiQpfTyeNX/VqTxsuOGGtGvX7gd/vS6nLc2aNYvi4uK6PH69HHEE3HUXnHoqbL89\nXHpp2h8pSZJU0G69FQ48EPbcM3YSpVq9x1cWLlzIxx9/zKeffkoymeSdd94hmUyy2Wab0aZNm1Rk\nrLdf/hI+/BAuuwy23TZsY5IkSVLqvfYaPPccPPJI7CRKh3oPTFdWVrLXXntxzDHHkEgkKC8vp2PH\njtx1112pyJcy11wT9t317QvPPx87jSRJUn667TbYckvo3Tt2EqVDIpnhyxmqq6spKSmhqqoqI9uW\n/tuKFXDkkfCvf8GLL8IOO2T08ZIkSXltwYJwr8NVV8GAAbHTKB0K6qjWxo3DAM+mm0L37uH8YUmS\nJKXGwIHQoAGcfnrsJEqXgioPAC1awPjxUFUVltOWL4+dSJIkKfetWAG33w79+nk8az4ruPIA0K4d\nVFbC9OnQvz/U1MROJEmSlNsqKmD+fLjggthJlE4FWR4AOneGIUNg+HC44orYaSRJknJXMgl//Sv0\n6AE77xw7jdKpYMsDwPHHw403wp//HK5QlyRJUt1NmgRvvAEXXRQ7idKt3vc85LqLLoIPPoCzzoKt\nt4auXWMnkiRJyi033wwdOsChh8ZOonQr6JUHgEQCbrkFunWDPn1g5szYiSRJknLH22+Hw2guuih8\nX6X8VvDlAaCoCIYNC/c+9OgBn34aO5EkSVJu+NvfYPPNoawsdhJlguXh/9toIxg7NjTmnj1h8eLY\niSRJkrLbggXwwANwzjmwwQax0ygTLA//ZYstwrLb++/DSSfBqlWxE0mSJGWvgQPDG69nnhk7iTLF\n8vA/fvITGDECnnoKzj03HD0mSZKk71pzKVz//tCyZew0yhTLw1oceWRo0gMHwk03xU4jSZKUfdZc\nCnf++bGTKJMK/qjWH3LqqeEI10svhW23DScxSZIk6T+XwvXs6aVwhcby8COuvRY+/BBOOQXatoUu\nXWInkiRJim/NpXC33BI7iTItkUxmdld/dXU1JSUlVFVVUVxcnMlHr5cVK8I2prffhhdeCMe5SpIk\nFbKjj4bPPoMZM7zbodA487AOjRvDyJHQokX4P8qXX8ZOJEmSFM/bb8Pjj3spXKGyPNRCy5bhCNeF\nC6F3b1gwK6mXAAAgAElEQVS+PHYiSZKkOG6+OVwKd9JJsZMoBstDLW2/PVRWwssvwy9+ATU1sRNJ\nkiRl1rx54VK4887zUrhCZXmog/32gyFDYNgw+P3vY6eRJEnKrFtvDaXhV7+KnUSxWB7q6IQT4MYb\n4brr4N57Y6eRJEnKjMWL4Y47wm3SG28cO41i8ajW9XDxxeEOiF/9CrbeGo46KnYiSZKk9LrnHli6\n1EvhCp1Hta6nVaugVy947jmYOhXat4+dSJIkKT2++QbatQvH1993X+w0iinatqWysjJKS0upqKiI\nFaFeiorg4YfDIHWPHjB3buxEkiRJ6VFRAZ9+CpdcEjuJYnPloZ4+/RQ6d4ZNN4Vnn4XmzWMnkiRJ\nSp2amrDDYrvtYMyY2GkUmwPT9dS2LYwbB++/D2VlYTuTJElSvnj8cXjrLbjssthJlA0sDynQvj08\n+ig8+WQ49zizazmSJEnpc/314bj6/fePnUTZwNOWUuSoo2DgQDj99DAHcfHFsRNJkiTVzwsvhMNh\nRo6ERCJ2GmUDy0MKnXZaOML1kktgm23CnRCSJEm56sYbYaedoLQ0dhJlC8tDiv3xj/Dhh3DKKWEe\nYr/9YieSJEmqu3ffhVGj4O67oYEb3fX/+VJIsQYNwvnHe+8dWvr778dOJEmSVHc33QRt2sDJJ8dO\nomxieUiDJk1CU99kEzj6aPjyy9iJJEmSau+zz+CBB+CCC8L3NdIaloc0adkSxo+Hr76CY4+FFSti\nJ5IkSaqdW2+Fxo3hzDNjJ1G2sTyk0Q47wOjR8NJL8ItfhEtWJEmSsllVFdx5J5xxBmy8cew0yjaW\nhzTr0gUefDBc637llbHTSJIk/bg77oBly+Cii2InUTayPGRAnz5www3wpz/BvffGTiNJkrR2S5fC\nzTeHHRNbbBE7jbKRR7VmyCWXhCNcf/Ur2HJL6NYtdiJJkqTvuueeMK952WWxkyhbJZLJZDKTD6yu\nrqakpISqqiqKi4sz+ejoVq0Kw9PPPAPPPgt77RU7kSRJUrBiBWy/PRx2WDhpSVobty1lUFERDBsG\nu+wSjnCdMyd2IkmSpODBB2HuXBgwIHYSZTNXHiKYPz/cPN20KUydGu6DkCRJimXVqvDmZocOMGJE\n7DTKZq48RNCmDTz+OMyb5x0QkiQpvuHD4f334Xe/i51E2c7yEMnOO4c7IF580TsgJElSPDU1cN11\n4TCXTp1ip1G2szxEdMABYX/hsGFw+eWx00iSpEJUWQlvveX3Iqodj2qNrE8f+OQTuPhi2GabcJSr\nJElSJiSTYdXhoIPCm5rSulgessCFF8JHH8Gvfx3ugOjZM3YiSZJUCCZOhJdfhieeiJ1EucLTlrLE\n6tVwwgkwYUK4B2KffWInkiRJ+e6QQ+Drr0OBSCRip1EucOYhSzRsCEOHQvv2YeXhww9jJ5IkSfls\n2jSYMiXMOlgcVFuuPGSZBQugS5dQJqZNg5YtYyeSJEn56Oij4eOP4fXXoYFvJ6uWfKlkmVatwh0Q\nX34JvXvD8uWxE0mSpHwzY0b4fmPAAIuD6saXSxbaYQcYMwZeeQX69fMOCEmSlFrXXBO+3zjppNhJ\nlGuilYeysjJKS0upqKiIFSGrde4MDz0EjzwCl14aO40kScoXr74aLqq94goo8txN1ZEzD1nu1lvh\n/PPDX889N3YaSZKU6449Ft54A955x/KguvMlk+XOOw/mzAkFYqutwhyEJEnS+njtNRg1Cu6/3+Kg\n9ePKQw6oqQl7EseOhcmTw5YmSZKkujruuHC6kqsOWl8OTOeABg3gwQehUyc45hiYPTt2IkmSlGtm\nzoSRI511UP248pBDvvwy3AFRUwPPPx+OdZUkSaqN448PBcJVB9WHKw85pGXLcCZzdXVYgVi6NHYi\nSZKUC2bOhMceC7dJWxxUH5aHHNOuHYwbB2++GeYgVq2KnUiSJGW7a64J30OcfHLsJMp1locctPfe\n8OijYRXi7LMhsxvPJElSLnn99f+sOjRqFDuNcp3lIUd16wb33gv33APXXhs7jSRJylbXXAPbbQen\nnBI7ifKBu95yWP/+MHdueCdhiy3gtNNiJ5IkSdnk9ddhxAj4xz9cdVBqWB5y3IAB8O9/w69+BZtt\nBj17xk4kSZKyxbXXuuqg1LI85LhEAm67DT77DE48MVwi99Ofxk4lSZJie+ONMCN5772uOih1vOch\nTyxbBkccAe+9F+6A2HHH2IkkSVJMxx4bjmh9913Lg1LHgek80bQpVFbCpptC164wf37sRJIkKZZX\nXoFRo+APf7A4KLVcecgzc+bAfvvB5pvDM89A8+axE0mSpEzr1i18T/Dmm9CwYew0yieuPOSZbbYJ\n9z/MmgV9+sDKlbETSZKkTHruOXjyyXBEq8VBqebKQ56aNAm6d4fycrj//jBYLUmS8lsyCYccAtXV\nMH06NPBtYqWYpy3lqcMPh8GD4Wc/gy23hD/9KXYiSZKUbhMnwrPPwpgxFgelh+Uhj5WXh0vkLrkE\n2raFs8+OnUiSJKVLMglXXAGdO0OPHrHTKF9ZHvLcRReFS+TOOScMUR97bOxEkiQpHcaMgZdeCqsP\nbldWujjzUABqasIqxOjR4RPKAQfETiRJklKppgb22gtatoSnn46dRvnM3XAFoEEDeOCBcIRraSm8\n/XbsRJIkKZUeeQRefx3++MfYSZTvXHkoIIsWwYEHhhMYXngBttgidiJJklRfq1bBHnvA9tvDuHGx\n0yjfufJQQDbeONwBUVMTjnGtqoqdSJIk1deQIfDuu3DttbGTqBBYHgrMllvCE0/Axx+H4ekVK2In\nkiRJ62vFCrj6ajj+eOjYMXYaFQLLQwHafXeorITnn4f+/cNKhCRJyj133RXeELzmmthJVCiilYey\nsjJKS0upqKiIFaGgHXggDB0KDz8Mv/lN7DSSJKmuqqvDVqVf/AJ22y12GhUKB6YL3N//DueeCzfc\nYImQJCmXXHVV+Po9a1bYlixlgpfEFbhzzoH58+HSS6F1a+jXL3YiSZK0LvPnw1/+AuedZ3FQZlke\nxDXXwOefw6mnhstlevaMnUiSJP2Ya6+FRo3gt7+NnUSFxoFpkUjAHXeEC+T69IFp02InkiRJP2T2\n7DAoPWAAbLJJ7DQqNM486FvLl0O3bjBzJjz3XLhwRpIkZZfy8vB1etYsaNo0dhoVGlce9K0mTWD0\naNhmG+jaFebMiZ1IkiT9t+nTYdiwcLeDxUExuPKg7/nsM9h/f9hgA5g6FTbdNHYiSZIEcNRR8Mkn\n8MYbUOTkqiJw5UHfs/nmMGECLFwIRx8NX38dO5EkSZo4EZ56Cv78Z4uD4nHlQT9oxgw45BDo3BnG\njg0rEZIkKfNqamDffcPX4mnTwmEnUgyuPOgHdewYZiCmTIH+/cMnLkmSlHnDh4d5h//7P4uD4rI8\n6EcdeigMHRqGsy68EDK7TiVJkpYvD8eyHnMMHHRQ7DQqdO6Y0zqdcEK4B+Kss6BNG/jd72InkiSp\ncNx2WxiSfvzx2Ekky4Nq6Ve/gvnz4fLLoXVrOO202IkkScp/X3wBf/pT+Dq8yy6x00iWB9XBlVfC\n55/DmWeG41t7946dSJKk/Hb11WHL8FVXxU4iBc48qNYSCbj1VjjuOCgrC4PUkiQpPd59FwYODKv+\nrVrFTiMFHtWqOluxItz/8Mor8Oyz0KFD7ESSJOWfXr1g5kx45x1o0iR2Gilw5UF11rgxjBoFO+wA\n3brBBx/ETiRJUn555hmorAxHs1oclE1cedB6+/xz2H//sBdz6lTYbLPYiSRJyn01NbDPPuEW6Rdf\n9F4HZRdXHrTeWreGCRNg2TI46ihYuDB2IkmSct+QITBjBvz1rxYHZR/Lg+plu+1Cgfj0U+jRA5Ys\niZ1IkqTctXRpuE/p+OPD6r6UbSwPqrfddw8X17zxRjiJacWK2IkkScpNN98ctgVff33sJNLaWR6U\nEvvuC6NHh+NbTz4ZVq+OnUiSpNwyb14YkD7nHNh++9hppLWzPChlDjsMHn4YRo4MF8lldhRfkqTc\n9rvfwQYbwBVXxE4i/TBvmFZK9eoFgwZBv36wySZwww0Oe0mStC4vvwz33Qe33w4tWsROI/0wy4NS\nrm9fWLQIzj8/fAIcMCB2IkmSsldNDZx3HvzkJ3DGGbHTSD/O8qC0OO+8cHTr734HG28MZ50VO5Ek\nSdlp6NBwn8PkyeFuBymbRbskrnv37hQVFVFeXk55eXkmIyhDkkm44AK47bbwidH/mSVJ+q7Fi2Hn\nncOxrI88EjuNtG7R+u2wYcO8YTrPJRLhyLlFi8JWpuLicBeEJEkKrrsurNTfdFPsJFLteNqS0qpB\nA/jHP0JpOOEEeO652IkkScoOs2eHW6QvvRS22SZ2Gql2om1bqqqqcuWhgCxfDkcfDdOnhz2dHTvG\nTiRJUly9esGrr8I770CzZrHTSLXjyoMyokmTcInczjtDt27w7ruxE0mSFM+ECVBZCTfeaHFQbnHl\nQRn15Zdw0EFhQGzqVNh669iJJEnKrJUroUMHaNUKnnnG+5CUW1x5UEa1bBnebWnYEI48EubPj51I\nkqTMuv32sAJ/yy0WB+Uey4Myrm1beOqpsPpw5JHw1VexE0mSlBkLFsAf/gCnnw577hk7jVR3lgdF\nscMOMHEizJ0bZiCqq2MnkiQp/QYMCKsN114bO4m0fiwPima33cIKxHvvQc+esGRJ7ESSJKXP88+H\n48uvuy7MO0i5yIFpRffCC2H7Upcu4eSJJk1iJ5IkKbVWrYJOnaBx4/B1r2HD2Imk9ePKg6Lbbz8Y\nOzZcIHfSSeEUCkmS8sltt8Ebb8Cdd1oclNssD8oKhxwCjz0Gjz8Op5wCq1fHTiRJUmp8+ilceSWc\nfXZYfZBymeVBWaN7dxg2DB59NJxCUVMTO5EkSfV30UWw4Ybwxz/GTiLVX1HsANJ/O+44uP9+6Ns3\nfKK99VbPwJYk5a4JE2D4cBgyBDbeOHYaqf4sD8o6J58MS5fCmWeGAvHnP1sgJEm5Z/ly+PWvw9bc\nn/0sdhopNSwPykpnnBEKxIUXhgLx+9/HTiRJUt3ccAN89FE4SdA3wZQvLA/KWhdcAF9/HYrDhhuG\nPaOSJOWC998P9zlccgnsumvsNFLqWB6U1S6/PFwed/HFoUCceWbsRJIk/bhkEs45BzbbzJVz5R/L\ng7JaIhHeuVmyBM46C5o1C0e5SpKUrUaOhCeegNGjw9ctKZ9YHpT1Egn4299CgejfH5o2hRNOiJ1K\nkqTv+/prOP98OOYYKC2NnUZKPcuDckKDBnD33bBsGZSXQ6NG0KtX7FSSJH3XH/4AX34ZjhqX8lG9\nLolbtWoVl112Ge3bt2ejjTaibdu29OvXj88++yxV+aRvNWwIDzwAxx4LffrAmDGxE0mS9B/Tp8PN\nN8NVV8G228ZOI6VHIplMJtf3N1dXV9OnTx/OOOMM2rdvz8KFCznvvPOoqanhpZde+sHfU1JSQlVV\nFcXFxesdXIVr5Uo46SQYNy7sKz366NiJJEmFbuVK2Hff8PcvvwxF7u1QnqpXeVibV155hZ/+9KfM\nmTOHLbfc8nu/bnlQKnzzDZx44n8G0rp2jZ1IklTIbrgBBgyAl16CTp1ip5HSp17bltZm0aJFJBIJ\nNvYOdqXRBhvA8OFw5JHQuzdMnBg7kSSpUM2eHbYqXXihxUH5L6UrDytWrGD//fdnt91244EHHljr\nv+PKg1JpxYowA/HMM2Eb06GHxk4kSSokySQccQR8+CG88Ua4k0jKZ3VaeXjooYdo3rw5zZs3p7i4\nmGnTpn37a6tWraJPnz4kEgnuuOOOlAeV1qZxY3jsMTjgAOjZE559NnYiSVIhuf9+ePppuOsui4MK\nQ51WHpYsWcL8+fO//Xnbtm1p3Ljxt8Xho48+4umnn2aTTTb5wY+xZuWhe/fuFP3PNFF5eTnl5eXr\n8Z+hQrdsWSgP//wnPPkk7L9/7ESSpHw3fz7sumu402Hw4NhppMyo97alNcXhgw8+YPLkybRo0eJH\n/323LSldli6FHj3glVdgwgTYb7/YiSRJ+aysDCZNgrffhk03jZ1Gyox6DUyvXr2a448/nhkzZjBk\nyBBWrlzJ/PnzmT9/PitXrkxVRqlWmjULdz/stRd06xZOvJAkKR3GjoWHH4ZbbrE4qLDUa+Vhzpw5\ntGvX7jv/LJlMkkgkmDx5MgcddND3fo8rD0q3xYtDeXjrrXAK0957x04kSconVVWwxx7hx/jxkEjE\nTiRlTsrveVgXy4MyoboajjoK3n03LCl37Bg7kSQpX5x+elh1eOMN2Gab2GmkzEr5PQ9SNiguDoPT\nO+4Y7oKYOTN2IklSPpgwAe69F266yeKgwuTKg/LaokXh/O2PPoLJk+EnP4mdSJKUq6qrw1alnXcO\nJcLtSipErjwor228cfgEv9VWcPjhYQ5CkqT1ccklsHBhWHmwOKhQWR6U91q0CIPTm28ebqB+443Y\niSRJuWbCBLjnHrcrSW5bUsH48suwhemTT8IQdYcOsRNJknJBdXXY9rrjjvDUU646qLC58qCC0bJl\nKA3bbguHHQYzZsROJEnKBb/5DXz1lduVJLA8qMCs2cK0ww5hBuKVV2InkiRls4kT4e674cYbw5tP\nUqFz25IKUlVVuEju7bfDka4//WnsRJKkbLNmu9IOO4QS4aqD5MqDClRJSSgNe+wR7oF4/vnYiSRJ\n2eb888PpSv/4h8VBWsPyoIJVXAxPPAF77QVdu8LUqbETSZKyxahRcP/9cMstbleS/pvlQQVto41g\n/HjYZ5+wjWnKlNiJJEmxzZ8Pp58OvXpB//6x00jZxfKggrfhhjB2LOy3H3TvDk8/HTuRJCmWZDIU\nhwYNwqC025Wk77I8SECzZlBZCQceCD16hHO8JUmFZ9AgGDMmHMvaunXsNFL2sTxI/1/TpjB6NBxy\nCBxzTNjOJEkqHB98ABdcAKeeGr4OSPo+y4P0X5o0CUNyXbtC794wYkTsRJKkTFi9Gvr2hVat4Oab\nY6eRspflQfofjRvDo4/CccfBSSfB0KGxE0mS0u3GG8Ox3YMHQ/PmsdNI2asodgApGzVqFEpD06Zw\nyimwdGkYoJMk5Z+XXoLf/x5++9sw+ybph1kepB/QsGG4GKhZMzjjjFAgzj8/dipJUipVV0N5OXTs\nCFdfHTuNlP0sD9KPaNAA/v73cJzrBRfAkiXwu9/FTiVJSpVzzoEFC2DChLDqLOnHWR6kdUgk4Prr\nQ4G4/PKwAnHttZ79LUm5buhQePDB8GP77WOnkXKD5UGqhUQCrroqFIjf/CasQPz1rxYIScpVH3wA\nZ50FJ58cfkiqHcuDVAeXXBJmIH7967ACceedYWuTJCl3rFwZ5hxatYLbb4+dRsot0cpDWVkZRUVF\nlJeXU15eHiuGVGdnnx0KxKmnhhWI++5zn6wk5ZKrroIZM2DqVCgujp1Gyi2JZDKZzOQDq6urKSkp\noaqqimL/H6scNnx4WOru2jX8fdOmsRNJktbl6afhiCPgT3+CAQNip5Fyj+VBqocnngiXye2zD1RW\nQklJ7ESSpB/yxRfQoQPstBNMnBiO5JZUN+7WluqhWzd46il4/XU49FD4/PPYiSRJa1NTEy79/OYb\nGDLE4iCtL8uDVE/77w9TpsDcueFm0o8/jp1IkvS/brghrBY/+CC0bRs7jZS7LA9SCrRvD9OmhRM8\n9t8f3nkndiJJ0hrPPQdXXBEu+ezWLXYaKbc58yCl0Ny5cNRRMH8+PP447L137ESSVNg+/xz22gt2\n2AEmTYIiD6mX6sWVBymFttgCnn023FR66KEweXLsRJJUuNbMOaxcCRUVFgcpFSwPUoq1aBFO8ejc\nGbp3h9GjYyeSpML05z+HQy2GDAlv7kiqP8uDlAYbbQRjx0LPnnD88TB4cOxEklRYnn4arrwSLr88\nbCeVlBqWBylNGjeGhx+G/v3Dj7/9LXYiSSoMn3wCZWVh++gf/hA7jZRf3P0npVHDhnDPPdCyJVx4\nYRikvu46SCRiJ5Ok/LRiBfTpA02ahDkH73OQUsvyIKVZIgHXXw+tW8Mll4QTme69Fxo1ip1MkvLP\nBRfAq6/C1KnQqlXsNFL+sTxIGXLxxbD55mEL0+efwyOPhNkISVJq3H8/DBwId98N++wTO42Un5x5\nkDLoZz+D8ePDO2KHHQYLFsROJEn54dVX4ayz4NRT4fTTY6eR8peXxEkRzJgBRx8NxcXwxBPQrl3s\nRJKUu776Cjp1CvNlU6eGeQdJ6eHKgxRBx47w/POQTEKXLuEdM0lS3a1eDT//OSxeDCNGWBykdLM8\nSJG0awfTpsHWW8PBB4eL5SRJdXP55TBhQjhZaZttYqeR8p/lQYqodetwkdH++4dtTA89FDuRJOWO\noUPDaXY33QRHHhk7jVQYLA9SZBttBJWVYZj65z8P90BkdhJJknLPyy+H4eh+/cLxrJIyw6NapSzQ\nqBHcdx9st11Ygv/wQ7jjDu+CkKS1mTsXeveGvfYKR7N68aaUOZYHKUskEnDVVbDttnDaafDJJzB8\neDiRSZIULF8Oxx4bPmc+9pgD0lKmuW1JyjL9+oXjW194AQ48EP7979iJJCk7JJNwxhnw+uswalS4\neFNSZlkepCx0+OHhJKaFC6FzZ5g5M3YiSYrvL3+BBx+EQYNg771jp5EKk+VBylJ77AEvvght2oQV\niCefjJ1IkuIZMwYuuwwGDIDy8thppMJleZCy2BZbwJQpoTz06AH/+EfsRJKUea++GgpDr17wxz/G\nTiMVNsuDlOU22ghGj4bTTw+D1AMGQE1N7FSSlBn//jf07Am77QZDhkADv3ORoop22lJZWRlFRUWU\nl5dT7vqj9KOKisLRrdtvD5deCu++G/b9brhh7GSSlD6LF4fiUFQU7sNp1ix2IkmJZDKz11FVV1dT\nUlJCVVUVxZ5BKdXZmgvldtop/P2WW8ZOJEmpt2pV2KY0dWo4QGKPPWInkgRuW5JyTmlp+EL6xRew\n777hllVJyjcXXhgOinjkEYuDlE0sD1IO6tABXnoJtt4aDjoofHGVpHzxt7/B3/8etmsedVTsNJL+\nm+VBylGbbQbPPBNuWj3xRLj22nCBkiTlsocfDqsOl10WLoSTlF2iDUxLqr8mTWDoUNh1V7jySnjn\nnXCca5MmsZNJUt1Nngx9+8LJJ8N118VOI2ltHJiW8sQjj0C/ftC+PYwaFVYmJClXvPEGHHBAmOUa\nNw422CB2Iklr47YlKU/06QPPPgsffxy++M6YETuRJNXOxx9Dt27Qrh2MGGFxkLKZ5UHKI3vvHQap\n27QJ7+A99FDsRJL04xYsgK5doVEjGD8e3JQgZTfLg5RnttwyrECccAL8/Ofwm9+E89IlKdtUV0P3\n7vDVV/DUU7D55rETSVoXB6alPNS0KQweDJ06wcUXw8yZMGwYtGgRO5kkBcuWwTHHwOzZMGUK7Lhj\n7ESSasOVBylPJRJw/vnhkqXp02GffcJAoiTFtnJlmNN65ZWwValDh9iJJNWW5UHKc4cfHr5Ab7QR\n7LdfGEaUpFhqaqB/f5gwAR57DLp0iZ1IUl1YHqQCsN128Pzz0KNHmIW44orwBVySMimZhHPOCdso\nH3ooDEpLyi3OPEgFYsMNwxfsjh1hwAB47TUYMgQ23jh2MkmF4oor4M474d57wxsZknKPKw9SAUkk\n4LLLwgVM06aFo11ffz12KkmF4IYbwq3RN90Ep54aO42k9WV5kApQ9+5hDqJ5c+jcGR54IHYiSfns\nL38Jb1xceWU4AU5S7rI8SAVq++3DHERZGfTrB2edBStWxE4lKd/cfDNccglcfjn84Q+x00iqL8uD\nVMCaNoV//APuvhsGDQq3Us+ZEzuVpHxxyy1w0UVhzuraa8PWyf/X3v3HRl3fcRx/FQp0A1poESyQ\nAGUDimJxpSC0/HIsDtg6f1DT6mRbCKsbmQoiuGQjY1tDUEY0gAJ1/oBJ69gMKD+3dUgHTn6JyA8B\nwy9hMAps67X0Wtve7Y93OiAwPNrrfe6uz0fyzV3To/cOAt7r+/68Px8AkY3wALRwMTHS1Kk2A3Hh\ngg1Ub97suioAkW7RIumpp2y5UkEBwQGIFoQHAJJsePrDD6Vhw2wmYs4cqb7edVUAItGSJdITT0jP\nPCPNm0dwAKIJ4QHA/yQmSuvW2fKCggI7YO7sWddVAYgkL71kZznMmCHNn09wAKIN4QHANVq1ssHG\nLVukTz+V0tKkTZtcVwUgEixZIk2bZsuVFiwgOADRiPAA4IZGjbKD5DIybBnTs89KtbWuqwIQrp5/\n3joO06dLCxcSHIBoRXgA8H/ddpstY3ruObuLOGaM9NlnrqsCEE78ftuCddYsO0H6N78hOADRjPAA\n4KZatbKhx7/9TTpzRho8WHrnHddVAQgHfr+FhrlzbTCa7ViB6Ed4ABCQ4cOlvXulkSOl73zHhiE/\n/9x1VQBc8flsvmHBAjvP4dlnXVcEIBSchYfc3FxlZ2erqKjIVQkAblFiorRmjfTCC9LixdKIEdKR\nI66rAhBqNTXSI49Iy5ZJhYW2LSuAliHG7/f7Q/mGHo9HCQkJKi8vV3x8fCjfGkAQ7d4tPfqodPq0\nDUfm57NcAWgJKiqkBx+0pYyrVtlzAC0Hy5YANErDoXKTJ0s/+pEtZSorc10VgOZ04YKd/7Jjh23h\nTHAAWh7CA4BGa99eWrpUWrtW+vvfpUGDpA0bXFcFoDmcOiVlZdmOa1u32u5rAFoewgOAJsvOlvbv\nl9LTpYkTbYiyqsp1VQCC5cABm3Gqq5O2b5fuvtt1RQBcITwACIrbb5fWr7dB6ldftSCxd6/rqgA0\n1Z//LGVmSl27WnDo29d1RQBcIjwACJqYGOs67NkjxcVJw4ZJ8+dL9fWuKwPQGIWFdsJ8VpZUWmo3\nCQC0bIQHAEE3cKANVE6fLv30p3Y2BFu6ApHD57NzG374Q9tJbe1aqWNH11UBCAeEBwDNom1b6zqU\nlgWnJFIAAAvuSURBVNoOLYMH22FSdCGA8Ob1Srm50nPP2TbMixdLsbGuqwIQLggPAJpVVpa0b59t\n5zprln19+LDrqgDcSFmZdO+90rp10ttvW/eQ81sAXI3wAKDZffnLdgeztFS6dIkuBBCOPv7Y5pRO\nnLC/q/ff77oiAOGI8AAgZLKypI8+kn78Y7oQQDj54x+l4cOlTp1sXmnIENcVAQhXhAcAIXWjLsTz\nz9OFAFzw+aQ5c6RJk6Rvf9u2Yu3Vy3VVAMIZ4QGAE1d3IWbPlu65h3MhgFDyeKQHHpB+/Wtp3jyp\nqMjCPQDcDOEBgDMNXYjt26XqaikjQ5o5U7p82XVlQHQ7eNDmG957T3r3XduWlcFoAIEgPABwbvhw\n6cMP7Q7okiV2TsT69a6rAqLTypXS0KG2/erOndLEia4rAhBJCA8AwkKbNnb388ABqX9/6Vvfkh5+\nWDp3znVlQHSorrZD3yZPlnJybDC6f3/XVQGINIQHAGGlb19p82bpzTdtSUVqqrR0qQ12AmicY8es\nw7dypfTKK9JrrzHfAKBxCA8Awk5MjPTII7aNa06OHTCXmSnt2eO6MiDyvP229LWvSZWV0gcfSFOm\nMN8AoPEIDwDCVmKiVFgobd1qQ9QZGdLUqXYKLoCbq66WnnxSeugh6RvfkHbvltLSXFcFINIRHgCE\nvVGjbKB60SI7zKpfP+nFF6XaWteVAeFp/34L28uW2d+V1aulhATXVQGIBoQHABEhNlaaNk06elTK\ny5OmT7cD5kpKXFcGhA+fz8JCRoZ9vWuX9MQTLFMCEDyEBwARpUsX6eWXbf4hMVEaN86WZZw86boy\nwK1z56QJE6SnnpIef9yCw6BBrqsCEG0IDwAi0t13S6Wl0qpVtuVkaqr0i19IVVWuKwNCy++X3nrL\ngsK+fdKmTdILL0hxca4rAxCNCA8AIlZMjC1hOnxYevppad48m4d49VWpvt51dUDzKyuTJk2ScnOl\ne++VPv5Yuu8+11UBiGaEBwARr0MHO536k0+kkSNtK8q0NOndd+2uLBBtGroNAwdaB+73v7frtttc\nVwYg2hEeAESNlBSpqMjWenfrJmVnS6NH2972QLQ4e/babsOhQ3YeCgCEgrPwkJubq+zsbBUVFbkq\nAUCUGjJE+stfpI0bpfJyO1l30iTpyBHXlQGNV18vLV4sDRggbdtmnQe6DQBCLcbvD21T3+PxKCEh\nQeXl5YqPjw/lWwNogerrpTfflH7+c+kf/7BD5ubMkZKTXVcGBG7vXik/37pq+fk239O5s+uqALRE\nLFsCENVat5YmT7auw/z5drc2JUWaMUM6f951dcDNVVban9UhQySvV9q+XVq6lOAAwB3CA4AWIS7O\ndmQ6flyaPVv67W+lPn2kmTNtxxognPh80sqVUv/+FhbmzbNT1keMcF0ZgJaO8ACgRenUyc6DOHnS\nwsTy5RYiZs2SLlxwXR0gvf++dM891jHLzLSB6FmzpDZtXFcGAIQHAC1U587Sr34lnThhJ/K+/LLU\nu7c0fbp05ozr6tASnTplOyhlZtqsTsMWrL17u64MAK4gPABo0ZKSpIICCxEzZkivv24zEVOmSEeP\nuq4OLUFlpfSzn9kuSlu3Sq+9ZoPRI0e6rgwArkd4AABJXbpYJ+LUKQsTGzbYh7mHH7a15kCw1dTY\n1qtf/aq0YIEto/v0U+n735da8X9nAGGKf54A4Crx8dIzz1gn4qWXpN27pfR0aexYae1aW04CNEVt\nrfTKKxYannxSuu8+6fBhOyW9QwfX1QHAzREeAOAG4uKkxx+3pUtvvSV9/rl0//1Sv37Siy9KHo/r\nChFpGs4cGTjQzhsZPlw6eNCWyjHXACBSEB4A4CZiY23p0vbt0o4dtgvOzJlSz542XH38uOsKEe5q\na6Xf/U666y7pu9+18PDRRxZKBwxwXR0A3BrCAwAEaOhQu3N88qT0k59IK1ZIX/mKdSQ2bmRJE65V\nVSUtWmR/Rh57TOrVywLo2rVSWprr6gCgcQgPAHCLevSwoerTp6Vly2w+YsIE26Xpl79kq9eW7l//\nsuH7Xr2sOzVypLRvnw3hDx3qujoAaJoYv9/vD+UbejweJSQkqLy8XPHx8aF8awBoFn6/ba25fLlU\nXCx5vRYmpk61x9hY1xUiFA4dkpYskd54w7pQU6bYDkp9+riuDACCh/AAAEHk8ViAWL5c2rNH6t5d\n+t73bNlKaqrr6hBsdXXSO+/Ylqtbtkjdukn5+dK0aVLXrq6rA4DgIzwAQDPZu1cqLJSKiqT//Me2\nfH3sMSkvjw+Wke7sWTvMbelSW6aWmWmB4aGHpLZtXVcHAM2H8AAAzaymRlq3Tlq50ta9+3zSmDFS\nTo70wAMEiUhRXW1dhtdflzZvltq1kx591ELD4MGuqwOA0CA8AEAIXbok/eEPdm3ZYvMSo0dfCRK3\n3+66QlytYZ7ljTesg/Tvf9v5DD/4gW3hm5DgukIACC3CAwA4cuGCtGaNBYmSEutIjBplIWLiRNvi\nE6HXEBhWr7b/NidP2g5bkyfb/Er//q4rBAB3CA8AEAYuXbL9/1evlv76VzvRul8/261pwgQLFe3a\nua4yevl80s6dVwLDZ5/ZcrIHH5QmTbJlZq1bu64SANwjPABAmKmstE7Ehg12nTkjtW8vjRsnff3r\n9kH2jjukVpzU0ySXLkl/+pMd8Ldpk3WCunWzoeecHDufgcAAANdqcniYO3euiouLdfr0abVt21bp\n6ekqKCjQ0P9zEg7hAQAC5/dLBw5cCRIffGBdiaQkm5UYPdrCxJ13Eia+SHW1dRe2bLGwsHOndRzS\n0qy7M368NGIEgQEAbqbJ4aG4uFhdu3ZVSkqKvF6vFi5cqNWrV+vYsWNKSkq67vWEBwBoPK/XAsR7\n79nVECYSE6WsLGnYMLuGDGGYt6JC2rFDKi2Vtm615zU19vsybpwFhm9+087iAAAEJujLlioqKpSQ\nkKCSkhKNHTv2uu8THgAgeLxe+1C8ZYv0/vt2N93jkWJipAED7GyJQYOuXD162PeijddrHZpdu+z3\nYNcu6ZNPrHOTlGQzI6NH2+Ndd9FdAIDGig3mD6utrdWyZcvUqVMnpaWlBfNHAwBu4EtfsmVLY8bY\n1z6fdPSoBYodO+ygujVrbI5Ckjp3tiVOgwbZY9++Up8+Uq9ekXG42eXL0rFj0qFDFhYOHrTHY8cs\nKMTGWjgYNUp6+mnrwqSmsqQLAIIlKJ2H9evXKzc3V1VVVerevbvWrFmj9PT0G76WzgMAhJbPJ506\nJe3ff+U6cEA6ckSqq7PXxMRIPXtakGi4eva0AeKuXe2xWzcpLq756vR6pX/+Uzp3zh4bnp84IR0/\nbtf581den5xsAeiOO+zxzjttfqE5awSAlu6WwsOqVauUn59vvzAmRhs3blRmZqa8Xq/OnTunixcv\nqrCwUCUlJdq5c6e6dOly3c8gPABAeKirk06ftg/nDR/Qr35eVnb9r+nYUerUSerQwZ536HDt1aaN\n3eVv3doeW7WSamttLqPhqq62pVUej1RefuVq6I40iI21Q/N695ZSUq69UlNtzgMAEFq3FB4uX76s\n81fd9unRo4fa3WDj8X79+mnKlCmaPXv2dd9rCA/jx49XbOy1q6by8vKUl5d3K/UDAJpJba1tX1pW\nZnf8Gx49HhtGrqy8/qqrk+rrrdvRcMXG2pKoq6/4eBtcbrji463DkZxsgSE52ZZYsdwIAMLLLc08\ntG/fXikpKV/4Op/Pp5qampu+pri4mM4DAISxNm1sJyJ2IwIANGjSwHRVVZUKCgqUnZ2t5ORkXbx4\nUYsXL9bZs2eVk5MTrBoBAAAAhIEmhYfWrVvr8OHDWrFihS5evKikpCRlZGRo27ZtSk1NDVaNAAAA\nAMJA0M95+CIMTAMAAACRiVE0AAAAAAEhPAAAAAAICOEBAAAAQEAIDwAAAAACQngAAAAAEBDCAwAA\nAICAEB4AAAAABCTk5zz4/X5VVFSoY8eOiomJCeVbAwAAAGiCkIcHAAAAAJGJZUsAAAAAAkJ4AAAA\nABAQwgMAAACAgBAeAAAAAASE8AAAAAAgIIQHAAAAAAEhPAAAAAAIyH8BnRDLTH1wXRsAAAAASUVO\nRK5CYII=\n", "text/plain": [ "Graphics object consisting of 1 graphics primitive" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "P.plot(1,3)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Determine the roots of this polynomial." ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[(-1/2*(I*sqrt(3) + 1)*(1/2*I*sqrt(3) + 1/2)^(1/3) - 1/2*(-I*sqrt(3) + 1)/(1/2*I*sqrt(3) + 1/2)^(1/3) + 1,\n", " 1),\n", " (-1/2*(1/2*I*sqrt(3) + 1/2)^(1/3)*(-I*sqrt(3) + 1) - 1/2*(I*sqrt(3) + 1)/(1/2*I*sqrt(3) + 1/2)^(1/3) + 1,\n", " 1),\n", " ((1/2*I*sqrt(3) + 1/2)^(1/3) + 1/(1/2*I*sqrt(3) + 1/2)^(1/3) + 1, 1)]" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "P.roots()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "What do you notice? Do you find this normal? Answer: We see radicals and complex numbers, but in the plot we saw that there is (at least) one real root." ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": false }, "outputs": [], "source": [] }, { "cell_type": "markdown", "metadata": {}, "source": [ "It is also possible to find the approximated values of roots when we specify the interval where to look for the root. Look for the function that allows to do this and determine an approximate value of the three roots by specifying the three intervals." ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "2.879385241571423" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "P.find_root(2.5,3)" ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "0.6527036446661393" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" } ], "source": [ "P.find_root(0,1)" ] }, { "cell_type": "code", "execution_count": 10, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "-0.532088886237956" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "P.find_root(-1,0)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " - multivariate polynomials\n", "\n", "Find in the documentation the way to define the ring of polynomials $R$ with two variables $y$ and $t$ with rationals coefficients." ] }, { "cell_type": "code", "execution_count": 12, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "Multivariate Polynomial Ring in y, t over Rational Field" ] }, "execution_count": 12, "metadata": {}, "output_type": "execute_result" } ], "source": [ "R. = PolynomialRing(QQ);R" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Define the polynomial $Q = y^2 - 1 - 2\\cdot t + t^3$." ] }, { "cell_type": "code", "execution_count": 19, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "t^3 + y^2 - 2*t - 1" ] }, "execution_count": 19, "metadata": {}, "output_type": "execute_result" } ], "source": [ "Q = y^2-1-2*t+t^3;Q" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Show it in 3 dimensions." ] }, { "cell_type": "code", "execution_count": 24, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/html": [ "\n", "\n" ], "text/plain": [ "Graphics3d Object" ] }, "execution_count": 24, "metadata": {}, "output_type": "execute_result" } ], "source": [ "plot3d(Q,(t,-20,20),(y,-20,20))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " - polynomial interpolation\n", " \n", " Find the (minimal degree) polynomial that passes through the following points: $[1,2]$, $[3,2]$, $[-1,-1]$ and $[4,3]$.\n" ] }, { "cell_type": "code", "execution_count": 31, "metadata": { "collapsed": false }, "outputs": [], "source": [ "T. = PolynomialRing(QQ);" ] }, { "cell_type": "code", "execution_count": 33, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "17/120*x^3 - 4/5*x^2 + 163/120*x + 13/10" ] }, "execution_count": 33, "metadata": {}, "output_type": "execute_result" } ], "source": [ "T.lagrange_polynomial([[1,2],[3,2],[-1,-1],[4,3]])" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": false }, "outputs": [], "source": [] }, { "cell_type": "markdown", "metadata": {}, "source": [ " - Polynomial factorization\n", "\n", "Find the factorizations of the polynomial $54\\cdot x^4+36\\cdot x^3-102\\cdot x^2-72\\cdot x-12$ over the following rings: $\\mathbb{Z}, \\mathbb{Q}$, and the algebraic numbers. With a bit more search, for $\\mathbb{Q}[\\sqrt{2}], GF(5)$." ] }, { "cell_type": "code", "execution_count": 47, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "(54) * (x - 1.414213562373095?) * (x + 1.414213562373095?) * (x + 0.3333333333333334?)^2" ] }, "execution_count": 47, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# A. = PolynomialRing(QQ)\n", "# A. = PolynomialRing(ZZ)\n", "A. = PolynomialRing(QQbar)\n", "F = 54*x^4 + 36*x^3 - 102*x^2 - 72*x - 12\n", "F.factor()" ] }, { "cell_type": "code", "execution_count": 51, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "(54) * (x - a) * (x + a) * (x + 1/3)^2" ] }, "execution_count": 51, "metadata": {}, "output_type": "execute_result" } ], "source": [ "t = var('t') # to reset the variable from the above 't'\n", "B. = QQ.extension(t^2 - 2); # 'a' is the new value added representing 'sqrt(2)'\n", "A. = PolynomialRing(B) # creating the polynomial ring\n", "F = A(54*x^4 + 36*x^3 - 102*x^2 - 72*x - 12) # making sure that the polynomial lives in 'A'\n", "F.factor()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Matrix calculus\n", "\n", "- Inverse of matrices" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Look in the documentation how to define a matrix with rational coefficients. Define the matrix:\n", "$$\n", "M = \n", "\\left(\n", "\\begin{array}{ccc}\n", "2 & 4 & 1 \\\\\n", "8 & 9 & 2 \\\\ \n", "3 & 7 & 3\n", "\\end{array}\n", "\\right)\n", "$$" ] }, { "cell_type": "code", "execution_count": 52, "metadata": { "collapsed": false }, "outputs": [], "source": [ "M = matrix(QQ,[[2,4,1],[8,9,2],[3,7,3]])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Do the same thing but in the Symbolic ring." ] }, { "cell_type": "code", "execution_count": 53, "metadata": { "collapsed": false }, "outputs": [], "source": [ "M = matrix(SR,[[2,4,1],[8,9,2],[3,7,3]])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Verify that it is invertible by computing the determinant:" ] }, { "cell_type": "code", "execution_count": 54, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "-17" ] }, "execution_count": 54, "metadata": {}, "output_type": "execute_result" } ], "source": [ "M.det() # or det(M), 'det' is one of the rare math. functions to be a 'universal' function" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Compute now the inverse of this matrix:" ] }, { "cell_type": "code", "execution_count": 55, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[-13/17 5/17 1/17]\n", "[ 18/17 -3/17 -4/17]\n", "[-29/17 2/17 14/17]" ] }, "execution_count": 55, "metadata": {}, "output_type": "execute_result" } ], "source": [ "M.inverse()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " - Kernel\n", " \n", " Now define the following matrix:\n", " $$\n", " N = \\left( \n", " \\begin{array}{ccc}\n", " 1 & 2 & 3 \\\\\n", " 4 & 5 & 6 \\\\\n", " 7 & 8 & 9\n", " \\end{array}\n", " \\right)\n", " $$" ] }, { "cell_type": "code", "execution_count": 56, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[1 2 3]\n", "[4 5 6]\n", "[7 8 9]" ] }, "execution_count": 56, "metadata": {}, "output_type": "execute_result" } ], "source": [ "N = matrix(3,3,range(1,10));N" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Determine now the kernel of the kernel:" ] }, { "cell_type": "code", "execution_count": 58, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "Free module of degree 3 and rank 1 over Integer Ring\n", "Echelon basis matrix:\n", "[ 1 -2 1]" ] }, "execution_count": 58, "metadata": {}, "output_type": "execute_result" } ], "source": [ "K = N.kernel();K" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Get a basis for the kernel." ] }, { "cell_type": "code", "execution_count": 59, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[\n", "(1, -2, 1)\n", "]" ] }, "execution_count": 59, "metadata": {}, "output_type": "execute_result" } ], "source": [ "K.basis()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " - Solution of linear equations\n", " \n", "Look up how to determine the solution(s) of the following linear equation:\n", " $$\n", " \\left(\\begin{array}{rrr}\n", "2 & 4 & 1 \\\\\n", "8 & 9 & 2 \\\\\n", "3 & 7 & 3\n", "\\end{array}\\right) \\times x = \n", "\\left(\\begin{array}{r}\n", "2 \\\\\n", "2 \\\\\n", "1\n", "\\end{array}\\right)\n", " $$" ] }, { "cell_type": "code", "execution_count": 62, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "(-15/17, 26/17, -40/17)" ] }, "execution_count": 62, "metadata": {}, "output_type": "execute_result" } ], "source": [ "A = matrix(3,3,[2,4,1,8,9,2,3,7,3])\n", "B = vector([2,2,1])\n", "A\\B" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " - Diagonalization of a matrix\n", " \n", " Let $M$ be the matrix\n", "$$\n", " M=\n", " \\left(\\begin{array}{rrr}\n", "0 & 1 & 3 \\\\\n", "1 & -1 & 0 \\\\\n", "1 & 0 & 1\n", "\\end{array}\\right)\n", "$$\n", " \n", " \n", " \n", " Find how to diagonalize $M$ by using the eigenmatrices and its diagonalization (if possible)." ] }, { "cell_type": "code", "execution_count": 84, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[ 0 1 3]\n", "[ 1 -1 0]\n", "[ 1 0 1]" ] }, "execution_count": 84, "metadata": {}, "output_type": "execute_result" } ], "source": [ "M = matrix(QQ,3,3,[0,1,3,1,-1,0,1,0,1]);M" ] }, { "cell_type": "code", "execution_count": 85, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[ -2 0 0]\n", "[ 0 -0.4142135623730951? 0]\n", "[ 0 0 2.414213562373095?]" ] }, "execution_count": 85, "metadata": {}, "output_type": "execute_result" } ], "source": [ "D,R = M.eigenmatrix_left();D" ] }, { "cell_type": "code", "execution_count": 88, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[ 0 1 3]\n", "[ 1 -1 0]\n", "[ 1 0 1]" ] }, "execution_count": 88, "metadata": {}, "output_type": "execute_result" } ], "source": [ "R^-1*D*R" ] }, { "cell_type": "code", "execution_count": 87, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "True" ] }, "execution_count": 87, "metadata": {}, "output_type": "execute_result" } ], "source": [ "R^-1*D*R == M" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " - Representation in Jordan form\n", " \n", "Define the following matrix $M=\n", "\\left(\\begin{array}{rrr}\n", "5 & 1 & -1 \\\\\n", "-1 & 7 & -4 \\\\\n", "0 & 0 & 3\n", "\\end{array}\\right)\n", "$.\n", "\n" ] }, { "cell_type": "code", "execution_count": 89, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[ 5 1 -1]\n", "[-1 7 -4]\n", "[ 0 0 3]" ] }, "execution_count": 89, "metadata": {}, "output_type": "execute_result" } ], "source": [ "M = matrix(QQ,3,3,[5,1,-1,-1,7,-4,0,0,3]);M" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "If the matrix is diagonalizable, we can use eigenmatrix to compute the diagonal matrix. Otherwise we should not do it!\n", "\n", "Try to diagonalize the matrix $M$:" ] }, { "cell_type": "code", "execution_count": 90, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[3 0 0]\n", "[0 6 0]\n", "[0 0 6]" ] }, "execution_count": 90, "metadata": {}, "output_type": "execute_result" } ], "source": [ "D,R = M.eigenmatrix_left();D" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Show the eigenmatrix, what do you notice?" ] }, { "cell_type": "code", "execution_count": 91, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[ 0 0 1]\n", "[ 1 -1 1]\n", "[ 0 0 0]" ] }, "execution_count": 91, "metadata": {}, "output_type": "execute_result" } ], "source": [ "R" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Here, we see that the eigenmatrix is not invertible. This is disapointing! In fact $M$ is not diagonalizable and we should not have used the method eigenmatrix_left.\n", "\n", "In fact, we must put the matrix in the Jordan form. Look for the function and give the eigenmatrix." ] }, { "cell_type": "code", "execution_count": 108, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "(\n", "[3|0 0] \n", "[-+---] [ 0 -1 1]\n", "[0|6 1] [ 1 -1 0]\n", "[0|0 6], [ 1 0 0]\n", ")" ] }, "execution_count": 108, "metadata": {}, "output_type": "execute_result" } ], "source": [ "JF,R = M.jordan_form(transformation=True);JF,R" ] }, { "cell_type": "code", "execution_count": 107, "metadata": { "collapsed": false }, "outputs": [], "source": [] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Verify that you really obtained the Jordan form." ] }, { "cell_type": "code", "execution_count": 110, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[ 5 1 -1]\n", "[-1 7 -4]\n", "[ 0 0 3]" ] }, "execution_count": 110, "metadata": {}, "output_type": "execute_result" } ], "source": [ "R*JF*R^-1" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**Exercise :** Compute $M^{11}$" ] }, { "cell_type": "code", "execution_count": 116, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[ -302330880 665127936 -665127936]\n", "[ -665127936 1027924992 -1027747845]\n", "[ 0 0 177147]" ] }, "execution_count": 116, "metadata": {}, "output_type": "execute_result" } ], "source": [ "R*JF^11*R^-1" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Derivation and Integration" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " - Define the following functions: $$sin(x)$$ \n", " $$fx^3+x-1$$\n", " $$\\ln(x^4-3)$$ \n", " Look in the documentation how to derivate a function and compute the first derivative of each functions." ] }, { "cell_type": "code", "execution_count": 120, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "(cos(x), 3*x^2 + 1, 4*x^3/(x^4 - 3))" ] }, "execution_count": 120, "metadata": {}, "output_type": "execute_result" } ], "source": [ "x = var('x')\n", "f = sin(x)\n", "g = x^3 + x - 1\n", "h = ln(x^4 - 3)\n", "\n", "f.derivative(),g.derivative(),h.derivative()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " - Compute the third derivative of the functions." ] }, { "cell_type": "code", "execution_count": 124, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "(-cos(x), 6, 128*x^9/(x^4 - 3)^3 - 144*x^5/(x^4 - 3)^2 + 24*x/(x^4 - 3))" ] }, "execution_count": 124, "metadata": {}, "output_type": "execute_result" } ], "source": [ "f.derivative(3),g.derivative(3),h.derivative(3)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " - Compute the following integral \n", " $$\n", " \\int x\\sin(x^3)dx\n", " $$" ] }, { "cell_type": "code", "execution_count": 126, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "-1/12*(x^3)^(1/3)*((sqrt(3) + I)*gamma(2/3, I*x^3) + (sqrt(3) - I)*gamma(2/3, -I*x^3))/x" ] }, "execution_count": 126, "metadata": {}, "output_type": "execute_result" } ], "source": [ "integrate(x*sin(x^3), x)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " - and this one :\n", " $$\n", " \\int_0^1 \\frac{x}{x^2+1} dx\n", " $$" ] }, { "cell_type": "code", "execution_count": 127, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "1/2*log(2)" ] }, "execution_count": 127, "metadata": {}, "output_type": "execute_result" } ], "source": [ "integrate(x/(x^2 + 1), x,0,1)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Differencial equations\n", "\n", "## Basic commands\n", "\n", " - We begin by defining a variable $t$ and a function $y$ that depends on this variable:\n" ] }, { "cell_type": "code", "execution_count": 154, "metadata": { "collapsed": false }, "outputs": [], "source": [ "t = var('t')" ] }, { "cell_type": "code", "execution_count": 155, "metadata": { "collapsed": true }, "outputs": [], "source": [ "y = function('y')(t)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Read in the documentation how the function 'desolve' works." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "For example, take the differential equation of the first order:\n", "\n", "$$y'+ 3y = e^t$$ " ] }, { "cell_type": "code", "execution_count": 156, "metadata": { "collapsed": false }, "outputs": [], "source": [ "de = diff(y,t) + 3*y == exp(t)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " - Vous pouvez obtenir la solution particulière vérifiant y(0) = 1 avec la commande :" ] }, { "cell_type": "code", "execution_count": 163, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "1/4*(e^(4*t) + 3)*e^(-3*t)" ] }, "execution_count": 163, "metadata": {}, "output_type": "execute_result" } ], "source": [ "Y = desolve(de,y,[0,1]);Y" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] }, { "cell_type": "markdown", "metadata": {}, "source": [ " **Exercise :** " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " - Resolve the following differential equation of the first order:\n", " $$\n", " y'-2y/(x-1)=(x+1)^3\n", " $$" ] }, { "cell_type": "code", "execution_count": 165, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "1/2*(x^2 + 2*_C + 10*x - 16/(x - 1) + 24*log(x - 1))*(x - 1)^2" ] }, "execution_count": 165, "metadata": {}, "output_type": "execute_result" } ], "source": [ "x = var('x')\n", "y = function('y')(x)\n", "de = diff(y,x) -2*y/(x - 1) == (x + 1)^3\n", "desolve(de,y)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " - Knowing that the second derivative of the function $x$ w.r.t. $t$ is obtained with *diff(x,t,2)*, resolve the following equation of the second order with constant coefficients:\n", " $$\n", " y\"+3y=t^2-7t+31\n", " $$" ] }, { "cell_type": "code", "execution_count": 166, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "1/3*t^2 + _K2*cos(sqrt(3)*t) + _K1*sin(sqrt(3)*t) - 7/3*t + 91/9" ] }, "execution_count": 166, "metadata": {}, "output_type": "execute_result" } ], "source": [ "t = var('t')\n", "y = function('y')(t)\n", "de = diff(y,t,2) + 3*y == t^2 - 7*t + 31\n", "desolve(de,y)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Graph theory\n", "\n", " - Digraph\n", " \n", " Look how to define the directed graph with vertices $V = \\{0, 1, 2, 3, 4\\}$ and arcs are $E = \\{ (0,1), (0,2), (0,3), (1,0), (2,4) \\}$.\n" ] }, { "cell_type": "code", "execution_count": 174, "metadata": { "collapsed": false }, "outputs": [], "source": [ "D = DiGraph([range(5),[(0,1),(0,2),(0,3),(1,0),(2,4)]])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " - Shortest path\n", "\n", "Determine the shortest path from 1 to 4." ] }, { "cell_type": "code", "execution_count": 175, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[1, 0, 2, 4]" ] }, "execution_count": 175, "metadata": {}, "output_type": "execute_result" } ], "source": [ "D.shortest_path(1,4)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Combinatorics" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " - Cartesian product\n", " \n", " Look in the documentation how to define the following cartesian products:\n", " $$\n", " \\{3,4,5\\} \\times \\{6,7,8\\} \\times \\{ 3 \\}\n", " $$\n", "\n", "List the elements of this product." ] }, { "cell_type": "code", "execution_count": 179, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[(3, 8, 3),\n", " (3, 6, 3),\n", " (3, 7, 3),\n", " (4, 8, 3),\n", " (4, 6, 3),\n", " (4, 7, 3),\n", " (5, 8, 3),\n", " (5, 6, 3),\n", " (5, 7, 3)]" ] }, "execution_count": 179, "metadata": {}, "output_type": "execute_result" } ], "source": [ "CP = cartesian_product([set([3,4,5]),set([6,7,8]),set([3])]);list(CP)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " - subsets\n", " \n", " Also, list all subsets of \\{1,2,3\\}." ] }, { "cell_type": "code", "execution_count": 181, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[[], [1], [2], [1, 2], [3], [1, 3], [2, 3], [1, 2, 3]]" ] }, "execution_count": 181, "metadata": {}, "output_type": "execute_result" } ], "source": [ "list(subsets(set([1,2,3])))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "List now the subsets of size 2 of $\\{1, 2, 3 \\}$." ] }, { "cell_type": "code", "execution_count": 185, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[{1, 2}, {1, 3}, {2, 3}]" ] }, "execution_count": 185, "metadata": {}, "output_type": "execute_result" } ], "source": [ "S = Subsets(set([1,2,3]), 2); S.list()" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "SageMath 7.5.beta3", "language": "", "name": "sagemath" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.10" } }, "nbformat": 4, "nbformat_minor": 0 }