{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Jupyter Notebook for turning in solutions to the problems in the Essentials of Paleomagnetism Textbook by L. Tauxe" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Problems in Chapter 3" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Problem 1a" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "To make a plot, we need to import the plotting package matplotlib.pyplot and tell ipython to plot inside the notebook:" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt # import the plotting module\n", "%matplotlib inline \n", "# This allows us to plot in the notebook environment" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Write your equation for magnetic energy in words (and Latex) here. " ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY8AAAEZCAYAAABvpam5AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XeYVOXZx/HvTRNQkGJDBLFGsKIICgKrqCDBBqKioKiJ\nJbFFY3xFjeRNoiHGGkt4VeyCKBZQBDGyCILYsSBGo1iwK1gABeR+/3jOyLDuLrOzc+bMzP4+1zXX\nzDlz5px7hmHufbq5OyIiIjVRL+kARESk+Ch5iIhIjSl5iIhIjSl5iIhIjSl5iIhIjSl5iIhIjSl5\nSNExswvM7Kak44jDut6bmQ03s5lZnrvczE7KPjqRNZQ8JK+iH8fJFfa9VcW+Iys7h7tf5u6/jo7r\nYGarzazovstmVmZmH6Tvi/m9eXQTqbWi+w8nRW8G0N3MDMDM2gANgN1SP5LRvm2Apyq+2MzqV3Fe\niyfc3DCzBrV5ec4CyQGLJB2HJEvJQ/LteaAhsFu03ROYDvynwr633f0TMxtpZveb2Z1m9jUwPNp3\nZ3RsKsEsMbNvzawbgJmdaGbzzewrM5tiZu0rC8bMHjWz0yvse8XMDo0edzez58xsiZk9a2Z7px23\n0Mz6pG3/FFdaqeFEM3sPeKLCNdYHHgM2j+L+xszaVPPevjGzvSqJfwczm2ZmX5rZAjMbXNn7TNPB\nzGZF55tqZq3TzrWXmc02s8Vm9rKZ9U57rtzM/mJmTwNLga2yuLaUECUPySt3XwHMBVI/TL2AmcCs\n6HFqX3qp4xDgPnffELibtateekb3G7p7M3efG/3wXwAcDmwUnX9sFSHdBgxNbZjZrsDmwKNm1gp4\nFLgaaAVcGe1vmXo7FWKprEqoF7AD0LfC57AU6Ad8FMXd3N0/rua9NXf3Z9LPESWgacBdwMbA0cAN\nZtaxivdqwDHAcGAToBHw++hcbYFHgP9195bR/gnpyYXwOf0K2AD4sobXlhKj5CFJmMGaRLEPIVHM\nTNvXMzomZba7TwRw9+9ZuxqnsuqTU4HL3P1Nd18NXEaoFmtXybGTgO3NbJtoexgwzt1XAb8E3nT3\nu919tbuPAxYAB1fxviqLZaS7L3f3HzI8fl3vLd0A4F13vz2K72XgAaCqEoADY9z97ehzHM+a0t5Q\nYLK7TwFw9ycIpcRfpr32Nnd/I/pM+9Xw2lJilDwkCU8B+0R/wW/s7v8F5hDaQloCO7J2yePDGp5/\nS+CaqPplMeGvZIC2FQ9M+xEdFtXjHw2kqo02B96v8JL3KjtPNT5Y9yFZ2xLolnqf0Xs9Bti0mtd8\nkvZ4OaEUkTrX4Arn6gFslnZ8+nvJ5tpSQmrTiCeSrWeADYFfA08DuPs3ZvYRcDKhKue96NjKegit\nq6rofeDP7l5VVVVFtwN3RLEsc/e50f5FwMAKx25JaKuAUPe/ftpzm/Fz1fVuquy5db23dO8DM9z9\nwHUcl4n3gTvd/eRqjkmPJ5fXliKkkofknbsvJ1SJnMPaJYxZ0b70Kqt1Ve18Dqwm9M5K+Rcwwsw6\nAZjZhtU15rr7HMIP4z8ISSRlMqFKa4iZNTCzowjtF49Ez78MHB091wUYRM26wn4KtDaz5jV4b+ke\njeIbamYNo9ueZrZDNdesqirsLuBgMzvQzOqbWeOoK3HbKl77SBbXlhKi5CFJmUFoaJ2Vtm8moYE7\nPaFUVfJwAHdfBvwVeDqqPunq7g8Bo4BxUQ+tV6nQYF2JO4CdCT+iROf+itCucC7wBaEReUC0H+Bi\nwg/7YmAkoTG/YpxVcvcFhIb8d6JeYW2qeW9fRT3J0p//FjiQUNW2CPiY0L7TqLrLVnicOteHwKHA\nCOAzQsniXNZOGD+91t2/y+LaUkIsycWgzGwMoUHuM3ffuYpjrgUOApYBw939pTyGKHWEmQ0Dfu3u\nvdZ5sIgkXvK4ldBro1Jm1h/Y1t23I9SF35ivwKTuMLOmwG+B/0s6FpFikWjycPeZhCJ/VQ4hNGYS\nNWK2MDP15pCcMbO+hGqaj4F7Eg5HpGgUem+rtqzdPfBDYAtCQ6NIrbn7VNZ0VxWRDCVdbZWJir1D\nNLGbiEjCCr3ksQhIHxW8RbRvLWamhCIikgV3z2qSy0IveUwEjoMwaRuwxN0rrbJy94xvq1c7X3zh\nzJ3r3HOPM3Kkc9hhzpZbOs2aOT16OKef7txxh7NwYebnTb9dcsklWb0uzptiKu64FJNiyvWtNhIt\neZjZWMIEeRtZWNfgEsKMq7j7aHefbGb9zextwmjeE3JzXWjdOty6dl37ucWL4eWX4aWXYOJE+P3v\noUkT6NULevcO99tuG84hIlJXJZo83H1IBsecvq5jcqllS9h333AL14c334QZM+DJJ+GPfwzJZMCA\ncOvVCxppWJSI1DGFXm2VODPYYQc45RS4+2748EO4/37YaCO46CLYZBMYPBjuugu+/XbN68rKyhKL\nuSqKKXOFGJdiyoxiyo9ER5jnipl5Uu/j00/hscdgwgR46ino2xeGDIGDDoLGjRMJSUQkI2aGZ9lg\nruSRQ19+GZLI2LEwbx4cdhgMHw49e6qNREQKj5JHgSSPdIsWwbhxcPPNIXGcfDIcdxy0apV0ZCIi\ngZJHASaPFHeYNQtGj4ZHHoGDD4ZTT4Xu3VUaEZFkKXkUcPJI9+WXcMcdcMMNoQRy7rkwcCA0KPSh\nmiJSkpQ8iiR5pPz4I0yaBFdcEXpvnX02nHgiNGuWdGQiUpfUJnmoq24C6tcPjekzZ4Z2kaefhq22\ngosvhq++WvfrRUSSpuSRsG7dYPx4ePbZ0O13++2VRESk8Cl5FIitt4b/+z94/nklEREpfEoeBaZD\nh58nkb/+FZYtSzoyEZE1lDwKVCqJPPMMvPJKSCK33BIa20VEkqbeVkVi7lw477xQjTVqFPTvr3Ei\nIlI76qpbB5IHhAGHjzwC558PbdrAtdfCjjsmHZWIFCt11a0jzMII9VdeCV19y8rgnHPg66+TjkxE\n6holjyLUoAGccQa8/npIHB07hpHrq1cnHZmI1BWqtioBc+fC6aeHRan+9S/YeeekIxKRYqBqqzqu\nW7eQQI4/Hvr0CYtUff990lGJSClT8igR9eqFad/nzYMFC2DXXcPSuSIicVC1VYl6+OFQldWvH/z9\n72FtdhGRdKq2kp859NDQoN6oUWgDmTw56YhEpJSo5FEHlJfDCSfA/vuHaeCbN086IhEpBCp5SLXK\nysLYEDPYZReYPj3piESk2KnkUcdMnhwa1gcOhL/9DZo2TToiEUmKSh6Ssf79Qynkiy+gS5fQO0tE\npKaUPOqgVq3gnnvgggtCO8i114Z5s0REMqVqqzruv/+FY46BjTaCW2+FTTZJOiIRyRdVW0nWttkG\nZs0Kgwp32w2mTk06IhEpBip5yE+mT4dhw8I0J3/6U5iAUURKl9bzUPLImc8+g2OPhVWrQrtImzZJ\nRyQicVG1leTMJpvAlClhbEiXLhoTIiKVU8lDqjRtGhx3HPz2tzBiRJh8UURKh6qtlDxis2gRHH10\nmNLk7ruhRYukIxKRXFG1lcSmbVt48snQK6trV5g/P+mIRKQQKHnIOjVsGAYSjhgBvXvDgw8mHZGI\nJE3VVlIjzz0HgwbB8OEwcqTaQUSKmdo8lDzy6tNP4YgjQvvHXXfBhhsmHZGIZENtHpJXm24K//43\ntG8f1k9fsCDpiEQk35Q8JCuNGsH118N550GvXmFsiIjUHetMHma2vpldbGY3RdvbmdmA+EOTYnDS\nSfDAA6EN5Prrk45GRPIlk5LHrcAKoHu0/RHw19gikqKzzz7w9NNw3XVw5plhahMRKW2ZJI9t3H0U\nIYHg7kvjDUmK0TbbwJw58MYbcMgh8M03SUckInHKJHn8YGZNUhtmtg3wQ3whSbFq0SIsc9u+PfTo\nAe+9l3REIhKXTJLHSGAKsIWZ3QM8CZwfZ1BSvBo2hBtvhBNPhL33hrlzk45IROJQ7TgPM6sHDAb+\nDewV7Z7r7p/nIbaMaZxHYZo0KTSo33BDGBciIoUl1kGCZvaCu++RVWR5ouRRuF5+GQYMCF16zzor\n6WhEJF3cyeNvwBfAvcBPjeXu/lU2F4yDkkdhe+896NcvJJFRozSliUihiDt5LAR+dpC7b5XNBSuc\nux9wNVAfuDnq1ZX+fBnwMPBOtGuCu/+lkvMoeRS4r74KvbDatYPbboP11ks6IhEpyrmtzKw+8Caw\nP7AIeA4Y4u5vpB1TBpzj7oes41xKHkVg+XIYOhQWLw4DC7U2iEiyYp3bKsYR5l2Bt919obuvBMYB\nh1YWQg6uJQWgSRMYPx522gl69oQPP0w6IhHJVpIjzNsCH6RtfxjtS+dAdzObZ2aTzaxTDq4rCapf\nH665Jixv26MHvP560hGJSDYaZHDMNu5+pJkdDWGEuVlOCgOZ1DO9CLRz92VmdhDwELB9ZQeOHDny\np8dlZWWUlZXlIESJg1nofdW2Ley3X6jC6tEj6ahESl95eTnl5eU5OVcmDeazgT7AbHfvHI0wH+vu\nXWt1YbO9gJHu3i/avgBYXbHRvMJr3gX2qNjTS20exWvqVBg2DG6/HQ46KOloROqWuNfzGEk8I8yf\nB7Yzsw5m1gg4CpiYfoCZbWpRMcfMuhKSXcF0EZba69sXHn44zMo7dmzS0YhIptZZbeXuj5vZi6wZ\nYX6mu39R2wu7+yozOx2YSuiqe4u7v2Fmp0TPjwaOAE4zs1XAMuDo2l5XCs/ee4fFpfr1Cz2xfvOb\npCMSkXXJqKuumbUFOhCSjQO4+1OxRlYDqrYqDe++CwccAMcfDxddFNpGRCQ+cQ8SHEWoUpoP/Jja\n7+4HZ3PBOCh5lI5PPgklkLIyuPJKjUYXiVPcyeM/wM7uXrDTsCt5lJYlS+Dgg6FDBxgzJszUKyK5\nF3eD+X+BRtmcXCQbLVqEXlhffQVHHgk/FOyfLSJ1V5UlDzP7Z/Rwc2A3wrTsqf/G7u5nxh9eZlTy\nKE0rVsAxx8DSpTBhAjRtmnREIqUllmorMxvOmoF8VvGxu9+ezQXjoORRulatghNOCFOZTJwIzZol\nHZFI6SjKiRFzScmjtK1eDaeeCq+9Fpa51YSKIrkRd5uHSKLq1YPRo2HPPaFPH/ii1qOMRKS2lDyk\nKJjB1VeHEellZaFLr4gkJ5OJEUUKghlcemloOO/VK4xKb9cu6ahE6qYqk4eZTarmdb6uBZpE4nLR\nRSGB9O4NTzwBW2+ddEQidU91JY8rqnlOrdOSqHPOCYtL9e4N06bBDjskHZFI3VJl8nD38tRjM2tK\nWFfjzXwEJZKJ004LCaRPn1AC6dgx6YhE6o5MlqE9BHiJMPstZtbZzCZW/yqR/Bg+HP72t5BA5s9P\nOhqRuiOTBvORQDdgOoC7v2RmqmWWgjFsWGhM33//UIW1445JRyRS+jJJHivdfUmFpWdXxxSPSFaG\nDg3jQQ44AB5/HHbaKemIREpbJsnjdTM7FmhgZtsBZwKz4w1LpOaOOSaUQFIJZOedk45IpHRlMkjw\nDGBHwqSIY4FvgLPjDEokW0OGhMGEBx4I8+YlHY1I6cp4bisz25AwvuObeEOqOc1tJRWNHw9nnglT\npsBuuyUdjUhhqs3cVuustjKzPYExQPNoewlwkrs/n80FRfLhyCNDFVa/fvDYY9C5c9IRiZSWTNo8\nxgC/cfeZAGa2T7RvlzgDE6mtwYNDI3oqgey+e9IRiZSOTJLHqlTiAHD3WWa2KsaYRHJm0KBQAjno\noFCFpRKISG5UN7fVHtHDGWY2mtBYDnAUMCPuwERyZeDAsCbIQQeFcSDqhSVSe9WtJFhO9SsJ7ht7\ndBlSg7lk4t574eyzw2y8nTolHY1I8mJpMHf3sqwjEilARx0FK1eGcSD//rcmUxSpjYzW8zCzAUAn\noHFqn7v/b1xBicRl6NCwLvr++8P06bDddklHJFKcMumqOxpoAuwH3AQMBubGHJdIbIYPDwmkTx8o\nL9d6ICLZWOcgQTN71d13NrNX3H0XM9sAmOLu++QnxHVTm4dk48YbYdSokEA6dEg6GpH8i3WQILA8\nul9mZm2BL4HNsrmYSCE57bRQAtlvP5gxQ0vaitREJsljkpm1BC4HXoj23RRfSCL5c8YZoRF9331D\nAmnbNumIRIpDxnNbAZhZY0Kj+Q7u/kxsUdWQqq2ktkaNgjFjQhVWmzZJRyOSH7WptqpR8ki74Afu\nXjCFfCUPyYW//AXuuSf0wtp006SjEYlf3G0eInXCRReFKqxUN96NNko6IpHCpeQhkmbkSFixAvr2\nDQMJW7RIOiKRwlTd9CSTqnldH3dvGk9INadqK8kl9zCNybPPhhUJmzVLOiKReMTS5mFmZdW8zt29\nYCZHVPKQXHOHU06Bt96CRx+FpgXzp5JI7uS9wbzQKHlIHH78MYxG//xzePhhWG+9pCMSyS0lDyUP\nicmqVXD00eH+vvugYcOkIxLJndokj3q5DkaklDRoELrvrlwJxx0XSiMiUsPkYWb1zax5XMGIFKJG\njWDChFB99atfhYWlROq6dSYPMxtrZs3NbH3gVeANM/tD/KGJFI7GjUO7x9tvhylNVEsqdV0mJY9O\n7v4NcBjwGNABGBZnUCKFaP31Q8+r556D885TApG6LZPk0cDMGhKSxyR3X8maJWlF6pTmzWHKlLAW\n+siRSUcjkpxMksdoYCGwAfCUmXUAvo4vJJHC1qpVSB7jx4cJFUXqohp31TUzA+q7+6p4Qqo5ddWV\nJCxaBL17w1lnhXYQkWITy8SIZjbM3e80s3NZU02VuogDV2ZzQZFS0bYtPPFESCBNmoSeWCJ1RXUT\nI6YmZGiG2jhEKtWhQ0gg++4bEsixxyYdkUh+ZLKG+c/qhMxsPXf/IdbIakDVVpK0+fOhTx+47joY\nNCjpaEQyE/cI83Iz2yrtYl2B57K5WEVm1s/MFpjZW2Z2fhXHXBs9P8/MOufiuiK51qkTTJ4Mv/lN\n6M4rUuoySR6XAo+Z2W/N7FJC76vhtb2wmdUHrgP6AZ2AIWbWscIx/YFt3X074GTgxtpeVyQunTvD\nxIlwwglhLRCRUrbO5OHuU4HTgGuAE4CD3P3FHFy7K/C2uy+Mxo6MAw6tcMwhwO1RHHOBFmamBUKl\nYHXrBvffHyZTnDUr6WhE4pPJ9CQXA/8EegIjgRlmNiAH124LfJC2/WG0b13HbJGDa4vEplcvuPtu\nGDgwjEYXKUWZLEPbGtjT3ZcDc8xsCnAz8Egtr51pC3fFxpxKXzcybbhvWVkZZWVlWQUlkgsHHgi3\n3AIDBoQBhbvsknREIlBeXk55eXlOzpXYeh5mthcw0t37RdsXAKvdfVTaMf8Cyt19XLS9AOjt7p9W\nOJd6W0lBGj8+LGn75JOwww5JRyOytlgGCaadfBPgD8COQONot7v7ftlcMM3zwHbRdCcfAUcBQyoc\nMxE4HRgXJZslFROHSCE78kj4/ns44AAoL4dttkk6IpHcyKTa6m7gXmAAcAqhp9Xntb2wu68ys9OB\nqUB94BZ3f8PMTomeH+3uk82sv5m9DSwlNNiLFJXjjoPly2H//WHGDGjfPumIRGovk0GCL7r77mb2\nirvvEu173t275CXCDKjaSorB1VfD9dfDU09BmzZJRyMSc7UVsCK6/yTqZfUR0DKbi4nUZWefDUuX\nhhJIeTlsvHHSEYlkL5Pk8VczawGcS+iy2xz4XaxRiZSoCy+EZctCb6wnn4SW+jNMilRiva1ySdVW\nUkzc4ZxzYM6c0I23WbOkI5K6Ku65rdIvlIuR5SJ1mhlceSXstlsYB7JsWdIRidRcjZIHPx+wJyJZ\nMIMbboAtt4TDDgvdeUWKSZXJw8weS59NNzI55nhE6ox69WDMGGjRIowHWbky6YhEMlddyWMMMNXM\nLjSzhgDufmF+whKpGxo0gLvuCu0gQ4fCjz8mHZFIZqptMDezDYA/An2BO1kzr5S7e8EsQ6sGcyl2\n338PBx8clrYdMyaUSkTiFmeD+UrgO8K0JM2ADaKb+oeI5FDjxvDQQ/DOO/Db34aSiEghq7LkYWb9\ngCuBScCf3L1g+4So5CGl4ptvwjxYPXrAFVeEhnWRuNSm5FFd8pgJnOrur9cmuHxQ8pBSsngx7Ldf\n6Mb75z8nHY2UsriSR9H8IhdRqCIZ+fxz6N07NKKPGJF0NFKqYpnbSr/GIsnZeOOwDnqvXtC0aZgX\nS6SQZDK3lYgkoE2btRPIyScnHZHIGtUmDzNrAExz933zFI+IpGnfHp54AsrKoEkTGDYs6YhEgmqT\nR7Rg02oza+HuS/IVlIisse22YQLF/fYLXXoHD046IpHMqq2WAq+a2eNAqruuu/uZ8YUlIuk6doQp\nU8JU7o0awaGHJh2R1HWZJI8Hols6NaaL5Nmuu8Kjj8Ivfwn164euvCJJ0XoeIkXm2WdD4rjjDujX\nL+lopJjFMs4j7eTvVrLb3X3rbC4YByUPqWvmzAlVV3ffHUaki2Qj7uSxUdpmY+AIoLW7X5zNBeOg\n5CF10axZMHAgjBsXGtNFairW5FHFBV90992zuWAclDykrnrqKTjiCLjvvjAiXaQmYhlhnnbyPVjT\nQF4P6ALUz+ZiIpJbvXrBvfeG7rsTJkDPnklHJHVFJtVW5axJHquAhcA/3P3NWCOrAZU8pK574gk4\n5pgwrXv37klHI8Ui79VWhUbJQwSmTg0j0CdNgm7dko5GikGci0FhZi3M7CozeyG6XWFmG2ZzMRGJ\nT9++cNttYUXC555LOhopdZksdjkG+AYYDBwJfAvcGmdQIpKd/v3hllvCOJAXX0w6GillmbR5zHP3\nXde1L0mqthJZ24MPwqmnhilNOndOOhopVLFWWwHLzeynPhxmtg9r5rgSkQJ0+OFw441hBPoLLyQd\njZSiTOa2OhW4I62dYzFwfHwhiUguDBwY5sDq3x8eeQT23DPpiKSUZFJtVd/df0wlD3f/Oi+R1YCq\nrUSqNmkSnHSSemHJz8VdbfWWmV0OtC3ExCEi1Tv44DW9sGbPTjoaKRWZJI/dgLeAm81srpmdYmbN\nY45LRHKof3+480447LAwJ5ZIbdVokKCZlQF3Ay2B+4A/u/vb8YSWOVVbiWRm2jQ49li4//4wtYnU\nbXEPEmxgZoea2UPA1cAVwNbAJGByNhcVkWQccECYhfeII6C8POlopJhl0mD+DlAO3Ozusys89093\nPyO+8DKjkodIzZSXw5FHwtix0KdP0tFIUuJez6OZu3+bVWR5ouQhUnMzZ8KgQXDXXWFtdKl7Yq22\nKvTEISLZ6dkzjEQfOhQefjjpaKTYZDJIUERKVI8e8Nhj8MtfwrJlMGRI0hFJsVDyEKnj9tgjrAfS\nt29IICedlHREUgwy6W21mZndYmZTou1OZqavl0gJ2Wmn0Ij+5z/DNdckHY0Ug0wGCd4GPA5sHm2/\nBfwuroBEJBnbbRfWRL/+erj00qSjkUKXSfLYyN3vBX4EcPeVhOVoRaTEtG8PM2bAPffABReAOjFK\nVTJJHt+ZWevUhpntBWiOK5ES1aZNqMJ6/HE46yxYvTrpiKQQZTLOYw/gn8COwOvAxsAR7j4v/vAy\no3EeIrn39dehF9b228NNN4Xp3aW0xD1IsDGhyuoXgAFvAvXc/ftsLhgHJQ+ReCxdGiZTbNUqTKzY\nqFHSEUkuxT0l+2x3X+nur7n7q+6+AtDEziJ1wPrrh3VAVq4MU7p/913SEUmhqDJ5mFmbqMqqqZnt\nbmZ7RPdlQNPaXNTMWpnZNDP7j5k9bmYtqjhuoZm9YmYvmdmztbmmiGSncWMYPx7atQvzYH3xRdIR\nSSGostrKzIYTlpvtAjyf9tS3wG3u/kDWFzX7O/CFu//dzM4HWrr7/1Ry3LvAHu7+1TrOp2orkZi5\nw4gR8NBDMHVq6JklxS3uNo9B7j4hq8iqPucCoLe7f2pmmwHl7r5DJce9C3Rx9y/XcT4lD5E8ufJK\nuPrqkEA6dkw6GqmNWJKHmQ1z9zvN7Fwg/SAD3N2vzOaC0bkXu3vL6LEBX6W2Kxz3DqFb8I/AaHe/\nqYrzKXmI5NEdd8Af/hAmVNS66MWrNsmjurmtUu0aG2RzYjObBmxWyVMXpm+4u5tZVb/8Pdz9YzPb\nGJhmZgvcfWZlB44cOfKnx2VlZZSVlWUTtohk4LjjoHVrGDAgTOnet2/SEUkmysvLKc/RKmDVVluZ\nWX3grNqUMqo47wKgzN0/MbM2wPTKqq0qvOYS4Dt3v6KS51TyEEnA00/DwIFwxRVhancpLrF11XX3\nH4E4JmmeSGiMJ7p/qOIBZtbUzJpFj9cHDgRejSEWEclSjx7w5JNw4YVw2WWazqQuyaTB/CqgIXAv\nsDS1391fzPqiZq2A8UB7YCFwpLsvMbPNgZvc/ZdmtjWQ6tHVALjb3S+r4nwqeYgk6KOPwmj0bt3g\nuuuggRZ7KApx97YqZ+0GcwDcfd9sLhgHJQ+R5H37LQweHBLHuHGwQVatpZJPcfW2+h3wNPBSNJNu\nwVLyECkMK1fCqafCvHnwyCOwWWVdZqRgxNXmsQVwNfCZmc0ws0vNbEBU5SQi8jMNG8LNN8Mhh0D3\n7rBgQdIRSVwyqbZajzDKfG+ge3S/xN0LZniQSh4ihee22+B//gfuuw969kw6GqlM3BMjNgGaAxtG\nt4+AZ7K5mIjUHcOHh5l4Bw0KgwqltFTX5nET0Ikwl9WzwBzgGXdfnL/wMqOSh0jhmj8/zMg7eHBY\n3rZeJn+ySl7EVfJoD6wHfAIsim5LsrmIiNRdnTrB3LkwZ04YUKhp3UvDukaY1yOsIJhq79gZ+JJQ\nAvljXiLMgEoeIoVvxQo47TR44QWYOFGz8haCWMd5RBdoR0gePYABQGt33zCbC8ZByUOkOLjDVVfB\nP/4BEybA3nsnHVHdFtc4j7NY07tqFWH1wKej+9eiqUsKgpKHSHF55BE44YQwtfuxxyYdTd0VV/K4\nCpgFzHH3j2oRX+yUPESKz6uvhvEgRxwR5sXSlCb5F3u1VaFT8hApTl9+CUOGwI8/hilNNt446Yjq\nlrjHeYiIxKJ1a3jsMejaFfbcE55/ft2vkcKg5CEiiapfP1RbXXEFHHQQ3Hpr0hFJJlRtJSIF4403\n4PDDYd/IqP8uAAAPT0lEQVR94ZproFGjpCMqbaq2EpGS0LEjPPssfPwx9O4NixYlHZFURclDRApK\n8+bwwANhffQuXeDxx5OOSCqjaisRKVjl5WFt9OOPhz/9Sd15c01ddZU8RErWZ5+FBPL99zB2LLRt\nm3REpUNtHiJSsjbZBKZMgX79QjXWlClJRySgkoeIFJGZM+GYY9ZM7964cdIRFTeVPESkTujZE15+\nGd57LwwsfO21pCOqu5Q8RKSotG4N998PZ5+9ZjzI6tVJR1X3qNpKRIrW22+HxvQWLcLI9DZtko6o\nuKjaSkTqpG23De0g3brBbrvBPfeENUMkfip5iEhJeP55GD4cttsObrwRNtss6YgKn0oeIlLndekS\nlrjt1Al23TWMCdHflPFRyUNESk56KeSGG9QWUhWVPERE0qRKITvuCLvsEqqx1CMrt1TyEJGS9vrr\ncMopsGoVjB4dqrQkUMlDRKQKO+4ITz0Fv/oVHHAAnHceLF2adFTFT8lDREpevXohebz2GnzySWhU\nHz9eDeq1oWorEalzZsyAs84Ka4dccw107px0RMlQtZWISA307h0a1IcODeumn3xymPpdMqfkISJ1\nUv36IWksWAAbbBDaRi6/HJYvTzqy4qDkISJ1WosWcOWVYZqTZ56B7beHm28OvbOkamrzEBFJM3cu\nXHABLFoEf/0rDBoEllWrQOHTMrRKHiKSQ+4wbVpIImZwySUwYEDpJRElDyUPEYnB6tXw4IPwl7+E\nhDJiRCiJ1K+fdGS5oeSh5CEiMXKHRx8NSWTJkpBEhgyBhg2Tjqx2lDyUPEQkD9zhySdDEnnnHTj9\n9DD4sGXLpCPLjsZ5iIjkgRn06QPTp8OECfDKK7D11nDaafDGG0lHl19KHiIiWejSBe68E+bPh002\nCeup9+0b2khWrEg6uvip2kpEJAd++AHuvRduuSUMPBw6FE48MQw+LFRq81DyEJEC8tZbcOutcPvt\n0K5dWJhq0CDYeOOkI1ubkoeSh4gUoFWr4PHH4Y47YMqUUNU1eDAcfnio6kqakoeSh4gUuGXLQgIZ\nPz7c77EHHHoo9OsXlstNYgBi0SUPMxsMjAR2APZ09xerOK4fcDVQH7jZ3UdVcZySh4gUjeXL4bHH\nYPLkkEgaNgxJpF8/2G8/aNYsP3EUY1fdV4HDgaeqOsDM6gPXAf2ATsAQM+uYn/Bqr7y8POkQfkYx\nZa4Q41JMmSmGmJo0gYEDwwSMH3wAkybBttvCdddBmzaheuvMM2HcOHj//cJctCqR5OHuC9z9P+s4\nrCvwtrsvdPeVwDjg0Pijy41i+AIXgkKMCQozLsWUmWKLyQx22gnOPTfMp/XFF3DttdC+faji2nPP\n0Oh+2GFw8cWhR9f8+bByZf7ir0yDZC9frbbAB2nbHwLdEopFRCQvGjeG7t3DDUKp45134MUXwzK6\n48fDH/8YSiTbbhsGKW65ZUg2qfu2bcOo96ZN42tLiS15mNk0YLNKnhrh7pMyOEUBFtRERPLLDLbZ\nJtwGD16zf/nyMJ5k4UJ4772QTObMCfcffRTm4FqxIqxX0rJluG/aFNZbL1SXtWtXy7iSbGg2s+nA\nuZU1mJvZXsBId+8XbV8ArK6s0dzMlGhERLKQbYN5IVRbVRX488B2ZtYB+Ag4ChhS2YHZvnkREclO\nIg3mZna4mX0A7AU8amaPRfs3N7NHAdx9FXA6MBWYD9zr7nVs6jERkcJUEoMERUQkv4p6Vl0z62dm\nC8zsLTM7P6EY2pnZdDN73cxeM7Mzo/2tzGyamf3HzB43sxYJxFbfzF4ys0kFFFMLM7vfzN4ws/lm\n1i3puMzsgujf71Uzu8fM1st3TGY2xsw+NbNX0/ZVGUMU81vR9//APMZ0efRvN8/MHjCzDZOOKe25\nc81stZm1ymdM1cVlZmdEn9drZjYqbX9S/35dzezZ6HfhOTPbM+uY3L0ob4RR528DHYCGwMtAxwTi\n2AzYLXq8AfAm0BH4O/CHaP/5wN8SiO0c4G5gYrRdCDHdDpwYPW4AbJhkXNH35x1gvWj7XuD4fMcE\n9AQ6A6+m7as0BsKg2Zej732H6P9BvTzFdEDqWsDfCiGmaH87YArwLtAqnzFV81ntC0wDGkbbGyf9\nWQHlQN/o8UHA9GxjKuaSR0EMInT3T9z95ejxd8AbhDEqhxB+KInuD8tnXGa2BdAfuJk1nRKSjmlD\noKe7j4HQruXuXycc1zfASqCpmTUAmhI6aOQ1JnefCSyusLuqGA4Fxrr7SndfSPiP3jUfMbn7NHdf\nHW3OBbZIOqbIlcAfKuzLS0zVxHUacFn0+4S7f57PuKqI6WPCH2wALYBF2cZUzMmjskGEbROKBYCo\nZ1hnwn+qTd390+ipT4FN8xzOVcB5wOq0fUnHtBXwuZndamYvmtlNZrZ+knG5+1fAFcD7hKSxxN2n\nJRlTmqpi2JzwfU9J6rt/IjA5epxYTGZ2KPChu79S4amkP6ftgF5m9oyZlZtZlwKI63+AK8zsfeBy\n4IJsYyrm5FFQLf1mtgEwATjL3b9Nf85DuTBv8ZrZAOAzd3+JKrpC5zumSANgd+AGd98dWEr4MicW\nl5ltA5xNKKpvDmxgZkOTjKkyGcSQ1/jM7EJghbvfU81hscdkZk2BEcAl6bureUk+P6cGQEt334vw\nh9z4ao7NV1y3AGe6e3vgd8CYao6tNqZiTh6LCPWcKe1YO3PmjZk1JCSOO939oWj3p2a2WfR8G+Cz\nPIbUHTjEzN4FxgL7mdmdCccE4d/nQ3d/Ltq+n5BMPkkwri7AbHf/0kP38AeAvROOKaWqf6+K3/0t\nWFP9EDszG06oEj02bXdSMW1DSPzzou/7FsALZrZpgjGlfEj4PhF951eb2UYJx9XV3R+MHt/Pmqqp\nGsdUzMnjp0GEZtaIMIhwYr6DMDMjZPP57n512lMTCQ2vRPcPVXxtXNx9hLu3c/etgKOBJ919WJIx\nRXF9AnxgZttHu/YHXgcmJRjXAmAvM2sS/VvuTxhXlGRMKVX9e00EjjazRma2FaF65Nl8BGRhmYTz\ngEPd/fsKseY9Jnd/1d03dfetou/7h8DuUXVfYp9T5CFgP4DoO9/I3b9IOK63zax39Hg/IDVBbc1j\niqPnQb5uhN4CbxIady5IKIZ9CO0KLwMvRbd+QCvgiegf53GgRULx9WZNb6vEYwJ2BZ4D5hH+Ktsw\n6bgIDa2vE5YKuJ3Q4ySvMRFKiB8BKwhteSdUFwOhquZtQvLrm6eYTgTeAt5L+67fkFBMP6Q+pwrP\nv0PU2ypfMVUVV/Q9ujP6Xr0AlCX873cCoaQ9N/q9mgN0zjYmDRIUEZEaK+ZqKxERSYiSh4iI1JiS\nh4iI1JiSh4iI1JiSh4iI1JiSh4iI1JiSh+SEmV1lZmelbU81s5vStq8ws99V8rpTzGxY9Hh4NJK6\noFiYdv8dM2sZbbeMttvn4NxP1+DYcjPbI8Njyyyair+25zezhenTnIuAkofkzizCtCiYWT2gNWGa\n55S9gbV+KM2svruPdvc7o13HE+aXSlQ0u+5P3P0D4EbCFORE96Pd/f3aXsvde9TkcOKdA6mq8zvV\nzxdVpYqfpZQOJQ/JlTmEBAGwI/Aa8K2FxZ/WI6xx8lL01+1VZvYccJaZXRIt4jOIMPr17mjG3cZm\ntkd0/PNmNiU1z1OKmTWLSgANou3m0XZ9M9stms00tWhRi+iYn/66NrONovmQUqWeiWb2b8IaDBVd\nRZjG5GxCkvxHZR+CmT0Yxfuamf062relhQWdWptZPTObaWb7R899F923MbOnLCzS86qZ7ZPJhx5N\nz/OUmb0Q3fZOe7q5mT1iYXGfG6PpVzCzA81sdnT8eAszG6/LGdHxr5jZL6LzrG9hwaG50b/ZIZV9\nlmbWtLLjpMjFNVxft7p3I0wN0Q44GTgF+F/CFDI9gBnRMdOB69JecwlwTtpzu0ePGwKzgdbR9lHA\nLZVccwxhniWi614ePX6FsHYIwJ+Aqyq5xkbAu9Hj4YQpHKqchgToS5iKpk81x7SM7psQpqVIbZ9E\nmFX1PODGtOO/je7PBUZEjw3YoJJz/xR72r4mrFnIajvguehxGbCcMGlgPcL0JoOi9zwDaBIddz5w\ncVXnj/a/C/w2enwacFP0+FLg2OhxC8JUQU0rfpZVHZf091W32t1UpJRcmk34q7w7YXGettHjr1m7\nyureas6Rqh75BaEE80T0B3N9wjw9Fd1MmJvqYcKP1q8sLDq1oYfFcCDMV3VfBvE/7u5Lqnn+oCiG\nnYF/V3HMWWaWWrRpC2B7YK6732JmRxKS6q6VvO5ZYIyFGZofcvd5GcQL0Ai4zsx2BX4kJJCfzulh\nYR/MbCxhHrbvCdWJs6PPtRHh321dHojuXwQGRo8PBA42s99H2+sB7QnVXNPSPsvKjmtHSCJSpJQ8\nJJeeJpQydib81f0B8HtC8khfN2BpNedI1bkb8Lq7d6/ugu4+O6q6KQPqu/t8S1tXO+1cKatYU13b\nuMJxy6q6jpntRphxd29glpmNI/z/STVK/4vwY9gH2Mvdvzez6YQfytS6E1tE768ZFT4Dd59pZj2B\nAcBtZnalr2kLqs7vgI/dfZiZ1Sckh59Om/4WWNN2Mc3dj8ng3Ol+iO5/ZO3fjYHu/lb6gWbWjZ//\nG//sOCluavOQXJpN+PH70oPFhGqKvan+r9vUj/u3QPPo8ZvAxma2F4Q1U8ysU2UvBu4grNWeWt72\na2BxWrvBMMLazQALCW0rAEdk8qaitoIbCQt9fUBYge0f7v6hu3eObqOj2BdHiWMHYK+004wizLB6\nCXBThUsQ9dz63N1vJpSmOlcVToXt5sAn0ePjCCW0lK5RYq0HHAnMBJ4BelhYBCvVbpFeWqmJqcCZ\nae8hFXPFGKs6ToqYkofk0muEXlbPpO17hbC061fVvC71F/JtwL/M7EXCd/MIYJSZpaa737vyl3MP\n0JIwBXXK8cDlZjYP2IXQ/gKhofu06Bqt065dXU+mXwML3T1VVXUD0DEqKaSbAjQws/nAZYROBFhY\nP2EPYJSHlfdWmFlqnY7UNfcFXo7iOhK4popYHjWzD6LbvVEsx0ef0S+A79LO+xxwHWF9knfc/UEP\n60kMB8ZGn83s6HXV8QqPU9t/BhpGjeivEdqWKh5T3XFSxDQluxQ9MzsCONjdj1/nwSKSE2rzkKJm\nZv8k9ILqn3QsInWJSh4iIlJjavMQEZEaU/IQEZEaU/IQEZEaU/IQEZEaU/IQEZEaU/IQEZEa+3+/\nDzbe+alQgQAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "thetas=np.arange(0,180,1) # makes an array of thetas from 0 to 180 at 1 degree increments. \n", "Es=np.cos(np.radians(thetas)) # replace this with YOUR equation - this is just and EXAMPLE. \n", "plt.plot(thetas,Es) # make a nice plot\n", "plt.title(\"Write your title here\")\n", "plt.xlabel(\"Write your X-axis Label here\")\n", "plt.ylabel(\"Write your Y-axis Label here\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Problem 1b" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "What are you trying to do. " ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "thermal energy is: 1\n" ] } ], "source": [ "# figure out thermal energy here and print it out\n", "YOUR_NUMBER=1 # obviously this is not the real number. what is it? \n", "print 'thermal energy is: ',YOUR_NUMBER" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "How does this compare with 1a? " ] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "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 }