1{
2 "metadata": {
3  "name": "",
4  "signature": "sha256:0462e55ea96e1efb52552045e241d0ed966342993dc50ef1bd4e9fe3b803c5b6"
5 },
6 "nbformat": 3,
7 "nbformat_minor": 0,
8 "worksheets": [
9  {
10   "cells": [
11    {
12     "cell_type": "code",
13     "collapsed": false,
14     "input": [
15      "from lcapy import Circuit, j, omega\n",
16      "cct = Circuit()\n",
17      "cct.add(\"\"\"\n",
18      "Vi 1 0_1; down\n",
19      "R 1 2; right, size=1.5\n",
20      "C 2 0; down\n",
21      "W 0_1 0; right\n",
22      "W 0 0_2; right, size=0.5\n",
23      "P1 2_2 0_2; down\n",
24      "W 2 2_2;right, size=0.5\"\"\")\n"
25     ],
26     "language": "python",
27     "metadata": {},
28     "outputs": [],
29     "prompt_number": 32
30    },
31    {
32     "cell_type": "code",
33     "collapsed": false,
34     "input": [
35      "cct.draw()"
36     ],
37     "language": "python",
38     "metadata": {},
39     "outputs": [
40      {
41       "metadata": {
42        "png": {
43         "height": 100,
44         "width": 200
45        }
46       },
47       "output_type": "display_data",
48       "png": "iVBORw0KGgoAAAANSUhEUgAAAZMAAADHCAQAAAB2fXa/AAAABGdBTUEAALGPC/xhBQAAAAFzUkdC\nAK7OHOkAAAAgY0hSTQAAeiYAAICEAAD6AAAAgOgAAHUwAADqYAAAOpgAABdwnLpRPAAAAAJiS0dE\nAP+Hj8y/AAAACXBIWXMAAADIAAAAyABj+uetAAASe0lEQVR42u2dz4ssWVbHP6e7F7MYaWKEcSU8\n46ELF68d4yG6kEaIQgU3M5CP/geM+gvaLHAxPbusxYDbTHAhsxipBN04MpIJs3DxRCqm25WNkjlj\nMwhPtJJGRRmn67iIqKzInxGZERn3Rtb53MXLjFcZ92RkfOPce8+954piGMZ+3nJtgGH4j8nEMEox\nmRhGKSYTwyjFZGIYpZhMDKMUk4lhlPKOawMyJNLUtQ1Gs0gMwFznri2pjwfeREIZcCuRazuM5pBE\nbomJ6TOTiYSu7amLU28iMX0CFsRA4PpSGE0hfUJ9mb9OGHIrF91uLYjLySoSEuo0N+JCp64vhtEE\nEjLR54X3QxKmeuHarjo4bXTp3KRxhiQgg8L7KRB3u1HtQd/EODtC+nkHHiDrwptMDKNACiw4g/Gt\nRzwZEDbOBx3L85VB4Gycq9NdePMmtZFYklqfH7qs/RSsRUpiYNrtkS7Ug4KixK6tONr6CbMan45R\nEle1t3B1QpQ7Itd21CvmTWoiITGh9I4+QQ842h/Urv30DIFXHfclJpPa9IGjb3QJSYCoMC7UYu2n\nRwZEZxERc+3OlC43ugiWL49qVjDIX9y4qP3kVyfpfnMrK+ZN6tFfvjriiS7B8lO9o+Y91ar91EhC\nn45PUlniWqdKp73JHcqMGxQlPPjTfRTNPcqw7dpPfGUiJgTLd3E3f9+HYt6kBpIQANeMgGOe6Akw\n1yumQE8OnPxZu/ZTXpmIAa90sTxwbN/LF1zrVOmuN2GGcqco3KLcPT49K306QVH6Cj0UZdBm7Se9\nLgG39IgLxetB6wrfyLUBy45oz7UdB9sdo+Q3d37LHxT/YMLy5maGHnYr1a39pL/n7cZBk0mNyhMm\n3C3f3jHx5aeuZP0ERR+e4Yfe6MXb/JgbvV7tJ70ugy0HJ66tqvmdXBvQ1UKIUuh45zdHZZ/IEGXZ\n8SbgDuW2rdqtHPhruzagq2X1Nl/e6BWfmuu3+fJGr9hDq1e7lYN/bdcGdLPk/amV2zK/dSuF0zZF\nkQunUpixbu1WDv69XRvQzbLt2b/pIXZ+euuzf91DnKp2K0f83q4NqGDi+7zmNX/t2o4Vm+7Y0mWu\nGujLw4prHXYilGoyq1W7L4UXfMpP+AkvXFuy1brsrvtO/s61OSXGPuOHKPcoyhcPRrsuu8alVkev\n9nx+x7hUcYj4dLX7Ufg0/1XvueeNX1LhmfwQlft3FZUv+I7/MvnntQPfdm2RUgzsbfxPhUDfY1hx\n438qhRnr1e5H4dO1A29cW1Qs8gZ9T/9JVb+p7yrKh85N2nsxN837mWub9j+1q8Q/9vmM8vhH3dp9\nKLzIPUmxfODaqqV1H6Lv6yNf0rd+6tyovQZ/tuWg84bXamBv439LbvT9TaPyG71e7X6UDV+i3G/3\nj06s++wdLfJN1waVGvzfzk3YXv51d0e7LNC3fzyrLMxYNp5VofZ/dH71PC9v6Xu6Cuo062MZ8hm/\nuHHw7/lj13aRQmH266rNATOCXVkOJWQGjPRy53ce0GdnDkwZkgDPd6WvrlT79/gTx9fvu/w8snJE\n+Sl/4NiqB/7yS1/+n8Lbj/iW3134LX2TL3jm2qpSq/cE+spj7fvCjNvCis3W3tIV8r5v8vWCL3lX\n5QvnRpWY/PnaBXXeM6lg886GUbUpJbubZVVu87q1t3SN3qy8vedz1xYVi/yn3Ged+O9rqCjfcW5S\nyeV8xo/I4ib33PPatT0Vrd4R6NseVtz4q51hxu1hxWZrb+kKveDN8ndVPvcsbvK+vEHf0XcVlXte\nex83URQ+5N9QXvvf3FpavGM0q+o41PYh46rDvXVrb+0qfcAdyv/609xase5DPuPfH+465+ZUMjlG\nXdtwoMVbAn27w4obn94aZtwdVmy29lavkwc9pSrF1sKfhhEQsJpmrgcs8pXre9Ex8/zvl0hMmJ/3\nxLUbW3Ct00rPnM55k80mzmEzrjYbWPvDis3W3uJVMm/yxBnDStrS7FXVp/mYBYVsKRISA6Nd0ZqG\nazfWca3TSs+cLnqTleHXw1eDrA7+Vl+L0kztrV0l8yZPG10w5nEztcwvjA84QSH7Vp4dclp96+ja\ntRuruNZppWdOB71J8Rl+XGDv0YMcEz2vW3tL18i8yVNH54yBREKy/IyHPs0zf9LnITvkQXnda9du\nFDCZnJKHhlN2mx/YgdY0T5raz5Ogtlq7sYJrd1bJNXey0aXkgb6sHBHYy8OMSsWwYrO1t3J9rNFl\n8DgEe1RgLw8zFs/TYu3GIyaTk6Kj/EavHvFY5Xrt33ZrN3JMJqcm6zof+zTPwozH3+b1ajdyTCan\n5poFo+oRj1U0ay4d5Uvq12488I5rA84dXci41mDsiOj427x27QYAXq+FXxoZM1Gpfx7DN0S7sc+v\nNboMoxSTiWGUYjIxjFJMJoZRisnEMEoxmRhGKSYTwyjFZGIYpVgU3hMkYrDlcEqqFkV3jsnEEzTl\nAuSWiGu9yo5JRI8bmfLKZvi6xWTiFxGwnLyhKalAnwkvXRv2tLG+iUdIDLA2x2kKRNn/GK4wmfhE\nTMGXrBC5Nu1pYzLxiQhItxzbJR6jJUwmPrHNmyTASNMjzmY0hsnEG7b1TGRIyJgr17Y9dUwm/rDm\nSySWCTGXasPBzrEBYX+IAGSQv45YcGWhRT8wmfhDDFw/NLok5oZEpuZJfMAaXZ6w3jPRKdfEWyew\nGK3juTeRkBD42o7Q2zmxOcqVAj0uXRtmeCoTiegREawE1SYAAjBlQcr4zLJPbcZMLPLuDZ7JRHrE\nZJtx7iYGegwkZcz0bCIKm95kW7DRcII3MpGYHj2C5YGUOXMWpMwzvyERATEQkTXFIiKQOWOuu97R\n3dqkNG/iD65T4isKYb5PbVZu6JftMkjEgNny7R2Dqnvc+lkYsLafFQHKcu9d//bgbeh7d2TjBvcG\nhPnmaYpyQ3LI7U5If7mHx52ve3hU+iYzNraj5u7hGLHJxLmdDisPGHCXvxlW36d27Szx0q/Mijup\nd6MwKGzWMyt6FBKUiULs576JjXz7jsjEYQ5hiRgu579e1xvqlYR+3vEf6dkMoUpED1jo8RnpPacr\nOYTdPUfi3I/MmnmeFDzTbbf7KU+p0BFv4igKLwkTAmCkz5t5muhCr3hJCkTcii1jMhrEiUxkyBCA\ny2YbSDrngjEQMpHExTczzhMHMpFbEmDBRfObMOtCX3ENBAxNKEZTtC4TGebR5QO6bhIc0ojSKy5Z\nAENLtGA0Q8sykYSETCSHTMOIDpspqyNeAXBjfRSjCVqVicQMgQWXp55colMuyZpeQe2TGU+eFmUi\nETcAvGpjuqKOGAEPdRpGDVqTiQQMCYDLtsJJeskUiMWWNhk1ac+b9ImAUfOjW3t4xRzoWw/FqEdL\nMpGABJi3O5FEF3lX3vyJUYu2vMkga3C1/fU0ZQTENjRs1KEVmUhIAkydTHLLpg32HdRsnA3trF7M\nGj2V57nKYC21dEAok7U/qjinWOdyTZ9YklZ7RcZZ0YJMJKYHjA/wJesDxgHBxrHqkZdrEgL6mEyM\nY2lhsnS2OvHIRVdK/YVJDFCUnuvp2FY2fpmOTKRvo9GV+RKX6YKu6QMxlmq0MTYawcdyLXf1T6IX\np/22J5eJxAQ43p9DFzKmZ6nhGqWpscNfd/1FqrBVJpLQK7xdiXZIv3CBLiv4iOyvXT/Hp/QIJDqb\nrF7u6cLS3ObY2WqMUXRbRg8ilNuqLX1mKLc1W7C1kyYQ7vo2Vpa/asKAPtHyiPXllmVn3CQfl9rm\nLRbMuai2pYBEhLj3Jeg8z8hrbCCBDGTGDSEQkMhEwjwJh5Gzr2+yINiapnRwwET47MZ0LhNgTEQo\n4ZllHq6N9OkzL05IlYiJXDKwAfRH9kXhU2BjtYb0mB8QAYmA1Itbc7q0x8iRQG4YMNKXxV9UU66Z\nbIlUPWX2tFZv0PU+AcFh6X2YbJ7DWetbsd7J+m+p2xMAMkMtjdNj2dfomrP59O0fmNY6xp+s6gsC\nbC3jIzdEjHdM4Umh6+nLm2Rfo2vBWqNLIqJDdgP0bIFtCiVbQjwhZEDMYmckaf7EBnxLKOubsLKk\naXBggO4h9akfzDGZ5EhIHxjt8Ri+tAG8oKzRVfAnMvCkM34sC0wmD2SDvbvnbKfePNy8YI83ySWR\nexMJifSq+P/SkztbPttRsnl2O32Jjq1nUmT/sqw5j95kyNXG/0/psnd5skjP/Ty7blEukxBAeqTr\n86F0rK/smdNJsjaA9T4qs3+GcO5NJKDPiacqG21j00Crs9+bLMgiHxvREolkIOVzfrLP+DIsbDvj\nVkSGri3wjfJGFxITrkZLJKKnV/TKhJI/r3zq5lsjsRSJ7SqtU0Em3Gx03hO9koCogp/YCFE6xKcZ\nAW7Z33lPqif3eCpUkcn1RrRkSj6kWHp+byLfns0IcIpOmYNsXVYgfVIbmNmgZHrczomL3FZZbMUQ\n5c71xDVFoYfSkQQFLV2NLb8fCUPXtvlYytLZXW2JlpDN7qq0HiEFAvFhOVS8tMdAx1wSycq2FhLI\nDeH57IPcJEdueC0DEp7rQm701d6/C7gDrvWq4olP90XvCBjvt/ZpIT0GwJgpEBIRMOrE5tMOOFYm\nE+Z6KX0WZbkUZULMXJ87/poRt8DV+e6wfhwS51524ThFlOccm4BoSk8GUOG2mxITOs9p4s9iY6/Q\n6bFTVo58vp7qe8hpz3/0t5VwOTmy7O9mOG92yYyQVF+6tOG8eFoyOTqdXVUXrXNJiejhUCYS+pHf\n5axoZvLShD/iY9dfpZw2kqNOiQidZoQf5HYYjdFMZ1/g4y4MG7TgOyVgRnBIN14CBiVByWn1zrjE\nTICRDXX6hygXXZBJC95EFzKiTyj9qre2LmReMq/okFGZspV6hlFCSz0xmRGy4Hn70yCkxw3mSzyl\nK96krb0Xr4HASbrNAbBwOYBgdJ+WZKIj5kDS9tp56ROyP4OIYZTS3r7wl0DAsM2ZuhLnvsT6JUYt\nWpOJTrkGIm7aqlEe6rIV+0ZN2vMm6BVTIG5nCakEDLO96LvQRTT8pkWZAK9IgUSSFuq6IQJGts21\nUZ9WZaILLlkAw1MLRYbEwNSGgY0maNeboCnZio9hhbwsRyKBTEiAh7oMoyYtywR0mqfrHshJRr0k\nYpInh6i+p5dh7KV1mYCOeMkCSJg0LRTpMSECplxYujajKRzIBDTlJSkQMZOm9hcHpM8NATDSC/Mk\nRnM4kQnonAtGQMAk2zm2LpLILJ8wf2kdd6NhnKbBGSxfDglrnCdmkr+cPe5rbsX/QkdSQrmuPmSY\nv7xjcMymmERLiRx5BisOf3+TSWUTHn2Bcku/ql8hZsBs+dYk0sFCR2Tiycp/iRkW1ivOGTPdNclE\nQiJieoXcxKMtCVyNDtCV9SaeyATynFG9tcW9UxaFlYoRwVp++5SxZZjqLiaTYw2K6BFX2OxhzBTb\nIbDjdEUmbWRWOQhNSfOUQTHr/iPzLQtSy5putIl3MsnQOfPHhEESMzl1wjLD2I2j8KJhdAmTiWGU\nYjIxjFJMJoZRisnEMEoxmRhGKSYTwyjFZGIYpZhMDKMUk4lhlGIyMYxSTCaGUYqnUx+N7iMRATA/\nh9VA5k2MEyCRTOiRbdXR6mYdp8FkYjSOhEyY6pVOdawXBExcW1QXk4nRPEMWhe1or4hk4NqkephM\njIaRkJjCwl2dM6fn2qp6mEyMpunB2nblc8Imk+C2j8nEaJrogKMdwWRiNM3muFa69WiHMJkYTbM9\ncbp5E8Mo0MD+Ar5hMjGMUkwmhlGKycRomu3JTju9xZ/JxDg9ne6+g8nEaJ70gKMdwWRiNM2Y9dGu\nkF1NsY5gMjEaRlPmxWaWhISMur2DgMnEaJ5LQkmW7xIWXNc4mweYTIzG0SmX9CUBCWRAj4uur2C0\nRb7GCdCRjElkAMz1uWtr6mMyMU6Cdr6hVcQaXYZRisnEMEoxmRhGKSYTwwnyQt7IPTCR/5AXrq3Z\nsO79t38kKir3b/8FmEwMJ8gHfMJXyXZn/gqfyNC1RSvWfVt+oM/e45v8gtx//a3/k2eirm2qYrZt\neH1myM94e+WAqjcPbHkm86/K3/LLAHzEt5CPvTHOeDrIcE0kIPLatVVLvqvyZ7lI4CPeQ7/WmbiJ\ndD5zoLHkt1HWWwe/4csvLC9/id8tvB/zKx3pm6T8uWsTjKfDz20c6UTfxDgvZEiycfDv9Ldc25Vb\n95rf/H7Bn/wa/2AyMVzQnS783/D7Kv/ijXHGk+JPKT6fle+5NqhgzI/1+2/4Vb7BRzzn9+D+/nfM\nmxhOkBf8gK8A8F/8oXrW95T35a/0yyDKJ/oN/fH/Azx7SPrUyoaOAAAAJXRFWHRkYXRlOmNyZWF0\nZQAyMDE1LTA0LTIwVDEwOjU5OjExKzEyOjAw4V25xQAAACV0RVh0ZGF0ZTptb2RpZnkAMjAxNS0w\nNC0yMFQxMDo1OToxMSsxMjowMJAAAXkAAAAUdEVYdHBkZjpWZXJzaW9uAFBERi0xLjUgBVwLOQAA\nAABJRU5ErkJggg==\n"
49      }
50     ],
51     "prompt_number": 33
52    },
53    {
54     "cell_type": "code",
55     "collapsed": false,
56     "input": [
57      "H = (cct.C.V / cct.Vi.V).simplify()"
58     ],
59     "language": "python",
60     "metadata": {},
61     "outputs": [],
62     "prompt_number": 34
63    },
64    {
65     "cell_type": "code",
66     "collapsed": false,
67     "input": [
68      "H(j * omega)"
69     ],
70     "language": "python",
71     "metadata": {},
72     "outputs": [
73      {
74       "latex": [
75        "$\\frac{1}{i C_{} R_{} \\omega + 1}$"
76       ],
77       "metadata": {},
78       "output_type": "pyout",
79       "prompt_number": 35,
80       "text": [
81        "     1     \n",
82        "\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n",
83        "\u2148\u22c5C\u22c5R\u22c5\u03c9 + 1"
84       ]
85      }
86     ],
87     "prompt_number": 35
88    },
89    {
90     "cell_type": "code",
91     "collapsed": false,
92     "input": [
93      "H(j * omega).rationalize_denominator()"
94     ],
95     "language": "python",
96     "metadata": {},
97     "outputs": [
98      {
99       "latex": [
100        "$\\frac{- i C_{} R_{} \\omega + 1}{C_{}^2 R_{}^2 \\omega^2 + 1}$"
101       ],
102       "metadata": {},
103       "output_type": "pyout",
104       "prompt_number": 36,
105       "text": [
106        "-\u2148\u22c5C\u22c5R\u22c5\u03c9 + 1\n",
107        "\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n",
108        " 2  2  2    \n",
109        "C \u22c5R \u22c5\u03c9  + 1"
110       ]
111      }
112     ],
113     "prompt_number": 36
114    },
115    {
116     "cell_type": "code",
117     "collapsed": false,
118     "input": [
119      "H(j * omega).real_imag"
120     ],
121     "language": "python",
122     "metadata": {},
123     "outputs": [
124      {
125       "latex": [
126        "$- \\frac{i C_{} R_{} \\omega}{C_{}^2 R_{}^2 \\omega^2 + 1} + \\frac{1}{C_{}^2 R_{}^2 \\omega^2 + 1}$"
127       ],
128       "metadata": {},
129       "output_type": "pyout",
130       "prompt_number": 42,
131       "text": [
132        "    \u2148\u22c5C\u22c5R\u22c5\u03c9           1      \n",
133        "- \u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500 + \u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n",
134        "   2  2  2        2  2  2    \n",
135        "  C \u22c5R \u22c5\u03c9  + 1   C \u22c5R \u22c5\u03c9  + 1"
136       ]
137      }
138     ],
139     "prompt_number": 42
140    },
141    {
142     "cell_type": "code",
143     "collapsed": false,
144     "input": [
145      "H(j * omega).magnitude\n"
146     ],
147     "language": "python",
148     "metadata": {},
149     "outputs": [
150      {
151       "latex": [
152        "$\\frac{1}{\\sqrt{C_{}^2 R_{}^2 \\omega^2 + 1}}$"
153       ],
154       "metadata": {},
155       "output_type": "pyout",
156       "prompt_number": 37,
157       "text": [
158        "        1        \n",
159        "\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n",
160        "   ______________\n",
161        "  \u2571  2  2  2     \n",
162        "\u2572\u2571  C \u22c5R \u22c5\u03c9  + 1 "
163       ]
164      }
165     ],
166     "prompt_number": 37
167    },
168    {
169     "cell_type": "code",
170     "collapsed": false,
171     "input": [
172      "H(j * omega).phase_degrees"
173     ],
174     "language": "python",
175     "metadata": {},
176     "outputs": [
177      {
178       "latex": [
179        "$\\frac{180.0}{\\pi} \\operatorname{atan2}{\\left (- C_{} R_{} \\omega,1 \\right )}$"
180       ],
181       "metadata": {},
182       "output_type": "pyout",
183       "prompt_number": 38,
184       "text": [
185        "180.0\u22c5atan2(-C\u22c5R\u22c5\u03c9, 1)\n",
186        "\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n",
187        "          \u03c0           "
188       ]
189      }
190     ],
191     "prompt_number": 38
192    },
193    {
194     "cell_type": "code",
195     "collapsed": false,
196     "input": [
197      "H1 = H.subs('C',1e-6).subs('R',1e3)\n",
198      "H1(j * omega)"
199     ],
200     "language": "python",
201     "metadata": {},
202     "outputs": [
203      {
204       "latex": [
205        "$\\frac{1}{0.001 i \\omega + 1}$"
206       ],
207       "metadata": {},
208       "output_type": "pyout",
209       "prompt_number": 39,
210       "text": [
211        "      1      \n",
212        "\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n",
213        "0.001\u22c5\u2148\u22c5\u03c9 + 1"
214       ]
215      }
216     ],
217     "prompt_number": 39
218    },
219    {
220     "cell_type": "code",
221     "collapsed": false,
222     "input": [
223      "from numpy import logspace\n",
224      "w = logspace(1, 6, 500)\n",
225      "%matplotlib inline\n",
226      "H1(j * omega).dB.plot(w, log_scale=True)"
227     ],
228     "language": "python",
229     "metadata": {},
230     "outputs": [
231      {
232       "metadata": {},
233       "output_type": "display_data",
234       "png": "iVBORw0KGgoAAAANSUhEUgAAAYwAAAEVCAYAAADzUNLBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmYFNXZ9/HvgIhhy4AaFVAHcX9ccCNiBAaDICCLRMSd\ncVfUJMYEIT4BjUaIW1xR4zYaA0oEIxpBQBjMg3GJMqIS34iCihqNBmUxCkK/f9zVnqaZ6anurq6q\n7v59rquv6arp6TrctnX63Kfq3CAiIiIiIiIiIiIiIiIiIiIiIiIiIiIxdQzwJvAWcFnEbRERkZhq\nDiwDqoAWQD2wT5QNEhERaBZ1AxrQHeswVgAbgIeBoVE2SERE4tlhdALeT9le6e0TEZEIxbHDSETd\nABER2dJWUTegAR8AO6ds74yNMr7Vvn3HxKpVH4baKBGREvA2sHuufxzHEcbfgT2wSe+tgZHAzNQX\nrFr1IYlEItDHhAkTAn1tY6/xuz/TdmPPCxGLb75JsGpVguXLEyxenGDBggR//nOC2toEN92UoGfP\nCZxzToLBgxMcdliCzp0TtGiRYLvtEnTvnmDkyARHHjmBu+9OMH9+go8/Lt5YhP258LNPsQg/Ftm+\nX5xiAXTN5+QcxxHGN8BFwNPYFVP3Av8o9EGrq6sDfW1jr/G7P9N26vMVK1Y02ZZspb5/8+ZQWWmP\nhhx4YDXp/6RNm+CTT2D5cnjnHZg/v5rnnoMHHoDXX4ett4YDDoDKymoefBCOOAK6doWKii2P39R2\nmLEI4rXZfC787FMsGt4uZCyyiYPf14cVi3xVBPZO4Up4vWXZq6mpoba2Nupm+JZIwAcfwGuvwZIl\n8PLL8NxzsGGDdRxHHAF9+sDBB0OzLMe/xRaLQlIsHMXCqbBvZTmf9+M4wpAs1NTURN2ErFRUQOfO\n9hgwwPYlEvD++7BokT1OPRVWrYL+/e01AwY0PsJJVWyxKCTFwlEsgqMRhsTS8uUwezbMmgULF0Lv\n3jByJAwZAm3bRt06keKU7wgjjpPekoW6urqom1AQXbrABRfAzJk2+jjhBJg61UYmJ51knUj6d4ZS\njUUuFAtHsQiOOgyJvXbtLE315JM28ujRwzqTffeFm26y9JWIFJ5SUlKUEgn461/hzjstdXXWWXDJ\nJdCxY9QtE4kvpaSkLFVUQK9eMGUKLF4M69fDfvvBuefapbwiEjx1GEVO+VnYdVe4+Wa4//46dtwR\nuneHiy6Cjz+OumXR0efCUSyCow5DSsZ3vwu//jX84x/QooXNcYwfD6tXR90ykdKgOQwpWStWWIcx\nbx5cey2ccoq7m1ykHOU7h1Gs//uowxDfnn8eLrwQWreG22+H/fePukUi0dCkd5lTftZpLBaHHw4v\nvggnnww//CH8/Ofw5Zfhti1s+lw4ikVw1GFIWWjeHM4/H5YuhY8+ggMPhGefjbpVIsVFKSkpSzNn\nwujRMGwYTJyo5UakPCglJZKDIUNsxdx166BbN3jhhahbJBJ/6jCKnPKzTraxaN8e7r8frrvOOpCr\nr4aNGwvTtrDpc+EoFsFRhyFlb/hwq8sxf77V4nj33ahbJBJPmsMQ8WzcCDfcANdfb2tUDR8edYtE\ngqX7MEQC9tJLMGIEHH+8TYi3aBF1i0SCoUnvMqf8rBNULA47zFJUb7xh92189FEgbxsqfS4cxSI4\n6jBEGrDttvCXv8DRR8Ohh1rBJpFyp5SUSBPmzIHTT4dLL7W7xLUelRQrzWGIhOC99+C442wF3Lvv\nhm22ibpFItnTHEaZU37WKWQsdtnFKvxt2AC9e8OHHxbsUIHQ58JRLIITVYcxAngD2AgcnPa7ccBb\nwJtAv5DbJdKoVq1g6lS7ye/737erqUTKid+hyT5AFbAJeBc7medjb++97gIuBV7x9u8LTAEOAzoB\n84A9vdemUkpKIvX443D22Vbp7+STo26NiD/5pqS2yvC7LsAlwEDgA+BD70A7AZ2BJ4HfAStyOG5j\nHc5QYCqwwXvfZUB34PkcjiFSMEOHwm672c/XX7dlRZopwSslLtNH/LfAE9joojdwEnCi93xv4C/A\ntQG3pyOwMmV7JTbSkEYoP+uEHYv997c6G88+CyedBF99FerhM9LnwlEsgpOpwzgBmIt920+3AZjj\nvaYxc4HXGngMzrKNyj1JbG23nZWAbdbMbvL79NOoWyRSOJlSUgC7AuuAT4EewJFYmugxH+99dA7t\n+QDYOWW7s7dvCzU1NVRVVQFQWVlJt27dqK6uBtw3inLYrq6ujlV7ynH7+efrOOccmDevmh49YMKE\nOjp3jr59SVHHJ+rt5L64tCfM7bq6OmprawG+PV/mI9Pkx3hglPd8KtAXqAO+DywBfpL30WEB8HPg\nZW87OendHTfpvTtbjjI06S2xdM898L//C48+CkceGXVrRDZXyPswTsJO4IdincMPgbHYyKFvrgf0\nHAe8DxyOzYXM8vYvBaZ5P2cBo1FKKqP0b5PlLA6xOPtsePBBW+l26tTo2hGHWMSFYhGcTCmpr4Cv\nvccyLDUF8A2wPs/jPkbjaa1rvIdIUerXD555Bo49FpYvh3HjtJyIlIZMH+N3sHRRBXCd95yU7d0K\n27SMlJKS2PvwQxg0CA4/HG67DZo3j7pFUu4KuZZULS4dVMGWqaEzcj1oANRhSFFYvdrSU+3awZQp\nWoNKolXIOYwarFM4I+158iExoPysE8dYtGtny6S3bGmpqs8/D+e4cYxFVBSL4GSaw7iUzBPONwbc\nFpGS1LIl/PGPtjx6z54wezZ00u2oUoQyDU2uwDqMvbC1nWZ6rz8WeBE4tdCNy0ApKSk6iQRcdx1M\nngyzZsE++0TdIik3YdTD+Cu2ntQab7st8BTQM9eDBkAdhhStBx+EMWPgscegR4+oWyPlJIx6GN9j\n8+VBNnj7JAaUn3WKJRannw73328LFz75ZGGOUSyxCINiERw/HcaDWArqCuBK4AXggQK2SaTkDRhg\nncU558B990XdGhF//A5NDsFSUAngWWBxwVrkj1JSUhL++U845hjrOMaO1Q1+UliFnMNoi5u3yOc1\nhaAOQ0rGRx/ZJbf9+9ukuDoNKZRCzmE8BtyOlUntkLK/A9AfuAN/q9ZKASk/6xRrLHbaCRYuhEWL\nbC2qb77J/z2LNRaFoFgEJ1OH0ReYjtW8WAR84T0WAccDj5D/IoQiAnToYHU1Vq6EE06IVzEmkaRi\nHfwqJSUl6euv4bTT4LPP4M9/hrZto26RlJIwLqsVkZC0bGnLonftqgp+Ej/qMIqc8rNOqcSieXO4\n6y7rMHr1sjRVtkolFkFQLILTVIlWEYlARQVMnAjbbmuV++bMgT33jLpVUu785rJ6YqVS7we2B9oA\nywvVKB80hyFl49574Ve/slVvDzoo6tZIMQtjLakrsBv39gL2xGptTwN+kOtBA6AOQ8rKjBlw/vkw\nfbqteCuSizAmvY8DhuJKtH6A3bAnMaD8rFPKsRg+3Aow/ehHNtJoSinHIluKRXD8dBhfA5tStlsX\nqC0ikkHfvvDEE3DWWVZfQyRsfoYmv8DmL/oBE4EzgSnALQVsV1OUkpKy9cYbtv7UuHEwenTUrZFi\nEsYcBlhn0c97/jQwN9cDBkQdhpS15cttxJFctFDEj7Bu3JsD/Nx7RN1ZSArlZ51yikWXLvDss/CH\nP9hII/37UznFoimKRXAydRhrsZVoG3qszvO41wH/AF4FZgDfTfndOOAt4E3cqEZE0nTqZIsWzp0L\nF10EmzY1/Tci+fAzNLka+BB4yNs+BegI/CqP4x4NPINNpk/y9o0F9sXmRw7DLt+dh13Km/6/glJS\nIp4vvoBjj7VRx333wVa6HVcaEUZKaggwGRtVrMaWNR+a6wE9c3GdwAtAZ+/5UGAqVgZ2BbAM6J7n\nsURK2ne/C08/DZ98Yivdfv111C2SUuWnw1gHnAo09x6nYOmqoJwJPOU97wikrpyzEhtpSCOUn3XK\nORatWsHjj0OzZjBkCMyaVRd1k2KjnD8XQfPTYZyM1cT42Huc4O1rylzgtQYeg1NeczmwHktDNUa5\nJxEfWraEhx+2gkxjxliqSiRIfrKdy7G0VLaObuL3NcBA4Icp+z4Adk7Z7uzt2/KPa2qoqqoCoLKy\nkm7dulFdXQ24bxTlsF1dXR2r9mg72u2ttoLTT69j9Wro08dSVW+8EZ/2RbGd3BeX9oS5XVdXR21t\nLcC358t8+Jn8uD9tO/mN/8w8jnsMcAPQG0hd8T856d0dN+m9O1uOMjTpLZJBImELFk6fbpX8Oimx\nK4Qz6f0X4Env8Qx2Cey6jH/RtFuxFW/nAouxSXWApdjChkuBWcBolJLKKPltQhSLVAsX1nH11XDG\nGbZY4TvvRN2i6OhzERw/KalH07anYHW987FHht9d4z1EJE9jxliZ1969LT21775Rt0iKWS5Dk72x\n0cbuAbclG0pJiWThoYfgF7+AJ5+EQw6JujUSlXxTUn5GGGtxaaEEdqXUZbkeUETCd+qp0KYNDBig\nmhqSOz9zGG2w+hdtgXZYOml6IRsl/ik/6ygWTkOxGDbMamoMHw6zZ4ffpqjocxEcPx3GMz73iUjM\n9e1rN/idfrqNNESykSmX9R2gFbAAqE7Z3w6Yjc1lREVzGCJ5WLwYBg6ESZNg1KioWyNhKeQcxnnA\nT7DlOl5O2b8GuC3XA4pI9A46CBYsgH79YM0aW+1WpCmZUlI3AV2wGhhdUh4HoA4jNpSfdRQLx08s\n9t7bamrcdBNcc82WNTVKhT4Xwck0wjgKmI8tbT68gd/PKEiLRCQ0VVXWafTrZ2tPTZoEFTknLKTU\nZfpoXAlMAGpp+G7rMwrRIJ80hyESoM8+szrh3bvDrbfaqrdSesKq6R036jBEArZ6NQwaBF27wj33\nqBBTKQpjLaltsBoYl2MjjgnA+FwPKMFSftZRLJxcYtGund2f8eGHcNJJsH598O2Kgj4XwfHTYTyO\nLW++Abvrey35Lz4oIjHUujU88QRs2GA3+v33v1G3SOLEz9DkdWC/QjckS0pJiRTQhg1QU2OjjZkz\nbQFDKX5hpKSewy6lFZEy0aIFPPgg7L47HH00rFoVdYskDvx0GD2xG/f+iSuzuqSQjRL/lJ91FAsn\niFg0bw6//z306GHV+z75JP92RUGfi+D4uQ5iQMFbISKxVFEBN94IEyZYTQ1V7ytvfnJZHRrYtwab\nBI+K5jBEQnbttXDnnfDMM9ClS9StkVyEUQ/jFWAXIJnFbA/8y3ucw+brTIlIiRozxmpq9OoFc+fa\n0iJSXvzMYczF0lLbeo9jsIp7FwJ3FK5p4ofys45i4RQqFqNHw9VX25xGfX1BDhE4fS6C46fD6AE8\nnbI9x9v3N2DrQjRKROJr1ChbPqR/f3j++ahbI2Hyk8uaC8wDHvZefwLQD+gPvAQcXLDWNU5zGCIR\ne+opu1dj2jSoro66NeJHGGtJbY8tB/IDb3sRtjDhF9jcxrJcD54HdRgiMbBgAYwcCbW1VpBJ4i2M\nG/f+DVwEHOQ9LvL2rSeazkJSKD/rKBZOWLHo08fuBD/jDHj00VAOmTV9LoLj5yqp7wFjgH2xsq1g\ny50flcdxr8LWp0oAnwE1wPve78YBZwIbgR9jcyYiElOHHw5PPw0DBsC6dSr5Wsr8zmE8glXeOw87\nuf8b60Ry1Ra7lwPgYuBA4GysU5oCHAZ0wuZO9gQ2pf29UlIiMfPmm7aMyLhxdjWVxE8YKaltgXuw\nFNRCrHBSPqMLcJ0FQBvgU+/5UGAqdlPgCizl1T3PY4lICPbeGxYuhOuvt5v8pPT46TCSq+L/CzgW\nuyqqfQDH/g3wHjZimejt6wisTHnNSmykIY1QftZRLJyoYrHbblby9b77YPz4eNQJ1+ciOH7mMH4D\nVAKXArcC7YBLfPzdXGDHBvb/EngCK8h0OTAWuInGS742+JGrqamhqqoKgMrKSrp160a1d21f8gOi\n7fLaTopLe6Lcrq+vj+z4y5bVMXEiXHllNWvXwuDBdVRURBePeu8Owzj99wlru66ujtraWoBvz5f5\niEOJ1l2Ap7CaG2O9fZO8n7OxS3pfSPsbzWGIxNyqVXap7f77wx132Oq3Eq0w7sPYDZuYrsKNSBLY\nVU652gN4y3t+MTZPcRpu0rs7btJ7d7YcZajDECkCa9bAkCHQsaPdq9GiRdQtKm9hTHr/GViOpaNu\nSHnkYyJWV6MeqMbSXQBLgWnez1nAaBpJSYlJT8eUM8XCiUss2ra1O8I//xxGjICvvw6/DXGJRSnw\nM4fxFXBLwMc9PsPvrvEeIlICvvMdeOwxOOUUG2089hi0ahV1qyQXfoYmpwFdsQUIU78fvFKQFvmj\nlJRIkfnmGzj7bHjnHXjySWjXLuoWlZ8w5jAmYZ3GMja/ga5PrgcNgDoMkSK0aRNcfDG8+CLMng3b\nbht1i8pLGHMYI4AuQG+sk0g+JAaUn3UUCyeusWjWDG67DY46yla4/de/Cn/MuMaiGPnpMF4jmBv1\nRESoqIBJk2yV21694L33om6R+OVnaLIQOACrfZGcw8j3stp8KSUlUgJ+9zu4+WYr+brHHlG3pvSF\nUdN7QgP7dLYWkbxdconVCa+uthVv99sv6hZJJn5SUnUNPBYWqD2SJeVnHcXCKaZYnHMOXHcd9O0L\nL78c/PsXUyzizs8IQ0SkoE4+2e7NGDAAZsyAI4+MukXSkDisJZULzWGIlKC5c+0GvylTbMQhwQrj\nsloRkVAcfTRMn24jjpkzo26NpPPTYQwGFgOrsMJHa4DVhWyU+Kf8rKNYOMUci549bf2pc8+Fhx/O\n//2KORZx42cO4ybgOOB1tiyVKiISuEMPhXnzoH9/qxN+1llRt0jA/30YRwEbC9yWbGgOQ6QMvPWW\nzWX87Gfwk59E3ZriF8Z9GJdhS40vwJVrTQA35npQERE/9tjDSr727Qtr18Ivf2l3iks0/MxhXAWs\nBbYB2niPtoVslPin/KyjWDilFItdd7VOY+pUGDcu+zrhpRSLqPkZYewEHF3ohoiINGannaCuzuY0\n1q6FW26xhQwlXH4Gd9cCz2D1MOJCcxgiZeiLL2DQIEtV3X03bKVbj7MSRj2MtUArbP5ig7cvAURZ\n/kQdhkiZWrcOhg2D9u3hoYdg662jblHxCOPGvTbe67bB5i7aEm1nISmUn3UUC6eUY9G6NTzxhNUH\nHz4c/vvfzK8v5ViEzU+H0auRh4hIJLbZBh591Mq8Dhpk8xpSeH6GJk/iljPfBugOvIzdmxEVpaRE\nhI0b4fzz4Y037O7wysqoWxRvYcxhpNsZuBkYnutBA6AOQ0QAu8z2kktg4UKYMwe23z7qFsVXFIsP\nrgT2yfWAEizlZx3FwimnWFRUWOW+Y4+F3r3hgw82/305xaLQ/FyUdmvK82ZANywlFYRLgeuA7YD/\nePvGAWdiS5H8GJgT0LFEpERVVMBVV1n1vl694JlnoKoq6laVHj9Dk1Epz78BVgCLAjj2zsDdwF7A\nIViHsS8wBTgM6ATMA/Zky0UPlZISkQbdfjv89reWntp776hbEy+FXktqK6A/cHKuB8jgRmAM8HjK\nvqHAVOx+jxXAMmyS/fkCHF9EStCFF9qlt0cdBbNmwYEHRt2i0tHUHMY3wC5Ay4CPOxSbC1mStr+j\ntz9pJTbSkEYoP+soFk65x6KmBm6+Gfr1g8mT66JuTsnwM4exHPg/YCbwpbfPz2q1c4EdG9h/OTZP\n0S9lX6YhUoO5p5qaGqq8JGVlZSXdunWjuroacP+zaLu8tpPi0p4ot+vr62PVnii2R4yoplUrOOGE\netavh5/+NF7tC2O7rq6O2tpagG/Pl/nwk8u6wvuZfuK+Msdj7oetTZXsfDoDHwDfB87w9k3yfs4G\nJgAvpL2H5jBExJf582HkSHjwQRgwIOrWRCuM+zBOAKb52Jer5Ww56d0dN+m9O1t2VuowRMS3v/0N\nhg6FO++05UTKVRj3YYzzuS9XqWf+pVhHtBQr2jSaRlJSYtLTMeVMsXAUC6euro4ePeDpp21C/KGH\nom5R8co0hzEAGIh9078F1yu1xa1aG4Td0rav8R4iIoE56CC7PyNZJ/y886JuUfHJNDQ5EDgI+DXw\nK++1CWANVq51VcFb1zilpEQkJ2+/bSVfL7oILr006taEK4w5jHbAOuzOa4Dm2GW2Xzb6F4WnDkNE\ncvb++9ZpnHwyjB9fPnXCw5jDmAN8J2W7FTYZLTGgXLWjWDiKhdNQLHbe2eqET58OY8ZkXye8XPnp\nMLbBqu4lrcE6DRGRorXDDlYnfOFCGD0aNqUvQCRb8DM0WYQtAphccPBQbEHCHoVqlA9KSYlIIFav\nhsGDYddd4b77SrtOeBhzGIcBjwAfets7ASOBv+d60ACowxCRwHz5pd2f0aYNTJlSunXCw5jDeAlb\nUfZ84AJgb6LtLCSFctWOYuEoFo6fWLRqBY8/bhX8hg1ruk54ufLTYbQGxgI/AV4DqoBjC9gmEZHQ\ntWwJ06ZB+/YwcCCsWRN1i+LHz9BkGjZ/cTrwP1gH8hx2n0ZUlJISkYLYuBEuuACWLLHl0du3j7pF\nwQkjJdUV+C2w3ttel+vBRETirnlzuOsuOOII6NMHPvkk6hbFh58O42s2vw+jq7dPYkC5akexcBQL\nJ5dYVFTADTfYgoUN1QkvV34uILsCW2a8M7aS7A+AmsI1SUQkehUVcOWVVr2vZ09bh6pLl6hbFS2/\nuaztsHoVFVi51E8L1iJ/NIchIqGZPBkmToS5c4u7Tniha3on37w3cCS2+GAL4LFcDygiUmxGj1ad\ncPA3hzEZOA+rv/2693xyIRsl/ilX7SgWjmLhBBWLUaNcnfAX0muAlgk/I4w+WCW85EortViBIxGR\nsjJihN3kN3iw3bPhldEuG35yWU8CFwErvO0q4DaivXlPcxgiEpkFC6xO+AMPFFed8DDWknoWW0/q\nRWwOozu2XMhqb3tIrgfPgzoMEYlUMdYJD2PSe3yG3+msHbG6ujqqy21c3AjFwlEsnELFIlknfOBA\nW7zw1FMDP0Ts+Okw/g78F6u4t5f3mEWwdb1FRIpOudUJ9zM0eQW7pLY9VhvjJWyZkFMK2K6mKCUl\nIrFRLHXCw1hLqgKr3z0cu5x2BLBfrgcUESk1Xbtaydff/97uDi/V77N+Ogyw6nqnAH/J8u+kwHS9\nvaNYOIqFE1YsyqFOuJ8T/0+Bcdjd3W9giw8uyPO4VwArgcXeI/XCtHHAW8CbQL88jyMiEppSrxOe\ncy4rTxOANcCNafv3xRY4PAzoBMwD9sTdNJikOQwRia241gkPYw6jUBpq9FBgKnYF1gpgGXbfh4hI\n0WjXztac+uQTOPFEWL++6b8pBlF2GBcDrwL3ApXevo5YqippJTbSkEYoV+0oFo5i4UQVi2Sd8E2b\nSqdOeCEHSnOBHRvYfzlwB/Brb/sq4AbgrEbep8HcU01NDVVVVQBUVlbSrVu3b2/OSX5AtF1e20lx\naU+U2/X19bFqT5Tb9fX1kR2/ZUu48MI6Jk2CgQOrmTkTXn45vOPX1dVRW1sL8O35Mh9+cll7YZfT\n7ojV9D4AWw7k6ryPbqqAJ4D9gbHevknez9nYfEf62pCawxCRohGXOuFhzGHcDfwSV9P7NeCkXA/o\n2Snl+XHeewLMBE4Etga6AHtga1iJiBStUqkT7qfDaMXm3/AT5L8syG+x+hqvYsWZLvH2LwWmeT9n\nAaPRelUZpadjypli4SgWTlxiUQp1wv3MYfwb2D1l+3jgozyPe3qG313jPURESkqyTnibNsVZJ9xP\nLqsr8HvgCGAVsBy763tF4ZrVJM1hiEhRi6JOeBj1MJJaYymsNbkeLEDqMESk6D3wAIwbF16d8DAm\nvS8FfobV8j7He34W0C3Xg0pw4pKfjQPFwlEsnDjHotjqhPuZwzgEOBS79LUCGIRd1XQ+8Cg2gS0i\nIjkopjrhfoYmf8UWB1zrbbcBngKOAV4G9ilM0zJSSkpESkoYdcLDSEltj7sHA+yS2h2wGhlf5Xpg\nERFx+vSxpURGjYIZM6JuTcP8dBh/xO7DmIAtS/4ctqJsa+x+CYlQnPOzYVMsHMXCKaZYJOuEX3gh\nPPRQ1K3Zkp85jKuwJTp+gN1Edx5W5xuiLdMqIlJy4lwnPJtc1g7ANrg7r98Lvjm+aQ5DREpaIeqE\n5zuH4WeEMQRbTbYj8AmwK/APbCFCEREpgGSd8L59Ye1aGD/e7hSPkp85jKuxmt7/xBYE/CFbrh4r\nESmm/GyhKRaOYuEUcyySdcJnzIhHnXA/HcYG4FPvtc2xet6HFrJRIiJidtjBLrmNQ51wPwOcedgS\n5BOB7bC01KHY2lJR0RyGiJSVIOqEh7GWVGvsfotm2FVR7bBLbT/L9aABUIchImXnyy9h+HBb7XbK\nFNh66+z+Powb98YDG7HUVC1wCzAm1wNKsIo5Pxs0xcJRLJxSikXUdcL9dBj9Gtg3MOiGiIhI01q2\ntDWnOnSAgQNhTYjrh2camlyAVbzrCrydsr8tsIhob9pTSkpEytqmTVYn/NVX/dcJL+QcxneB9sAk\n4LKU164h2vkLUIchIkIiYTf1zZ8Pc+bA976X+fWFnMNoDqwGLsQ6idXeIwF0yPWAEqxSys/mS7Fw\nFAunlGMRdp3wTBdmvYJbBiRdAtgt+OaIiEg2wqwTHvGN5jlTSkpEJE1TdcLDWEsKYCjQCxtZLMSq\n74mISIyMHg2tW8NRRxWmTrify2onAT8G3sAWHfwxdte3xEAp52ezpVg4ioVTbrEoZJ1wPx3GIOxe\njPuAe7HSrMcGcOyLsQ7odTavCz4OeAt4k4bvARERkQxGjLDlQwYPhiD7Sz+5rCVAH9yltNtiCxAe\nkMdx+wC/xG4A3ICVgf03sC9Wze8woBO2jtWeQPpyW5rDEBFpQnqd8DDmMCZiV0zVedu9gbG5HtBz\ngfe+G7ztf3s/hwJTvf0rgGVAd+D5PI8nIlJ2knXChw2DO+7I//0ypaQmA0diJ/AewAxguvf84TyP\nuwc2if481hEll0vvCKxMed1KbKQhjSi3/GwmioWjWDjlHosePWD2bKsTnq9MI4x/AtdhJ/FHsI5j\ncRbvPRdwWwJHAAAKCklEQVTYsYH9l3vHbQ8cjqWfptH4fR0N5p5qamqoqqoCoLKykm7dulFdXQ24\nD4i2y2s7KS7tiXK7vr4+Vu2Jcru+vj5W7Qlzu66ujtraWgCGD69i8mTy4ieXVQWcCIwEWmFzDFOx\nDiVXs7CrrxZ628uwzuNsb3uS93M2MIEtK/xpDkNEJEth1MNIdRBwP7A/tnRIrs7DRi4TsEntecAu\nuEnv7rhJ793ZcpShDkNEJEth1MPYChiCnchnY5e7Ds/1gJ77sBTUa9ho5XRv/1IsPbUUG4WMpvHl\nSQTlZ1MpFo5i4SgWwck0h9EPS0UNAl7ETuznAmsDOO4G4LRGfneN9xARkRjJNDSZj3US04H/hNMc\n35SSEhHJUthzGHGhDkNEJEthzGFIjCk/6ygWjmLhKBbBUYchIiK+KCUlIlImlJISEZFQqMMocsrP\nOoqFo1g4ikVw1GGIiIgvmsMQESkTmsMQEZFQqMMocsrPOoqFo1g4ikVw1GGIiIgvmsMQESkTmsMQ\nEZFQqMMocsrPOoqFo1g4ikVw1GGIiIgvmsMQESkTmsMQEZFQqMMocsrPOoqFo1g4ikVw1GGIiIgv\nmsMQESkTmsMQEZFQRNVhPAws9h7LvZ9J44C3gDeBfuE3rbgoP+soFo5i4SgWwYmqwzgROMh7TPce\nAPsCI72fxwCT0Sgoo/r6+qibEBuKhaNYOIpFcKI+GVcAJwBTve2h3vMNwApgGdA9kpYVic8//zzq\nJsSGYuEoFo5iEZyoO4yewMfA2952R2Blyu9XAp3CaEg2w1Y/r23sNX73Z9ou9BBbsWj82Pm+NptY\n+NmnWDS8XchYZPvepRSLQnYYc4HXGngMTnnNScCUJt4nlMuhivUkuWLFiibbki3FovFj5/vasE4M\nioV7HnQsyrnDiPKy2q2wEcTBwIfevrHez0nez9nABOCFtL9dBnQtdANFRErM28DuUTciF8cAC9L2\n7QvUA1sDXbB/XLHeKyIiUlK2ivDYI3GT3UlLgWnez2+A0YSUkhIRERERERERERERkVhpHnUDAtIF\nuB44DfhTxG2J2lDgUuAU4AvgnWibE6m9gauAUUA7Nl+Cphy1Bp7Drkp8K+K2RKka+AN2U/Ba4N1I\nWxOtCuA3wDBge+DVTC+O+sa9oCwHzo66ETHxOHAucD52YUE5exO4AFuKpn/EbYmDMcAjUTciBjYB\na4CWbH6jcDkaht0cvZ4yjEW5jy5SXQ90i7oRMTAYmAUMj7ohETsa+wIxChgUcVuilrxU/3vAQ1E2\nJAYuA87xnjd5/ozzCOM+bNmQ19L2H4N9c3wL+8eWg2xiUQH8FjtJluKqa9l+Lp4ABmAnylKTTSx6\nA4cDJ2MniFK7vymbWCQv1f8cG2WUmmxisRKLA9jIq2j1xFazTf1HN8fu8q4CWmAnxH2ADsCdlG4n\nkk0sLgb+DtwBnBdqK8ORTSx6AzcDdwE/DbWV4cgmFkmjgIEhtS9M2cTiOOx88TDQK9RWhiObWHwH\nuAe4BUvfFrUqNv9H98CWC0kai1tOpNRVoVgkVaFYJFWhWCRVoVgkVVGAWMQ5JdWQTsD7KduhrWYb\nQ4qFo1g4ioWjWDiBxKLYOgwtE+IoFo5i4SgWjmLhBBKLYuswPgB2TtnemTK8FMyjWDiKhaNYOIqF\nUxaxqGLzPNxW2Aq2VdiKtukTeqWsCsUiqQrFIqkKxSKpCsUiqYoyi8VU7I7Ur7Hc2xne/gHA/8Nm\n/MdF07TQKRaOYuEoFo5i4SgWIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiVsGLbG/l4FPMbaPP++\nJ/AG8AqlVRvhWOCKPN+jDjgkZXssVk+jIQcA9+Z5PBEpY48AM8n/xJXJmixeW8GWhYPuxOqfp9sq\n5xbFwwJghwb2N8/yPQ5O2Z4PbJvh9XVYFTsRkay0AVYAuwD/SNlfjZ1Y/uTtTy2ROdDb93esoMsT\n3v4rgEtTXve6977gOow2wDzgZWAJMMTbX4Uth/CA93epC7CdDXwGvOO1ozfwV6wu+pvY4pzXAS8C\nr2K10sE6ndu818wF/gL8yPvdCqygF8Ch2EkXoDVWCe0FbDSTbF8NMAOrkPhPrFpi0jHev6feO06F\n95rtvN83w4qGpZ/EdwYWpWzXYh3j81jp3sOA57x2LAL29F73Hayo0FKvTc/jRhjtgP/zno/A1imq\nBxamHOcy4EJERLJ0CnaSAngW9021GisJ2RE7AT4HHAFsA7wH7Oq9bgo2OgGYwOYdxmts2WE0B9p6\nz7fDTqRgHcZGoHsj7bwfV++7GktxJdtwLnC597wl8JL3fsOBOV77dwJWpbzHchruMK7BjWQqsU6s\nFdZhvO21vSXW4XQCtmfzeFR6P8cDP/Ge96PhOswnArem/Rtn4kZXbXEjjb7Ao97zn2HV1wD2Bzbg\n/rsNx40Ul3j/brCOJKkPNqqUElXsw26Jr5OA33nP/+Rtv+Jtv4gtjgb2LbUL8CX2Tf9db/9U3Dd6\nP5oBE7E5iU1Yh5RMj7zrHbMxqWmqF1Pa0A87cR7vbbcD9vCOMQWrMfARlqppSj9gMPBzb7sl1ukl\ngGdwHd9SrFPqgHW0ybYk6y7fh42AbgbOxDqDdLt47Ur1J1xNhErgQWB3b1/yPNDTe1+wTnlJyt/3\n944NNip5AJiGjUSSPvLaLiVKHYYUQgfs2+Z+2AmpuffzF97vv0557Ubsc5he4CX1JP4Nm9du2aaB\nY56CjSwO9t5zecrr1jXR3tRjp7/2IiwdlGpgWvsaa2t6O4fjRj5J38dfPJJWAh8DR2GppZMaeV36\nXM2XKc+vwjqp47AT/IKU36X/XVJ34Hzv+QXe9iAsZXYI8B/vb1W0qIQVWwElKQ7HY99gq7DRwy7Y\nCbxnI69PYCma3XApmJG4k88KXGrkYO8907UDPsFOuH1S3sePxk6STwOjcV+s9sTSSM967WuGpWaq\nU/5mBZaKAjevkXyvH6dsH5Th2Als/qAX7ht7h5Tf34PNuUyj4RP0u8CODf2DPO1wI7yalP3P4q6C\n2g+78gngf7D5muSxumIjsQnAv4HO3v6dcCMiKUHqMKQQTgQeS9s3Hfs2nKDhk9xX2Ml5Njbpvdp7\nJP+2AzZpfSHWuSQl3+uP2Il6CXAam0+0N/WtN5HyM/W192ApolewFM0d2GjpMWyksBRLzfwNd+K/\nEkvrvISNNpLvdxXQwmvf697rGjpm0qdYSm4GlrabmvK7J7BJ9IbSUWApo4PT9qUe41osffcKbvSH\n9+9r4/27rsT+O1RgdRRmpf39Eiwmi3Cpq+5YpyMiUnCtU57fjpvcjbv72Xw0UWiHsvnVSQ2Zj5uY\nztccGr5EN10duqxWRELyU2AxdiPdH2h4riKOUq+0KrSxWNrriCZeNxA3ignDAbgrrERERERERERE\nRERERERERERERERERERK1f8HmcQCFnRNfeEAAAAASUVORK5CYII=\n",
235       "text": [
236        "<matplotlib.figure.Figure at 0x7f637db2d350>"
237       ]
238      }
239     ],
240     "prompt_number": 40
241    },
242    {
243     "cell_type": "code",
244     "collapsed": false,
245     "input": [
246      "H1(j * omega).phase_degrees.plot(w, log_scale=True)"
247     ],
248     "language": "python",
249     "metadata": {},
250     "outputs": [
251      {
252       "metadata": {},
253       "output_type": "display_data",
254       "png": "iVBORw0KGgoAAAANSUhEUgAAAYwAAAEVCAYAAADzUNLBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xu8VXP+x/HXqZQoc2ri16CcririKKWY0QnlXrmVWzoT\ngwqlRGUo18l1jEuuw8Eo90uhlMvOoKJy6EJXIcwwklJStH9/fNZu7U7nss7ea6+1L+/n47Efe6+1\n197r08e2vuf7/a7v9wsiIiIiIiIiIiIiIiIiIiIiIiIiIiKSpo4BPgWWAVeEHIuIiKSpmsByoADY\nCSgF2oYZkIiIQI2wAyhHZ6zAWAVsAZ4EeocZkIiIpGeBsRfwZdz2amefiIiEKB0LjGjYAYiIyI5q\nhR1AOb4CmsRtN8FqGds0bLhndM2arwMNSkQkC6wAWib64XSsYcwFWmGd3rWBfsDk+APWrPmazp2j\n7LtvlMaNo9StG6VmzSgNGkRp1ixKYWGUbt2i9OoVZcCAKJddFuWmm6I8/HCUKVOizJkTZeXKKOvX\nR9m6NUo0GmXs2LFEo1FPDy/HVnSM1/2VbVf02q+HcpEeufCyT7kIPhfV/b50ygXQIpmLczrWMH4F\nLgJew+6Y+ifwSdmD5szZfnvzZli3Dn78Edautecff4QffoDvvoNvv4VFi+x1/CMahcaNoV69IpYu\nhb333vHRuDHUiCtai4qKqvxHVHSM1/2Vbce/XrVqVZWxVJeXf191jlUuqj6mvP1e9ikX5W+nMhfV\nyYPX44PKRbLyfPumYEWd0jJpGzbAf/4Dq1fb48sv3dex7XXroHlzaNnSfbRqZc9Nm25fmAStuLiY\nkpKS8AJII8qFS7lwKReuvLw8SOK6n441jEDtuiu0aGGPimzYACtXwrJlsHw5fPghPP20bf/4I+y3\nH7Rvv/2jUaNg4i8uLg7mRBlAuXApFy7lwj85X8NI1tq1sHAhLFiw/SM/Hzp3hkMOseeOHWGXXcKO\nVkRyWbI1DBUYKRCNWu1jzhx4/317XrgQ2rWDI46wxx//CPXqJX+uSCTiaxtlJlMuXMqFS7lwqUkq\nDeXlQevW9ujf3/Zt2gQffABvvgnjx8PcuXDggXDkkXDCCXDwweH2hYiIVEU1jJBs3AjvvQczZsCU\nKXY314knQq9eVojUrRt2hCKSbdQklSWWLbOC46WX4OOPoU8fOOss6N4datYMOzoRyQbJFhhqBEkT\nrVrB8OEwcyYsXgwHHACjRtltu5ddZvvKE4lEAo0znSkXLuXCpVz4RwVGGvrDH+DSS62f4/XXoU4d\nOOooq2088wxs2RJ2hCKSi9QklSE2b4YXX4QJE2DpUrjoIhg82G7fFRHxQk1SOaJ2bejbFyIRmD4d\nliyxkeajR8N//xt2dCKSC1RgZKD994dHH7Umq08+idC2LYwYAd9/H3Zk4VJbtUu5cCkX/lGBkcEK\nCmDYMJtUceNGaNMGbr4Zfv457MhEJBupDyOLLFliTVRz59rgwDPOsEGEIiKgcRhSjnfegSFDYPfd\nrZO8deuwIxKRdKBO7xxXXvvsH/8I8+bBccfBoYfCuHG5cSuu2qpdyoVLufCPCowsVauWDQQsLbUJ\nELt0qXjwn4iIF2E1SZ0GjAPaAJ2A+XHvjQYGAr8BlwDTy/m8mqSqIRqFhx6CMWPsMXSoJjoUyUWZ\n2ofRBtgK3A+MwC0w2gETsUJkL+B1oLVzbDwVGAlYscJmz23YEB57zJ5FJHdkah/Gp8DScvb3BiYB\nW4BVwHKgc3BhZZ7qtM+2aGFzVbVubQs6zZ2burjCoLZql3LhUi78k24NE3sCq+O2V2M1DfHJTjvB\n7bfDLbfAscfCww+HHZGIZIpUNknNABqXs38MMMV5/RbbN0ndBcwGnnC2HwJeBZ4v8x1qkvLBkiW2\neNOpp8INN6hfQyTbpfOKez0S+MxXQJO47b2dfTsoLi6moKAAgPz8fAoLC7ctwxirgmq76u1Zs6B7\n9wizZsHUqUXUrZte8Wlb29pOfDsSiVBSUgKw7XqZDK8lTVugAOt8/hzrg/DDW8BlwDxnO9bp3Rm3\n07slULY6oRqGI+LDesW//AIDB8LKlfDqq9CggT+xBc2PXGQL5cKlXLhS2endDLgT63i+DzgH+DN2\nZ9MK4B9YIZKIk4AvgS7AK8BUZ/9i4GnneSowmB0LC/FZnTrwr3/ZIL/u3eHbb8OOSETSUWUlzdPA\ng0AEu2sp3k5Ad+A8oG9KIqucahgpEI3CNdfAk0/CG2/AXrrdQCSrZOo4jGSpwEihW26Be++Ft96C\nffYJOxoR8UsQ4zD6Ars5r68CXgA6JHpC8Vesg8tPI0faaPCjjoL//Mf3r0+ZVOQiUykXLuXCP14K\njKuAdcAfgSOBfwL3pjIoCd/QoTBgAPToAWvWhB2NiKQDL1WTUqAQGA8swMZIfAgclMK4qqImqQBE\no3DFFbYs7BtvQP36YUckIskIog/jFWwsRA+skNgEzAEOTPSkPlCBEZBoFAYNsnmoXn3VRoqLSGYK\nqg9jGtATWAs0AEYmekLxV6rbZ/Py4J577NbbIUOsAElXaqt2KRcu5cI/XgqMDcB3WB8GwK/Y2AzJ\nETVrwqRJMGeOzUMlIrnJS9VkHNAR2BebanwvbIzGYakLq0pqkgrBl19C165w993Qp0/Y0YhIdQXR\nJHUSNu34Bmf7K0DdnzmoSRN48UX4y19g0aKwoxGRoHkpMH5h+wWMdk1RLJKAoNtnDz4Ybr0VTjkF\n1q8P9NRVUlu1S7lwKRf+8VJgPIPNH5UPnA+8gU07LjlqwADo1g3OPTe9O8FFxF9e27J6Og+A17C1\nLsKkPoyQbdoEhx1mS74OGxZ2NCLiRVBzSRUArbCCYhegJhBmg4QKjDTw2WfQpYv1a3TtGnY0IlKV\nIDq9z8eape5ztvcGXkz0hOKvMNtnmzWD++6zWsZPP4UWxjZqq3YpFy7lwj9eCowh2BiMdc72UmCP\nlEUkGeWkk+Dww+HSS8OORERSzUvV5H1sBbzY/FG1sDW4D0hhXFVRk1QaWb8eCgvhtts0PkMknQXR\nJDUTuBLru+iBNU9NSfSEjluAT4CPgOeB38W9NxpYhi0D23PHj0q6qV8fHn8cLrwws6ZDF5Hq8VJg\nXIFNDbIAuAB4FfhrkuedDuyHTWC4FCskwNb07uc8HwNM8BhjzkqX9tlDD4XzzoPzzw/vVtt0yUU6\nUC5cyoV/qroY18LW134AONV5PEjy62zPwB0MOAfrSAcbUT4JWxJ2FTZnVeckzyUBueoqm9X22WfD\njkREUsFLW9ZLwCXA5ymKYQpWSEwE7gJmY2tugA0QnAo8V+Yz6sNIU7Nm2SjwhQuhYcOwoxGReMn2\nYdTycExDYBHW+R2bTyoK9KriczOAxuXsH4PbB3IlsBkrLCpSbslQXFxMQUEBAPn5+RQWFlJUVAS4\nVVBtB7/dtSscckiEs8+GV18NPx5tazuXtyORCCUlJQDbrpfJ8FLSFFWwP5LkuYuBv2DLvm5y9o1y\nnsc7z9OAsVizVTzVMByRSGTbDyVdrF8P++0HjzwCRx4Z3HnTMRdhUS5cyoUriBpGJNEvr8Qx2CJM\n3XALC4DJWG3jdmwa9VZYzUYySP36MGECXHABLFgAdeuGHZGI+MFLSVPeFCA/Ah8AI4CVCZx3GVAb\nWONszwIGO6/HAAOxhZqGYnNXlaUaRgY4+WTo0AH+muw9dSLiiyDmkroe+BLrmAY4HWiBDeS7kIqb\nrFJJBUYG+Owzmw69tNTW0hCRcAUxcK8XNr35OufxAHA08CS2vreEKNbBlY6aNYOLLoLLLgvmfOmc\ni6ApFy7lwj9eCoyN2GC6Gs6jL26/g/7Ml0pdcQXMng36f1Yk83mpmrQA/gF0cbZnA8OwpVo7Au+k\nJrRKqUkqgzz7LFx7LcyfD7W83GYhIikR1HoY6UYFRgaJRu322tNOg0GDwo5GJHcF0YexL7Ys6yJn\n+wCSn0tKfJIJ7bN5eXDLLVbLSOW6GZmQi6AoFy7lwj9eCowHsVtdNzvbC4AzUhaRZKWOHaF7d7j9\n9rAjEZFEeamazAUOxl0PA6AUKExVUB6oSSoDffYZdOoEixfDHlqCSyRwQTRJfQe0jNs+Ffgm0RNK\n7mrWDM4+25qmRCTzeCkwLsLGYbQBvgYuBdR1mSYyrX32r3+FJ5+EZcv8/+5My0UqKRcu5cI/XgqM\nFdgEgY2wDvDDsLUqRKqtUSMYPtzWzhCRzFJZW9aIuNfldRiE2X2pPowMtmEDtGgBM2ZA+/ZhRyOS\nO1LZh1EfqIcNzhuEzR67NzZ/VIdETyiy664wciRcc03YkYhIdVRWYIwDrgGaYAXECGA4VoDsk/LI\nxJNMbZ8dNAjefRc++si/78zUXKSCcuFSLvzjpQ9jD2yN7Zgtzj6RhO2yC1x+OYwbF3YkIuKVl7as\nK7HJB593ju8DPAXcmMK4qqI+jCzw88/QsiW8/DIcdFDVx4tIcoKaS6oj8Ces8/ttbBBfmFRgZIk7\n74Q33oCXXgo7EpHsl+pO75h5wB3YrLUfVnBMdVwHfISNGH8D6yeJGY2tyPcp0DPB788Zmd4+e/75\nMHeuLbKUrEzPhZ+UC5dy4Z/KCowXgHuwi3bDuP0NsQWU7nWOScTNwIHY9CIvAmOd/e2w5q922Lrf\nE6qIUTLczjvDpZfC+PFhRyIiVamqanIEcCY2WG9PZ9/X2BoYTwARH2IYDfwOGOW83grc5Lw3Dbtb\na3aZz6hJKousXw/Nm8OsWdanISKpkWyTVFXL2bzpPFLhBqA/8DPQ2dm3J9sXDqux8R+SxerXt9ts\nb74ZHngg7GhEpCKpXP9sBtC4nP1jgCnY3VdXYjWLO4A/V/A95VYliouLKSgoACA/P5/CwkKKiooA\nt80yF7bj22fTIZ5Etzt0gIEDixg3DpYuTez7YvvS4d8T9nZpaSnDhg1Lm3jC3L7jjjty+vpQUlIC\nsO16mYx0WHGvKfAqsD9WeADEWrSnYf0bc8p8Rk1Sjkgksu2HkumGDoXatW2xpURkUy6SpVy4lAtX\npi7R2gq7EwrgYqxJqj/W2T3R2d4LeB2bWr1s6aACIwt98YWNx1i+HBo0CDsakewTxHoYYGMwYk1G\nuwPNEj2h42/Yyn2lQBHuRIeLgaed56nAYCpokpLs07Qp9OoFd98ddiQiUh4vBcY44HLsDiaA2sC/\nkjzvqUB77LbaU4Bv4967EatVtAFeS/I8WS++/T4bXH453HUXbNxY/c9mWy6SoVy4lAv/eCkwTgJ6\nAxuc7a9IfMCeSKXatoWuXeGxx8KORETK8tKW9T7WpxBb03tXYBZwQArjqor6MLLYzJlwwQW29ncN\nDdsU8U0QfRjPYEu05gPnY1N5PJToCUWqcvjhtmbG1KlhRyIi8bwUGLcAzzmP1sBVwJ2pDEq8y8b2\n2bw8my7k73+v3ueyMReJUi5cyoV/vBQYu2K1isuAB4G6wE6pDEqkb1/45BN/F1gSkeR4acuaD/wR\naAC8C3wAbAbOSmFcVVEfRg4YPx6WLIFHHgk7EpHsEMTAvVhn98VY7eJmbGryAxM9qQ9UYOSANWts\nMsLFi6FxeZPMiEi1BDVwrytWo3ilmp+TFMvm9tmGDeGMM2DCBG/HZ3Muqku5cCkX/vFy4R+GDdp7\nAVgEtADeSmVQIjFDh8L999tyriISrnSYfDARapLKISecAKecAn+uaD5jEfEkiD6MPbCpQdphfRhg\n8zsdkehJfaACI4dMmwZjxsC8eXbLrYgkJog+jCew9bWbY/NKrQLmJnpC8VcutM/27Gmr8s0uu+5i\nGbmQC6+UC5dy4R8vBcbvsZHdm4GZ2Ky1YdYuJMfUqAFDhmgWW5GweamazAa6ANOxEd5fY9OFtEhh\nXFVRk1SOWbsWmjWzwXy6xVYkMUE0Sd2AzSM1Ahvt/RBwaaInFElEfj706wcPPhh2JCK5y0uBMQVY\niy14VAR0ACb7dP4RwFagYdy+0dhqfJ8CPX06T9bKpfbZIUPgvvtgy5by38+lXFRFuXApF/6p5eGY\nPYC/AAVxx0eBgUmeuwnQA/g8bl87oJ/zHFuitTVWqEiOa98eWrWCF16wuaZEJFhe2rJmAW8D83Av\n3FFs9tpkPANcB7wEdATWYLWLrcBNzjHTsDuzyt4foz6MHPXss3DnnfD222FHIpJ5ku3D8FLDqAtc\nkegJKtAbWA18XGb/nmxfOKzGahoiAPTuDcOG2Sy2B4Y5m5lIDvLSh/EycHwC3z0D6/co++iF1STG\nxh1bWYmnqkQlcq19dqed4MIL4Z57dnwv13JRGeXCpVz4p7Iaxk+4F+sx2DiMWHdjFNitiu/uUcH+\n/YFm2Iy3AHtjzV2HYOuFN4k7dm9n3w6Ki4spKCgAID8/n8LCQoqKigD3B6Lt7Nxu1y7C+PFw001F\nNGiw4wUh7PjSYbu0tDSt4glzu7S0NK3iCXI7EolQUlICsO16mYx0mGjhM9w+jHbARGwN8Vind0t2\nrGWoDyPHnX02HHQQjBgRdiQimSOIcRh5wCnA34HbgJMSPVkF4q/8i4GnneepwGDUJCXlGDIE7r0X\ntur+OZHAeCkwJgAXYB3Ui4ALnX1+aY7VLmJuxGoVbYDXfDxPVirbHJMrunSB3XaD6dPdfbmai/Io\nFy7lwj9e7pLqjjUVxf6WK8FqACKhycuDwYNtcaVjjgk7GpHc4KUt62XgImyWWrABfHcDJ6QmJE/U\nhyFs3AhNm8LcueBDf55I1guiD2M34BNsptoIVruoj00Z4tcUISLVtssucM45Nl2IiKSel5KmqJx9\nUeezUawgCZpqGI5IJLLtdrpctGwZHHYYfPEFzJ6d27mIl+u/i3jKhSuIkd6RRL9cJNVatbLba595\nBpo0qfp4EUlcOozDSIRqGLLN5Mlw441Vr8gnkuuC6MMQSWvHHw/ffGNrfotI6qjAyHC6xxxq1rT5\npa66KhJ2KGlDvwuXcuEfLwXGicCHwA/AeuexLpVBiVTXuefalOdr1lR9rIgkxktb1gpsOpCFpM9C\nRurDkB2cfTZ06ADDh4cdiUh6CqIPYzU2JUi6FBYi5RoyxEZ+a34pkdTwUmBcgU0EOBpbg3sEoL/h\n0oTaZ12bNkV2mF8qV+l34VIu/OOlwLgOWxtjZ6Ce86ifyqBEEhE/v5SI+M9LW9ZCbNGjdKI+DCnX\nxo02gG/ePM0vJVJWEH0YrwJHJ3oCkSDF5pe6//6wIxHJPl4KjMFYH8YmdFtt2lH7rCuWi8GD4eGH\nYdOmcOMJk34XLuXCP14KjHrOcTtjfRf1qXo976qMw+6++tB5HBv33mhgGfAp0DPJ80gOatUKCgtt\nfikR8Y+XtqzDK9j/dhLnHYvVVG4vsz+2pncn3DW9W7PjLb3qw5BKvfQSjB8Ps2aFHYlI+ghittrL\ncdfV3hnoDMwDjkj0pI7ygu4NTAK2YAs2LXfOp2nlpFpOOAEuuQTmz7fBfCKSPC9NUidg04OcCPTA\n7pha68O5LwY+Av4J5Dv79sSaqmJWYzUNqYDaZ13xuYjNL3XPPeHFEyb9LlzKhX8SmXxwNdDWw3Ez\ngAXlPHoB9wLNgELgG+C2Sr5HbU+SkHPPheefhx9+CDsSkezgpUnqrrjXNbCLvJeJpHt4jOEhbLlX\ngK+A+GVw9nb27aC4uJgC50b7/Px8CgsLt62qFfuLIhe2i4qK0iqedNs+/ni48soIffumRzxBbsek\nSzxhbcf2pUs8QW5HIhFKSkoAtl0vk+Gl82NA3Otfsb6Fd5M87x+wmgXApVgn95m4nd6dcTu9W7Jj\nLUOd3uLJrFnQvz8sXQo1NJm/5LhUD9yrhQ3ae9R5PEHyhQXATcDHWB9GN6zQAFgMPO08T8XGgKhk\nqETZvyZzWXm56NIF6teHGTOCjydM+l24lAv/VNUk9SvQFKgD/OLjec+p5L0bnYdI0vLybBbbe+6B\nozVfgUhSvFRNHgfaAJOBjc6+KDuOoQiSmqTEM80vJWKCmEtqBfCKc6xmq5WMo/mlRPzhpcBYjE3l\ncU3c45MUxiTVoPZZV2W5yLX5pfS7cCkX/vFSYIz2uE8kbWl+KZHkVdaWdSxwHNAPeDLu2PrY7a+d\nUxtapdSHIdWm+aUk16WyD+NrbIDeJud5HjAX6/zW/SaScU44Ab7+2uaXEpHqq6zA+AgoweaO+pfz\n+lHgJfy9xVaSoPZZV1W5qFkTLrggN+aX0u/CpVz4x0sfxnSgbtz2LtgIbJGMc955ml9KJFFe2rJK\nsfmjqtoXJPVhSMLOPtumPB8+POxIRIIVxDiMDUDHuO2DgZ8TPaFI2AYPhgkT4Lffwo5EJLN4KTCG\nAc8A7ziPp7C1LCQNqH3W5TUXXbvC738PU6ZUfWym0u/CpVz4x8v05h8A+zqPPGyt7S2pDEoklfLy\nYMQIuPVW6NMn7GhEMoeXtqxdgeHYJIR/AVphhcfLKYyrKurDkKT8+qsN5ps0yWa0FckFQfRhPAJs\nBg51tr8Gbkj0hCLpoFYtGDYMbqtsrUcR2Y6XAqMFtn7FZmd7Q+rCkepS+6yrurkYOBDeegtWrkxN\nPGHS78KlXPjHS4HxC9uPw2iBBu5JFqhf38Zl3HFH2JGIZAYvbVk9gSux+aNmAIcBxcBbSZ77YmxF\nvd+w6dOvcPaPBgY6+y/BBg6WpT4M8cVXX0H79rB8OTRsGHY0IqmVbB+G1w82Ag5xjp8N/C/REzq6\nA2OwyQ23ALsD3+Gu6d0Jd03v1sDWMp9XgSG+GTAA2rSB0ZqDWbJcEJ3eedi620cBRwB/SvRkcQYB\nf8O9Pfc757k3MMnZvwpYTriz4qY9tc+6Es3FiBFw113ZtVaGfhcu5cI/XgqMCcAFwMfAQuf1hCTP\n2wo4HKutRLDR4wB7AqvjjluN1TREUuaAA2yqkEceCTsSkfTmpWryKdZUFGsWqoGtwtemis/NABqX\ns/9K7LbcN4GhWPPTU0Bz4C6sEHnCOfYh4FXg+TLfoSYp8dWsWXDGGbBsGey0U9jRiKRGsk1SXkZ6\nL8cG7a1ytps6+6rSo5L3BuEWAh9ghVEj4CugSdxxezv7dlBcXExBQQEA+fn5FBYWUlRUBLhVUG1r\nuzrbLVoUMXEi7LNPesSjbW0nux2JRCgpKQHYdr1MhpeS5m2sFvA+EMX6FD4A1jnbvRI47wVY89NY\nrFP7dawginV6d8bt9G7pnCeeahiOSCSy7YeS65LNxZtvwqBBsHixrZ2RyfS7cCkXriBqGFdX8l6i\nV+2HnccCbEDgOc7+xcDTzvOv2G23KhkkEN272621zz0HffuGHY1I+vFS0tTDpjP/DXcSwqmEOwGh\nahiSEq+8AmPGQGmpTVIokk2CuK32baAO1kT0GtAfW65VJOscdxzUqAEvhzm1pkia8joOYyNwMnY7\n7WnYOt+SBmIdXOJPLvLyrIZx/fWQyZVY/S5cyoV/vBQYAF2Bs7ApPKrzOZGMc/LJsGEDvPpq2JGI\npBcvbVndgBHAu9istS2w8ROXpDCuqqgPQ1Lq+eetljF3rjVRiWSDoOaSSjcqMCSlolHo1AlGjYJT\nTw07GhF/BNHpLWlM7bMuP3ORl2c1jKuvht9+8+1rA6PfhUu58I8KDJEKHH20jcuYODHsSETSg5qk\nRCoxc6atzPfpp5pjSjJfEE1S+wJvAIuc7QOAvyZ6QpFM0q0btGgBDz4YdiQi4fNSYDyILXYUW9N7\nAXBGyiKSalH7rCtVubjlFrj2Wvjxx5R8fUrod+FSLvzjpcDYBZgTtx0l3GlBRAJ14IE2Anz8+LAj\nEQmXl7asqdj6288ABwGnAucCx6YwrqqoD0MC9dVXttDS/Pmwzz5hRyOSmCDGYbQAHgAOBX4APsNG\nfa9K9KQ+UIEhgRs7FpYvhyeeqPpYkXQURKf3CuBIbIGjfYHDCLewkDhqn3WlOhcjR0IkAu+/n9LT\n+EK/C5dy4R8v62GMYMc1KX4E5gGlvkckkqbq1YPrroNhw+CddzRliOQeL1WTicDBwBTn+OOxO6X2\nAZ7F5pcKmpqkJBRbt8Khh8L559v4DJFMEkQfxr+xDu6fnO16wKvAMVgto20C530Sa94CyAfWYh3q\nAKOBgdiCTZcA08v5vAoMCc38+XbX1OLFNhJcJFME0YexO+4YDLBbav8PWyNjU4LnPR0rIA4CnnMe\nYGt693Oej8HW31DFvxJqn3UFlYsOHWxCwjFjAjldQvS7cCkX/vHSh/EENg7jRaxkOhFrptoVW3s7\nGXlAX6C7s90bmIQVSquA5UBnYHaS5xHx1fXXQ9u2cO65NqutSC7wWjXphN0dFcXWxZjr0/kPB25z\nvh/gLqxwiN24+BA2DuS5Mp9Tk5SE7rHH4M47Yc4cqFkz7GhEqpZsk5SXGgbAB8AXwM5YodHU2a7M\nDKBxOfvHYB3oYFOMVDUXaLklQ3FxMQUFBQDk5+dTWFhIUVER4FZBta3tVG73719ESQlcdFGEfv3C\nj0fb2i67HYlEKCkpAdh2vUyGl5KmF1YL2BP4Frs76hNgvyTPXQtYDXQAvnb2jXKeY5MwTAPGsv3U\nJKAaxjaRSGTbDyXXhZGLlSuhc2d47z1o3TrQU1dKvwuXcuEKotP7emxN76VAM2wQX9kLeCKOwgqe\nr+P2TcY6xGs752oFZMAwKclVzZvbCPCBAzNzoSWR6vBS0swDOgIfYbWB34CPsWnOk/EIMAubdiTe\nGOy22l+xtcNfK+ezqmFI2ti6FYqK4JRTYOjQsKMRqVgQ4zBeB04C/oZND/ItNpDv0ERP6gMVGJJW\nli2Drl1h9mxo2TLsaETKF0STVG9szMWlWJ/CcuzWWkkDsQ4uCTcXrVrZ+t9nngmbN1d9fKrpd+FS\nLvzjpcC4GmuG2gKUAHcCl6cwJpGMdPHFsMceVnCIZCMvVZMPcaftiFkAtPc/HM/UJCVp6bvvoLAQ\nHn0Ujjoq7GhEtpfKJqlBWMGwr/Mce6zCOr1FpIzdd7cBfQMGWOEhkk0qKzAmYn0Vk4ETnNcnYndM\nnZX60MQLtc+60iUXRx4J55wD/fuHd6ttuuQiHSgX/qmswKgJrAOGAOud1+uwkdeao1OkEtdeC5s2\nqT9DskvN9W9uAAAPlklEQVRlbVmrqGBaDmd/c9+j8U59GJL2vv3WJia8/XYboyEStiDGYaQjFRiS\nEebNg2OOsaVd90t2Mh2RJAUxDgNsLMZtwK1oDEZaUfusKx1z0bEj3HYb9OkD338f3HnTMRdhUS78\n46XAGI+tfLcIm/vpEmzUt4h4cM45cNJJ0Lu39WuIZCovVZMFQCE2eA+sM7wUjcMQ8WzrVjjrLNiy\nBZ56SutnSDiCaJKKYutux+RTcWe4iJSjRg0oKYH//Q+GDwf9vSOZyEuB8TdgPvCo85gH3JjKoMQ7\ntc+60j0XderAiy/CG2/A+PFVH5+MdM9FkJQL/1S24t4EbPDeJGAmtoxqFFvk6JvUhyaSffLzYfp0\n6NbNCpDhw8OOSMS7ytqyhgH9sJX2nsIKjg+DCMoD9WFIRvvyS1tDY+hQuOSSsKORXBHEOIwCbBW8\nfsAuuLWOpYme1AcqMCTjff65FRojR8LgwWFHI7kgiE7vVdittQdhBcdJ2O21yeiMLb36IfAB1twV\nMxpYBnwK9EzyPFlP7bOuTMvFPvtYf8att8KNN/rbEZ5puUgl5cI/XgqMWkAvrGYxDbuQn5zkeW8G\nrsIKoaudbYB2WE2mHXAM1o/idXChSMZp3hzeeQcmTbL+jK1bw45IpGKVVU16YjWK47HawCRs5tqf\nfDjvJOAF4GngDOccZ2O1i63ATc5x04BxwOwyn1eTlGSVH36AE0+EZs3gn/+E2rXDjkiyUSqbpEYB\ns4C22HQgE/GnsIh9923AF8AtWEEB1sG+Ou641cBePp1TJG01aGB3T61bBz172ngNkXRT2W21RyT5\n3TOAxuXsvxKbXuQSrJZxGvAw0KOC7ym3KlFcXExBQQEA+fn5FBYWUlRUBLhtlrmwHd8+mw7xhLkd\n25cu8SSy/cIL0L9/hPbtYfr0Itq3T+z7SktLGTZsWOj/nnTYvuOOO3L6+lBSUgKw7XqZjLBmq10H\n7BYXw1rgd1jNA6yTHaxJaiwwp8zn1STliEQi234ouS6bcjFxot1ye999iU2Nnk25SJZy4crU6c3n\nA5diAwKPxAqITlhn90TsLqq9gNeBluxYy1CBIVlv7lzo18+mR7/tNth557AjkkwX1PTmfjsfuzOq\nFLje2QZYjHWELwamAoPRvFWSow4+GObPt7XBu3SBJUvCjkhyXVgFxlzgEGwW3K5sP4L8RqxW0QZ4\nLfjQMkt8+32uy8Zc/O53NrvtoEFw2GHwj394u/U2G3ORKOXCPxrjIJLm8vLgggvgvffg2Wfh8MNh\naZjzLEjO0hKtIhlk61a4+2649lrrFB85Un0b4l2m9mGISAJq1LDJCufOhdJSaNcOJk/W+hoSDBUY\nGU7ts65cykVBATz3HNx/P1xxBfToAe+/776fS7moinLhHxUYIhmsRw/4+GPo2xdOPtnWDl+4MOyo\nJFupD0MkS/z8M9x7L9x0E3TvDpdfDh06hB2VpBP1YYgIAHXr2oy3y5dDp07QuzcceSRMm6Y+DvGH\nCowMp/ZZl3Jh6teHjh0jrFwJxcVW09hvP7jjDvj++7CjC55+F/5RgSGSpXbaCfr3h48+ss7xefOg\nRQs4+2x480347bewI5RMoz4MkRyyZg089pg9vvnGOstPP92mHsnL1KuBeJapkw8mSwWGSJKWLLFp\nRyZNsg7z3r3h+OOhWzeoUyfs6CQV1Omd49Q+61IuXF5yse++cPXVsHixDf7bYw8YN86e+/SBBx6w\nDvRM/9tMvwv/qMAQyXF5eXDAAXDllTZf1YoVcOqpMHOm1Tb23hvOPNP6QZYsyfwCRBKnJikRqVA0\nagXIzJnu46efbOr1Tp3c5z33VB9IJlAfhogE6ptvbC6ruXPhgw/sUaMG7L+/3b4b/2jQIOxoJV6m\nFhgHAvcBuwKrgLOA9c57o4GBwG/Yut/Ty/m8CgyHlp90KReuIHMRjcJXX8GiRe5j4ULrG6lXz27l\njT2aN3df7757MLUS/S5cyRYYtfwLpVoeAoYD/wb+DIwErsaWaO3nPMeWaG0NeFgyJjeVlpbqfwaH\ncuEKMhd5edbPsffecPTR7v6tW60gWbnSmrVWrICXX3Zf//wz7LVXxY/Gja1QqVcvuYJFvwv/hFVg\ntMIKC7BCYRpWYPQGJgFbsJrHcmx979nBh5gZ1q5dG3YIaUO5cKVDLmrUgCZN7NGt247v//STFSjx\nj5Ur4d//ttf/+Q/873+wZQs0amSFR6NG27/Oz7dVCXfbzX2Of123bnrkIluEdZfUIqxwADgNaOK8\n3hNYHXfcaqymkXLVufXOy7EVHeN1f2Xbqb5NULmo+NzJHludXHjZl8m5mDs3wr77whFH2Ij0UaPg\nzjvh4osjzJ4Nq1ZZobJ2LdxxR4RHHrGp3Hv3htatYcWKCCtWQCQCjz8OI0dGGDQITjwR2raN8Pvf\nQ+3aNhljs2bWx9KpkxVexx4Lp5xio97PPx+GDYMxY+C66+C222wSx5ISmDjRVjmcPNnm5HrzTbjr\nrgjvv2/rkSxaBMuWweefW9/O99/D+vXwyy/uUrrZ9LtIZQ1jBtC4nP1jsD6KO4GrgMnA5kq+J5DO\niuq0c3o5tqJjvO6vbDv+9apVqzzFXB3KRdUxJnpsdXLhZV8u5KJuXVi0KMJpp21/7LhxEcaNKyp3\nO/Z682Y455xV3HCDNYFt3Og+Ktr+4Qf39ZYtsHnz9o/lyyM0bly0w/7YI/4z1pQWoXbtImrWZIdH\nrVr2vG5dhEaNdjzmv/+N0LTp9vs//zxCy5ZF5OVZLa5GDVi2LEKbNu6+JUsitGtXRI0asHhxhPbt\n7XWy0uEuqdbA48AhwChn33jneRowFphT5jPLgRaBRCcikj1WAC3DDqK6dneeawCPAcXOdjugFKgN\nNMP+celQqImISEguAZY4jxvLvDcGq0F8ChyNiIiIiIiIiIiIiIjkopphB+CTZsCtQH/gmZBjCVtv\nYAQ23cqPwMpwwwlVG+A6YACwG/BhuOGEblfgPeBrYFnIsYSpCLszszPwE/B5qNGEKw+4AeiD3Yz0\nUWUHZ8v05p8B54UdRJp4CTgfuBCbZiWXfQoMAk5HN1AAXA48FXYQaWArNnddHbYfKJyL+mCDozeT\ng7nI9dpFvFuBwrCDSAMnAlOBk8MOJGQ9sD8gBgDHhxxL2GK36u8B/CvMQNLAFcBfnNdVXj/TuYbx\nMPBfYEGZ/cdgfzkuw/6xuaA6ucgDbsIukqVBBRig6v4upgDHYhfKbFOdXHQDugBnYheIbBvfVJ1c\nxGaPWIvVMrJNdXKxGssDZPgkr38CDmL7f3RNbIxGAbATdkFsCzTEpkvP1kKkOrm4GJgL3AtcEGiU\nwahOLroB/wDuB4YFGmUwqpOLmAHAcQHFF6Tq5OIk7HrxJHB4oFEGozq5qIvNHn4n1nyb0QrY/h/d\nFZsuJGYU7nQi2a4A5SKmAOUipgDlIqYA5SKmgBTkIp2bpMqzF/Bl3HZgs9mmIeXCpVy4lAuXcuHy\nJReZVmBomT2XcuFSLlzKhUu5cPmSi0wrML7CXTsD53XO3QrmUC5cyoVLuXApF66cyEUB27fD1cJm\nsC3AZrQt26GXzQpQLmIKUC5iClAuYgpQLmIKyLFcTMJGpP6Ctb392dl/LDbL7XJgdDihBU65cCkX\nLuXCpVy4lAsRERERERERERERERERERERERERERERERGRLNYHm2N/3xSe46ckP/8nYBEwn+xaG+EE\nYFyS3xEBOsZtj8LW0yjPAcA/kzyfiOSwp4DJJH/hqsz6ahybx44LB92HrX9eVq2EI0oPbwH/V87+\nmtX8jg5x228Cv6/k+Ai2ip2ISLXUA1YBTYFP4vYXYReWZ5z98UtkHufsm4st6DLF2T8OGBF33ELn\ne8EtMOoBrwPzgI+BXs7+Amw6hEedz8VPwHYe8D2w0omjG/BvbF30T7HJOW8B3gc+wtZKByt07naO\nmQG8ApzivLcKW9AL4GDsoguwK7YS2hysNhOLrxh4HlshcSm2WmLMMc6/p9Q5T55zTCPn/RrYomFl\nL+JNgHfjtkuwgnE2tnRvJ+A9J453gdbOcXWxRYUWOzHNxq1h7Aa847w+DZunqBSYGXeeK4AhiIhU\n01nYRQrgbdy/VIuwJSH3xC6A7wGHAjsDXwD7OMdNxGonAGPZvsBYwI4FRk2gvvO6EXYhBSswfgM6\nVxDnI7jrfRdhTVyxGM4HrnRe1wE+cL7vZGC6E/8fgB/ivuMzyi8wbsStyeRjhdguWIGxwom9Dlbg\n7AXszvb5yHeerwaGOq97Uv46zKcDd5X5N07GrV3Vx61pHAU867wejq2+BtAe2IL73+1k3Jrix86/\nG6wgiemO1SolS2V6tVvS1xnA353Xzzjb853t97HJ0cD+Sm0GbMT+0v/c2T8J9y96L2oAf8P6JLZi\nBVKseeRz55wViW+mej8uhp7YhfNUZ3s3oJVzjonYGgPfYE01VekJnAhc5mzXwQq9KPAGbsG3GCuU\nGmIFbSyW2LrLD2M1oH8AA7HCoKymTlzxnsFdEyEfeAxo6eyLXQf+5HwvWKH8cdznj3bODVYreRR4\nGquJxHzjxC5ZSgWGpEJD7K/N/bELUk3neaTz/i9xx/6G/Q7LLvASfxH/le3Xbtm5nHOehdUsOjjf\n+VnccRuqiDf+3GWPvQhrDop3XJn4Koq1bJwn49Z8Yg7BWz5iVgP/BY7AmpbOqOC4sn01G+NeX4cV\nUidhF/i34t4r+7mYzsCFzutBzvbxWJNZR2CN81ktWpTFMm0BJckMp2J/wRZgtYem2AX8TxUcH8Wa\naJrjNsH0w734rMJtGungfGdZuwHfYhfc7nHf40VFF8nXgMG4f1i1xpqR3nbiq4E1zRTFfWYV1hQF\nbr9G7Lsuids+qJJzR7H+g8Nx/2JvGPf+Q1ify9OUf4H+HGhc3j/IsRtuDa84bv/buHdB7Y/d+QSw\nH9ZfEztXC6wmNhb4Dtjb2f8H3BqRZCEVGJIKpwMvlNn3HPbXcJTyL3KbsIvzNKzTe53ziH22IdZp\nPQQrXGJi3/UEdqH+GOjP9h3tVf3VG417jj/2IayJaD7WRHMvVlt6AaspLMaaZmbhXvivwZp1PsBq\nG7Hvuw7YyYlvoXNceeeM+R/WJPc81mw3Ke69KVgnennNUWBNRh3K7Is/x81Y89183Nofzr+vnvPv\nugb775CHraMwtcznP8Zy8i5u01VnrNAREUm5XeNe34PbuZvuHmH72kSqHcz2dyeV503cjulkTaf8\nW3TLiqDbakUkIMOAD7GBdI9Tfl9FOoq/0yrVRmHNXodWcdxxuLWYIByAe4eViIiIiIiIiIiIiIiI\niIiIiIiIiIiISLb6f5zIgihRf0DAAAAAAElFTkSuQmCC\n",
255       "text": [
256        "<matplotlib.figure.Figure at 0x7f637db2db10>"
257       ]
258      }
259     ],
260     "prompt_number": 41
261    },
262    {
263     "cell_type": "code",
264     "collapsed": false,
265     "input": [],
266     "language": "python",
267     "metadata": {},
268     "outputs": [],
269     "prompt_number": 41
270    }
271   ],
272   "metadata": {}
273  }
274 ]
275}