{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "\n", "\n", "\n", "# 4. Graphing" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Graphing your data shows relationships much more clearly and quickly than presenting the same information in a table. A picture is worth a thousand words, so a well-made graph can save a lot of writing when you report your results." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 4.1 Rough Graphs" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Rough graphs are made for your own benefit. They don't need to be as polished as graphs that will appear in a report, so they can hand-drawn. Graphing each point as you take it is not a good idea, because that is inefficient. It can prejudice you about the value of the next data point. So take a set data points and then graph them all. \n", "\n", "It is very useful to make rough graphs during lab when you have the chance to do something about any problems. You'll know if you need to check any data point, since any strange measurement will stand out in a graph. You'll also be able to spot regions where you should take more data. Typically, people take approximately evenly spaced data points over the entire range of the controllable (or “independent”) variable, which is certainly a good way to start. A graph of that inital data will tell you if there are regions where you should look more closely: regions where your graph is changing rapidly, going through a minimum or maximum, or changing curvature, for example. The graph helps you identify interesting sections where you should get more data, and saves you from taking lots of data in regions where little is happening. \n", "\n", "For example, Figure 4.1 below shows the initial data taken for a phenomenon called mechanical resonance. All you need to know about resonance is that the amplitude (a measure of the response of an oscillating system) depends on the frequency of oscillating force driving the system. Notice that the driving frequencies were approximately evenly spaced. For the apparatus used, the highest and lowest frequencies attainable were easy to find, and the experimenter chose to space the frequencies evenly to get roughly 10 different frequencies between the extremes. From the graph of the initial data, you can see that the response doesn't change very much at either very high or very low frequencies, but something interesting happens near 5 cycles/s. \n", " \n", "
\n", " \n", "Figure 4.1: Initial set of data (filled circles) for amplitude response versus driving frequency for a resonant system. Notice that the data points are evenly spaced.\n", "
\n", "\n", "Suppose that the experimenter noticed this, too, and went back to take more data in the interesting range of frequencies. The frequency spacing used the second time is smaller than used in the first set by about a factor of ten, yielding 15 more measurements in the interest region. The result of adding the second set of measurements is shown in Figure 4.2. As you can see, the shape of the graph is now much better defined. Also, the additional data show that the high amplitude at 4.5 cycles/s was not due to a mistake. The experimenter could, of course, have taken data with the closer spacing over the entire frequency range, but that would waste time on measurements at both low and high frequencies where nothing much is happening. The strategy of taking coarsely spaced data and then backing up to take more data in interesting regions is a good compromise between completeness and efficiency, but you usually can’t identify the interesting regions without making a graph. \n", " \n", "
\n", " \n", "Figure 4.2: Amplitude vs. driving frequency in a resonant system, after adding more data points (triangles) around 4 to 6 cycles per second.\n", "
" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 4.2 Polished Graphs" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "For reports, you should produce more polished graphs using a computer. \n", "\n", "In particular, data points plotted on any graph should include uncertainty bars (sometimes called error bars) showing the uncertainty range associated with each data point. You should show both vertical and horizontal uncertainty bars, if they are both known. If the uncertainties are too small to be visible on the graph, you should mention that in your report so it's clear that you haven't forgotten them. Draw uncertainty bars by indicating the “best guess” value (typically either a single measured value or mean of a set of measurements) with a dot, and then drawing \"bars\" through the dot, whose lengths span the 95% confidence range of that value. An example of uncertainty bars is shown in Figure 4.3. The single data point plotted corresponds to a measured pendulum period $T$ of 1.93 s $\\pm$ 0.03 s for an initial release angle $\\theta$ of 20° $\\pm$ 2°. \n", "\n", "
\n", " \n", "Figure 4.3: A graph of pendulum period T versus initial angle θ, showing how uncertainty bars indicate the uncertainty ranges associated with the displayed measurement value.\n", "
\n", "\n", "There are several other details that are important for polished graphs appearing in reports. You should make sure that each graph is labeled well. Use the following checklist to make a graph has the correct features and format. \n", "\n", "   ◻ The graph has a figure number which you should refer to in the text of the report. \n", "   ◻ The graph has an appropriate title which give more information than axis labels. \n", "   ◻ The axes are clearly labeled, including units. \n", "   ◻ The measured data points are clearly plotted, including uncertainty bars, but no lines connecting data points. \n", "   ◻ If appropriate, a curve fit is included. Also, the fitting formula used, the values of the parameters, and the uncertainties of the parameters are given. \n", "   ◻ If there is more than one data set on a graph, use different symbols and include a legend. " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 4.3 Making Graphs with Python" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "It is difficult to make graphs that follow all of the guidelines above by using Excel. In addition, Excel isn't very useful for fitting data with uncertainties. The matplotlib plotting library for Python (part of pylab) makes publication quality figures. Python can also be used for fitting data (see Ch. 5). On the computational tutorial webpage, there is extensive documentation on making graphs with Python. An example of making a scatter plot of data with uncertainty bars is shown below." ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAy4AAAIqCAYAAADCVR4JAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAWJQAAFiUBSVIk8AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3Xm4JGV5uP/7Yd/FAQUXkEE2NzSggKAgIERjXMKiRkV+Jm6JCqL5uqIS44LRBEGNRkUxEgXELcYVhx0JIIi4z0wcRAyLOOz7zDy/P94qTk9Pd58+feqcU+ec+3NdfXV3be/bVdVV9VS9S2QmkiRJktRma810BiRJkiRpPAYukiRJklrPwEWSJElS6xm4SJIkSWo9AxdJkiRJrWfgIkmSJKn1DFwkSZIktZ6BiyRJkqTWM3CRJEmS1HoGLpIkSZJaz8BFkiRJUusZuEiSJElqPQMXSZIkSa03qwOXiLg2IjIinjbTedHsEBHPrPaZpTOdF81eEfG+aj/67EznpUkR8ayIOCcibql+X0bEy2Y6X6OIiB2q/K/oMW7k7TfOck+txh07ar7bxmOmenG/0ExpReASEad0nCQHvd4403mdjSJinYi4smM9NnZSncpld6Vzao/94f6I+FNELImIr0fE2yNiu6lIv0+eFkTEcRHx7ulKU1MrItaqtulxEbHZTOdnOkXE/sC3gWcAGwM3AjcAd89gtqQ5KyKeGBF/FxGfi4ifRcSK6tx26hDzPiYi3hkRP4iI6yLivoi4NSIuiYh3RMSDpuM3SNNtnZnOQJf7geUDxt/Z9X0pcAdw15TlaG54I/DEWbjsXu4Dbq4+B7AZsADYAXgB8L6I+Arwusz8U4/57wR+A/yugbwsAN4DrATe28DyNPPWomxTgM8Ct/WZ7o+U/ei66cjUNDmK8vu/DPxNZt4zw/mZSlO1/f6vWu5NDS9Xc9N/Ao+b6EwRsR9wbsegBG4FHgTsUb1eGxHPysxfNpBPqTXaFrj8KDOfMezEE5l2voqIbYHjgN8CmwIPmQ3LHuCCzHxmVz4eDOwFvAI4DHgRsHdE7JmZq12YZObFwC7TkE/NYZl5InDiTOejYfUF1H/M8aBlyrZfZr4FeEvTy9WcdR/wE+DHwGXAC4FnDpyjWJdyo/crwKnA+Zl5Z0RsBBwCfBTYBvhWRDxurv+fNb+0oqiYptTHKcU+3kA5SM6WZQ8tM2/OzO9m5guB5wL3Ug7aX5mpPEmz0IbV+x0zmgtp/nhKZu6Wma/OzM9QimYO4zfAzpn50urcdydAZt6VmadSbt4BbA8c2ny2pZkzqwOX8SrnR8TjI+IrEfHHiLg7In4VEe+JiPX7Vc4cpsJZRLyymuaHXcPX6ah/8ciIeFxE/EeVz/sj4swey3peRPxXRFxflVG9ofp+0KjrpWPZh1Iu5L+Rmd+Z7PKma9mTkZnfBt5afd0nIp7dOX7Q9q3qN7wiSuXkP1Xb7I8R8YuIODkiDu6Y9kJgSfV17R71b47tmHaTiHhRRHypWtYt1f64JCI+FRE79PotPfan7ap8/CEi7o2I30bEhyNi00HrpNoPP12ld1dE3BwRV0XEiRGxW595HhoRx1flru+IiDurz++rnnANraovkhHxP+NMd0Q13XURsXbH8M2q/+0VEXF79dv/LyIui4h/jojHTiQ/A9I/lXIXs/b7rm362Y5p+x0/urfZY6rtfl21Dq+IiJd2TB8R8dqIuLxaz3+qpn/kOHldGBEfj4jF1Ta9LSJ+HBFviXLXddjf/EB+gTrNCzp+Q/cxbu0ox7/zq/3onmo//PeI2L5PGqv95yLiORHxvYi4MSJWRcTrh8zr0P/PIZY1sHJ+RGxQ7XO/qX7jddV2Gfi0NvpUzo+uCv0RsWtEnBHleH9PlHPTOyNivQHL3igi3ltt8848PbZ7+RMVEZtHxL9GxNXV/+v31TZ9+JDzHxYR34+ImzrmPzUintRj2oUxVkdx4x7jf12NvyUi1rhGqfab1c773duz2k8urf5Tt0bEoog4cGJr5YFlj3T8HkZmrhxxvt9n5rIB4xcB11Zfdx8lDWhkv9g1Ij4fEcuqffaWiLgwIl4dEX1L/ETEFlHOT7/rSPfTUY6pg87hD1wTRsQ2EfHJKu17I+LHPabfNyJOr+a7tzqunBURL+qetmu+tSLi5RHxw2qfvy/Kefm0iHjKMOtm1LxHxCMi4nUR8d2IWBpjx/8rohyzetZt6l5vEfHs6n9xc/X6QUTs2TH95hHxwWo/vycirqm+b9Bj2d3nvfr4dn2MXXu/IwYc3yYkM2f8BZxCKaN57gTnu7aa72k9xv05cE81vi7/eW/1+ULgQ9Xnz3bN98xq+NIB6b6ymuaHXcPX6UjvZZS6N3Xa9wBndky7HqUseXblsfP7ByaxTjet1s+dwLZd6+vYSW6vkZbdMc1nR0jz1F7rvM+061PKsCfwxWG3L3B61/q/hfIkqf5+Yce03+xII4Hru15v7Jj2jV3Lva1jX0zgdmD/Hvnp3J+eR6n/Ve8n93eM+x9gnT7r4hhKHZx62rpOWP19jfUJ7EupR1RPc2/XPFcDO05g2+3YMe/2A6b7TjXNRzuGPRj4dcf8K6v10Pmb3jeZ/bkjrY9X265e7o1d2/RfOqZ9X699uWubHV5t23pfWtUx7mjKjaMzOtbxHR3jlwEP7pPPw1n92HZX1/50JfCQIX/z2h2/r16nf+oYdkbHtBsDP+xI577qd3Xm4y97pPHAf45yUyGrdXEzsAJ4/ZB5Hfr/WU2/QzV8RY9l9dx+1bjNgEu79v/62Hw78JIBy62PU8f2ywvwbEqDB/Vv6NyXz+zz2zenFCnql6eX9svTEOv1EcD/dm3Hel+8nrHzXa9j5todv7n+fZ3HjpXAq3rMd001/qCu4Vt1beM/6xr/2I48rtdrewKfrz7fTznWdubl+SOsn5GO3yMeg+p1eWoDy7qyWtaJI84/8n7Rsd5Wdq2nFR3ffwhs2GO+bSl1UHulewPwqgH7Y32N8SrKcSwp1yl3Aj/umC6Aj/TYrp3H6FOBtXqk8SDg7I7pVrH6tdtK4O9GWN/D5v0bXfm+uWs9LwYe3mP5ncfho6p8r+zK+12UYvcPBX7B2DVD53H2Gz2W3Xne++sqz1kte7VjNLDxpPftJv5sDfzBTql+1LkjbuindQ1/KGMXeT8CHlcNX5cSUNzB2MF1qgKX26udu047gEd3TPuxarr/pTzW3aQavgnwWsYOuIePuE5PquZ/W4/1NdnAZaRlM02BSzV9fUF49TDbFziAsRPvG4BNO7bbwyn1Zz7UNU/fC6Ou6V5GKU+/F/CgathawGMYC16vo+sg3rU/3Qyc1bE/bVDth/UJ9NU90n1xx/ynAbt0/aYjgA93zbM9YweyT1GCjrWqeR4PfL8adxWw9gS234+r+d7ZZ/wWjAVje3QMfy9jJ6xnUwVolP/yTsDbgL+dzP48YJ0/csB0wwQut1BOMttV4x4EfIax48P7qnX9EsqNjKAEjXXwtMaNi2ofur96fQDYphq+NvBUSjn5BL49wm/veyOoGv/Zavzd1b63fjV8F+B8xk5yj+6ar/7P3UU5UZ4EPLQatyHwiCHy1uj/s9/2q8Z9nrGLhpcD61bDnwRczliwNkrgsqqa/0vAo6pxmwDv7NhvDu6x3P/s2G9e0vE/eAJwCWPns1ECl/oi7AbgL6ku1oD9KPUX69/b60Lx2I7f9XbGzmOPBL7K2IXcPn3W0/u6hh/O2EVk0nEDqBr/d9Xws/tsz5ur/ezVwEbVuO2BC6rxv2cCx61q/pGO36O8aChwodQ3rS8Y1wgcp2G/OLRjf30r1Y0UynHuzykX1wl8omu+6NhW/0c55kc17mnVfDcPSPfajnSvBPbs/A92fH5zx297DbB5NXxDyvXYddX4/9cjjW9V434CPKve7pSbbO+s1vsKYK8Jru9h8/4h4B2U4+4G1bB1gf0ZO89+s8fy6+NwHYi8t2N/Xki5AZrAxZQbs7+k+t9W2+3VjAWeB3ctu/u89z+MXausB/wtYzdr/m3S/5PJLqCJF2OBy32seee6fn1+wIbuDlzez9jB5EE95qvvmPW68GgqcFlc71Q95t+FcqD/E9XJa0AerxxhfT6ZcrL4JdVJt2t9jRy4TGbZTG/g8i7GTqhrdQzvF7i8oxr+rQnkaajAZZxlBHBOtZyXDtiffkrHHcaOaT5Zjf9B1/D1KQf+pFS2HjY/p1Xz/Euf8esDP6umecEElvumap5f9Bn/2j7b5QfV8DePuo4nuD2aDFx+RddFEuWC57cd07ykx7JfUR9Deoy7uBr3hj752oKxk+6TJvjbBz3BfjRjdyPXCBQpT2Pq3/W5rnHP7Pi9Q++LXcto9P85YPtt3/E7X9Zjvi0pLYb1W+54gUsC3+mT3+9W4z/dNbzzieULe8z3YMoF2ISPRZSLnXrZT+8xfmfGbo50/zc3ZeyJ4j/1mHdtyo3DZM1Ao75rfkHX8I9Xw+tz+Ne7xtfHp/f02Z4JvKhHXh7J2IX83qPsg33WX9/j94jLaypwqW+M3gpsMcL8k9kv1qEEiAk8t8/yd6AEmPdR3cSohh/E2Hl7zx7zbc/YBfCgwOVP9HnqTGkN9M4q/7v3mWafajk30VGagRKoJOWpxeZ95q2PVWs8mRhnnY+b9yGWsWU1/yqqm1od4zqPw5/ps27r8fcCC3tM8wV6H6M6z3vX9Vo3jF03r2CIm1WDXm2r47Iu5VFxr9dEytUfUr1/KjNv7TH+y5RH1VPpY9m/JY8jKQe8r2Vmv2Z5z6TcVX1iRAzdWleUugH/TrlAel1m3j/OLEOb7LIz85GZGZn5yqbyNEBnk8nD7Dt1s7db9SpXPVWy/KO/XX3dZ8Ck/5KZvRpA+Eb1/viu4QcBD6McJIZq5SgiNmGsIucJffJ7L+VOap3GsE6jHEwfGxG79hj/19X7l7qG19vlYRNIqy0+nF1l2DNzFeVCB0pxiC/3mG9R9b5DZ3niiNiZctf3TuDTvRLM0gT496uvk64n1+EQyn/pD5QnEt3p3kkpegFw2ID/0IdHTH+6/p+HUn7n7ylPOVaTmTdRjoGTcXyf4f3+y/X57OrMPKNHnm6eRJ4Oq94vzMwLeiz7N5RzUS/Pojwtupexbd8570pKQAGwf0Rs2TH6/Op9j4jYsGP4fpSLm5Mo23zfiIiu8QDn9cnTbzPz9B55uZbytAzWXL8jm8Dxe9pExF8Ar6u+Hpu9uwUYz2T2iwMpgeLSzPxWrwkycymlOOa6jG1TGNvXz8vMS3rM91uGa3TnlMz8Y59xhwMbUVpiu7zXBJl5EeX4vAXwZx2jjqzeP5+Zt/RZfn0OO3DEY9WgvA9UHZ8uphzDnjpg0g/2mPe3lJtPAKdl7zpU9blp0H/o3/qsm89Rgpq1gb8aMP+42tYc8nk5ySaOo1RM3bn6emGvaTIzo1Sufslk0hrHxQPG7V29vyginjtgurqC8jaU+hTDOArYDfhSZp4z3sQTNJXLnko5xDRnUQLFpwBnR8RnKHcJG+nnIUrT0a+nHNQfTblb2X1QG1Th8bI+w/9QvXcHZ3tV71dk5vVDZvMpjN05+fHq1wurqS80thlyuWTm/0XEeZQ7eS+hFDUDIEpF9KdXX7sDl+9QLiaPqQL4LwMXZebtw6Y9g37WZ/iN1fsvqgufbnXLQkEpXlbfAKmPG+sDvxuwfTap3ofePkOoG3E4vwq+ejm7et+Uckd1cdf4O4Cfj5j+lP4/O3T+zn7HjfMod1VHtUYl4Uq//3J94dTzfFZZ4+JySPXv7RcI1ON6nSvrea/oc4MQSl8jqyjHut0oT1DJzN9ExPXA1pRj1TkRsQWlSe5fZOYNEXERpajQ44GfRcRO1fT3UYqi9NJv3UL/9TuuBo7f0yIiHk8JuINSXOhjIy5qMvtFfZzaptrG/dSVyDuPU8Pu60cMGA/DXX/tPU7+FlTv2zB2/q3n/YeIeMM4ediEUjdtUN+EvQzKOwARsReliNvelLpIazRyQf/98c4qSOnlRsqTl37H6frcNOg/dG6vgZm5qrruPpyx/WskbQtcmrCA8qeFwZ2L/d8U52NQoFHfPd60eo1nqFaCImIbSrnF2yhlOBszlcueIvUfKyllLgeqTqSvo9zp2696ERHLgO9RHo1eOUpGIuIA4L8YO7gkZT3WF6QbUioE9zr41PpdqNfLWLdr+FbV+0SeLNb7ZXTMP8jQrVdVvkQJXF4cEW/vuDB8cZXmTzLz113zfJ5ycP5bSn2DlwOrIuKnlHX6qQkEZtOt3/Fn5ZDjYfXtWm+fdZia7TNI/dT3DwOmubbj80NYM3C5aUAwMNBU/j+71L9z0Plh0DoYz8rM7Ndhcr//cv2kYirOZ5P5vePuE5l5V0TcTLlz3V1y4ALKRcx+lKeQ+1KOA+dW48+jBC77UW4C1HfmLxlQmmHQDY1+63egho7fUy4iHk0JDDenXPj3vDEbEd8E9uwx6j8zsz63T2a/qI9T6zPx41RT+/ow118bMdwxsnOarav3zYeYr553ooHLwJvUEfE2Sv3G+jq3bhCjLpHxIEod2H7746Dz5bDnpkH/oUHHx3rcpPr8a1tRsSb0vQ05zQY1c1iv99dVRafGew26+9DpeEqU/0HgjijNOD7wYmzdrFcNm8iBdiqXPRWeUL3/bsAd4tVkaUd/IaUlrv+ilBVdSKkQekVEvHXA7D1FxPrAFykHkR9QKhhumJmbZ+bWmbk1Y0W5mtx3R1lWvV/+ccj9cpiO0jqdSTm4PoqxO1fQv5gYWbwS2BX4J8rFzH2UO3PvAZZUFxbzQb19Lhty+0xFkcz1JzHvSE2/1qbi/zmi6T7HzPQ5bbz0R90n6rv5+3W9nzfO+LqY2ZSbweP3hFRPhBZRLsovp7Tu1y9A3oLexfF7NqM7KNk+w+vj1FeGPE69r2PeptbhMNdfHx4yf6f2mPc5Q857LRPXN+8R8UTGgpYTKa3sbZCZCzr2x7rI6UwfN3ppJE9zMXCpm5GDweXi+42r28Ffo63qDhP9g3erH7c10v9Eh0dV7x+k3HnqftWPDt9Vfe/3eH+6l92o6mSzf/V1QkUoMvP6zPxoZj6fcldgL8pFUgDvj4n3GbIPZd3cRGmK86Is9UQ6DXNXaqLquyqPGjjV6ur98iETqVc1rKrc63err38NUBX/2I3ynz1twLw/y8x3V0VJN6c0D/0LSjD9H9HR78scVm+fnWfg99Z3AQftT519z4xURns8U/D/7Fbne1Cxn+mub1XnaZTz2bDLHuX3jrtPVEW366ff3ftEHZjsVR2zuwOXH1Pqc/ULbKbDTB2/hxYRD6MELY+iHBP/fEDRPTLzaUPc6JjMfjGZ65up3Ndrk8lfXcy36Wu3YdV18L6dmW/MzF/lmn0Bzej+yHD7zKTOD3MucKnuMvym+tqvY8roN46xYkVbD7g4GKmDoQ51GcbnxYBOmDSy1zL2yHmNCrbDqu72X0I5WNSVyjr3m/pJzqC7CPXF3K8HFG+Y6JOLYdRlwP8sIrYeOOWYSxm72zOpynMD1E9VXljt+/XTlvOHvTuVmfdmqfRZdxL2CEq58yZ0Pp1r2x2r+rixGaWs/XS6onrfq1cHZJX6yddtlFZ3ptQQ/89R1L/z6V2Vwjvt12f4VPlJ9T6o8vfTB4wbpP69+w6Ypt/vrefdZcAx5hmMXWdc0TXuF5QbjRtQmsjdFfhlVhWTM3MFZZ9/SET8JeVYuoLSUtl0manj91CqG0yLKHXKlgLPzNEq43ebzH5RH6ceVzUoMhH1vj7ofzzqvl6r87d/TLAz5Y55Dxk41dSp98ef9BoZpTPqPaYvOz313C+qhgrq7dp9LJiQORe4VL5evb86evcq/kL63yX6NaUS6NqUtstXU/0RXzDJ/J1CucO8DeO0+DSRP9aAOymRmXWLQADvqoYNHTRN5bKbVLWo8qHq6wWZ+f1B03fM17dH1+oEWj+J6ywWUbd0tFaf/QzGnjzt3CuNKr+TPRD38gPKxdw6wD8PM0N1l65+zPzuQU9dqp5yN+k3foBvUZ7IPYRy8d23mFiVzqCedu/u+DyZIkwPyFKs8M7q67DlmKdFZv6cscrH/1zdze4pSi/rzfRSXNR9cjyU0qxld3obA/9QT5tDFs8c1oj/z1HUv3NbSt2r7nwsoPRnMJ3q89nCiDi0e2SUnrJHzVPdQtPTImKNwCgidmSshalu36M0uLAeY9u+c961Kf28AJyTpcWjB2RmMvZE/FjK9ci5XYupn668u3r/cZYW7KbLTB2/x1VdG5xF6U/md8AB2Vx9v8nsFz9g7Hrgo4OeDve4vqn39f2iRw/0EbEdpV7UZJxOaYp5I8auFYbN3ynV+1MjYmDjTiMERcOo98cn9Bn/Lma4rhXw99UxqduRlKcxKxnbziOZq4HLiZQnJw8HvhsRjwGIiHUj4qWUjtR6Vtiu7qrUTfidFBF7R8RaEbF2RDyL8qfsV3Z0KNUFSN3ax/sj4qSIWFiPr+qIHBwRp9K7udQpUV2MZvU6dvw5Jrz8a6tlf7bh5W4eEc+KiNOB/6ZqdYkSoA7rnyPijIh4fnVxUi9764j4BCXIXEXp7Rd4oOnB+tHxK/os90LKBfZDgC/UdyYjYsOIeBWlo8wm7pCtJkvTyf+v+npERHy5vvsVxcMi4jUR8dGuWd9Cqej3CODian08cDEYETtGxJsoAf6TRsjX3YwFRx+gtAB4P/2b1jwnIk6MiKfH6k0DP57SvCKUSuG/7Bg32f34F9X7y1tYBO31lDo+TwTOj4gD6jxWx6nHR8S7KR3bPrSpRLO0QnNy9fXDEfHK+kKu2q++Q6lvcielD46mTfj/OYrqd36h+vqZiHhZ/VQ8SjPe36dcqE+bzFzMWDHKz0XEizvy9HhKADGoaPOgZZ9LqRgfwFcj4i/qJ00R8XRK0c6eTxqytOxXN+18TES8rb6ZEaWlwDMoTbKuolxQ9VLXV6kvUruLgZ03zvipNqXH7+oGw5b1i7HAe/3O4d03iarv36EcB/5ACVp+P2o+uk1yv7iP0gIplCazvxcRe3TMv25EPDki/hlY0jX7WYw15/uNiPjzjvn2puzr3UX1Jvrb/shYQP2q6tz4uHp8RGxQnW8+SVd9qsz8b0rnjFD2h/dEx9PGiFgQES+IiG8x5A3DCTqren9+RLw1qqbEI+KhEfEvlHN+49cTE7Qx5br7sVBuOkXEK4BPVOM/nZmTaeCEkTuAafLF2BOIcyc436AO0/6CsQ6S6pal6u/nUXaqpKvn1mreHRirK5OUk3Hd6dHlwBurz4M6oOzbeV3HtJ/umL5uqeRmxjpAS+CsBtfzwE4iu/I/oU4qx1t21zST6YDyXsY6Jb2BEkR2rsOVlLv3C/osp18HlB/vWs6tjPXeXL/e0mN57+8YfztwdfV6fcc0b+pazi2M9RI/qf2JcTrBpAQiq7ry2LnO1ujQk9LizHUd09xPKeN9T9fv2KdXmkNsy2d1Lee/Bkz7865tu7wrH3cAz2hqP67mf1XH/HdRguCrgeM7phmmA8p+26xvr+3DLIPyJLhz37y32j73sfp6nVAnXww4nlbjN6YUS+lM95aO73dTKgUP9Z+bYN4m/P8c9N8YtA0oRfEu7VjuPVV69f/nJQOWO14HlH07iBy0niitZf60T55u68jTnSOs20dQAt3Ofb7uWPL6jv9Dr3ytTSmO23msWN7xfSXw6gFp78bq23CrrvHrM3buTeDZo/ynBm2bIdbPSMfvIZf9vq5l93t1H2f+pmPc7fTvuPt64IwR/3Mj7xfV/K9k9WPS3ZTrqhUdw+7vMd92jHVgmZTrrzrd6zrSXaMzY8Y5hnVN+x5WPzfeWe27KzuGLekx3yaUenWd2+dmxv6P9WuNTh7Hyc9QeacETnUaq6o817/j3/vt5wzXufqF1TRrdL47aBmsfs76a8auMTqvuxO4CNh4lP2x8zVXn7iQmd+h3KX5GuXPsj7lT/guSsdsdTGLNZ68ZOkcaS/KXa6bKAfn31NaNXoa5WJpsvlbkZmvppQh/U/KxdF6lKYVr6nyfQRjHQKqWI+xVlC2oJzA/5fyZ34HsH1mviQzJ9oE4UeAoykHpMWUOz7rU7bFaZSDSa87KO8B3k5prnNtShHER9FRzCgz/5XyePtiysF7HcoTi3dRyq1Pen/qp8rzbpSbA7+jNGN4N6UflY/So2nrLPUGdqb8rour/G1ezfdjSgMNu2fppGsUZzH2pAr6FBOrvAI4jlKE5BrK/2MVpVf6k4DHZ7k72JgsrVe9htJ2/0pKsaFHMVZvakZlueu3I+WJ1U8o/4HNKRewF1H2q11ysne11kz3TuBgSrGk+k70hpSg7jOUbfHfTabZYdT/54Rl5m2UctrHUe4IB+VEfBrlnHJpE+lMME/LKS3xvZ9Sl6HO05eqPNVNT4/b9HuPZf8BeDLleHAN5Th2C2Wb7gYsGzDvysx8KeXp9g8p++AmlCZrvwQ8OTN7dpZauZKx4i+/zswbOkdmqQhf19dbSdm/p9VMHr8H6Lx224T+HXdvxVh/JBMymf2imv+zwC6UEjC/pAQsm1Guqc6mBITb95jvakqrkR+nXHetQwkMPl2lW3cuPeF9vSudf6zS+Sxj/6mNKcHRdyl1ZffuMd8dmfk8SgMxX6fs6xtRzq1LKNdzhzL21Klph1GudepqDVCOx0dk5mumKM2JuIBy/XwmY0HLrylPuQ7IBop6RhUtzTsRcTFl5R6Rqzd3J0nSrBERrwE+BSzKiTdTLs0aEfFB4G3AyTk1Tb5rgqqiq3UQtU2O1gz00ObsE5dBqjKae1Hu4Jw9zuSSJLVSVQetvrt71qBppdksIragFJMD9/V5a84GLhHx2qqy4PYdlVc3iYj/j7HKVV/OzFF7HJYkacpFxMKI+FxEPC2qzn2rxhj2pFzAPZZShOZzg5YjtV2UBpFOjIjd6gZZqgZXnkkpMvxQSvHwbwxYjOawOVtULCKOB+qelFdQyt5uzuptyh80Ql0ISZKmTUTsQqnXVbuF0pJY3ZrY3cChmfnd7nml2SRK6631flxXfN+Esdb8bgKelZmXz0D21MN0FxWby50ffomyo+9HaR1jC0plXImVAAAgAElEQVTw8gtKG+X/nv07lJIkqS1+T+kr5WBKwxkPoWr1iFIp/l+rRmWk2e4KxhpR2p6yr99HaYDiu5R9van+ajQLzdknLpIkSZLmjjlbx0WSJEnS3GHgIkmSJKn1DFwkSZIktZ6BiyRJkqTWM3CRJEmS1HpzuTnkaRERy4DNgKtnOCuSJEma27YDbsvMhTOdkZlg4DJ5m2244YYLHvOYxyyY6YxIkiRp7vrVr37F3XffPdPZmDEGLpN39WMe85gFl19uJ66SJEmaOrvvvjtXXHHF1TOdj5liHRdJkiRJrWfgIkmSJKn1DFwkSZIktZ6BiyRJkqTWM3CRJEmS1HoGLpIkSZJaz8BFkiRJUusZuEiSJElqPQMXSZIkSa1n4CJJkiSp9QxcJEmSJLWegYskSZKk1jNwkSRJktR6Bi6SJEmSWs/ARZIkSVLrGbhIkiRJaj0DF0mSJEmtZ+AiSZIkqfXWmekMSJIkSbPVCWctfuDzMQftNIM5mfsMXCRJkqQRnbhoyQOfDVymlkXFJEmSJLWegYskSZKk1jNwkSRJktR6Bi6SJEmSWs/ARZIkSVLrGbhIkiRJaj0DF0mSJEmtZ+AiSZIkqfUMXCRJkiS1noGLJEmSpNabVYFLRGwREa+MiK9HxNKIuDsibo2ICyPibyNira7pt4uIHPA6baZ+iyRJkma3xTfcPvC7mrXOTGdggg4HPglcB5wDXANsBRwCfBZ4dkQcnpnZNd9PgW/0WN7PpzCvkiRJmoMuWnoTJy5awqXLlq82/OATzmePhQs4+sAd2WeHLWcod3PXbAtcFgPPA76dmavqgRHxDuBS4FBKEPPVrvmuzMzjpiuTkiRJmptOv+wa3v61n7Gq+zZ55dJlyzni5Es4/pBdeeFTtpnezM1xs6qoWGaenZnf6gxaquHXA5+qvj5j2jMmSZKkOe+ipTcNDFpqqxLe9rWruGjpTdOTsXlitj1xGeT+6n1Fj3EPj4jXAFsAfwIuzsyrpi1nkiRJmvVOXLRk3KCltirhpEVLLDLWoDkRuETEOsDLq6/f6zHJQdWrc55zgSMz85oh07i8z6hdhsymJEmSZqnFN9y+Rp2W8VyybDmLb7idnbbadIpyNb/MqqJiAxwPPB74TmZ+v2P4XcA/AbsDD65e+1Eq9j8DWBQRG09vViVJkjTbjFrsy+JizZn1T1wi4ijgzcCvgSM6x2XmjcC7u2Y5PyIOBi4E9gReCZw4XjqZuXuf9C8Hdpt4ziVJkjRb3HFPr9oIUzef1jSrn7hExOsoQccvgf0zc6jnd5m5gtJ8MsC+U5Q9SZIkzRGbbDDa/f5R59OaZm3gEhFvBD5O6Ytl/6plsYn4Y/VuUTFJkiQNNGoleyvnN2dWBi4R8VbgBOBKStBy4wiL2at6/21jGZMkSdKctNNWm7LHwgUTmmfPhQusmN+gWRe4RMS7KJXxLwcOzMy+NZ4iYs+IWK/H8AOAY6qvp05JRiVJkjSnHH3gjqwVw027VsBRB+44tRmaZ2ZVobuIOBJ4L7ASuAA4KmKNvefqzDyl+vwh4HFV08fXVsN2BQ6oPr8rM380lXmWJEnS3LDPDlvywUOeMG4nlGsFHH/IrhYTa9isClyAhdX72sAb+0xzHnBK9fmLwF8BTwGeDawL3ACcAXw8My+YspxKkiRpznnRU7blkQ/eiJMWLeGSHv267LlwAUcduKNByxSYVYFLZh4HHDeB6U8GTp6q/EiSJGn+2WeHLdlnhy1ZfMPtHHzC+Q8M/8Ex+1qnZQrNujoukiRJUht0BykGLVPLwEWSJElS6xm4SJIkSWo9AxdJkiRJrWfgIkmSJKn1DFwkSZIktZ6BiyRJkqTWM3CRJEmS1HoGLpIkSZJaz8BFkiRJUusZuEiSJElqvXVmOgOSJEnSbHX0gTvOdBbmDQMXSZIkaUTHHLTTTGdh3rComCRJkqTWM3CRJEmS1HoGLpIkSZJaz8BFkiRJUusZuEiSJElqPQMXSZIkSa1n4CJJkiSp9QxcJEmSJLWegYskSZKk1jNwkSRJktR6Bi6SJEmSWs/ARZIkSVLrGbhIkiRJaj0DF0mSJEmtZ+AiSZIkqfUMXCRJkiS1noGLJEmSpNYzcJEkSZLUegYukiRJklrPwEWSJElS6xm4SJIkSWo9AxdJkiRJrWfgIkmSJKn1DFwkSZIktZ6BiyRJkqTWM3CRJEmS1HoGLpIkSZJaz8BFkiRJUusZuEiSJElqPQMXSZIkSa1n4CJJkiSp9QxcJEmSJLWegYskSZKk1jNwkSRJktR6Bi6SJEmSWs/ARZIkSVLrGbhIkiRJaj0DF0mSJEmtZ+AiSZIkqfUMXCRJkiS1noGLJEmSpNYzcJEkSZLUegYukiRJklrPwEWSJElS6xm4SJIkSWo9AxdJkiRJrWfgIkmSJKn1DFwkSZIktZ6BiyRJkqTWM3CRJEmS1HoGLpIkSZJab1YFLhGxRUS8MiK+HhFLI+LuiLg1Ii6MiL+NiJ6/JyL2jojvRMTyiLgrIq6KiDdGxNrT/RskSZIkTdw6M52BCToc+CRwHXAOcA2wFXAI8Fng2RFxeGZmPUNEPB/4KnAPcDqwHHgucAKwT7VMSZIkSS022wKXxcDzgG9n5qp6YES8A7gUOJQSxHy1Gr4Z8BlgJfCMzPxxNfxdwNnAYRHx4sw8bVp/hSRJkqQJmVVFxTLz7Mz8VmfQUg2/HvhU9fUZHaMOAx4CnFYHLdX09wDHVl//bupyLEmSJKkJsypwGcf91fuKjmEHVO/f6zH9+cBdwN4Rsf5UZkySJEnS5My2omI9RcQ6wMurr51Bys7V++LueTJzRUQsAx4HbA/8apw0Lu8zapeJ5VaSJEnSRM2VJy7HA48HvpOZ3+8Y/qDq/dY+89XDN5+qjEmSJEmavFn/xCUijgLeDPwaOGKis1fvOXAqIDN375P+5cBuE0xXkiRJ0gTM6icuEfE64ETgl8D+mbm8a5L6icqD6G2zrukkSZIktdCsDVwi4o3Ax4GfU4KW63tM9pvqface868DLKRU5v/tVOVTkiRJ0uTNysAlIt5K6UDySkrQcmOfSc+u3p/VY9y+wEbAjzLz3uZzKUmSJKkpsy5wqTqPPB64HDgwM28aMPmZwE3AiyPiyR3L2AB4X/X1k1OVV0mSJEnNmFWV8yPiSOC9wErgAuCoiOie7OrMPAUgM2+LiFdRAphzI+I0YDnwPEpTyWcCp09P7iVJkiSNalYFLpQ6KQBrA2/sM815wCn1l8z8RkTsB7wTOBTYAFgKvAk4KTPHbVFMkiRJ0syaVYFLZh4HHDfCfBcBf9F0fiRJkiRNj1lXx0WSJEnS/GPgIkmSJKn1DFwkSZIktZ6BiyRJkqTWM3CRJEmS1HoGLpIkSZJaz8BFkiRJUusZuEiSJElqPQMXSZIkSa1n4CJJkiSp9QxcJEmSJLWegYskSZKk1jNwkSRJktR6Bi6SJEmSWs/ARZIkSVLrGbhIkiRJaj0DF0mSJEmtZ+AiSZIkqfUMXCRJkiS1noGLJEmSpNYzcJEkSZLUegYukiRJklrPwEWSJElS6xm4SJIkSWo9AxdJkiRJrWfgIkmSJKn1DFwkSZIktZ6BiyRJkqTWM3CRJEmS1HoGLpIkSZJaz8BFkiRJUuut09SCImIn4JnAvsC2wJbA3cCNwJXAOcDZmXlPU2lKkiRJmh8mHbhExIuBvwf2qQf1mOxA4E3ALRFxCvCxzLx6smlLkiTNpBPOWvzA52MO2mkGcyLNfSMHLhGxP/CvwBOBW4AvABcBlwHXA8uBDYEtgF2AvYCDgWOAv4+Ik4D3Z+Ztk/kBkiRJM+XERUse+GzgIk2tyTxxWQRcAbwI+GZm3tdjmtur19XA94DjImJH4LXA64E7gH+aRB4kSZIkzQOTCVwOy8yvTXSmzFwCvDkiPgxsN4n0JUmSJM0TIwcuowQtXfNfTylSJkmSJEkD2RyyJEmSpNZrrDnkThERwNbAur3GZ+Y1U5GuJEmSpLmp0cAlIg4H3gY8AVi7z2TZdLqSJEmS5rYmO6B8HXASsAK4EPhD9VmSJGnOWXzD7Wt832mrTWcoN9Lc1+STj2OAG4G9M3NZg8uVJElqjYuW3sSJi5Zw6bLlqw0/+ITz2WPhAo4+cEf22WHLGcqdNHc1WTn/EcBXDFokSdJcdfpl13DEyZesEbTULl22nCNOvoQzLvv9NOdMmvuaDFx+D6zf4PIkSZJa46KlN/H2r/2MVTl4ulUJb/vaVVy09KbpyZg0TzQZuJwCPDsiLNwpSZLmnBMXLRk3aKmtSjhp0ZKpzZA0zzQZuHwIuAz4YUTsZwAjSZLmisU33N63eFg/lyxbvkYFfkmjayxwycyVwCeAHYCzgVsiYmWPly2NSZKkWWXUYl8WF5Oa02RzyM8HzqT037IM+D9sDlmSJM0Bd9wz2iXNqPNJWlOTzSEfB9wFPCczL2xwuZIkSTNqkw1Gu2QadT5Ja2qyjsvOwJcNWiRJ0lwzar8s9uciNafJwOUm4L4GlydJktQKO221KXssXDChefZcuICdtrKtIqkpTQYuXwUOioh1G1ymJElSKxx94I6sFcNNu1bAUQfuOLUZkuaZJgOXY4Gbga9ExHYNLleSJGnG7bPDlnzwkCeMG7ysFXD8IbtaTExqWJM1xn4GrAvsCTw3Im4Bbu0xXWbmoxtMV5IkaVq86Cnb8sgHb8RJi5ZwSY9+XfZcuICjDtzRoEWaAk0GLmtRmj++pmNYr3sSQz5klSRJap99dtiSfXbYksU33M7BJ5z/wPAfHLOvdVqkKdRY4JKZ2zW1LEmSpLbrDlIMWqSp1WQdF0mSJEmaEgYukiRJklpv5MAlIg6dTMIR8bCIeOpkliFJkiRpfpjME5evRMTlEfGiiFh/2JkiYueIOAFYCjxzEulLkiRJmicmUzn/QOAE4MvArRHxTeAi4MfAdZQ+XTYAtgB2AfYC/hx4MnAfcBLw0UmkL0mSJGmeGDlwycxzIuLPgL8GXge8HDhiwCwB3AKcCJyYmb8bNW1JkqQ2OPrAHWc6C9K8ManmkDMzgS8BX4qInSlFv54GbEt50nI3cCNwFXAucHZm3j2ZNCVJktrimIN2muksSPNGk/24/Ab4DfCJppYpSZIkSWBzyJIkSZJmgVkXuETEYRHxsYi4ICJui4iMiFP7TLtdNb7f67Tpzr8kSZKkiWusqNg0OhZ4InAHcC2lxbLx/BT4Ro/hP28wX5IkSZKmyGwMXI6hBCxLgf2Ac4aY58rMPG4qMyVJkiRp6sy6wCUzHwhUImImsyJJkiRpmsy6wGVED4+I11CaaP4TcHFmXjWRBUTE5X1GDVNUTZIkSdIkzJfA5aDq9YCIOBc4MjOvmZEcSZIkSRraXA9c7gL+iVIx/7fVsF2B44D9gUUR8aTMvHO8BWXm7r2GV09idmskt5IkSZJ6arw55Ih4bkScFhE/jYilHcMfExFviYhHNJ1mP5l5Y2a+OzOvyMxbqtf5wMHAJcAOwCunKz+SJEmSRtNY4BLFFyhPNw4HHg0s7JjkZuADwMuaSnNUmbkC+Gz1dd+ZzIskSZKk8TX5xOXvgSOAzwMLgI90jszM64GLgOc0mOZk/LF633hGcyFJkiRpXE0GLn9L6ejxVZl5K5A9plnC6k9hZtJe1ftvB04lSZIkacY1GbjsDJyTmb0CltqNwEMaTHOgiNgzItbrMfwASkeWAKdOV34kSZIkjabJVsVWABuMM80jgDsmk0hEvAB4QfV16+r9qRFxSvX5psz8h+rzh4DHVU0fX1sN2xU4oPr8rsz80WTyI0mSJGnqNRm4/BJ4RkREr6cuEbEBJWD4ySTTeRJwZNew7asXwO+AOnD5IvBXwFOAZwPrAjcAZwAfz8wLJpkXSZIkSdOgyaJiX6T0In9CRKy23IhYG/hX4OHAKZNJJDOPy8wY8NquY9qTM/MvM3O7zNwkM9fPzG0z80UGLZIkSdLs0eQTl38HngccRWkO+XaAiDiTUhH+4cA3M/M/G0xTkiRJ0jzQ2BOXzFwJ/CXwXmA9YCcggEOAjSg92B/eVHqSJEmS5o8mn7jUHTseFxH/SAlctgBuBX5dBTaSJEmSNGGNBi61qnL+b6Zi2ZIkSZLmn8aKikXEoyPi5RGxRZ/xW1bjt+81XpIkSZL6abJVsbcB/wLc1mf8rcBHgP/XYJqSJEmS5oEmA5dnAD/MzPt7jayGn8VY54+SJEmSNJQmA5dHAFePM801lGaRJUmSJGloTQYu9wGbjTPNpkA2mKYkSZKkeaDJwOXnwHMiYt1eIyNiPUo/L79sME1JkiRJ80CTgcupwLbAGRGxdeeI6vsZwDbAfzSYpiRJkqR5oMl+XD4NHAo8HzgoIq4C/kCp+7IrsBHwQ+BTDaYpSZIkaR5o7IlLZq4C/gI4Hrgf2IsSyOxFqf/yAeA51XSSJEmSNLQmn7jUTR6/IyKOBXYBNgduAX5twCJJkiRpVI0GLrUqSLESviRJkqRGNFk5X5IkSZKmRKNPXCJiR+BoYA/gwcDaPSbLzHx0k+lKkiRJmtsaC1wi4qmUVsM2BFYAN1Tva0zaVJqSJEmS5ocmn7h8EFgfeC3wuczsFbRIkiRJ0oQ1Gbg8BTgzMz/d4DIlSZIkqdHK+fcB1zS4PEmSJEkCmg1cfgT8WYPLkyRJkiSg2cDlHcDeEXFEg8uUJEmSpEbruDwfOBs4JSJeCVwO3NJjuszMf2owXUmSJElzXJOBy3Edn59evXpJwMBFkiRJ0tCaDFz2b3BZkiRJkvSAxgKXzDyvqWVJkiRJUqcmK+dLkiRJ0pQwcJEkSZLUeo0GLhHxsIj4REQsjYi7I2Jlj9eKJtOUJEmSNPc1VsclIh4BXApsBfwCWB/4HXAvsH2V1pXArU2lKUmSJGl+aPKJy7uBrYFnZeYTq2Gfz8xdKIHL94ENgUMaTFOSJEnSPNBk4PLnwPcy84fdIzLzWuBwSuDyjw2mKUmSJGkeaDJw2ZpSRKy2khKoAJCZdwBnAc9vME1JkiRJ80CTgcttwHod328GHtE1za3AQxpMU5IkSdI80GTg8jtgm47vPwUOiIiNACJiLeBg4NoG05QkSZI0DzQZuCwC9o+IdavvXwAeDvwoIj4MXAQ8Dji9wTQlSZIkzQONNYcMnEwpHrYlcF1mnhoRuwNvAHatpjkNeH+DaUqSJEmaBxoLXDJzCfChrmHHRMQHKM0hX52ZNzSVniRJkqT5o8knLj1l5h+BP051OpIkSZLmrsbquETEyoh41zjTvDMiVjSVpiRJkqT5ocnK+VG9hplOkiRJkobWZOAyjAcD90xzmpIkSZJmuUnVcYmIfbsGbddjGMDawLbAS4HfTCZNSZIkSfPPZCvnnwtk9TmBI6tXLwGsAt48yTQlSZIkzTOTDVzeSwlYAng3JZA5r8d0K4E/Aedk5q8nmaYkSZKkeWZSgUtmHld/jogjgW9k5kmTzZQkSZIkdWqyA8qFTS1LkiRJkjo1FrhExNrA+pl5V9fwA4DnA3cBn87MZU2lKUmSJGl+aLI55I8AyyPiQfWAiHgxcBbwBuCtwKURsU2DaUqSJEmaB5oMXPalVL6/tWPYe4BbgJcDbwE2B97UYJqSJEmS5oEmA5dtgKX1l4jYHtgZ+FhmnpqZHwG+CzyrwTQlSZIkzQNNBi6bAbd1fN+H0lTy9zqG/QJ4ZINpSpIkSZoHmgxcrgM6WxZ7JnA3cHnHsE2AFQ2mKUmSJGkeaKxVMeB/gOdFxF8C9wCHAYsy8/6OabYH/tBgmpIkSZLmgSafuHygWt43ge8D6wHvr0dGxGbAM4BLGkxTkiRJ0jzQZAeUP4uIPYEjq0GnZ+ZlHZPsCvwA+HJTaUqSJEmaH5osKkZm/gz4hz7jLgQubDI9SZIkSfNDk0XFJEmSJGlKjPzEJSJeXn38embe3vF9XJn5H6OmK0mSJGn+mUxRsVMo/bT8D3B7x/dBoprGwEWSJEnS0CYTuPwNJQi5rvr+islnR5IkSZLWNHLgkpmndH3/wqRzI0mSJEk9WDlfkiRJUuvNusAlIg6LiI9FxAURcVtEZEScOs48e0fEdyJieUTcFRFXRcQbI2Lt6cq3JEmSpNFNplWx3444a2bmo0dNFzgWeCJwB3AtsMugiSPi+cBXgXuA04HlwHOBE4B9gMMnkRdJkiRJ02AylfPXYs1WxNYDHlZ9XgH8CdiiI53rgPsmkSbAMZSAZSmwH3BOvwkjYjPgM8BK4BmZ+eNq+LuAs4HDIuLFmXnaJPMkSZIkaQqNXFQsM7fLzIX1i/IU5A+U5pH3BzbIzIcBGwAHAJdQAo5dJ5PhzDwnM5dk5nhNLwMcBjwEOK0OWqpl3EN5cgPwd5PJjyRJkqSp12Qdl/cDm1OebJyXmasAMnNVZp5LCWYWVNNNlwOq9+/1GHc+cBewd0SsP31ZkiRJkjRRkykq1u2vgC9nZs+iYJl5T0R8E3gxcFSD6Q6yc/W+uEd+VkTEMuBxwPbArwYtKCIu7zNqYB0bSZIkSZPX5BOXLYB1x5lm3Wq66fKg6v3WPuPr4ZtPQ14kSZIkjajJJy7/S6ns/p7MXCNQiIgHU+qcjNoa2VSI6n3c+jKZuXvPBZQnMbs1mSlJkiRJq2vyicungIcDl0bEyyNiu4jYsHo/klI5f2vgEw2mOZ46gHpQn/GbdU0nSZIkqYUae+KSmR+PiB2BNwCf7zFJAB/LzH9rKs0h/AZ4MrATsFodlYhYB1hIaba5TU+BJEmSJHVp8okLmXk0pVPHzwE/oQQEPwFOBp5WjZ9OZ1fvz+oxbl9gI+BHmXnv9GVJkiRJ0kQ1WccFgMy8GLi46eWO6EzgQ8CLI+JjHR1QbgC8r5rmkzOVOUmSJEnDaTxwmWoR8QLgBdXXrav3p0bEKdXnmzLzHwAy87aIeBUlgDk3Ik4DlgPPozSVfCZw+nTlXZIkSdJoZl3gAjwJOLJr2PbVC+B3wD/UIzLzGxGxH/BO4FBgA2Ap8CbgpMwct0UxSZIkSTNr1gUumXkccNwE57kI+IupyI8kSZKkqTfrAhdJktrkhLMWP/D5mIN2msGcSNLcZuAiSdIknLhoyQOfDVwkaeo02hyyJEmSJE0FAxdJkiRJrWfgIknSiBbfcPvA75Kk5ljHRZKkCbpo6U2cuGgJly5bvtrwg084nz0WLuDoA3dknx22nKHcSdLc5BMXSZIm4PTLruGIky9ZI2ipXbpsOUecfAlnXPb7ac6ZJM1tBi6SJA3poqU38fav/YxV43RdvCrhbV+7iouW3jQ9GZOkecDARZKkIZ24aMm4QUttVcJJHU0lS5Imx8BFkqQhLL7h9r7Fw/q5ZNlyK+xLUkMMXCRJGsKoxb4sLiZJzTBwkSRpCHfcs2Ja55Mkrc7ARZKkIWyywWg9CIw6nyRpdQYukiQNYdR+WezPRZKaYeAiSdIQdtpqU/ZYuGBC8+y5cAE7bbXpFOVIkuYXAxdJkoZ09IE7slYMN+1aAUcduOPUZkiS5hEDF0mShrTPDlvywUOeMG7wslbA8YfsajExSWqQNQYlSZqAFz1lWx754I04adESLunRr8ueCxdw1IE7GrRIUsMMXCRJmqB9dtiSfXbYksU33M7BJ5z/wPAfHLOvdVokaYpYVEySpBF1BykGLZI0dQxcJEmSJLWegYskSZKk1jNwkSRJktR6Vs6XJGkSjravFkmaFgYukiRNwjEH7TTTWZCkecGiYpIkSZJaz8BFkiRJUusZuEiSJElqPQMXSZIkSa1n4CJJkiSp9QxcJEmSJLWegYskSZKk1jNwkSRJktR6Bi6SJEmSWs/ARZIkSVLrGbhIkiRJaj0DF0mSJEmtZ+AiSZIkqfUMXCRJkiS1noGLJEmSpNYzcJEkSZLUegYukiRJklrPwEWSJElS6xm4SJIkSWo9AxdJkiRJrWfgIkmSJKn1DFwkSZIktZ6BiyRJkqTWM3CRJEmS1HoGLpIkSZJaz8BFkiRJUusZuEiSJElqPQMXSZIkSa1n4CJJkiSp9QxcJEmSJLWegYskSZKk1jNwkSRJktR6Bi6SJEmSWs/ARZIkSVLrGbhIkiRJaj0DF0mSJEmtZ+AiSZIkqfUMXCRJkiS1noGLJEmSpNabF4FLRFwdEdnndf1M50+SJEnSYOvMdAam0a3AR3sMv2O6MyJJkiRpYuZT4HJLZh4305mQJEmSNHHzoqiYJEmSpNltPj1xWT8iXgZsC9wJXAWcn5krZzZbkiRJksYznwKXrYEvdg1bFhGvyMzzxps5Ii7vM2qXSedMkiRJ0kDzJXD5PHAB8AvgdmB74PXAq4HvRsRTM/OnM5g/SbPUCWctfuDzMQftNIM5kSRpbpsXgUtm/mPXoJ8Dr42IO4A3A8cBfzXOMnbvNbx6ErNbA9mUNAuduGjJA58NXCRJmjrzvXL+p6r3fWc0F5IkSZIGmu+By43V+8YzmgtJkiRJA833wOWp1ftvZzQXkiRJkgaa84FLRDwuIhb0GP4o4OPV11OnN1eSJEmSJmI+VM4/HHhbRJwDLKO0KvZo4DnABsB3gI/MXPYkSZIkjWc+BC7nADsDf0YpGrYxcAtwIaVfly9mZs5c9iRJkiSNZ84HLlXnkuN2MClJE7X4htvX+L7TVpvOUG4kSZrb5nzgIklNu2jpTZy4aAmXLlu+2vCDTzifPRYu4OgDd2SfHbacodxJkjQ3zfnK+ZLUpNMvu4YjTr5kjaCldumy5Rxx8iWccdnvpzlnkiTNbQYukjSki5bexNu/9jNWjVMrblXC2752FRctvWl6MiZJ0jxg4CJJQzpx0ZJxg5baqoSTFi2Z2gxJkjSPGLhI0hAW33B73+Jh/VyybPkaFfglSdJoDFwkaQijFvuyuJgkSc0wcJGkIdxxz4ppnU+SJK3OwEWShrDJBqO1Hj/qfJIkaXUGLpI0hFH7ZbE/F0mSmmHgIklD2KkmkCAAABBwSURBVGmrTdlj4YIJzbPnwgXstNWmU5QjSZLmFwMXSRrS0QfuyFrx/7d350F2lXUax78PiBtgWDJCOYgEScTCZRRlmSgIGVLUTOk4jJbllIiOyzDqBFMypWCJoGOBNSoSXMYFRdEqtVTEFXEiCoKCG5uCCQiKo6DIYoAAAr/545yGpu1Oupvue07u/X6qut7u85577q/TN8l9+l3O9M7dLLBi2eL5LUiSpBFicJGkaVq620KOP+TJGw0vmwVOOOQpThOTJGkOuWpUkmbgRc/cmZ22fSSrVq/lgknu67L3ou1YsWyxoUWSpDlmcJGkGVq620KW7raQNdevY/mJ59x3/KyV+7mmRZKkeeJUMUmapYkhxdAiSdL8MbhIkiRJ6j2DiyRJkqTeM7hIkiRJ6j2DiyRJkqTeM7hIkiRJ6j23Q5akB+GIZYu7LkGSpJFgcJGkB2HlQUu6LkGSpJHgVDFJkiRJvWdwkSRJktR7BhdJkiRJvWdwkSRJktR7BhdJkiRJvWdwkSRJktR7BhdJkiRJvWdwkSRJktR7BhdJkiRJvWdwkSRJktR7BhdJkiRJvWdwkSRJktR7BhdJkiRJvWdwkSRJktR7BhdJkiRJvWdwkSRJktR7BhdJkiRJvWdwkSRJktR7BhdJkiRJvWdwkSRJktR7BhdJkiRJvWdwkSRJktR7BhdJkiRJvWdwkSRJktR7BhdJkiRJvfeQrguQNnUnfmvNfZ+vPGhJh5VIkiQNL4OL9CCdtHrtfZ8bXCRJkuaHU8WkB2HN9es2+LUkSZLmhiMu0iycd+UNnLR6LRdefeMDji8/8Rz2WrQdRyxbzNLdFnZUnSRJ0vBxxEWaoc/+8NccesoFfxFaxlx49Y0cesoFfO6H1w64MkmSpOFlcJFm4Lwrb+CoL17KvbXh8+4teNMXL+G8K28YTGGSJElDzuAizcBJq9duNLSMubdg1biF+5IkSZo9g4s0TWuuXzfl9LCpXHD1jS7YlyRJmgMGF2maZjvty+likiRJD57BRZqmW++4e6CPkyRJ0v0MLtI0bfXw2e0ePtvHSZIk6X4GF2maZntfFu/nIkmS9OAZXKRpWrLD1uy1aLsZPWbvRduxZIet56kiSZKk0WFwkWbgiGWL2SzTO3ezwIpli+e3IEmSpBFhcJFmYOluCzn+kCdvNLxsFjjhkKc4TUySJGmOuGpYmqEXPXNndtr2kaxavZYLJrmvy96LtmPFssWGFkmSpDlkcJFmYeluC1m620LWXL+O5Seec9/xs1bu55oWSZKkeeBUMelBmBhSDC2SJEnzY2SCS5KdknwsyW+T3JnkmiTvTbJt17VJkiRJ2rCRmCqW5PHA+cCjgTOAK4C9gCOAg5Msrao/dliiNmFHuHOYJEnSvBuJ4AJ8gCa0rKiqk8cOJnkPsBJ4B3B4R7VpE7fyoCVdlyBJkjT0hn6qWJJdgeXANcD7J3S/FbgNODTJlgMuTZIkSdI0DX1wAQ5s27Oq6t7xHVW1DjgPeCSwz6ALkyRJkjQ9ozBV7Altu2aK/rU0IzJLgNVTXSTJj6foeurll1/OnnvuOfsKJUmSpI24/PLLAXbpuIzOjEJwWdC2t0zRP3Z8m1lef7P169ff85Of/OTiWT5em77d2/aKTqtQV/z5y9fAaPPnr0G+BnYB/jSA5+mlUQguG5O2rQ2dVFWTDqmMjcRM1a/h52tgtPnzl6+B0ebPX74GBmcU1riMjagsmKL/URPOkyRJktQzoxBcftG2U+1ZO3YTjqnWwEiSJEnq2CgEl7PbdnmSB3y/SbYGlgLrgR8MujBJkiRJ0zP0waWqrgLOolnM9NoJ3ccBWwKfrKrbBlyaJEmSpGkalcX5rwHOB1YlWQZcDuwNHEAzRezNHdYmSZIkaSNStcHNtIZGkscCbwMOBrYHfgd8CTiuqm7ssjZJkiRJGzYywUWSJEnSpmvo17hIkiRJ2vQZXCRJkiT1nsFFkiRJUu8ZXCRJkiT1nsFFkiRJUu8ZXCRJkiT1nsFllpLslORjSX6b5M4k1yR5b5Jtu65N8yvJC5KcnOTcJH9KUkk+1XVdGowk2yd5ZZLTk1yZZH2SW5J8L8krkvjv6ghI8s4kq5Nc274Gbkzy0yRvTbJ91/Vp8JIc2v5/UEle2XU9ml/t+76a4uO6rusbVt7HZRaSPB44H3g0cAZwBbAXcADwC2BpVf2xuwo1n5JcBDwVuBX4DbA78OmqekmnhWkgkhwOfJDmJrZnA78GdgAOARYAXwBeWP7jOtSS3AX8BPg58HtgS2Af4BnAb4F9qura7irUILU3ub4U2BzYCnhVVX2026o0n5JcA2wDvHeS7lur6l2DrWg0PKTrAjZRH6AJLSuq6uSxg0neA6wE3gEc3lFtmn8raQLLlcD+NG9eNTrWAM8DvlZV944dTHI0cCHwzzQh5gvdlKcBeVRV3THxYJJ3AEcDRwGvGXhVGrgkAT4O/BH4InBktxVpgG6uqmO7LmKUOKVhhpLsCiwHrgHeP6H7rcBtwKFJthxwaRqQqjq7qtb6G/XRVFXfrqqvjA8t7fHrgP9pv3zOwAvTQE0WWlqfa9vFg6pFnVsBHAi8nOY9gKR5YnCZuQPb9qxJ3risA84DHkkzZUDSaPlz297daRXq0nPb9pJOq9BAJHkicAJwUlWd03U9GriHJXlJkqOTHJHkgCSbd13UMHOq2Mw9oW3XTNG/lmZEZgmweiAVSepckocAL22/PLPLWjQ4SY6kWdOwgGZ9y7NoQssJXdal+df+nT+NZp3b0R2Xo27sSPMaGO/qJC+vqu92UdCwM7jM3IK2vWWK/rHj2wygFkn9cQLwJODrVfXNrovRwBxJsznDmDOBl1XVHzqqR4NzDPA04FlVtb7rYjRwHwfOBX4GrAN2BV4HvBr4RpJ9q+riDusbSk4Vm3tpW9c/SCMiyQrgDTQ7DB7acTkaoKrasapC85vXQ2jevPw0ydO7rUzzKcleNKMs766q73ddjwavqo5r1zxeX1W3V9VlVXU48B7gEcCx3VY4nAwuMzc2orJgiv5HTThP0hBL8lrgJJptcQ+oqhs7LkkdaN+8nE4zVXh74JMdl6R5Mm6K2BrgLR2Xo/4Z26Rlv06rGFIGl5n7RdsumaJ/bCeZqdbASBoSSV4PvA+4jCa0eNOxEVdVv6IJsXskWdh1PZoXW9G8B3gicMf4Gw/S7C4K8JH22GT3+NBw+33burvsPHCNy8yN3bNjeZLNJtzHYWtgKbAe+EEXxUkajCRvpFnXchFwUFXd0HFJ6o/HtO09nVah+XIncMoUfU+nWffyPZpfdDqNbPTs27a/7LSKIWVwmaGquirJWTTTAV4LnDyu+ziahP2hqnIvd2lIJXkL8Dbgx8Byp4eNliS709x47roJxzcD3k5zg+Lzq+qmLurT/GoX4r9ysr4kx9IEl09U1UcHWZcGJ8kewO8m/tuf5HE0o/AAnxp4YSPA4DI7rwHOB1YlWQZcDuwNHEAzRezNHdameZbk+cDz2y93bNt9k5zafn5DVXnn5CGV5DCa0HIPzY4yK5obZz/ANVV16oBL0+AcDPx3knOAq2jumL4DsD/N4vzrgFd1V56kefZC4E1JzgauptlV7PHAPwAPB74OvKu78oaXwWUW2lGXZ9C8eTkY+Hvgd8Aq4Dh/+zr0/gY4bMKxXdsPgF/RbJGq4bSobTcHXj/FOd8FTh1INerC/wIfppka/FSa7e9vo/nF1WnAKv8fkIba2TT39XsazdSwLYGbaaYIngacVlXuLjsP4p+rJEmSpL5zVzFJkiRJvWdwkSRJktR7BhdJkiRJvWdwkSRJktR7BhdJkiRJvWdwkSRJktR7BhdJkiRJvWdwkSRJktR7BhdJkiRJvWdwkSRJktR7BhdJkiRJvWdwkSRJktR7BhdJGnFJnpOkkhzbdS0zkeShSdYm+dosH79n+32/Yq5rkyTNPYOLJA25JLu0b9BP7bqWObYC2A04ZjYPrqofA18C/ivJVnNZmCRp7hlcJEkXAk8E3td1IdOVZEvgzcC32gAyW8cDO9KEIElSjxlcJGnEVdXtVXVFVd3QdS0z8C/ANsCpD+YiVXUhcAXwb0k2n4O6JEnzxOAiSUOsXbdydfvlYe2UsbGPl7XnTLrGJcl32uNbJDkmyVVJ7khyRZJXjTvv8CSXJlmf5DdJjksy6f8vSfZO8vkk1yW5K8m1ST6U5DEz/NZeAdxFM9Vr4nNsneQtSS5L8qck69raP5tkz0mu9RlgZ+DvZliDJGmAHtJ1AZKkefUdmpGJI4CLeeAb/YumeY3PAHsDXwf+DLwA+HCSPwNPAQ4DvgqsBp5Hs+bkduCd4y+S5OXAR4A7gS8D1wKLgVcCz02yT1X9emPFJFkAPAP4YVXdPqEvwJnA3wLfBz4K3A08FngOcC4wcWrZeW17EPDNjT2/JKkbBhdJGmJV9Z0k19AEl4uq6thZXGZn4ElVdTNAknfTTK86EbgZeEpV/V/bdyxwJXBkkndX1d3t8SXAh4BrgP3Hzm/7DgS+BZwE/NM06tkX2Bz40SR9T6IJLV+qqgdcqx0FWjDJY37YtvtN47klSR1xqpgkaWPeNBZaAKrql8D3aEZy3j4+hLTnfQVYCPz1uGv8O7AFcMT489vHfJtmBOa5SbaeRj07t+3vNnDO+okHqureqrppkuO3AHeMu64kqYcccZEkbcxkIxu/bdvJdvQaCyY7Ab9qP9+3bfdP8sxJHvNomlGUJVNcc7zt2/YvQgjwc5opcC9O8jjgDJqQ9aOqumsD17wR2GEjzytJ6pDBRZK0Qe2IxER3t+2G+rYYd2wsbPznRp5uOvdTGRtNefjEjqq6p516dgzNWpyxdTbrknwCOKqqbp3kmo9gklEaSVJ/OFVMkjQIYwFnQVVlAx/fnca1ft+220/WWVU3VdXKqnos9y/+vwJ4HfDBiee3a1+2GXddSVIPGVwkafjd07Zd3qfkB2377Dm41iVtu/vGTqyqK6vqFGB/4FbgHyc57QlAmP4ua5KkDhhcJGn43QQU3S4+fx/NVsontjuMPUCShyaZbqj5GfAHYJ9JrrMoyR6TPGZb4GFMPh1s7DpnT/P5JUkdcI2LJA25qro1yQXAs5N8GlhDMwrz5aq6ZMOPnrMarkjyr8DHgJ8lObOtYwuaQPVsmjAynVGUSnI68Ooke1TVz8Z1PxU4PcmPgctoNhH4K5qRli2YcG+Z1nKaP48zZvv9SZLmn8FFkkbDoTT3XTkYeDHN1KjfcP+0q3lXVZ9KcjHwBuAAmsBwG024+Dzw2Rlc7gPAq4GXAm8cd/xHwPE0U8MOphlp+QPNTmWrquob4y/S3szy+cBXq+raWXxbkqQBSVV1XYMkSTOW5Js0IyyLqmpWO4Il+Q9gFbBfVZ07l/VJkuaWa1wkSZuqI2ludPma2Tw4ySOAo4AvGFokqf+cKiZJ2iRV1aXtupmtZ3mJXYAPA6fOVU2SpPnjVDFJkiRJvedUMUmSJEm9Z3CRJEmS1HsGF0mSJEm9Z3CRJEmS1HsGF0mSJEm9Z3CRJEmS1HsGF0mSJEm9Z3CRJEmS1HsGF0mSJEm9Z3CRJEmS1HsGF0mSJEm9Z3CRJEmS1HsGF0mSJEm99/9ksq2jp58wFQAAAABJRU5ErkJggg==" }, "execution_count": 8, "metadata": { "image/png": { "height": 277, "width": 407 } }, "output_type": "execute_result" } ], "source": [ "from pylab import *\n", "t = array([1,2,3,4,5])\n", "d = array([0.9,4.1,8.7,16.5,24.9])\n", "derr = array([0.6,0.9,0.75,0.9,1.2])\n", "figure()\n", "scatter(t,d)\n", "errorbar(t, d, derr, ls='None')\n", "title('Figure 4.4: Distance vs. time for sliding down a 12-degree ramp')\n", "xlim(0)\n", "ylim(0)\n", "xlabel('time (s)')\n", "ylabel('distance (m)')\n", "show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 4.4 Graphical Analysis" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Analyzing data on a graph often involves finding the slope and the intercept of a line. You will frequently manipulate your data so that the resulting graph is approximately a straight line (more about this in Ch. 6). The slope and/or intercept of such a line can usually be used to extract information about the physical system under investigation. \n", "Determining the slope and intercept of a linear graph is such a common and important task that you will learn to use a program to help you do it (see Ch. 5). However, it is good to be able to estimate roughly the slope and intercept of a rough graph by hand, so that you can see if your measurements are at least roughly correct before you enter them all into the computer. Estimating the slope and the intercept also allows you to check the computer's results, which will help you to catch simple errors like entering the data incorrectly. \n", "\n", "Start by drawing by eye the line that you think best matches your data. This line doesn't need to pass directly through any data points, but should be the best fit to all of the data. Next, find the slope and the intercept of that line you've drawn.\n", "The slope of a line is defined as the “rise over run” or the change in the vertical coordinate value divided by the change in the horizontal coordinate value. To determine the slope, you must first choose two points on your line. They need not be actual data points but they must lie on your line. They should be about as far apart on the graph as possible to minimize the effects of estimating their locations. Mark each of those points by a dot with a circle around it to distinguish them from data points. Read the coordinates $(x_1,y_1)$ and $(x_2,y_2)$ of the points off the graph. Calculate the slope $m$ using\n", "\n", "\\begin{equation}\n", "m = \\frac{\\rm rise}{\\rm run} = \\frac{\\Delta y}{\\Delta x} = \\frac{y_2 - y_1}{x_2 - x_1}. \\tag{4.1}\n", "\\end{equation}\n", "\n", "Using the slope, the $y$-intercept $b$ is\n", "\n", "\\begin{equation}\n", "b = y -mx. \\tag{4.2}\n", "\\end{equation}\n", "\n", "Since the $y$-intercept is defined as the value of $y$ where a line intersects the $y$-axis (where $x = 0$), you can also read the intercept directly off the graph as long as the graph shows the $y$-axis." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 4.5 Exercises" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**4.1** Make a polished graph of the data provided in the table below) Use the checklist in section 4.2 to make sure that you have included everything.\n", "\n", "\n", "| Distance fallen (m) | Kinetic energy (J) |\n", "| --- | --- |\n", "| 0.200 ± 0.003 | 2.00 ± 0.10 |\n", "| 0.400 ± 0.003 | 3.78 ± 0.15 |\n", "| 0.500 ± 0.003 | 5.65 ± 0.20 |\n", "| 0.800 ± 0.003 | 8.01 ± 0.25 |\n", "| 1.000 ± 0.003 | 9.82 ± 0.30 | \n", "\n", "**4.2** Draw what you think is the best possible line through the data points in the graph you just created in the previous exercise, and find the slope and intercept of this line." ] } ], "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.8.5" } }, "nbformat": 4, "nbformat_minor": 1 }