{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# A sample Jupyter notebook" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import sympy as sy\n", "from sympy import init_printing\n", "from IPython.display import display\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAD4CAYAAADhNOGaAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO3deXhU5dnH8e+dnR1CQgIhIUAChH2JIG7sgljBBQVXbFXUFq11qVtbfdtXi0sLtloRxbrUVoFaxapEVlFZg4YlQEgIgYQtgZAACdmf94858Y0xIYGZycnM3J/rmitztpnfYck995zlEWMMSimlfJef3QGUUkrZSwuBUkr5OC0ESinl47QQKKWUj9NCoJRSPi7A7gDnIywszMTGxtodQymlPMqWLVuOGWPCa8/3yEIQGxtLcnKy3TGUUsqjiMj+uubrV0NKKeXjtBAopZSP00KglFI+TguBUkr5OC0ESinl41xSCETkTRHJFZEd9SwXEfmLiGSIyDYRGVpj2UwRSbceM12RRymlVOO5qiN4C5h0luVXAPHWYxbwKoCIhAJPASOA4cBTItLBRZmUUko1gksKgTFmLZB/llWmAu8Yhw1AexHpDEwElhtj8o0xJ4DlnL2gqCZyqqScpVsP8e6G/RwsOGN3HKWUGzXVBWVRQHaN6RxrXn3zf0REZuHoJoiJiXFPSh937HQpy3ceJSn1COsyjlNWWQXAb4GBXdsxsV8kE/tFENepjb1BlVIu1VSFQOqYZ84y/8czjVkALABITEzU0XRcJOdEMUmpR0nacYTk/flUGYgObcFtI7sxqX8k7VsGfV8cXkhK44WkNHqGt2Jiv0gm9Y9kQFQ7ROr6a1RKeYqmKgQ5QHSN6a7AIWv+6Frz1zRRJp9kjCE99zRJO46wLPUIqYdOAtAnsg2zx8YzqV8kCZ3b/OCXe1yn1tw7uieHC8+wfOdRlu04wmtrM/nbmr10aRfC5f0imdgvkgtiOxDgryeiKeVpxFVDVYpILPBfY0z/OpZdCcwGJuM4MPwXY8xw62DxFqD6LKJvgWHGmLMdbyAxMdHovYYar6rKsDWngKTUo3yReoTMY0UADI1pb33dE0lsWKtzes0TRWWs2HWUpNSjfJWeR2lFFR1aBjKhbwQT+0VycVwYIYH+7tgdpdR5EpEtxpjEH813RSEQkX/h+GQfBhzFcSZQIIAxZr44Pl6+jONAcDHwU2NMsrXtz4AnrJd6xhjz94beTwtB41RVGV5ZncF7Gw9w5GQJAX7CyJ4dubxfJJf3jSCibYhL3qeotIIv9+SRlHqEVbtyOVVaQasgfyb2i+SpKf1o1yLQJe+jlHKOWwtBU9NC0LDyyioeXbKND787yOje4UwZ1IVxfSJo19K9v5TLKqpYt/cYSalHWLIlh57hrXnnZ8Pp5KKio5Q6f1oIfEhJeSW/eO9bVu7O5aEJvZg9Ns6WA7pfpedx97tbCGsdzD/uGEFMx5ZNnkEp9f/qKwR6ZM/LFJ4p57aFm1iVlssfru7PfePibTur59L4cN67cwQnS8q5bv46dh0+aUsOpdTZaSHwIrmnSpixYAPfZZ/gLzOGcOuF3eyOxJCYDiy+eyT+Ikx/bT3JWWc9D0ApZQMtBF4iO7+Y6+evJ+tYEW/MvICrBnWxO9L34iPasOTekYS1DuaWhRtZvTvX7khKqRq0EHiBtCOnuO7VdRQUl/PeXSMY1etHQ5LarmuHliy6ZyRxnVpz1zvJfJxy0O5ISimLFgIPt2X/Ca6fvw4RWHzPSIbGNN979oW1DuZfd13IsG4d+OX7Kby9LsvuSEoptBB4tDVpudzyxkZCWwWx5J6L6BXR/O8B1CYkkLd/NpwJfSN4amkqc5fvwRPPXFPKm2gh8FBLtx7izreT6R7WisX3XER0qOecmhkS6M+rNw9l2rCuvLQynaeXplJVpcVAKbs01b2GlAu9uz6L3y1N5YLYUN6YmUjbEM+7cjfA348Xpg2kQ8tAXv9qHyeKy3nx+kEEBehnE6WamhYCD2KM4S8rM5i7Yg/jEzrx8k1DPfp+PiLCE5MTCG0VzHPLdnOypJxXbx5GiyDP3SelPJF+/PIQVVWG//lkJ3NX7OHaoVHMv2WYRxeBaiLCvaN78sdrB7B2Tx63LNxIYXG53bGU8ilaCDxAeWUVDy5K4a11WdxxSXdenDbI6273fOPwGF65aSjbcwq54bX15J4ssTuSUj7Du36beKGKyirueXcLH6Uc4pGJvfnNlQn4+XnnQDBXDOjMm7dfQPaJYq6bv45DOkSmUk1CC0Ez9+qavazcncv/TOnHL8bYc/O4pnRJfBj/vOtC8k+X8dCirXo2kVJNQAtBM7Ytp4CXVqYzZVAXZl4Ua3ecJjM4uj2/u6ov6zOP8+Y3++yOo5TX00LQTJ0pq+RXH6QQ1jqYP0z90aBvXu+GxGjGJ0Tw/LI00o6csjuOUl7NJYVARCaJSJqIZIjIY3UsnysiKdZjj4gU1FhWWWPZUlfk8QbPLdvN3rwiXrx+kNsHk2mORIQ51w2gTUgAD3yQQmlFpd2RlPJaThcCEfEHXgGuAPoCN4pI35rrGGN+ZYwZbIwZDPwV+LDG4jPVy4wxU5zN4w3W7snjrXVZ/PTiWC6JD7M7jm3CWgfz3HUD2XX4JHOXp9sdRymv5YqOYDiQYYzJNMaUAe8DU8+y/o3Av1zwvl6poLiMR5ZsJb5Tax6d1MfuOLYb3zeCG4dH89ravWzap2MZKOUOrigEUUB2jekca96PiEg3oDuwqsbsEBFJFpENInJ1fW8iIrOs9ZLz8vJcELv5Mcbw5Ec7yC8qY+70wV5xwZgr/ObKvkR3aMmDi1I4VaIXmynlaq4oBHWdz1jfOX8zgCXGmJpf+MZYY2jeBMwTkZ51bWiMWWCMSTTGJIaHN7/77bvCxymH+HTbYR4Y34v+Ue3sjtNstAoOYO70QRwqOMPvP9lpdxylvI4rCkEOEF1juitwqJ51Z1DrayFjzCHrZyawBhjigkwe52DBGX778Q6GdevAPaPqrIU+bVi3UH4+Oo7FW3JYtuOw3XGU8iquKASbgXgR6S4iQTh+2f/o7B8R6Q10ANbXmNdBRIKt52HAxYDPfeSrqjI8bF08NfeGwfh76ZXDzrp/XDz9o9ry+IfbyT2lt6BQylWcLgTGmApgNpAE7AIWGWNSReT3IlLzLKAbgffND0chSQCSRWQrsBqYY4zxuULw5jf7WJ95nN9d1ZeYjp4zrkBTCwrwY970wRSXVfLokm06oI1SLiKe+J8pMTHRJCcn2x3DJdKOnOKql7/msvhwXr9tmNffQsIV3vpmH09/spNnrunPzSO62R1HKY8hIlusY7I/oFcW26i0opIHPkihbUgAc64boEWgkW4bGcul8WH87393se9Ykd1xlPJ4WghsNHd5OrsOn2TOtQMJax1sdxyP4ecnvDDNMZrZrz5IoaKyyu5ISnk0LQQ22bQvn9fW7uXG4dGM7xthdxyPE9kuhGeu6U9KdgF/W7PX7jhKeTQtBDY4VVLOg4tSiO7Qkt9c2bfhDVSdfjKwC1cP7sJLK9PZml3Q8AZKqTppIbDB7z/ZyaGCM8ydPohWwTpstDP+Z2p/OrUJ5leLUjhTpjemU+p8aCFoYst2HGHxlhx+PjqOYd1C7Y7j8dq1CORP1w8iM6+IP36+y+44SnkkLQRNKPdUCU/8ZzsDotrxy/HxdsfxGhfFhXHHJd15Z/1+1qTl2h1HKY+jhaCJGGN4dMk2ikormDt9EIFeNvi83R6Z2JteEa359ZJtnCgqszuOUh5Ffxs1kX9uOsDqtDwev6IPcZ3a2B3H64QE+jN3+mBOFJfx5Efb9apjpc6BFoImcOB4Mf/7311cGh/GbSNj7Y7jtfp1aceDE3rz2fYjLN1a330PlVK1aSFoAs8t2w3A89MG4qc3lHOrWZf1YFB0e579bJeeRaRUI2khcLMt+0/w6fbDzLqsB53btbA7jtfz9xN+c2UCR0+WsvDrTLvjKOURtBC4kTGGZz/bRXibYGZd1sPuOD7jgthQJvaL4NU1e8k7VWp3HKWaPS0EbpSUeoQt+0/w4IReeuFYE3t0Uh9KK6p4aeUeu6Mo1expIXCTsooq5ny+m/hOrbl+WFe74/icHuGtuWlEDP/alE1G7mm74yjVrGkhcJN/btxP1vFinpicQIBeM2CLX46Lp0WgP3M+3213FKWaNZf8hhKRSSKSJiIZIvJYHctvF5E8EUmxHnfWWDZTRNKtx0xX5LFb4ZlyXlqZzkU9OzK6d7jdcXxWx9bB3Du6Jyt2HWVD5nG74yjVbDldCETEH3gFuALoC9woInXdUvMDY8xg6/GGtW0o8BQwAhgOPCUiHZzNZLdX1+yl4Ew5T0xO0MFmbHbHJd3p3C6EZz/bRVWVXmSmVF1c0REMBzKMMZnGmDLgfWBqI7edCCw3xuQbY04Ay4FJLshkm4MFZ3jzm31cMziK/lHt7I7j80IC/Xn48t5syynkk216kZlSdXFFIYgCsmtM51jzartORLaJyBIRiT7HbRGRWSKSLCLJeXl5LojtHi8mpQHw0MTeNidR1a4ZEkXfzm15flkaJeV6kZlStbmiENT13UftHvwTINYYMxBYAbx9Dts6ZhqzwBiTaIxJDA9vnt+77zhYyH++O8gdl3Qnqr1ePNZc+PkJT16ZwMGCM7yzPsvuOEo1O64oBDlAdI3prsAPenBjzHFjTPWVPa8Dwxq7racwxvDMp7sIbRXEvaN72h1H1XJxXBije4fz11UZendSpWpxRSHYDMSLSHcRCQJmAEtrriAinWtMTgGqRxBJAi4XkQ7WQeLLrXkeZ3VaLuszj/PLcfG0DQm0O46qw+NXJFBUWsFfV2XYHUWpZsXpQmCMqQBm4/gFvgtYZIxJFZHfi8gUa7X7RSRVRLYC9wO3W9vmA3/AUUw2A7+35nmUisoq/vjZbrqHteKmETF2x1H16B3ZhhsSo3l3Qxb7jxfZHUepZkM88b7tiYmJJjk52e4Y3/vnxgM88Z/tzL9lKJP6d254A2Wb3JMljHphDWP7dOKVm4faHUepJiUiW4wxibXn6yWvTioqreDPy/eQ2K0DE/tF2h1HNaBT2xBmXdaDT7cfZsv+E3bHUapZ0ELgpNfWZnLsdClPXKkXj3mKWZf1ILxNMM9+tktHMlMKLQROOXqyhNfXZnLlwM4MjfH4C6J9RqvgAB6c0Ist+0+QlHrE7jhK2U4LgRPmLt9DRVUVj07sY3cUdY6uH9aV+E6tmfP5bsoqquyOo5SttBCcp7Qjp1iUnM2tF8YS07Gl3XHUOQrw9+PxyX3IOl7MPzfutzuOUrbSQnCe/vj5LloHB3Df2Di7o6jzNKZ3Jy7q2ZGXVqZzsqTc7jhK2UYLwXn4Ov0Ya9LymD02jg6tguyOo86TiPDE5AROFJfzt9V77Y6jlG20EJyjqirHOMRdO7TgtpGxdsdRTuof1Y5rh0Tx5jf7OFhwxu44StlCC8E5+s93B9l5+CSPTOxNSKC/3XGUC1TfKfZP1p1jlfI1WgjOQUl5JS9+kcbAru24amAXu+MoF4lq34KfXdydD787yI6DhXbHUarJaSE4Bwu/3sfhwhKemJyAn59ePOZNfj6mJx1aBvLMp3qRmfI9Wgga6djpUl5ds5fxCRFc2KOj3XGUi7UNCeSX4+JZn3mc1Wm5dsdRqklpIWikv63ey5nySh67Qi8e81Y3jehG97BWPL8sTcc3Vj5FC0Ej5J4s4b2N+7l2SBRxnVrbHUe5SVCAHw+Mj2f3kVN66wnlU7QQNML8LzOpqDLM1ovHvN5PBnahR3grXlqZrl2B8hkuKQQiMklE0kQkQ0Qeq2P5gyKy0xq8fqWIdKuxrFJEUqzH0trb2i33lKMbuGZIFN06trI7jnIzfz/h/rGOruCLndoVKN/gdCEQEX/gFeAKoC9wo4j0rbXad0CiNXj9EuD5GsvOGGMGW48pNDOvVXcDY7Qb8BVXDepCj7BWzFuhXYHyDa7oCIYDGcaYTGNMGfA+MLXmCsaY1caYYmtyA45B6pu96m7g6sFRxIZpN+Ar/P2E+8bFWV3BUbvjKOV2rigEUUB2jekca1597gA+rzEdIiLJIrJBRK6ubyMRmWWtl5yXl+dc4kZa8GUmZRVVemzAB101sAvdw/RYgfINrigEdV1ZVef/HBG5BUgEXqgxO8YaQ/MmYJ6I9KxrW2PMAmNMojEmMTw83NnMDco7Vco/Nu7n6iFRdNduwOcE+Ptx39g4dh0+yfJd2hUo7+aKQpADRNeY7gocqr2SiIwHngSmGGNKq+cbYw5ZPzOBNcAQF2Ry2oK1eymrqOK+sfF2R1E2mTKoC7EdW/LSinS92lh5NVcUgs1AvIh0F5EgYAbwg7N/RGQI8BqOIpBbY34HEQm2nocBFwM7XZDJKcdOl/LuBsexAe0GfJejK4hn5+GTLNdjBcqLOV0IjDEVwGwgCdgFLDLGpIrI70Wk+iygF4DWwOJap4kmAMkishVYDcwxxtheCBas1WMDymHqYKsrWKldgfJeAa54EWPMZ8Bnteb9rsbz8fVstw4Y4IoMrnLsdCnvrt/P1MFR9AjXq4h9XYC/H7PHxvPw4q2s2JXLhL4RdkdSyuX0yuJaXl+bSWlFpXYD6ntXD+5Ct44tmbdij3YFyitpIajh2OlS3lm/nymDutBTuwFlCfD3Y/aYOFIPnWTlLr0zqfI+WghqeP2r6m5AzxRSP+S4xUhL5q3UrkB5Hy0EluOnS3ln3X6uGtRF7zCqfiTA349fjIljx8GTrNqtXYHyLloILK9/tY+Sikq9bkDV65ohUcSEtmSeXlegvIwWAiC/qIx31mcxRbsBdRaB1rGC7QcLtStQXkULAY5jA2fKK7lPzxRSDbhmaBTRoS30ugLlVXy+EOQXlfH2uiyuGtiFuE5t7I6jmrlAfz/uGxPPtpxCHdtYeQ2fLwRvWN3A/eO0G1CN831XoMcKlJfw6UJwwuoGfqLdgDoH1ccKtuYUsiataW6JrpQ7+XQheOPrTIrLK7lfjw2oc3Tt0K507dCCeXqsQHkBny0EJ4rKeOubLK4c0Jn4CO0G1Ln5vivILmDNHu0KlGfz2UKw8Ot9jm5gnF43oM7PtUO7EtVejxUoz+eTheBEURlvrcti8oDO9NJuQJ2noAA/Zo+NIyW7gC+1K1AezCcLwcKv91FUVsH9ehWxctJ11V2BHitQHsznCkFBsdUN9O9M70jtBpRzggIc9yD67kABa9OP2R1HqfPikkIgIpNEJE1EMkTksTqWB4vIB9byjSISW2PZ49b8NBGZ6Io8Z7Pw632cLq3QYwPKZaYNqz5WoHcmVZ7J6UIgIv7AK8AVQF/gRhHpW2u1O4ATxpg4YC7wnLVtXxxjHPcDJgF/s17PLQqKHWcKTR4Qqd2AcpmgAD9+PqYn3x4o4CvtCpQHckVHMBzIMMZkGmPKgPeBqbXWmQq8bT1fAowTEbHmv2+MKTXG7AMyrNdzize/3scp7QaUG1w/LJou7UJ0FDPlNhm5p/np3zex/3iRy1/bFYUgCsiuMZ1jzatzHWuw+0KgYyO3BUBEZolIsogk5+Wd3xkax4vKuHJgZ/pEtj2v7ZWqj6MriOPbAwV8naFdgXK9v65KZ0NmPq2DXTLU/A+4ohBIHfNqfySqb53GbOuYacwCY0yiMSYxPDz8HCM6PHPNAP4yY8h5batUQ65P7EqXdiF6XYFyuYzc03yy9RC3XdSNjq2DXf76rigEOUB0jemuwKH61hGRAKAdkN/IbV3K36+u2qOU84ID/Ll3TBzJ+0/wTcZxu+MoL/LyqnSCA/yZdWkPt7y+KwrBZiBeRLqLSBCOg79La62zFJhpPZ8GrDKOj0xLgRnWWUXdgXhgkwsyKWWLGxK70lmPFSgX2pt3mqVbD3HbSPd0A+CCQmB95z8bSAJ2AYuMMaki8nsRmWKtthDoKCIZwIPAY9a2qcAiYCewDPiFMabS2UxK2SU4wJ+fj+5J8v4TrNurXYFy3surMggO8Oeuy9zTDQCIJ35qSUxMNMnJyXbHUKpOpRWVjHp+DdGhLVh090gcJ8gpde4y804z/s9fcuelPXhicoLTryciW4wxibXn+9yVxUq5W3CAPz8f05PNWSdYr12BcsLLqzIICvDjLjcdG6imhUApN7ghMZrItiHM0zOI1HnKzDvNRykHufXCboS3cc+xgWpaCJRyg5BAf+4d3ZNNWfmsz9SuQJ27l1c7uoFZl/V0+3tpIVDKTaZfEE1E22DmrUi3O4ryMPuOFfHRdwe5ZYT7uwHQQqCU24QE+nPvqJ5s2pevxwrUOak+NjBrlHuPDVTTQqCUG80YHkOnNsHMW7HH7ijKQ2QdK+KjlIPcPKIbndqENMl7aiFQyo1CAh3XFWzUrkA10surMwjwE+5uom4AtBAo5XbVXcFLK7UrUGeXdayI/3zXtN0AaCFQyu2qzyDakJnPBj2DSJ1FdTdwTxN2A6CFQKkmcePwGMLbBPOSnkGk6rH/uKMbuGlEDJ3aNl03AFoIlGoS1WcQrc88zkbtClQdXl7l6AbuHeX+6wZq00KgVBO5aYTVFazUrkD90IHjxXxoUzcAWgiUajIhgf7cM6on6/YeZ9O+fLvjqGbk5dXp+PsJ99jQDYAWAqWa1M0jYghrrWcQqf934HgxH357kJuGxxBhQzcAWgiUalKOrqAH32RoV6AcXlmdgZ+fcO9oe7oB0EKgVJO7eUQ37QoUANn5xfz72xxbuwFwshCISKiILBeRdOtnhzrWGSwi60UkVUS2icj0GsveEpF9IpJiPQY7k0cpT9Ai6P+7gs1Z2hX4sldWZ+An9h0bqOZsR/AYsNIYEw+stKZrKwZuM8b0AyYB80SkfY3ljxhjBluPFCfzKOURHF1BkF5X4MOy84tZsiWHG4dHE9nOvm4AnC8EU4G3redvA1fXXsEYs8cYk249PwTkAuFOvq9SHq1FkD93X9aTrzOOkaxdgU/62xpHN3Dv6Di7ozhdCCKMMYcBrJ+dzrayiAwHgoC9NWY/Y31lNFdE6r3xtojMEpFkEUnOy8tzMrZS9rv5whg6tgrS6wp8UHZ+MYuTc5jRDLoBaEQhEJEVIrKjjsfUc3kjEekMvAv81BhTZc1+HOgDXACEAo/Wt70xZoExJtEYkxgerg2F8nwtgwK4e1QPvko/xpb92hX4kr+t2Wt1A/YeG6jWYCEwxow3xvSv4/ExcNT6BV/9iz63rtcQkbbAp8BvjDEbarz2YeNQCvwdGO6KnVLKU9xyYTc6tgrSUcx8SM6JYhYnZzP9gmg6t2thdxzA+a+GlgIzreczgY9rryAiQcB/gHeMMYtrLasuIoLj+MIOJ/Mo5VFaBgUw67LqruCE3XFUE2hu3QA4XwjmABNEJB2YYE0jIoki8oa1zg3AZcDtdZwm+p6IbAe2A2HA/zqZRymPc+vIboTqsQKfUN0N3HBBV7q0bx7dAECAMxsbY44D4+qYnwzcaT3/B/CPerYf68z7K+UNqruCOZ/vZnNWPhfEhtodSbnJy6syAPh5MzhTqCa9slipZuC2kd2IaBvMM5/uwhhjdxzlBmlHTrEoOZtbL4xtVt0AaCFQqlloGRTAQxN6k5JdwKfbD9sdR7nBHz/fRavgAO4b27y6AdBCoFSzcd2wrvSJbMPzy9Iorai0O45yoW8yjrEmLY/ZY+Lo0CrI7jg/ooVAqWbC3094fHICB/KLeXf9frvjKBepqjI88+kuotq3YOZFsXbHqZMWAqWakVG9wrk0Poy/rsqgsLjc7jjKBf7z3UF2Hj7Jryf1JiTQ3+44ddJCoFQz8/gVCZwsKeeVNRl2R1FOKimv5MUv0hjYtR1XDexid5x6aSFQqpnp26Ut1w3tylvfZJGdX2x3HOWEhV/v43BhCU9MTsDPT+yOUy8tBEo1Qw9d3gs/P3ghKc3uKOo8HT9dyqtr9jI+oRMX9uhod5yz0kKgVDPUuV0L7rykB0u3HmJrdoHdcdR5+MvKdM6UV/LYFX3sjtIgLQRKNVN3j+pBx1ZBPPOZXmTmaTLzTvPexgPMuCCauE5t7I7TIC0ESjVTbUICeWBCLzbty2f5zqN2x1Hn4LlluwkO8OOB8b3sjtIoWgiUasZmXBBNj/BWzFm2m/LKqoY3ULbbtC+fpNSj3DOqJ+Ft6h1rq1nRQqBUMxbo78fjVySQmVfE+5uz7Y6jGmCM4dnPdhHRNpg7L+1hd5xG00KgVDM3PqETw7uHMm/5Hk6V6EVmzdmn2w+Tkl3AQ5f3pkVQ87x4rC5aCJRq5kSEJycncLyojNe+zLQ7jqpHaUUlzy3bTZ/INlw3tKvdcc6JU4VAREJFZLmIpFs/O9SzXmWNQWmW1pjfXUQ2Wtt/YI1mppSqZVB0e6YM6sLrX2VyuPCM3XFUHd5dv5/s/DM8PjkB/2Z88VhdnO0IHgNWGmPigZXWdF3OGGMGW48pNeY/B8y1tj8B3OFkHqW81iMTe2MM/OmLPXZHUbUUFpfz11UZXBofxqhe4XbHOWfOFoKpwNvW87dxjDvcKNY4xWOBJeezvVK+Jjq0JbdfHMu/v81h56GTdsdRNby8Op2TJeU8MTnB7ijnxdlCEGGMOQxg/exUz3ohIpIsIhtEpPqXfUegwBhTYU3nAFH1vZGIzLJeIzkvL8/J2Ep5pl+MjqNtSCDP6kVmzUZ2fjFvr9vPdUO7ktC5rd1xzkuDhUBEVojIjjoeU8/hfWKMMYnATcA8EekJ1PUlWr3/so0xC4wxicaYxPBwz2u9lHKFdi0DuW9sHF9nHOPLPfqBqDl4PikNPz/H/aE8VYOFwBgz3hjTv47Hx8BREekMYP3Mrec1Dlk/M4E1wBDgGNBeRAKs1boCh5zeI6W83K0juxET2pI/frabyirtCuyUkl3AJ1sPceclPejcrnmNQ3wunP1qaCkw03o+E/i49goi0kFEgq3nYcDFwE7j6GtXA9POtr1S6oeCA/z59aTepB09xb+35Ngdx2cZY3j20110bBXE3aM85+KxujhbCOYAE0QkHZhgTSMiiSLyhrVOApAsIltx/OKfY4zZaS17FHhQRDJwHDNY6GQepXzClQM6Mzi6PTg6hbQAAA82SURBVC9+kUZxWUXDGyiXW77zKJuy8nlgQi/ahATaHccp4okHnBITE01ycrLdMZSyVXJWPtPmr+fBCb24f1y83XF8SnllFRPnrQUg6YHLCPT3jGtzRWSLdbz2BzwjvVLqRxJjQ5nYL4L5X+4l91SJ3XF8yvubDpCZV8TjVyR4TBE4G8/fA6V82KOT+lBWUcW8Fel2R/EZp0rKmbcineHdQxmfUN8Z855FC4FSHqxHeGtuHhHD+5sOsOfoKbvj+IRX1+zleFEZT05OwHFdrOfTQqCUh7t/XDxtWwTyyOKtOmaBm+04WMiCtZlcMySKQdHt7Y7jMloIlPJwHVsH88zVA9iaU8jLqzLsjuO1SsoreeCDFMJaB/PUVX3tjuNSWgiU8gJXDuzMtUOieHl1Bt8dOGF3HK805/PdZOSe5oXrB9K+pXfdKFkLgVJe4ump/YhsG8KDi7bqtQUu9lV6Hm+ty+L2i2K5NN77bnGjhUApL9E2JJAXrx9E1vEinv1sl91xvEZBcRkPL95Kz/BWPHZFH7vjuIUWAqW8yMieHbnzku78Y8MBVqfVeesvdY5++3Eqx0+XMW/6EEICPWf4yXOhhUApL/PQ5b3pHdGGXy/ZRn5Rmd1xPNrHKQf5ZOshHhgfz4Cu7eyO4zZaCJTyMiGB/sydPpjC4nKe+HC7jltwng4VnOE3H+1gaEx77hnV0+44bqWFQCkv1LdLWx66vBfLUo/w728P2h3H41RVGR5evJXKKsPc6YMJ8ILbSJyNd++dUj7szkt7MLx7KE8vTSU7v9juOB7l7+uyWLf3OL/7SV+6dWxldxy300KglJfy9xP+dP0gAB5atFUHsWmkPUdP8dyy3YxP6MT0C6LtjtMktBAo5cWiQ1vy9JR+bMrK542vMu2O0+yVVVTxwPsptAkO4I/XDvSaewk1RAuBUl7uuqFRTOoXyYtfpLHz0Em74zRr81bsYefhk/zx2gGEtwm2O06TcaoQiEioiCwXkXTrZ4c61hkjIik1HiUicrW17C0R2Vdj2WBn8iilfkxEePbaAbRvGcSvPkihpLzS7kjN0uasfOZ/uZfpidFc3i/S7jhNytmO4DFgpTEmHlhpTf+AMWa1MWawMWYwMBYoBr6oscoj1cuNMSlO5lFK1SG0VRDPTxtI2tFT/OmLNLvjNDunSsr51QcpdO3Qkt962Q3lGsPZQjAVeNt6/jZwdQPrTwM+N8boKQxKNbExvTtxy4UxvPH1PtbtPWZ3nGblD//dyaGCM/z5hkG0Dg6wO06Tc7YQRBhjDgNYPxsarmcG8K9a854RkW0iMldE6v1STkRmiUiyiCTn5eU5l1opH/XE5ARiO7bi4UVbOVlSbnecZiEp9QiLknO4d3RPEmND7Y5jiwYLgYisEJEddTymnssbiUhnYACQVGP240Af4AIgFHi0vu2NMQuMMYnGmMTwcO+7+59STaFlUAB/vmEQR0+V8vTHqXbHsV3eqVIe/3A7/bq05ZfjetkdxzYN9kDGmPH1LRORoyLS2Rhz2PpFf7a7XN0A/McY8/3HkOpuAigVkb8DDzcyt1LqPA2J6cDsMXG8tDKdcQkRXDmws92RbGGM4bF/b+N0aQXzpg8mKMB3T6J0ds+XAjOt5zOBj8+y7o3U+lrIKh6I42Tdq4EdTuZRSjXC7LFxDOrajic/2s7RkyV2x7HFvzZls3J3Lo9N6kN8RBu749jK2UIwB5ggIunABGsaEUkUkTeqVxKRWCAa+LLW9u+JyHZgOxAG/K+TeZRSjRDo78fc6YMpKa/k4cVbfe7GdFnHivjDf3dySVwYt18Ua3cc2zl1eNwYcxwYV8f8ZODOGtNZQFQd64115v2VUuevR3hrnryyL7/9aAevrc30+jtsVisqreCXH6QQ6C+8cP1A/Px84+rhs/G986SUUt+7ZUQM6/ceY87nuymrqOK+sXFefVuFE0Vl3P7WZnYcLOSVm4bQuV0LuyM1C1oIlPJhIsJLM4YQEriNPy/fQ35RGb/7SV+v/JR8uPAMty7cxIH8Yl69eajPXT18NloIlPJxgf5+vDhtEO1bBPHmN/soPFPO89MGEuhF9+DPzDvNrQs3UXimnHd+NpwLe3S0O1KzooVAKYWfn/DbnyTQsXUQLySlUXimnFduGkqLIM8fo3fHwUJmvrkJgPdnXUj/KO8dcvJ8eU/JV0o5RUT4xZg4nrmmP6vTcrntzY0UnvHsq4/X7z3OjAUbCAn0Z/E9I7UI1EMLgVLqB24e0Y2/3jiElOwCpr+2ntxTnnmdQVLqEWb+fROd24Ww5N6R9AhvbXekZksLgVLqR34ysAsLZ17A/uPFXD9/PQeOe9Z9IhclZ3PvP7bQt3NbFt09Us8OaoAWAqVUnS7rFc57d42g8Ew5181fx+4jnjGozetrM/n1km1cHBfGe3eOoEOrILsjNXtaCJRS9Roa04FFd4/ET+CG+evZsj/f7kj1Msbw3LLdPPPZLq4c0Jk3ZibSygdvKX0+tBAopc6qV0QbltxzER1bB3PzGxtZnXa2e0vao7LK8PiH23l1zV5uGhHDX24cQnCA55/x1FS0ECilGhQd2pLF94ykZ3hr7no7mY9TDtod6XulFZXM/ue3vL85m9lj4njm6v74e+EFce6khUAp1ShhrYP516wLGdatAw98kMI767PsjsTp0gp+9tZmPt9xhN9cmcDDE3t79S0y3EULgVKq0dqGBPL2z4Yzrk8Ev/s4lXkr9th259L8ojJufn0DGzLz+dP1g7jz0h625PAGeiRFKXVOQgL9mX/LUB7993bmrUhnb14R0xOjGdEjtEluS5FfVMaKnUeZv3YvB0+c4bVbhjG+b4Tb39ebaSFQSp2zAH8/Xpg2kM7tQlj49T4+2XqIdi0CGZfQiYn9IrksPtylt6c4WHCGL1KPkJR6hE378qkyEB3agnd+NpwRet8gp4knDkiRmJhokpOT7Y6hlALOlFWyNj2PpNQjrNyVS+GZcloE+jOqVzgT+0cwtk8E7VoEnvPrZuSeIin1KEmpR9iWUwhAr4jWTOwXycR+kfTr0laPB5wjEdlijEn80XxnCoGIXA88DSQAw60BaepabxLwEuAPvGGMqR7JrDvwPo6B678FbjXGlDX0vloIlGqeyiur2JiZT5L16T33VCkBfsLInh2Z2C+Sy/tF0KlNSJ3bGmPYfrCQZTsc2+7NKwJgcHR765d/hN4mwknuKgQJQBXwGvBwXYVARPyBPTiGsswBNgM3GmN2isgi4ENjzPsiMh/Yaox5taH31UKgVPNXVWVIySlwFIUdR8g6XoyI4yK1if0imNgvkqj2LdicdYKk1CN8kXqEQ4Ul+PsJI7qHMql/JBP6RujtIVzILYWgxouvof5CMBJ42hgz0Zp+3Fo0B8gDIo0xFbXXOxstBEp5FmMMe46e/v7T/s7DjttVtAzyp7iskuAAPy6ND2divwjGJ0TobSHcpL5C0BQHi6OA7BrTOcAIoCNQYIypqDH/R+MaVxORWcAsgJiYGPckVUq5hYjQO7INvSPb8Mvx8WTnF3//9c+l8WGM6hWut4OwUYN/8iKyAqhrTLcnjTEfN+I96jqaY84yv07GmAXAAnB0BI14X6VUMxUd2lLP+29GGiwExpjxTr5HDhBdY7orcAg4BrQXkQCrK6ier5RSqgk1xZXFm4F4EekuIkHADGCpcRycWA1Ms9abCTSmw1BKKeVCThUCEblGRHKAkcCnIpJkze8iIp8BWJ/2ZwNJwC5gkTEm1XqJR4EHRSQDxzGDhc7kUUopde70gjKllPIR9Z01pDedU0opH6eFQCmlfJwWAqWU8nFaCJRSysd55MFiEckD9p/n5mE4rmHwJbrPvkH32fs5u7/djDHhtWd6ZCFwhogk13XU3JvpPvsG3Wfv56791a+GlFLKx2khUEopH+eLhWCB3QFsoPvsG3SfvZ9b9tfnjhEopZT6IV/sCJRSStWghUAppXyc1xYCEZkkImkikiEij9WxPFhEPrCWbxSR2KZP6VqN2OcHRWSniGwTkZUi0s2OnK7U0D7XWG+aiBgR8ehTDRuzvyJyg/X3nCoi/2zqjK7WiH/XMSKyWkS+s/5tT7YjpyuJyJsikisiO+pZLiLyF+vPZJuIDHXqDY0xXvcA/IG9QA8gCNgK9K21zs+B+dbzGcAHdudugn0eA7S0nt/rC/tsrdcGWAtsABLtzu3mv+N44DuggzXdye7cTbDPC4B7red9gSy7c7tgvy8DhgI76lk+Gfgcx0iPFwIbnXk/b+0IhgMZxphMY0wZ8D4wtdY6U4G3redLgHEiUtfwmZ6iwX02xqw2xhRbkxtwjArnyRrz9wzwB+B5oKQpw7lBY/b3LuAVY8wJAGNMbhNndLXG7LMB2lrP2+EFIx0aY9YC+WdZZSrwjnHYgGO0x87n+37eWgiigOwa0znWvDrXMY7BcwpxDI7jqRqzzzXdgeMThSdrcJ9FZAgQbYz5b1MGc5PG/B33AnqJyDciskFEJjVZOvdozD4/DdxiDZL1GXBf00Sz1bn+fz+rBscs9lB1fbKvfZ5sY9bxJI3eHxG5BUgERrk1kfuddZ9FxA+YC9zeVIHcrDF/xwE4vh4ajaPj+0pE+htjCtyczV0as883Am8ZY/4kIiOBd619rnJ/PNu49PeXt3YEOUB0jemu/Lhd/H4dEQnA0VKerRVr7hqzz4jIeOBJYIoxprSJsrlLQ/vcBugPrBGRLBzfpS714APGjf13/bExptwYsw9Iw1EYPFVj9vkOYBGAMWY9EILj5mzerFH/3xvLWwvBZiBeRLqLSBCOg8FLa62zFJhpPZ8GrDLWURgP1eA+W1+TvIajCHj6d8fQwD4bYwqNMWHGmFhjTCyO4yJTjDGeOs5pY/5df4TjpABEJAzHV0WZTZrStRqzzweAcQAikoCjEOQ1acqmtxS4zTp76EKg0Bhz+HxfzCu/GjLGVIjIbCAJx1kHbxpjUkXk90CyMWYpsBBHC5mBoxOYYV9i5zVyn18AWgOLrePiB4wxU2wL7aRG7rPXaOT+JgGXi8hOoBJ4xBhz3L7UzmnkPj8EvC4iv8Lx9cjtHv6hDhH5F46v98KsYx9PAYEAxpj5OI6FTAYygGLgp069n4f/eSmllHKSt341pJRSqpG0ECillI/TQqCUUj5OC4FSSvk4LQRKKeXjtBAopZSP00KglFI+7v8A5RRxFEl0mAgAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "xn = np.linspace(0., 1., 21)\n", "yn = np.sin(2.0*np.pi*xn)\n", "plt.plot(xn,yn);" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/latex": [ "$\\displaystyle x^{2} - x y^{3} + y$" ], "text/plain": [ "x**2 - x*y**3 + y" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "xs, ys = sy.var('x y', real=True)\n", "mex = xs**2 + ys - ys**3 * xs\n", "mex" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "math expression $$\\int_0^\\infty f(x)dx$$." ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/latex": [ "$\\displaystyle x^{2} - x y^{3} + y$" ], "text/plain": [ "x**2 - x*y**3 + y" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "display(mex)" ] } ], "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.6" } }, "nbformat": 4, "nbformat_minor": 4 }