')) 68 | .append( 69 | textarea 70 | ) 71 | ) 72 | ); 73 | 74 | IPython.dialog.modal({ 75 | title: "Please, talk and when ready click the button", 76 | body: dialogform, 77 | buttons: { 78 | "OK": { class : "btn-primary", 79 | click: function() { 80 | talkstop(); 81 | var corr_input = $.trim($(textarea).val()); 82 | console.log(corr_input); 83 | IPython.notebook.get_selected_cell().set_text(corr_input); 84 | } 85 | }, 86 | Cancel: { click: function() { 87 | talkstop(); 88 | } 89 | } 90 | } 91 | }); 92 | 93 | var talkstart = function (){ 94 | console.log('entrostart'); 95 | // Start recognising 96 | recognizer.onresult = function(event) { 97 | console.log('onresult'); 98 | var transcription = $('#talkarea'); 99 | transcription.val(''); 100 | for (var i = event.resultIndex; i < event.results.length; i++) { 101 | if (event.results[i].isFinal) { 102 | transcription.val( 103 | transcription.val() + event.results[i][0].transcript + '\n(Confidence: ' + event.results[i][0].confidence + ')' 104 | ); 105 | } else { 106 | transcription.val( 107 | transcription.val() + event.results[i][0].transcript 108 | ); 109 | }; 110 | }; 111 | }; 112 | // Error on recognising 113 | recognizer.onerror = function(event) { 114 | console.log('onerror') 115 | transcription.val('Recognition error: ' + event.message); 116 | }; 117 | // Start listening 118 | recognizer.start(); 119 | }; 120 | 121 | var talkstop = function () { 122 | // Stop listening 123 | console.log('entrostop') 124 | recognizer.stop(); 125 | }; 126 | 127 | talkstart(); 128 | 129 | }; 130 | 131 | } else { 132 | 133 | var talk_func = function(){ 134 | 135 | var dialogupdate = $('') 136 | .append( 137 | $('') 138 | .text("Sorry, it seems that your browser doesn't support this functionality") 139 | ); 140 | 141 | IPython.dialog.modal({ 142 | title: "Not recognised", 143 | body: dialogupdate, 144 | }); 145 | 146 | }; 147 | 148 | }; 149 | 150 | var talk_button = function () { 151 | if (!IPython.toolbar) { 152 | $([IPython.events]).on("app_initialized.NotebookApp", talk_button); 153 | return; 154 | } 155 | if ($("#talk_button").length === 0) { 156 | IPython.toolbar.add_buttons_group([ 157 | { 158 | 'label' : 'Talk to IPython', 159 | 'icon' : 'fa-pencil-square-o', 160 | 'callback': toggle_talk, 161 | 'id' : 'talk_button' 162 | }, 163 | ]); 164 | }; 165 | }; 166 | 167 | var toggle_talk = function(){ 168 | talk_func() 169 | }; 170 | 171 | var load_css = function () { 172 | var link = document.createElement("link"); 173 | link.type = "text/css"; 174 | link.rel = "stylesheet"; 175 | link.href = require.toUrl("./talk.css"); 176 | console.log(link); 177 | document.getElementsByTagName("head")[0].appendChild(link); 178 | }; 179 | 180 | var load_ipython_extension = function () { 181 | load_css(); 182 | talk_button(); 183 | }; 184 | 185 | return { 186 | load_ipython_extension : load_ipython_extension, 187 | toggle_talk : toggle_talk, 188 | talk_func : talk_func, 189 | 190 | }; 191 | 192 | }); 193 | -------------------------------------------------------------------------------- /IPython-201501-pythonmadrid-kikocorreoso/03-Extensiones_JS.ipynb: -------------------------------------------------------------------------------- 1 | { 2 | "metadata": { 3 | "name": "", 4 | "signature": "sha256:164f84e40d61b1915a7bd08c82774e35dd7b020c5f40c8be5e042865cf64e40f" 5 | }, 6 | "nbformat": 3, 7 | "nbformat_minor": 0, 8 | "worksheets": [ 9 | { 10 | "cells": [ 11 | { 12 | "cell_type": "heading", 13 | "level": 1, 14 | "metadata": {}, 15 | "source": [ 16 | "Creando una extensi\u00f3n js simple" 17 | ] 18 | }, 19 | { 20 | "cell_type": "markdown", 21 | "metadata": {}, 22 | "source": [ 23 | "
\n",
24 | "La creaci\u00f3n de extensiones js no est\u00e1 bien cubierta en [los documentos oficiales de IPython](http://ipython.org/ipython-doc/dev/config/extensions/) (o quiz\u00e1 no he sido capaz de encontrar la informaci\u00f3n \u00fatil) y te toca navegar para encontrar informaci\u00f3n \u00fatil.\n",
25 | "
"
26 | ]
27 | },
28 | {
29 | "cell_type": "markdown",
30 | "metadata": {},
31 | "source": [
32 | "Vamos a crear una extensi\u00f3n que nos permitir\u00e1 hablarle al notebook de IPython y el texto se escribir\u00e1 en una celda del mismo. Ha sido probado con IPython 2.3.0 y 3.0.0.dev en chrome/chromium v40 (necesitar\u00e1s que la [web speech API funcione en tu navegador](https://dvcs.w3.org/hg/speech-api/raw-file/tip/speechapi.html))."
33 | ]
34 | },
35 | {
36 | "cell_type": "markdown",
37 | "metadata": {},
38 | "source": [
39 | "Para hacerlo funcionar debes crear la carpeta 'talk' en ***IPYTHONDIR/nbextensions*** y copiar los ficheros que se encuentran [aqu\u00ed](-------------------------------) en la carpeta 'talk' que acabas de crear. "
40 | ]
41 | },
42 | {
43 | "cell_type": "markdown",
44 | "metadata": {},
45 | "source": [
46 | "Ahora estamos listos para cargar la extensi\u00f3n:"
47 | ]
48 | },
49 | {
50 | "cell_type": "code",
51 | "collapsed": false,
52 | "input": [
53 | "%%javascript\n",
54 | "IPython.load_extensions('talk/talk');"
55 | ],
56 | "language": "python",
57 | "metadata": {},
58 | "outputs": [
59 | {
60 | "javascript": [
61 | "IPython.load_extensions('talk/talk');"
62 | ],
63 | "metadata": {},
64 | "output_type": "display_data",
65 | "text": [
66 | "'))\n", 173 | " .append(\n", 174 | " textarea\n", 175 | " )\n", 176 | " )\n", 177 | " );\n", 178 | " \n", 179 | " IPython.dialog.modal({\n", 180 | " title: \"Please, talk and when ready click the button\",\n", 181 | " body: dialogform,\n", 182 | " buttons: {\n", 183 | " \"OK\": { class : \"btn-primary\",\n", 184 | " click: function() {\n", 185 | " talkstop();\n", 186 | " var corr_input = $.trim($(textarea).val());\n", 187 | " console.log(corr_input);\n", 188 | " IPython.notebook.get_selected_cell().set_text(corr_input);\n", 189 | " }\n", 190 | " },\n", 191 | " Cancel: { click: function() {\n", 192 | " talkstop();\n", 193 | " }\n", 194 | " }\n", 195 | " }\n", 196 | " });\n", 197 | " \n", 198 | " var talkstart = function (){\n", 199 | " console.log('entrostart');\n", 200 | " // Start recognising\n", 201 | " recognizer.onresult = function(event) {\n", 202 | " console.log('onresult');\n", 203 | " var transcription = $('#talkarea');\n", 204 | " transcription.val('');\n", 205 | " for (var i = event.resultIndex; i < event.results.length; i++) {\n", 206 | " if (event.results[i].isFinal) {\n", 207 | " transcription.val(\n", 208 | " transcription.val() + event.results[i][0].transcript + '\\n(Confidence: ' + event.results[i][0].confidence + ')'\n", 209 | " );\n", 210 | " } else {\n", 211 | " transcription.val(\n", 212 | " transcription.val() + event.results[i][0].transcript\n", 213 | " );\n", 214 | " };\n", 215 | " };\n", 216 | " };\n", 217 | " // Error on recognising\n", 218 | " recognizer.onerror = function(event) {\n", 219 | " console.log('onerror')\n", 220 | " transcription.val('Recognition error: ' + event.message);\n", 221 | " };\n", 222 | " // Start listening\n", 223 | " recognizer.start();\n", 224 | " };\n", 225 | " \n", 226 | " var talkstop = function () {\n", 227 | " // Stop listening\n", 228 | " console.log('entrostop')\n", 229 | " recognizer.stop();\n", 230 | " };\n", 231 | " \n", 232 | " talkstart();\n", 233 | " \n", 234 | " };\n", 235 | " \n", 236 | " } else {\n", 237 | " \n", 238 | " var talk_func = function(){\n", 239 | " \n", 240 | " var dialogupdate = $('')\n", 241 | " .append(\n", 242 | " $('')\n", 243 | " .text(\"Sorry, it seems that your browser doesn't support this functionality\")\n", 244 | " );\n", 245 | " \n", 246 | " IPython.dialog.modal({\n", 247 | " title: \"Not recognised\",\n", 248 | " body: dialogupdate,\n", 249 | " });\n", 250 | " \n", 251 | " }; \n", 252 | " \n", 253 | " };\n", 254 | " \n", 255 | " var talk_button = function () {\n", 256 | " if (!IPython.toolbar) {\n", 257 | " $([IPython.events]).on(\"app_initialized.NotebookApp\", talk_button);\n", 258 | " return;\n", 259 | " }\n", 260 | " if ($(\"#talk_button\").length === 0) {\n", 261 | " IPython.toolbar.add_buttons_group([\n", 262 | " {\n", 263 | " 'label' : 'Talk to IPython',\n", 264 | " 'icon' : 'fa-pencil-square-o',\n", 265 | " 'callback': toggle_talk,\n", 266 | " 'id' : 'talk_button'\n", 267 | " },\n", 268 | " ]);\n", 269 | " };\n", 270 | " };\n", 271 | " \n", 272 | " var toggle_talk = function(){\n", 273 | " talk_func()\n", 274 | " };\n", 275 | " \n", 276 | " var load_css = function () {\n", 277 | " var link = document.createElement(\"link\");\n", 278 | " link.type = \"text/css\";\n", 279 | " link.rel = \"stylesheet\";\n", 280 | " link.href = require.toUrl(\"./talk.css\");\n", 281 | " console.log(link);\n", 282 | " document.getElementsByTagName(\"head\")[0].appendChild(link);\n", 283 | " };\n", 284 | " \n", 285 | " var load_ipython_extension = function () {\n", 286 | " load_css();\n", 287 | " talk_button();\n", 288 | " };\n", 289 | "\n", 290 | " return {\n", 291 | " load_ipython_extension : load_ipython_extension,\n", 292 | " toggle_talk : toggle_talk,\n", 293 | " talk_func : talk_func,\n", 294 | " \n", 295 | " };\n", 296 | " \n", 297 | "});" 298 | ], 299 | "language": "python", 300 | "metadata": {}, 301 | "outputs": [] 302 | }, 303 | { 304 | "cell_type": "markdown", 305 | "metadata": {}, 306 | "source": [ 307 | "Si deseas que la extensi\u00f3n est\u00e9 disponible siempre que abras un nuevo notebook deber\u00e1s a\u00f1adir algo de c\u00f3digo en ***custom.js*** del perfil de usuario que quieras usar. Por ejemplo, en los notebooks previos hemos creado el perfil pymad por lo que para incluirlo en el ***custom.js*** de ese perfil deberemos incluir el siguiente c\u00f3digo:" 308 | ] 309 | }, 310 | { 311 | "cell_type": "markdown", 312 | "metadata": {}, 313 | "source": [ 314 | " // we want strict javascript that fails\n", 315 | " // on ambiguous syntax\n", 316 | " \"using strict\";\n", 317 | "\n", 318 | " $([IPython.events]).on('app_initialized.NotebookApp', function(){\n", 319 | " require(['custom/talk/main'],function(talk){\n", 320 | " talk.parameters('dummy');\n", 321 | " console.log('Talk extension loaded correctly');\n", 322 | " });\n", 323 | " });" 324 | ] 325 | }, 326 | { 327 | "cell_type": "heading", 328 | "level": 1, 329 | "metadata": {}, 330 | "source": [ 331 | "M\u00e1s info" 332 | ] 333 | }, 334 | { 335 | "cell_type": "markdown", 336 | "metadata": {}, 337 | "source": [ 338 | "http://nbviewer.ipython.org/github/ipython-books/cookbook-code/blob/master/notebooks/chapter03_notebook/03_controls.ipynb\n", 339 | "\n", 340 | "https://github.com/ipython/ipython/wiki/Extensions-Index\n", 341 | "\n", 342 | "https://github.com/ipython-contrib/IPython-notebook-extensions\n", 343 | "\n", 344 | "http://www.damian.oquanta.info/posts/loader-and-writer-ipython-nbextensions-to-easily-edit-your-text.html\n", 345 | "\n", 346 | "http://www.damian.oquanta.info/posts/a-poor-man-spell-checker-for-the-ipython-notebook.html\n", 347 | "\n", 348 | "http://www.damian.oquanta.info/posts/tweet-me-said-an-ipython-notebook-cell-and-i-did-it.html\n", 349 | "\n", 350 | "https://github.com/minrk/ipython_extensions\n", 351 | "\n", 352 | "https://www.pfenninger.org/posts/ipython-notebook-extensions-to-ease-day-to-day-work/" 353 | ] 354 | } 355 | ], 356 | "metadata": {} 357 | } 358 | ] 359 | } -------------------------------------------------------------------------------- /IPython-201501-pythonmadrid-kikocorreoso/02-Manipulando_con_js.ipynb: -------------------------------------------------------------------------------- 1 | { 2 | "metadata": { 3 | "name": "", 4 | "signature": "sha256:1ddb191faaab8c9365bd3b6183f165b3448d26a1c9fb33bf4b0c3e05dc47be7c" 5 | }, 6 | "nbformat": 3, 7 | "nbformat_minor": 0, 8 | "worksheets": [ 9 | { 10 | "cells": [ 11 | { 12 | "cell_type": "heading", 13 | "level": 1, 14 | "metadata": {}, 15 | "source": [ 16 | "Notebook javascript API" 17 | ] 18 | }, 19 | { 20 | "cell_type": "markdown", 21 | "metadata": {}, 22 | "source": [ 23 | "El notebook dispone de una API javascript 'pyth\u00f3nica'. Desde esta API puedes acceder a la configuraci\u00f3n de code mirror, puedes a\u00f1adir botones a la barra de herramientas, habilitar o deshabilitar ciertas celdas de c\u00f3digo,..." 24 | ] 25 | }, 26 | { 27 | "cell_type": "heading", 28 | "level": 2, 29 | "metadata": {}, 30 | "source": [ 31 | "Algunos ejemplos" 32 | ] 33 | }, 34 | { 35 | "cell_type": "markdown", 36 | "metadata": {}, 37 | "source": [ 38 | "--------------------------------------" 39 | ] 40 | }, 41 | { 42 | "cell_type": "heading", 43 | "level": 3, 44 | "metadata": {}, 45 | "source": [ 46 | "A\u00f1adir n\u00fameros de l\u00ednea en las celdas" 47 | ] 48 | }, 49 | { 50 | "cell_type": "markdown", 51 | "metadata": {}, 52 | "source": [ 53 | "A\u00f1adir n\u00fameros de l\u00ednea a todas las celdas nuevas:\n", 54 | "\n", 55 | " IPython.Cell.options_default.cm_config.lineNumbers = true;\n", 56 | "\n", 57 | "o solo a las celdas de ***codigo***:\n", 58 | "\n", 59 | " IPython.CodeCell.options_default.cm_config.lineNumbers = true;" 60 | ] 61 | }, 62 | { 63 | "cell_type": "code", 64 | "collapsed": false, 65 | "input": [ 66 | "%%javascript\n", 67 | "IPython.CodeCell.options_default.cm_config.lineNumbers = false;" 68 | ], 69 | "language": "python", 70 | "metadata": {}, 71 | "outputs": [ 72 | { 73 | "javascript": [ 74 | "IPython.CodeCell.options_default.cm_config.lineNumbers = false;" 75 | ], 76 | "metadata": {}, 77 | "output_type": "display_data", 78 | "text": [ 79 | "
\n",
145 | "Esta funcionalidad la he probado solo en versiones modernas de Chrome. En firefox v35 no lo he conseguido hacer funcionar.\n",
146 | "
"
147 | ]
148 | },
149 | {
150 | "cell_type": "markdown",
151 | "metadata": {},
152 | "source": [
153 | "Se pueden a\u00f1adir atajos de teclado de forma muy sencilla usando la js IPython API. "
154 | ]
155 | },
156 | {
157 | "cell_type": "code",
158 | "collapsed": false,
159 | "input": [
160 | "%%javascript\n",
161 | "IPython.keyboard_manager.command_shortcuts.add_shortcut(\"Ctrl-Shift-Alt-h\", {\n",
162 | " help : 'Talking cells',\n",
163 | " help_index : 'zzz',\n",
164 | " handler : function() { \n",
165 | " var input = IPython.notebook.get_selected_cell().get_text();\n",
166 | " var msg = new SpeechSynthesisUtterance();\n",
167 | " msg.text = input;\n",
168 | " msg.lang = \"es-ES\";\n",
169 | " msg.voiceURI = 'native';\n",
170 | " msg.volume = 1; // 0 to 1\n",
171 | " msg.rate = 0.5; // 0.1 to 10\n",
172 | " msg.pitch = 2; //0 to 2\n",
173 | " window.speechSynthesis.speak(msg);\n",
174 | " return false;\n",
175 | " }}\n",
176 | ");"
177 | ],
178 | "language": "python",
179 | "metadata": {},
180 | "outputs": [
181 | {
182 | "javascript": [
183 | "IPython.keyboard_manager.command_shortcuts.add_shortcut(\"Ctrl-Shift-Alt-h\", {\n",
184 | " help : 'Talking cells',\n",
185 | " help_index : 'zzz',\n",
186 | " handler : function() { \n",
187 | " var input = IPython.notebook.get_selected_cell().get_text();\n",
188 | " var msg = new SpeechSynthesisUtterance();\n",
189 | " msg.text = input;\n",
190 | " msg.lang = \"es-ES\";\n",
191 | " msg.voiceURI = 'native';\n",
192 | " msg.volume = 1; // 0 to 1\n",
193 | " msg.rate = 0.5; // 0.1 to 10\n",
194 | " msg.pitch = 2; //0 to 2\n",
195 | " window.speechSynthesis.speak(msg);\n",
196 | " return false;\n",
197 | " }}\n",
198 | ");"
199 | ],
200 | "metadata": {},
201 | "output_type": "display_data",
202 | "text": [
203 | "En el notebook existen dos modos: ***edit*** y ***command***. Se pueden a\u00f1adir atajos de teclado para trabajar tanto en un modo como en el otro por lo que si un atajo se ha creado para funcionar en un modo en el otro modo no funcionar\u00e1 a no ser que tambi\u00e9n lo creemos de forma expl\u00edcita.
"
355 | ]
356 | },
357 | {
358 | "cell_type": "code",
359 | "collapsed": false,
360 | "input": [
361 | "%%javascript\n",
362 | "\n",
363 | "IPython.keyboard_manager.command_shortcuts.add_shortcut('Ctrl-Alt-c', {\n",
364 | " help : 'clear all outputs',\n",
365 | " help_index : 'zzz',\n",
366 | " handler : function (event) {\n",
367 | " IPython.notebook.clear_all_output(); \n",
368 | " return false;\n",
369 | " }}\n",
370 | ");\n",
371 | "IPython.keyboard_manager.edit_shortcuts.add_shortcut('Ctrl-Alt-c', {\n",
372 | " help : 'clear all outputs',\n",
373 | " help_index : 'zzz',\n",
374 | " handler : function (event) {\n",
375 | " IPython.notebook.clear_all_output(); \n",
376 | " return false;\n",
377 | " }}\n",
378 | ");"
379 | ],
380 | "language": "python",
381 | "metadata": {},
382 | "outputs": []
383 | },
384 | {
385 | "cell_type": "code",
386 | "collapsed": false,
387 | "input": [
388 | "%%javascript\n",
389 | "\n",
390 | "IPython.keyboard_manager.command_shortcuts.remove_shortcut('Ctrl-Alt-c');\n",
391 | "IPython.keyboard_manager.edit_shortcuts.remove_shortcut('Ctrl-Alt-c');"
392 | ],
393 | "language": "python",
394 | "metadata": {},
395 | "outputs": []
396 | },
397 | {
398 | "cell_type": "heading",
399 | "level": 4,
400 | "metadata": {},
401 | "source": [
402 | "M\u00e1s info sobre atajos de teclado"
403 | ]
404 | },
405 | {
406 | "cell_type": "markdown",
407 | "metadata": {},
408 | "source": [
409 | "http://nbviewer.ipython.org/github/ipython/ipython/blob/2.x/examples/Notebook/User%20Interface.ipynb#Keyboard-shortcut-customization"
410 | ]
411 | },
412 | {
413 | "cell_type": "markdown",
414 | "metadata": {},
415 | "source": [
416 | "------------------------------"
417 | ]
418 | },
419 | {
420 | "cell_type": "heading",
421 | "level": 1,
422 | "metadata": {},
423 | "source": [
424 | "Custom.js y Custom.css"
425 | ]
426 | },
427 | {
428 | "cell_type": "markdown",
429 | "metadata": {},
430 | "source": [
431 | "Si deseas que algunas de las modificaciones que has hecho est\u00e9s disponibles cuando vuelvas a abrir un notebook puedes especificarlo en el fichero ***custom.js*** disponible en cada uno de los perfiles de usuario. Tambi\u00e9n, si quieres modificar la apariencia, puedes usar el fichero ***custom.css***"
432 | ]
433 | },
434 | {
435 | "cell_type": "code",
436 | "collapsed": false,
437 | "input": [
438 | "!ipython profile create pymad\n",
439 | "profile_dir = !ipython locate profile pymad\n",
440 | "profile_dir = profile_dir[0]\n",
441 | "profile_dir"
442 | ],
443 | "language": "python",
444 | "metadata": {},
445 | "outputs": [
446 | {
447 | "output_type": "stream",
448 | "stream": "stdout",
449 | "text": [
450 | "[ProfileCreate] Generating default config file: '/home/kiko/.ipython/profile_pymad/ipython_config.py'\n",
451 | "[ProfileCreate] Generating default config file: '/home/kiko/.ipython/profile_pymad/ipython_kernel_config.py'\n",
452 | "[ProfileCreate] Generating default config file: '/home/kiko/.ipython/profile_pymad/ipython_notebook_config.py'\n",
453 | "[ProfileCreate] Generating default config file: '/home/kiko/.ipython/profile_pymad/ipython_nbconvert_config.py'\n"
454 | ]
455 | },
456 | {
457 | "metadata": {},
458 | "output_type": "pyout",
459 | "prompt_number": 11,
460 | "text": [
461 | "'/home/kiko/.ipython/profile_pymad'"
462 | ]
463 | }
464 | ],
465 | "prompt_number": 11
466 | },
467 | {
468 | "cell_type": "markdown",
469 | "metadata": {},
470 | "source": [
471 | "Vamos a cargar el fichero custom.js del perfil *pymad*, El fichero cargado lo guardaremos usando `%%writefile` y con el fichero *custom.js* actualizado abriremos un nuevo notebook para ese perfil y veremos como se han aplicado los cambios."
472 | ]
473 | },
474 | {
475 | "cell_type": "code",
476 | "collapsed": false,
477 | "input": [
478 | "%load $profile_dir/static/custom/custom.js"
479 | ],
480 | "language": "python",
481 | "metadata": {},
482 | "outputs": [],
483 | "prompt_number": 12
484 | },
485 | {
486 | "cell_type": "code",
487 | "collapsed": false,
488 | "input": [
489 | "%%writefile $profile_dir/static/custom/custom.js\n",
490 | "// leave at least 2 line with only a star on it below, or doc generation fails\n",
491 | "/**\n",
492 | " *\n",
493 | " *\n",
494 | " * Placeholder for custom user javascript\n",
495 | " * mainly to be overridden in profile/static/custom/custom.js\n",
496 | " * This will always be an empty file in IPython\n",
497 | " *\n",
498 | " * User could add any javascript in the `profile/static/custom/custom.js` file\n",
499 | " * (and should create it if it does not exist).\n",
500 | " * It will be executed by the ipython notebook at load time.\n",
501 | " *\n",
502 | " * Same thing with `profile/static/custom/custom.css` to inject custom css into the notebook.\n",
503 | " * \n",
504 | " * Classes and functions are available at load time and may be accessed plainly: \n",
505 | " * \n",
506 | " * IPython.Cell.options_default.cm_config.extraKeys['Home'] = 'goLineLeft';\n",
507 | " * IPython.Cell.options_default.cm_config.extraKeys['End'] = 'goLineRight';\n",
508 | " * \n",
509 | " * Instances are created later however and must be accessed using events:\n",
510 | " * \n",
511 | " * $([IPython.events]).on(\"app_initialized.NotebookApp\", function () {\n",
512 | " * IPython.keyboard_manager....\n",
513 | " * });\n",
514 | " *\n",
515 | " * __Example 1:__\n",
516 | " *\n",
517 | " * Create a custom button in toolbar that execute `%qtconsole` in kernel\n",
518 | " * and hence open a qtconsole attached to the same kernel as the current notebook\n",
519 | " *\n",
520 | " * IPython.events.on('app_initialized.NotebookApp', function(){\n",
521 | " * IPython.toolbar.add_buttons_group([\n",
522 | " * {\n",
523 | " * 'label' : 'run qtconsole',\n",
524 | " * 'icon' : 'icon-terminal', // select your icon from http://fortawesome.github.io/Font-Awesome/icons\n",
525 | " * 'callback': function () {\n",
526 | " * IPython.notebook.kernel.execute('%qtconsole')\n",
527 | " * }\n",
528 | " * }\n",
529 | " * // add more button here if needed.\n",
530 | " * ]);\n",
531 | " * });\n",
532 | " *\n",
533 | " * __Example 2:__\n",
534 | " *\n",
535 | " * At the completion of the dashboard loading, load an unofficial javascript extension\n",
536 | " * that is installed in profile/static/custom/ \n",
537 | " *\n",
538 | " * IPython.events.on('app_initialized.DashboardApp', function(){\n",
539 | " * require(['custom/unofficial_extension.js'])\n",
540 | " * });\n",
541 | " *\n",
542 | " * __Example 3:__\n",
543 | " *\n",
544 | " * Use `jQuery.getScript(url [, success(script, textStatus, jqXHR)] );`\n",
545 | " * to load custom script into the notebook.\n",
546 | " *\n",
547 | " * // to load the metadata ui extension example.\n",
548 | " * $.getScript('/static/notebook/js/celltoolbarpresets/example.js');\n",
549 | " * // or\n",
550 | " * // to load the metadata ui extension to control slideshow mode / reveal js for nbconvert\n",
551 | " * $.getScript('/static/notebook/js/celltoolbarpresets/slideshow.js');\n",
552 | " *\n",
553 | " *\n",
554 | " * @module IPython\n",
555 | " * @namespace IPython\n",
556 | " * @class customjs\n",
557 | " * @static\n",
558 | " */\n",
559 | "\n",
560 | "$([IPython.events]).on(\"app_initialized.NotebookApp\", function () { \n",
561 | " \n",
562 | " var colors = [[\"Red\", \"PaleVioletRed\"],\n",
563 | " [\"Orange\", \"DarkOrange\"], \n",
564 | " [\"Green\", \"LightGreen\"], \n",
565 | " [\"None\", \"#F7F7F7\"]];\n",
566 | " \n",
567 | " var cell_flag_init = IPython.CellToolbar.utils.select_ui_generator(\n",
568 | " \n",
569 | " // List\n",
570 | " colors, \n",
571 | " // Setter, called when the select state has been changed. \n",
572 | " // Sets the color of the cell background\n",
573 | " // add the color value in the cell's metadata.\n",
574 | " function(cell, value){\n",
575 | " if (cell.metadata.bgcolor == undefined){\n",
576 | " cell.metadata.bgcolor = {};\n",
577 | " }\n",
578 | " cell.metadata.bgcolor = value;\n",
579 | " var prompt = cell.element.find('div.input_area');\n",
580 | " prompt.css(\"background-color\", value);\n",
581 | " },\n",
582 | " \n",
583 | " // Getter, called when the control is rendered. Determines the initial\n",
584 | " // state of the checkbox. If the flag doesn't in the metadata, default\n",
585 | " // to false.\n",
586 | " function(cell){\n",
587 | " if (cell.metadata.bgcolor === undefined || cell.metadata.bgcolor == false) {\n",
588 | " return \"default\";\n",
589 | " } else {\n",
590 | " return cell.metadata.bgcolor;\n",
591 | " }\n",
592 | " },\n",
593 | " // Label\n",
594 | " \"Choose bg color\"\n",
595 | " );\n",
596 | "\n",
597 | " IPython.CellToolbar.register_callback(\"choose_bg_color\", cell_flag_init);\n",
598 | " \n",
599 | " // Create and register the toolbar with IPython.\n",
600 | " IPython.CellToolbar.register_preset('Choose bg color', [\"choose_bg_color\"]);\n",
601 | " console.log(\"Background colors modifications loaded\")\n",
602 | " \n",
603 | "});\n",
604 | "\n",
605 | "$([IPython.events]).on('notebook_loaded.Notebook', function(){\n",
606 | " \n",
607 | " var cells = IPython.notebook.get_cells();\n",
608 | " for(var i in cells){\n",
609 | " var cell = cells[i];\n",
610 | " if (cell.metadata.bgcolor){\n",
611 | " var prompt = cell.element.find('div.input_area');\n",
612 | " prompt.css(\"background-color\", cell.metadata.bgcolor);\n",
613 | " }\n",
614 | " }\n",
615 | " \n",
616 | "});"
617 | ],
618 | "language": "python",
619 | "metadata": {},
620 | "outputs": [
621 | {
622 | "output_type": "stream",
623 | "stream": "stdout",
624 | "text": [
625 | "Overwriting /home/kiko/.ipython/profile_pymad/static/custom/custom.js\n"
626 | ]
627 | }
628 | ],
629 | "prompt_number": 13
630 | },
631 | {
632 | "cell_type": "markdown",
633 | "metadata": {},
634 | "source": [
635 | "... things happened in the live talk ..."
636 | ]
637 | },
638 | {
639 | "cell_type": "heading",
640 | "level": 1,
641 | "metadata": {},
642 | "source": [
643 | "M\u00e1s informaci\u00f3n."
644 | ]
645 | },
646 | {
647 | "cell_type": "markdown",
648 | "metadata": {},
649 | "source": [
650 | "http://nbviewer.ipython.org/github/ipython/ipython/blob/master/examples/Notebook/JavaScript%20Notebook%20Extensions.ipynb\n",
651 | "\n",
652 | "https://blog.safaribooksonline.com/2013/12/18/ipython-notebook-plugins/\n",
653 | "\n",
654 | "http://nbviewer.ipython.org/github/payne92/notebooks/blob/master/00%20Javascript%20In%20Notebooks.ipynb"
655 | ]
656 | }
657 | ],
658 | "metadata": {}
659 | }
660 | ]
661 | }
--------------------------------------------------------------------------------
/PythonCientifico-201411-PyConEs-CachemePybonacci/03-Intro-SciPy.ipynb:
--------------------------------------------------------------------------------
1 | {
2 | "metadata": {
3 | "name": "",
4 | "signature": "sha256:7cf591a9240964e4db513cade5897e7fba9b3c501b92ece91b4edf5a704c6970"
5 | },
6 | "nbformat": 3,
7 | "nbformat_minor": 0,
8 | "worksheets": [
9 | {
10 | "cells": [
11 | {
12 | "cell_type": "markdown",
13 | "metadata": {},
14 | "source": [
15 | "# Taller Python Cient\u00edfico - PyConES 2014\n",
16 | "## Parte III: SciPy\n",
17 | "\n",
18 | "#### Juan Luis Cano, Kiko Correoso y \u00c1lex S\u00e1ez ([@Pybonacci](https://twitter.com/Pybonacci))\n",
19 | "#### Fran Navarro, Zuria Bauer y Daniel Domene ([@CAChemEorg](http://twitter.com/CAChemE))"
20 | ]
21 | },
22 | {
23 | "cell_type": "markdown",
24 | "metadata": {},
25 | "source": [
26 | "![\"Licencia](\"http://i.creativecommons.org/l/by/4.0/88x31.png\")
Taller de por Pybonacci y CAChemE se distribuye bajo una Licencia Creative Commons Atribuci\u00f3n 4.0 Internacional." 27 | ] 28 | }, 29 | { 30 | "cell_type": "heading", 31 | "level": 1, 32 | "metadata": {}, 33 | "source": [ 34 | "\u00bfQu\u00e9 es SciPy?" 35 | ] 36 | }, 37 | { 38 | "cell_type": "markdown", 39 | "metadata": {}, 40 | "source": [ 41 | "[SciPy](http://www.scipy.org/) es un conjunto de paquetes open-source para computaci\u00f3n cient\u00edfica general (mat\u00e9maticas, ciencia e ingenier\u00eda). La mayor\u00eda hacen uso NumPy (biblioteca de Python para matrices multidimensionales), otros ofrecen interfaces a programas muy conocidos y utilizados escritos en Fortran, o C++. SciPy ofrece as\u00ed un gran n\u00famero de algoritmos cient\u00edficos de alto nivel, algunos de los temas cubiertos son:\n", 42 | "\n", 43 | "* Funciones especiales ([scipy.special](http://docs.scipy.org/doc/scipy/reference/special.html))\n", 44 | "* Integraci\u00f3n ([scipy.integrate](http://docs.scipy.org/doc/scipy/reference/integrate.html))\n", 45 | "* Estad\u00edstica ([scipy.stats](http://docs.scipy.org/doc/scipy/reference/stats.html))\n", 46 | "* Algebra lineal ([scipy.linalg](http://docs.scipy.org/doc/scipy/reference/linalg.html))\n", 47 | "* Optimizaci\u00f3n ([scipy.optimize](http://docs.scipy.org/doc/scipy/reference/optimize.html))\n", 48 | "* Interpolaci\u00f3n ([scipy.interpolate](http://docs.scipy.org/doc/scipy/reference/interpolate.html))\n", 49 | "* Transformadas de Fourier ([scipy.fftpack](http://docs.scipy.org/doc/scipy/reference/fftpack.html))\n", 50 | "* Procesamiento de se\u00f1ales ([scipy.signal](http://docs.scipy.org/doc/scipy/reference/signal.html))\n", 51 | "* Matrices dispersas ([scipy.sparse](http://docs.scipy.org/doc/scipy/reference/sparse.html))\n", 52 | "* Tratamiento de im\u00e1genes multi-dimensionales ([scipy.ndimage](http://docs.scipy.org/doc/scipy/reference/ndimage.html))\n", 53 | "* Lectura/escritura de archivos (en ingl\u00e9s, I/O) ([scipy.io](http://docs.scipy.org/doc/scipy/reference/io.html))" 54 | ] 55 | }, 56 | { 57 | "cell_type": "markdown", 58 | "metadata": {}, 59 | "source": [ 60 | "Antes de ponernos a programar, vale la pena comprobar si la funcionalidad que deseamos utilizar no est\u00e1 ya implementada en SciPy. Como programadores no profesionales, los cient\u00edficos a menudo tendemos a reinventar la rueda, lo que conduce a errores, c\u00f3digo no optimizado, dif\u00edcil de compartir e imposible de mantener. Por el contrario, los programas de SciPy est\u00e1n optimizados y probados, y por lo tanto deben ser utilizados cuando sea posible.\n", 61 | "\n", 62 | "Como siempre, para hacer el uso de estas funciones, __importaremos__ lo paquetes que vayamos a necesitar:" 63 | ] 64 | }, 65 | { 66 | "cell_type": "markdown", 67 | "metadata": {}, 68 | "source": [ 69 | " import scipy as ...\n", 70 | "\n", 71 | "\n", 72 | "Bueno, en realidad `SciPy` nos ofrece muchas funciones y rara vez querremos usarlo todo en una misma sesi\u00f3n. __T\u00edpicamente usaremos:__\n", 73 | "\n", 74 | " from scipy import ..." 75 | ] 76 | }, 77 | { 78 | "cell_type": "markdown", 79 | "metadata": {}, 80 | "source": [ 81 | "Debido a limitaci\u00f3n de tiempo en este taller, vamos a centrarnos solo en el paquete de integraci\u00f3n. Si quieres, puedes darle un vistazo r\u00e1pido a [Scientific Lectures: SciPy, high-level scientific computing](https://scipy-lectures.github.io/intro/scipy.html) donde puedes encontrar ejemplos m\u00e1s variados." 82 | ] 83 | }, 84 | { 85 | "cell_type": "heading", 86 | "level": 1, 87 | "metadata": {}, 88 | "source": [ 89 | "Integraci\u00f3n mediante m\u00e9todos num\u00e9ricos con SciPy/Python" 90 | ] 91 | }, 92 | { 93 | "cell_type": "markdown", 94 | "metadata": {}, 95 | "source": [ 96 | "SciPy proporciona el paquete `integrate` que proporciona algunas __t\u00e9cnicas de integraci\u00f3n tanto de funciones como de ecuaciones diferenciales__. En primer lugar import\u00e9mos el paquete y el resto de librer\u00edas que vayamos a necesitar. " 97 | ] 98 | }, 99 | { 100 | "cell_type": "code", 101 | "collapsed": false, 102 | "input": [ 103 | "import numpy as np\n", 104 | "import matplotlib.pyplot as plt\n", 105 | "\n", 106 | "%matplotlib inline" 107 | ], 108 | "language": "python", 109 | "metadata": {}, 110 | "outputs": [], 111 | "prompt_number": 1 112 | }, 113 | { 114 | "cell_type": "code", 115 | "collapsed": false, 116 | "input": [ 117 | "from scipy import integrate" 118 | ], 119 | "language": "python", 120 | "metadata": {}, 121 | "outputs": [], 122 | "prompt_number": 2 123 | }, 124 | { 125 | "cell_type": "heading", 126 | "level": 2, 127 | "metadata": {}, 128 | "source": [ 129 | "Cuadratura de Gauss" 130 | ] 131 | }, 132 | { 133 | "cell_type": "markdown", 134 | "metadata": {}, 135 | "source": [ 136 | "Como se puede ver en la [ayuda](http://docs.scipy.org/doc/scipy/reference/tutorial/integrate.html), si queremos realizar una _integraci\u00f3n num\u00e9rica_ de una funci\u00f3n de una variable, podemos utilizar el [m\u00e9todo de cuadratura de Gauss](http://es.wikipedia.org/wiki/Cuadratura_de_Gauss) [`quad`](http://docs.scipy.org/doc/scipy-0.14.0/reference/generated/scipy.integrate.quad.html) (aunque tambi\u00e9n podemos usar `trapz`, `simps`...). La forma de acceder a ella tal y como hemos importado el paquete ser\u00eda ejecutando `integrate.quad`. Sin embargo, es m\u00e1s com\u00fan importarla del siguiete modo:" 137 | ] 138 | }, 139 | { 140 | "cell_type": "code", 141 | "collapsed": false, 142 | "input": [ 143 | "from scipy.integrate import quad" 144 | ], 145 | "language": "python", 146 | "metadata": {}, 147 | "outputs": [], 148 | "prompt_number": 3 149 | }, 150 | { 151 | "cell_type": "markdown", 152 | "metadata": {}, 153 | "source": [ 154 | "\u00bfQu\u00e9 es lo primero que necesitamos hacer para integrar una funci\u00f3n? Pues s\u00ed, la funci\u00f3n... definamos una:\n", 155 | "\n", 156 | "$$f(x) = x sin(x)$$" 157 | ] 158 | }, 159 | { 160 | "cell_type": "code", 161 | "collapsed": false, 162 | "input": [ 163 | "def fun(x):\n", 164 | " return x * np.sin(x)" 165 | ], 166 | "language": "python", 167 | "metadata": {}, 168 | "outputs": [], 169 | "prompt_number": 4 170 | }, 171 | { 172 | "cell_type": "markdown", 173 | "metadata": {}, 174 | "source": [ 175 | "Integraremos la funci\u00f3n en el intervalo $[2, 9]$:\n", 176 | "$$ \\int_{2}^{9} x sin(x) dx$$" 177 | ] 178 | }, 179 | { 180 | "cell_type": "code", 181 | "collapsed": false, 182 | "input": [ 183 | "x_inferior = 2\n", 184 | "x_superior = 9\n", 185 | "\n", 186 | "value, err = quad(fun, x_inferior, x_superior)\n", 187 | "\n", 188 | "print('El resultado es: %f con un error de %e' %(value, err))" 189 | ], 190 | "language": "python", 191 | "metadata": {}, 192 | "outputs": [ 193 | { 194 | "output_type": "stream", 195 | "stream": "stdout", 196 | "text": [ 197 | "El resultado es: 6.870700 con un error de 2.864870e-13\n" 198 | ] 199 | } 200 | ], 201 | "prompt_number": 5 202 | }, 203 | { 204 | "cell_type": "markdown", 205 | "metadata": {}, 206 | "source": [ 207 | "Si lo deseamos, podemos hacer uso de lo que acabamos de aprender en la secci\u00f3n previa de este taller y representar la funci\u00f3n e intervalo de integraci\u00f3n que hemos tomado. F\u00edjate en como hemos hecho uso de $\\LaTeX$ en el t\u00edtulo del gr\u00e1fico." 208 | ] 209 | }, 210 | { 211 | "cell_type": "code", 212 | "collapsed": false, 213 | "input": [ 214 | "# Vector con valores de la variable independiente\n", 215 | "x = np.linspace(0,10,100)\n", 216 | "\n", 217 | "# Hacemos uso de la funci\u00f3n para obtener los valores\n", 218 | "# de la variable dependiente\n", 219 | "y = fun(x)\n", 220 | "\n", 221 | "# Representamos la figura\n", 222 | "plt.plot(x,y, linewidth = 2)\n", 223 | "\n", 224 | "# A\u00f1adimos un t\u00edtulo\n", 225 | "plt.title(u'$y = x sin(x)$', fontsize = 25)\n", 226 | "\n", 227 | "# Creamos el intervalo para sombrear el \u00e1rea de integraci\u00f3n\n", 228 | "x_fill = np.linspace(2,9,100)\n", 229 | "# Obtenemos los puntos sobre la curva\n", 230 | "y_fill = fun(x_fill)\n", 231 | "\n", 232 | "# Dibujamos el \u00e1rea sombreada\n", 233 | "plt.fill_between(x_fill,y_fill, color = 'gray', alpha = 0.5)\n", 234 | "\n", 235 | "# Representamos los resultados con una gradilla\n", 236 | "plt.grid()" 237 | ], 238 | "language": "python", 239 | "metadata": {}, 240 | "outputs": [ 241 | { 242 | "metadata": {}, 243 | "output_type": "display_data", 244 | "png": "iVBORw0KGgoAAAANSUhEUgAAAXMAAAEVCAYAAAASFwXVAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd4VFX6wPHvOzMJJZTQkRpQkCIQehNERH7ogq4NxYq9\nrauuu7qWVVfXuhYWwY4i7loQUVFxUdCwgoAgBJCiIKKC9B5akpnz++PMMCEmkGTKvTPzfp4nz8yd\n3Ln3nZOZN3fee865YoxBKaVUYvM4HYBSSqnIaTJXSqkkoMlcKaWSgCZzpZRKAprMlVIqCWgyV0qp\nJKDJXCmlkoAmc6WUSgKazJVymIjUFpFJIpIZhW2dKyJ/jUZcKrFoMlfKQSKSAXwA/N0YszPS7Rlj\nJgEZInJvxMGphKLJXKlSiEiGiMwSkZ9EpGqMdjMOeMMYszSK27wPGCwig6O4TeVymsyVKl1ToBew\nAzgQ7Y2LyBnA8caY56K5XWNMALgFGC8i1aO5beVeohNtKVU6EakHbDfG+KO8XQ+wDLjbGDM5mtsu\nso+pwCxjzMOx2L5yFz0yV+oIjDFbop3Ig4YBtbH18lh5CbhVRNJjuA/lEprMlXLGJcDHMfpHEfIp\nUA0YFMN9KJfQZK5UnIlIGjAE+CKW+zHG7AXmA2fGcj/KHTSZqxKJyP0iMk9EVojIkBJ+f4mIbBCR\nhnGO63ci8qmIzBSR2SLSQ0S8IvJE8LGZInL5EZ5fR0SeFpGvRGSGiLwrIrcWW+dor722iEwVkQUi\n8pqIVArGcLOIfBaM4TsRubGUMDoDVYHFsX69wELsSVyV7Iwx+qM/h/0AlwN3Bu+/A6wrYZ2vgHyg\nWhzjugx4AagcXL4F29PkBeBcIB2YCviBzBKeXx/4AfhDkcdOAn4FTi3Ha38e6I0tYQSA+4F3geFF\n1vlT8HcnlfD8q4K/qxHL1xt8zs3AQcDr9PtKf2L744vg/4BKQiLiA0YYY0J9lNtik0XRdTKAbsAi\nY0xeCdsYDlwfQRibjTHnF9vmMcDZxpiiJYPlQE2gnjFmkoh0wpYvVgG7S9jubcAeY8yYIo91Awyw\npoyvPR04zhgzR0SaBh++HRhkjPmqyKrfBm97ADOLxZEFGGNMSTFG8/UCrAfSgCbAT6XtTyU+Teaq\nuH7AFAARyQbaAc8UW6cv9r2TU9IGjDETgYlRjusi4Ilij7UP3r4V3O9iETkB+MXYvtbFnQA0FZF2\nxpjlwec8CTwJICIDgQ+D9ztS8mvvDswI3g+VL8YUS+Rg+6gDbCshjlrA0UZ7RuP1AuwN3tY8yv5U\ngtNkrg5jjPmC8Im5q4O3rxRb7aTgbfEjzpgxxhRPbAAnY8sVnxdZb/kRNjMXOA1YLCLLgP8BL5ng\n6EtjzOdFtnVl8Paw126MmQ3MDi4ODN6+XcK+Qr+bVcLvKgO/+UZTbD/ReL1gSyxgS0IqiekJUFWi\nYMlhOLDcGFP8RN0AbPnhy3jHFRKM7yTgW2PM1jI+7RHgAWyNvAPwB2ChiJxVbNvp2CPj3BJee1Gn\nYGvYC4s9vxIwFNt235fwvAPYeneZVfD1gv3HEdqnSmKazFVpTgTqAB8VfTA4R0l3bKLb40RgQT2B\n6oRLHkdljCk0xtxvjGmOrSH/OfirO4ut+nvsgJ5XS9uWiDQBjgNyjDHFh1GfHozt36U8fTvlP1Iu\n9+sNCu1nRzmfpxKMlllUadoFb+cXe/yI9XIAETkfuC6CfW8xxgwvts362N4fq4MPnRq8zSm23vXA\nemPMlCKPDQPGAw8aY0YBGGM2AE+JSH+gUbH9X4ktT/w7+Py7gRnGmLlF1jkleFtScr0CKAzuExG5\nJ/j8OcHf/4id2bB6af8QI3m9xTTGlmZ+KeX3KkloMlelyQjeFu8lMSJ4m1PaE40xb1NyHblCRKQW\ntndIbRGpE4zpvOCvlxVZrwZwPuFEG3Id9gTgxmLbFaAl8FqRx6oGn/+RMWZHcIDPScaYh4ptM1QT\nPyyZi0gV4P+Az4wxG4PPH2CM+UeR1ZYEb1tSQl/zKLzeoloC3xljCo+wjkoCWmZRpfkMe0Q3CCA4\nUOUeYCTxr5dnAXWB17GJ7U5gOlAAdAzGVxd4E9tHvPgQ+YXAvyjyD0ZEagMvAxuAUcXWLwS+C95/\nGBjDbw0EfjXGfFfsccG2WyjpPgI8W2ydXGwvk64lvVgif71FdcaOCVBJTpO5KpExJhe4EDhdROZi\nk3vo6G7xkfpIxyCWRdg5ultge9D4jTE3A+cA94jIbGyifqhIKaOo+7FlhmnBEZrTgEnYf0inF02G\nxph92IFDZ4jI15RQwgieIA3w238CRZ8/LPj8X02xWRGDR8kfY3unxOL1huKshu1H/2Fp66jkEfEU\nuCJyJ3Ax9s29FLjcGHPwyM9SiUjs5cgeBv5WQtlBlYOI/A5b3mkYqxKIiJwLPAc0Nsbkx2Ifyj0i\nOjIXkSxsX+QuxpgOgBe4IPKwlJNE5AoR2R5MBqHHPMClwC7s1KoqAsaYj4FN2GH5sXINMFYTeWqI\ntMyyG1vHqxrsB1sVO3xYJbZbsf2TNxR57AagDXCzMWazI1Eln78CdwVPxEaViPTEjhh9MtrbVu4U\nUTI3xmzHvll+xg7E2GmMmR6NwJSjfgAuDI52REQuAh4F/miMmeBoZEnEGPMhdr6V0mZXrBAR8WLr\n+Tc6PBZAxVFENXMRORZ7cqUf9uv3O8AkY8x/ohOeckJwkqdnsQNnKmGT+2PGmCVHfKIqt+CkZZ9h\nE++iKG3zIaDQGHNfNLanEkOkyfx87NShVwWXLwF6GWNuLLKOXmRUKaUqwBhT5hJcpDXzlUAvEakS\nrPsNwn5tLB6Q/hjDfffd53gMbvnRtkjdtliyxDBgQD525uHwj8cTwPbItMtpaX7uvLOQvXudj9mJ\nn/KKaASosVNwTgAWYLsmLgRejGSbyWzt2rVOh+Aa2hZhqdIWeXlw221+Xn7ZQyCQRkZGPr///a+0\nbbuHVq3yqFMnn3vvXUC/fiuZNy+TnJyGPPII/PvfhYwd62XYsKifJ04qEQ/nN8Y8DjwehViUUklq\n+3YYMqSQ+fN9eDwBzjxzPZdfvpaaNQ/vYl+9eiFDhmxkyJCNnH32r4wa1Yo1a6pzxhkwZkyAG2/U\ncY6l0ZaJo5EjRzodgmtoW4Qle1ts3AgnnljA/Pk+6tffx0svfcMtt6z+TSIHGDIkfMnVDh128+KL\nC7niih8A+MMfPIwde6SZC1JbxCNAj7oDERPrfSil3GntWjj55ELWrvXRtGkeTz65lHr1yj9A/N13\nGzFmTGsAnnsuwHXXJf9xqIhg4ngCVJVDTk6O0yG4hrZFWLK2xa5dcMopftau9XHccbsYPXrxURN5\nbm5uiY+fc86v3HDDKgCuv97Dm2/qAWJxmsyVUlFnDFx8cSFr1nhp0WI3Tz+9hMzMgoi2ed5567n2\nWju9+9VXG378MRqRJg8tsyilou6RR/zcdZeXjIwCXnjhGxo3js5V64yB++5rx5df1qdHj0Jmz/bh\nS9KrMmiZRSnlqBkzDPfcY1PLXXetjFoiBxCB2277njp1DvD11z4efFBPiIZoMo+jZK2NVoS2RVgy\ntcXmzXD++QECAeGii9bSp8+2cj2/tJp5UTVrFnLnnSsBeOghD3NKndE9tWgyV0pFza23+tm2zUun\nTtu5/PK1MdtP1647GT78J/x+4ZJL/OTrJL9aM1dKRcf06YZTTxXS0/288sr8qJZXSpKfL1x5ZVfW\nravGqFEBbr45uY5NtWaulIq7Awfg2mtt/fqSS36KeSIHSE83XHfdGgDuv9+wY0fMd+lqmszjKJlq\no5HStghLhrZ4+OEAa9b4aNYsj/PP/6XC2ylLzbyoPn2206nTdnbu9PLggzG5+l7C0GSulIrId9/B\no4/a+7fdtoq0tPiVVUXg+uvt0fmYMZ6U7nuuNXOlVESGDvXz8cdehgxZzx13rHIkhoceOp7p04/h\nvPP8TJzodSSGaNOauVIqbubMgY8/9lKpUiFXX73WsTiuvHItaWl+3nnHyzffOBaGozSZx1Ey1Eaj\nRdsiLFHbwhi4/XZbpz733HXUrh3ZcH0of808pGHDg5x55joAHnkkNWvnmsyVUhUyYwbMmuUjIyOf\nCy5Y53Q4nHfer3i9Ad57z8uaNU5HE39aM1dKlZsx0K1bIQsX+rjqqh+46KKK92CJplDt/IYb/Iwd\nm9i1c62ZK6Vi7r33DAsX+qhV6wBnn73e6XAOCX1DePVVYVv5ZhJIeJrM4yhRa6OxoG0RlmhtEQjA\n3XeHBgj9TJUqgahtu6I185Bjj91Lt25b2b/fw9ix0YsrEWgyV0qVy9SpsHKlj7p19zN06Aanw/mN\n0NH5M88YDsR+IKpraM1cKVUuffoUMmeOj+uvX8Xw4e4psYQYA1df3ZkffqjJ888brr22zGVnV9Ga\nuVIqZubMgTlzfGRkFDB06EanwymRSNGjcz+pciypyTyOEq02GkvaFmGJ1BahPtxnnLGeqlWjf2GI\nSGvmIf37b6VGjYMsW+ZLmUFEESdzEckUkUkiskJElotIr2gEppRyl++/h48+8uLz+TnnnF+dDueI\n0tMNgwdvAuCFF1LjakQR18xF5DVgpjHmFRHxARnGmF1Ffq81c6WSwFVX+Rk3zstpp63n9tudmYOl\nPNaurcrll/egWrUAGzd6yMhwOqLyiWvNXERqAv2MMa8AGGMKiyZypVRy2LQJXn9dEDGuGO1ZFllZ\n+2jXbgd5eR4mTkz+A8pIyywtgC0i8qqILBSRl0SkajQCS0aJVBuNNW2LsERoixdeCJCf76Fnzy00\na7Y/ZvuJVs08ZOhQW2p5/vnkn6/FF4XndwH+YIyZLyKjgL8C9xZdaeTIkWRlZQGQmZlJdnY2AwYM\nAMJvZF1OreUQt8Tj5HJubq6r4im+7PfD88/3AyA7ewq5uXlkZ2cD4eQbreXVq1dHdXu1a39Kevo6\nvv56EMuXw+bN8W+/si7n5OQwfvx4gEP5sjwiqpmLSENgjjGmRXD5ROCvxpihRdbRmrlSCWzyZMM5\n5wiNG+cxYcICPAnWB+6JJ1rx8ceNueWWAE8/nTjBx7VmbozZCPwiIq2DDw0ClkWyTaWUuzzzjC1R\nnHnmhoRL5MCh/vATJhjy8x0OJoai8ae5CfiPiCwGOgIPR2GbSal4iSGVaVuEubktvvsOcnLSSE8v\nZMiQTTHfX7Rr5gDHH7+HZs32sH27l+nTo75514g4mRtjFhtjuhtjOhljztbeLEolj2eftX20Tzll\nE9WrJ+ZJRBEYNGgLAP/5T/L2Ode5WZRSJdq7Fxo1CrB7t4cXXlhA69Z5TodUYevXV+Hii3tSrVqA\nzZs9VKnidERHp3OzKKWi4s03Dbt3e2jTZmdCJ3KAxo3306rVLvLyPHzyidPRxIYm8zhyc2003rQt\nwtzaFs8/b0sSv/99/Ibux6JmHjJwYKjUkpjloqPRZK6U+o1vv4VvvrGzIw4YsNXpcKLi5JNtMp86\n1UNeYn/RKJEm8zgKDRRQ2hZFubEtxo2zR+Unn7yJSpXid8We0GCfWGjQ4CDt2u3gwAEPU6Yk33k8\nTeZKqcMUFMDrr9v7p58e++6I8XTKKcnbq0WTeRy5tTbqBG2LMLe1xccfw7ZtXpo330ObNnviuu9Y\n1swBBgzYisdj+OwzLzt2xHRXcafJXCl1mBdftCcITzttE5KYV1wrVe3a+XTsuJ2CAuH995Or1KL9\nzJVSh2zcCE2aGMAwceIcatcucDqkqPvgg2MYNep4hg4t5MMPI51rMHa0n7lSqsImTAjg9ws9emxN\nykQO0LfvNgCmT/ewb5/DwUSRJvM4cltt1EnaFmFuaQtj4OWXbc+V3/3OmROfsa6ZA9Stm0/r1js5\ncMCTVHO1aDJXSgGwYAGsWuWjZs2D9Oy53elwYurEE+3rmzw5eQYQac1cKQXATTf5GTPGy9ln/8xN\nN61xOpyY+vHHDK64ojt16vjZtMmL1+t0RL+lNXOlVLkVFsJbb9n7gwdvcTaYOMjK2kvDhnvZts3L\n3LlORxMdmszjyC21UTfQtghzQ1tMnw5bt3pp3DiP1q3j27e8qHjUzMFOixs6ETp5cnIMINJkrpRi\nwgRbOx48eHPS9S0vzYkn2mT+wQfJUQbWmrlSKS4vD+rXD7B/v4f//GcujRodcDqkuPD7hbPO6s2e\nPemsWAFt2jgd0eG0Zq6UKpf33zfs3++hbdsdKZPIAbxeQ8+e9ug8GUaDajKPIzfURt1C2yLM6bZ4\n/XVbMx48eLOjcUD8auYhoVLLe+8lft1ck7lSKWzzZpgxw4vXG+Dkk5Nj3vLy6N59Bz5fgAULvGxP\n8K71mszjyI3zVjtF2yLMybZ4+22D3y90776VmjWdH74fy/nMS1K1qp/27XcQCEjCjwbVZK5UCgvN\n6z1oUOodlYf07LkTgI8+SuzRoJrM48jp2qibaFuEOdUWv/wC8+b5SE8vpE+fbY7EUFy8a+bAoakL\npk0TErnjXVSSuYh4RWSRiHwYje0ppWLv7bdt5urVaytVqiT+CcCKatFiL7VrH2DzZi9LljgdTcVF\n68j8ZmA5kMD/12JP68Rh2hZhTrXFG2/YBD5woHtKLPGumYMdDdqjhz06nzo1ftc7jbaIk7mINAFO\nB14GUmTsmFKJ7YcfYNEiH5UrF9KrV4J344iCXr3sNeQ+/jiFkznwNPAXIHFbIU60ThymbRHmRFu8\n/bb9uPbps4VKldzz0XWiZg7QtesOPJ4A8+Z52b3bkRAiFtE1k0RkKLDZGLNIRAaUtt7IkSPJysoC\nIDMzk+zs7ENfLUNvZF1OreUQt8Tj5HJubm7c9//mmycCHlq0+ITc3N2HyhuhZOrU8urVqx3bf5s2\nu1i+fDH/+hf87W/la89oLOfk5DB+/HiAQ/myPCKam0VEHgYuAQqBykAN4F1jzKVF1tG5WZRykZUr\noW1byMgoYPLkr0hP188nwOuvN+OVV1py5ZV+Xn7Z+QnO4zo3izHmLmNMU2NMC+AC4POiiVwp5T5v\nvWXLKn37btFEXkToJOh//0tCdlGMdj/zBGyC+CleYkhl2hZh8W6LN9+0yfyUU9zTiyXEqZo5QKtW\nedSseZD1672sWOFYGBUWtWRujJlpjDkjWttTSkXf8uXw/fc+qlXLp0uXHU6H4yoeD3TpYo/Op09P\nvONSHQEaR071J3YjbYuweLbFO++ESixb8fncl7Cc6GdeVLduuwCYMcPRMCpEk7lSKSTUJTEVZ0gs\ni9C3lZwce13URKLJPI60ThymbREWr7ZYuRJWrPCRkeHeEouTNXOAhg0P0rBhHrt3C99842go5abJ\nXKkUMWlSaKDQVtLS3FdicYuOHbcAiVdq0WQeR1onDtO2CItXW4RLLO6YIbEkTtfMQZO5UsrFVq2C\nb7/1UbVqAV276lwsR9Khgz2fMHs27N/vcDDloMk8jrROHKZtERaPtgiVWHr33urqgUJO18wBatTI\nJzvbcPAgfPWV09GUnSZzpVJAqMQyYID2YimLk0+2//ASqdQS0dwsZdqBzs2ilKPWrIFjj4XKlQt5\n//2vXDVLohvt3LkTn28Yw4Z56dED5s1zJo64zs2ilHK/d98NX1FIE3nZ9OsHPh8sWAA7dzodTdlo\nMo8jrROHaVuExbot3nnHXlFowAD39mIJcUPNHKBaNejVCwIBmDnT6WjKRpO5Ukls/XqYP99Herr/\n0KyAqmwGDbK3iVI312QeR9q3OkzbIiyWbfHee7bE0q3btoS4aLMb+pmHhP4semSulHLcxImJU2Jx\nm549oVIlWLIEtiVA82kyjyOtE4dpW4TFqi22bIHZs734fAF6906AbIR7auYAlSvbujnAl186G0tZ\naDJXKkl98IEhEBCys7dTrVqCTQHoEiedZG8TodSiyTyOtE4cpm0RFqu2CPVi6d8/cQYKualmDuG6\neSJ8kdRkrlQS2rULvvjCi8dj6NcvMUosbtSrF6Snw+LFsMOdswYfosk8jrROHKZtERaLtvjoIygo\nEE44YQeZmQVR336suKlmDlClCvToYS/wPGuW09EcmSZzpZLQO+/YGnkilVjcKlQ3d/vxh87NolSS\n2bcP6tYNsH+/h4kT51Cv3kGnQ0ooO3fuZNiwYXi9XgCmT4dTT4WuXe3w/njRuVmUSnGffgr793to\n3XqnJvIo6N3bztOyaJE9F+FWmszjSOvEYdoWYdFui1AvlpNOSrwTn26rmQNkZED37naeFjfXzSNO\n5iLSVES+EJFlIvKtiPwxGoEppcqvoAA++sh+M+/XT+vl0ZIIQ/ujcWReANxqjGkP9AJuFJG2Udhu\n0tG+1WHaFmHRbIucHNi920Pz5nto2jSBrnkW5LZ+5iGJMHgo4mRujNlojMkN3s8DVgCNIt2uUqr8\nwgOFEq/E4mZ9+4LXC998A3l5TkdTsqjWzEUkC+gMOHRtDnfTOnGYtkVYtNrC74f337f3E7XE4saa\nOdj5zTt3tm08d67T0ZTMF60NiUg1YBJwc/AI/ZCRI0eSlZUFQGZmJtnZ2Ye+WobeyLqcWsshbonH\nyeXc3NyobG/uXNiy5Utq1TrAccdVBsLJMVS+cPvy6tWrHY8nLy+PYcOG/aZ9TzwRFizI4d//hkGD\nBvzm95Eu5+TkMH78eIBD+bI8otLPXETSgI+AT4wxo4r9TvuZKxUHf/pTgKef9nDuuT9z441rnA4n\nYRXvZx4yeTKccw6ccortex5r5e1nHvGRuYgIMA5YXjyRp6odO2yf1EWLYNkyOxXp1q2GrVv9eDxQ\nvbqhRg2oX1/o3NlD164eOneGOnWcjlwlKmPg3XcDgEfr5THSt6+9nTPH9hpKS3M2nuKiUWbpC1wM\nLBGRRcHH7jTG/DcK204IoTrahx/aOTGWLStpLQFmAQMOe/Ttt8P3s7MLGTHCw/DhHirwLSuh5OTk\naI+WoGi0xeLF8PPPPjIzD9CunYtHthxFbm6ua3u0NGgArVrBqlWQm2v7nrtJxMncGDOLFB189P33\n8NJLMGECbN4cfrxyZejUCbp0gY4d4ZhjwOPZwfTpC8nOrsq+fT727vWyZUslVq2qxqpVGaxeXZ3c\nXB+5uXDHHdC3byF33eXltNMEKfMXLZWqJk2yR+UnnriVYtUBFUX9+tlkPmuW+5K5zs1SToEATJkC\n//rX4RPvtGwJw4bZn3797LSZRW3ZsoWvvvqKmjVrlrjdgwc9zJtXm5ycusyZU5cDB+z/2U6dCrn3\nXi9nnaVJXZWubdtCVq708fjji+ne3eVztbpcaTVzgFdfhSuugLPOsjX0WIp7zTxVFBbaksgjj4TL\nKFWrwgUXwDXX2GkyI0m2lSoF6N9/K/37b2XfPi9TphzD2283YfHiypxzDpx0UiEvveSjVavovB6V\nPL77Dlau9JGRkU929k6nw0lq/frZ21mz7HkKNx1gpWR5pDyMgXfegbZt4eKLbSJv0gRGjYJff4Vx\n4+yFX8vyR126dGmZ9lm1qp8LLljHW299zU03fU+NGgeZOdPHCScEePDBAPn5Eb4oFyjeRTGVRdoW\nkycHAOjdextpaYn9Ldit/cxDjj3W1s63bLHlFjfRZH4E//ufnTFt+HBYvRqOOw5efhl++AFuvhlK\nqZhETaVKAc4++1cmTFjA4MG/kp/v4d57PfTq5eenn2K7b5U4bL08MSfWSjQicOKJ9r7bLvKsybwE\n69bZBH7SSTBvnv1P/PzzsGIFXHnlb+vhZdWhQ4cKPa9mzQLuvPN7nngil/r197FokZfsbD/TpiXu\nUZj2ZAmLpC1++QUWLvRRqVIh3bptj15QDnFrT5aiipZa3ESTeRH5+fDYY9CmjS2tVKkC995rv05d\ne62d09hJXbvu5KWXFtK9+1Z27vRy2mnw4IMBkuj8siqn996zf/zu3bdRuXLA4WhSgx6Zu9zs2ZCd\nDX/9K+zda0d6rVwJf/87VK8enX2UtWZ+JDVqFPLoo99y2WV2hN+993q4+mo/fn/Em44rrZmHRdIW\nyTaxlttr5mC7HVerZsutGzY4HU1YyifzXbvg+uvtf9sVK+yggGnTYNIkaNbM6ehK5vHAyJE/8+CD\nS0lP9zNunJfhwws5qBeVSSlbtsBXX3nx+QL07p34JZZE4fPZc2lgDwLdIqWT+YcfQrt2th7u88E9\n98CSJTB4cGz2V9GaeWn69t3O448voWrVAiZP9nH66X7XTs9ZnNbMwyraFu+/bwgEhE6dtlOtWmF0\ng3JIItTMITy0/6uvnI2jqJRM5lu3wkUXwRln2O6FvXrZeVQefNCO3kwknTrtYtSoxWRmHuTzz70M\nHernwAGno1LxMGlS6PJwiTndbSLr08fe6pG5gyZNgvbt4Y037AnOUaPsWekTToj9vqNRMy9Jq1Z5\njB6dS61aB5g508t55xVSUBCTXUWN1szDKtIWO3fCF1948XgM/folR70cEqNmDnZsiccDCxfCfpdc\n0CllkvnmzXDeefZn82Z7Tb+lS21/8WSYy6Jp0/088cRSqlXL56OPfFx2mZ+Adm5IWh99BAUFwgkn\nbCcz0+X/uZNQjRrQoYMdGb5ggdPRWEmfzI2BN9+0tfFJk+xZ6GefhRkz7GiueIp2zby4li338thj\nS6lcuZA33/Tyxz/6XdttUWvmYRVpi4kTbY082QYKJUrNHMKlFrfUzZM6ma9fD2eeCRdeCNu2waBB\n8O23tveKJ0lfebt2e/jHP5aSluZn7FgvY8bo4XmyycuDTz+1b+BEvTxcMnBb3TwpU5oxdth9+/a2\nx0qNGnaq2k8/hebNnYsrVjXz4rp23cVf/vIdALfeKnz6qfsOz7VmHlbetvjkEzvLZtu2O6hXLwkm\n6ikiUWrmcHiPFjd8A066ZL5yJZx8Mlx9te1DPmwYLF8OV13lrhnOYu3UUzdz4YVr8fuFc881rFzp\ndEQqWpK1xJJosrKgYUP7rd8Nk24lTTI/cADuv9+Ozpo5E+rVsz1WPvgAGjd2Ojor1jXz4q68ci19\n+25mzx4Pp5/uZ4eLprnWmnlYedriwAH45BP7se3fP/lKLIlUMxdxV908KZL5xx/bM8t//7udX+XK\nK+0R+ogH5+0WAAAacklEQVQRqXU0XpzHA3ff/R0tW+7mxx+9XHqpe0+IqrL57DPYu9fDccft4phj\ndECB00KlFjfUzRM6ma9eDUOH2p/Vq+2c4zNn2np57dpOR/db8aqZF1Wlip9//GM5GRkFfPSRlyef\ndMcJUa2Zh5WnLSZODA0USs4SSyLVzEGPzCO2eTP88Y+2u+HHH9sTnE89ZS9q27+/09G5zzHHHOCO\nO2zR/M47hTlzHA5IVUh+PkyZYr9q9u+/xeFoFEDnzlCpkj0v53QZM6GS+a5dtpRy7LHwzDO2w/7l\nl9sLK996K6SlOR3hkcW7Zl5Uv37bOPfcnyksFM47z882hw/stGYeVta2mDEDdu/2kJW1m2bNXDLs\nMMoSqWYONpF362bvz53rbCwJkcy3bLGTYDVvbk9y5uXZ0srixfDKK/biEerorr32R9q02cn69Vo/\nT0RvvWVLLCefnHwnPhOZW+rmrk7mS5bAjTfaJP7QQ/bI/KSTbF38ww/tSc9E4kTNvCifz3D//SvI\nyChg6lQvL73kXP1ca+ZhZWmLggL44ANbYknmibUSrWYO7qmbR5zMRWSIiKwUkVUickek29uxA159\n1c4X3KmTHXq/f789Ep89G3JytC4eiQYNDnLLLd8DtjS1erXDAakymTEDdu3y0Lz5Hpo33+d0OKqI\n0NzmX39tS79OiSiZi4gXGAMMAdoBI0SkbXm2YQysXQsvvghDhkD9+nDFFbb+VKOGPTJfutQeiYf+\nAyYqJ2vmRQ0atIUBAzayb5+Hiy4qdOQNqDXzsLK0xdtv+4PrJveJz0SrmYPNWS1b2iuUffutc3FE\nelXLHsBqY8xaABF5CzgTWFHSygcOwJo19mhw1SqYP99OP7t+fXgdrxcGDrTzjZ9/PmRkRBihKtGf\n/rSapUtr8vXXVXj44QD33uvqiltKKyiA99+39wcMSN4SSyLr3dvmtjlz7OUnnRBpMm8M/FJkeR3Q\ns/hKdevak5alXdasdm17xeszzrA/detGGJVL5ebm0rVrV6fDAOxgqltuWczf/taLBx6AE0/cRvv2\n8TtEnz17Nn1DZ45S3NHaIicnnZ07a9G48S7q1t2cMFeTqoilS5c6/g3WVKBnQK9e8J//2IrC9dfH\nIKgyiDSZl+lVb9s2EshCBGrXzqRly2x69hzACSdApUo5NGsGAwcOAMIng0JfPZNluUePHjRu3JgN\nwSvAdgv2Z1oQnAzZieXGjWH69MnMnFmbq67qyBNPzGP58iVAuCQUOmkb7WUAr9cbs+0n0vKPP/6I\nNzipfkm/f++9Y4ERnH76Ple9f2KxnJOTw4YNGxyNJy0tDU9wWtWyfr5797bLM2bkkJNTsXyRk5PD\n+PHjAcjKyqK8pCL/hQ49WaQXcL8xZkhw+U4gYIx5rMg6ZvNmQ7Vq9pJsqTy83o327IF27fysW+fl\nqqt+4KKLfjn6k1Tc+P3CWWf1Ys+eSixdGp8rYqnyKyiAmjVtZ40tW6JTXRARjDFlzpiRFkoXAK1E\nJEtE0oHzgSnFV6pXz16iTRO5+1SvDuPG2aPC117L4uefqzgbkDrMwoWZ7NlTiVatCmnf3uloVGnS\n0qB7d3vfqcFDESVzY0wh8AdgGrAceNsYU+LJT+XevtWDB8Oll/opKPDy2GOt43K5uUTsTxwrR2qL\nGTPqAXDhhZ6UOBhy62ekLEJdFJ2aLiPiLgzGmE+MMccbY44zxjwSjaBU/I0a5aV+fT/Ll9figw8a\nOR2OAvLzhVmzbDK/4ALtbeR2vXrZW6eOzCOqmZdpByIm1vtQ0fHuu4ZzzxWqVi1gwoT51KmTXFex\nSTSzZ9fhnns60KFDIUuWRNpXQcXapk32YhXVqsHOnZFfKD7eNXOVRM4+WzjtND/79qXxzDMtnQ4n\n5YVKLCNG6Mc0ETRoAC1a2G7YTgwe0ndJHLm9HigCzz3npUqVADNnNuTrr2vFbF9aMw8rqS327/fw\n1Ve2S0QqlVjc/hk5Gifr5qnzLlFl0rw53Hefvf/008dx8KC+RZwwd24dDh700bVrIS1aOB2NKitN\n5ikiUeYj+dOfPLRt62fjxgxef71ZTPaRiHNwxEpJbfH55/ao/KKLIiy8JphE+YyUJnQSVJO5coW0\nNHj5ZZtE3n67Gb/8on3P42nvXi/z5tVFxDB8eAr0R0winTrZMTWrVsHWOE+jo8k8jhKpHtinD1x2\nmZ/CQg+jRx8b9QtZaM08rHhbzJpVh4ICL336+Gnc2KGgHJJIn5GSpKWFrzz09dfx3bcmc1Wqxx/3\nUqNGgAUL6jJ7dh2nw0kZn31WH0i9EkuycKq/ufYzV0c0erTh5puFBg328dprC6hUybmrE6WC7dvT\nOe+83ni9sGGDUEf/hyacyZPhnHPg1FPh008rvh3tZ66i6oYbhPbt/WzaVJU33mjqdDhJb8aMegQC\nwv/9n18TeYIKHZnPm0dcpsYI0WQeR4lYD/T5bN9zgDffbMaGDZWjsl2tmYcVbYtPP7UllssuS80R\nn4n4GSmuUSNo2hR274aVK+O3X03m6qj69YMLLrATcT37rHZ6jpWff67K6tU1qV49wNChTkejIuFE\nF0VN5nGUyH1on3jCS9WqAWbNasDChZkRb0/7mYeF2iJ04vPsswNUjs4XoISTyJ+Ropw4CarJXJVJ\n48Zw1132/ujRx+L3a//naDImnMxTtcSSTEIjQTWZJ6lErwfedpuHrCw/P/1UnSlTjoloW1ozD8vN\nzeXbb2uwaVNVGjXyc9JJTkfknET/jIR07mz7nC9bZmvn8aDJXJVZ5crw1FP2LfPKK1ns2qVHkNES\nOvF54YWCRz+VCa9yZZvQjYH58+OzT33bxFEy1AN//3thwAA/eXnpvPpq8wpvR2vmYe3adSEnxybz\nSy5J7Y9kMnxGQuJdN0/td44qNxEYM8aLx2P48MPG/PhjVadDSnizZtUhLy+dTp0K6djR6WhUtGgy\nT2LJUg9s3x6uuSZAIOBh7NiWFZq3RWvmYW+/vRqAK6/U4fvJ8hmBw5N5PAbBazJXFfLgg3belm++\nqcucObWdDidhbd5cie+/r0FaWoCLLtIeQskkKwvq17ezJ65ZE/v9aTKPo2SqB9atC/ffb++PHXss\nBQXlS0RaM7emTWsAnMwZZwSorf8Tk+ozIhLfwUOazFWF/eEPHo47zs+vv2bw3nuNnA4n4RgDn3zS\nAICrr9aeQcmo6DwtsabJPI6SqR4Ith/t00/bt9CECVns3JlW5udqzRyWLKnJhg0Z1Kkzg0GDnI7G\nHZLtMxLPk6ARJXMR+aeIrBCRxSIyWURqRiswlRh+9zth0KBC9u5NY9y4indVTEUff2yPyv/v/wSv\nnvtMSt26gccDubmwf39s9xXRfOYiciowwxgTEJFHAYwxfy22js5nnuSWL4eOHQ3GwEsvLaBly71O\nh+R6+/Z5Ofvs3hw86GPVKjjuOKcjUrHSqRMsWQKzZkHfvmV/XlznMzfGfGaMCc3YOw9oEsn2VGJq\n1w6uvTZAICA880zFuiqmmhkz6nHwoI/evQs1kSe5eJVaolkzvwKYGsXtJZ1kqwcW9cADXmrW9JOb\nW4fZs4/eLSOVa+bGwPvv2xPGN97oTer3RXklY1v07GlvY30S9Kin0EXkM6BhCb+6yxjzYXCdu4F8\nY8wbJW1j5MiRZGVlAZCZmUl2dvahLkihP54uJ/7y3/8u3HJLDqNG7adHj6qkp5tDSTvUFbF4Ei/t\n98m8vHZtVdasGUDt2n7q1fuS3NxcV/z93LAcai+3xBONZY8HYABz5x55/ZycHMaPHw9wKF+WR8TX\nABWRkcDVwCnGmAMl/F5r5imioADat/ezapWXa65ZzYgR65wOyZUefvh4PvvsGP785wD//Kd2KEt2\ngQDUqmVnT1y3zk4nXRZxrZmLyBDgL8CZJSVylVrS0mD0aNst4/XXm7N9e9m7KqaKXbt8hCbVuv56\nTeSpwOOJT6kl0nfTM0A14DMRWSQiz0YhpqSVjPXA4oYMgSFD/Ozfn8bLL2eVul6q1synTWtIQYGX\nU04ppGVL+1gqvC/KKlnbIh4nQSPtzdLKGNPcGNM5+HNDtAJTietf//KSlmb4738b8d131ZwOxzUC\nAfjgA3tRj5tu0hGfqSQeyTzimvlRd6A185R0220BnnrKQ7t2OxgzZjGic0jxzTe1+POfO9GokZ+f\nfvLi03yeMrZts/MZVakCu3bZkuTRxLVmrlRp7r3XQ926fpYvr8WMGfWcDscVJk+2R+XXXCOayFNM\nnTp2YNj+/fDtt7HZhybzOErWemBJataERx+1b6/nn2/J/v2Hj1dPtZr5+vVVmDOnHmlphmuvPfxj\nl0rvi6NJ5raIdalFk7mKmcsvF7p0KWTbtipMmNDM6XAcNWlSI4wRRowI0LCkURsq6cU6mWvNXMXU\n/PnQs6fB6zWMGzefZs1iPNuQC+3Z4+O883px8KCPJUugQwenI1JO+OYbO/FWq1bw/fdHX19r5spV\nuneHyy8PUFjoYfToY1Ny3pYpU47h4EEfAwcWaiJPYR07QuXKsGqVPSEabZrM4yiZ64FH8uij4UvM\nffllHSB1auYFBcLkyXbI3+23l3zWM1XfFyVJ5rZIS7NH5gBffx397WsyVzFXrx489JC9P3bssRw4\nkDpvuy++qMf27ZVp06aQwYOdjkY5rXdvexuLy8ilzqfKBUKT66Si667z0KGDn82bq/Laa81T4hqg\nxsDEifao/M9/9pba1z6V3xfFJXtbxPIkqCZzFRc+H7z0khcRwzvvNOXHHzOcDinm5s2rzQ8/1KRe\nPT8XXaSjptTh1wQNBI68bnlpMo+jZK4HlkXPnnDVVQH8fg8PPLAx6m9mNzEGxo+33TFvv12oXLn0\ndVP9fVFUsrdFo0bQtKmdQXHlyuhuW5O5iqvHH/dSr56ftWurM3Vq8na4nj+/Nt99l0mdOn6dHVEd\nJlZ1c32XxVGy1wPLIjMTRo3yAAN44YWW7NiRfNPkGgOvvho+Ks84SkVJ3xdhqdAWsaqbazJXcTdi\nhDBwoJ+8vHRGjz7W6XCibsGCWqxcmUnt2n5uuEE/YupwmsyTQLLXA8tKBK688kuqVAmQk9OQWbPq\nOB1S1BStld92m1CtDDMA6/siLBXaonNn2+d82TJbO48WTebKEY0awcMP2/tPPdWKPXuSYxrBb76p\nxfLltcjM9HPTTfrxUr9VuTJ06WL/8c+fH73t6rstjlKhHlhWAwYM4KabPPTsWciOHZUZO7al0yFF\nzO+H555rAdhaefXqZXuevi/CUqUtQqWWaJ4E1WSuHOP1wmuv+UhPDzBtWiPmzq3ldEgRmTatAWvW\n1KBxYz+33KIfLVW6WNTN9R0XR6lQDyyrUFscfzw88IB97IknWrN7d2KWW/bv9zJunD0qf+wxD1Wq\nlP25+r4IS5W2KJrMozX5nCZz5bjbbvPQvbufbduq8MQTrRJyZsW33mrC9u2V6dy5kBEjdLSnOrLm\nzaFhQzt74urV0dmmzmeuXGHNGujUKUBenoe//GUFp5++yemQymzLlnQuvrgH+fk+vvwSTjzR6YhU\nIjjrLHj/fZgwAS655Le/1/nMVUJq2RLGjrXv22eeacX69eWoUzjsxRdbkJ/v46yz/JrIVZmFRoJ+\n9VV0tqfJPI5SpR5YFiW1xSWXCOed5+fAAR8PPNCGwkL3lyu+/roW06cfQ3p6gH/+03v0J5RA3xdh\nqdQW0R7WH3EyF5HbRCQgIrWjEZBKXSLw4otemjTx8/33NXn2WXd3V9y3z8sTT7QG4P774djkG8yq\nYqhbNzub6NKlsGdP5NuLqGYuIk2Bl4Djga7GmO0lrKM1c1Uuc+dCv36GwkLhzjuXM3jwZqdDKtGo\nUcfxwQdNyM4uZP58H77E7IijHNSjhx04NH06nHLK4b+Ld838KeD2CLeh1GF69YLRo+0BwJNPtmb1\n6jKMiY+zJUtq8sEHTfD5DOPHayJXFRPNUkuFk7mInAmsM8YsiTyM1JBK9cCjOVpbXHedh8su85Of\n7+Oee9q5qv/5/v0eHnvMllfuuMPQqVNk29P3RViqtUWfPvY2Gsn8iJ8QEfkMKGnS6buBO4GiVzUs\n9evAyJEjycrKAiAzM5Ps7OxDw3ZDfzxdTq3lkCOt//zzXmbPnsHq1V7+9reO/POf37J8+UKAQ5ed\nC10YOl7LixblMmFCM379NYM2bfz07/8lOTmRtUdubq7jfw+3LIfa2y3xxHoZ7PKcOQP4/PMcJkwY\nD3AoX5ZHhWrmInICMAPYF3yoCbAe6GGM2VxsXa2Zqwr75Rfo3t3Ppk1e+vXbxH33rcBbsU4jUTFp\nUmPGjm1F1aoB5s/30K6dc7GoxGcMNGkCv/4KK1ZAmzbh38WlZm6M+dYY08AY08IY0wJYB3QpnsiV\nilTTpvDpp16qV/fz5ZcNGD36OMdGiC5ZUpPnn7ddVl59VTSRq4iJRK9uHq1+5nroXQbFSwyprDxt\n0bEjTJniIT09wJQpTZgwoVnsAivF1q3p3HdfW/x+D7fe6mf48Oj1gdf3RVgqtkWobh7p4KGoJHNj\nTMuSuiUqFS0DBghvvCF4PIbx41vy8stZcTtC37Ejjdtu68DOnZXp39/P4487WOdRSSdaR+Y6N4tK\nKBMmGK64Avx+Ydiwddx88+qY1tB37kzj1ls7snZtddq08fO//3mpVy92+1Op5+BBqFEDCgpg+3Z7\nnVzQuVlUkrv0UuHddyE9PcCHHzbhwQfbkJ8fm2H/u3aFE3nr1oXk5GgiV9FXqRJ07WpPhs6bV/Ht\naDKPo1SsB5YmkrY480xh2jShevUAM2c25IYbOvPzz9GdmGvt2qrcdFMn1q6tTqtWhcyc6aNBg6ju\n4hB9X4SlaltEo9SiyVwlpAEDhP/9z0Pz5n5++KEG11zTlalTG0Sljv7JJw247rou/PJLNVq39jNz\npo+GJY22UCpKonESVGvmKqHt3g3XXlvIW2/Z8W89emzhqqt+olWrvHJva/v2dJ57rgXTpx8DwIUX\n+nnhBS/V3DebgEoyGzbYi5xXrw47dthLKpa3Zq7JXCWFCRMMN9xg2LvXftns23czl176M61a5SFH\n+Ths3FiJN99syiefHENBgZcqVQKMGQOXX+456nOVipaWLeHHH2HRIsjO1mTuajk5OUWG8aa2WLTF\n5s3w8MOFPP+8h4MHbVKvV28/3bvvoGvXHdSunU+lSgHS0wNs2FCZ5ctrsGxZdZYty8Tvt+ufcUYh\njz7qo23bqIZ2RPq+CEvltrjkEvj3v2HMGLjxxvInc/fMXqRUhOrXh1GjfNxxBzz8sJ+33oItW6ow\ndWoVpk5tVOrzvF7DiBGF3H23j/bt9SOhnNG3r03ms2fbZF5eemSuklYgALm58MknAf73vwC7dhn2\n74eDB4X69Q19+njo08dL795ol0PluKVL7WjnZs3gp5+0zKKUUgkpEIDatWHXLjvBXNOmOmjItVK1\nD21JtC3CtC3CUrktPJ5wf/PZsyvw/OiGo5RSqqL69rW3FUnmWmZRSimX+OILGDgQunSBhQu1Zq6U\nUglp716oWdPO0xIIaM3ctVK5HlictkWYtkVYqrdFRgZ07mxPhpaXJnOllHKRUN28vLTMopRSLvLO\nOzB8OICWWZRSKmFV9Mhck3kcpXo9sChtizBtizBtCzt7YqtW5X+eJnOllHKZ5cvL/xytmSullAvp\nNUCVUioFRZTMReQmEVkhIt+KyGPRCipZaT0wTNsiTNsiTNui4iqczEXkZOAMoKMx5gTgiahFlaRy\nc3OdDsE1tC3CtC3CtC0qLpIj8+uBR4wxBQDGmC3RCSl57dy50+kQXEPbIkzbIkzbouIiSeatgP4i\nMldEckSkW7SCUkopVT5HvEaWiHwGNCzhV3cHn1vLGNNLRLoDE4GW0Q8xeaxdu9bpEFxD2yJM2yJM\n26LiKtw1UUQ+AR41xswMLq8GehpjthVbT/slKqVUBcTrgs7vAwOBmSLSGkgvnsjLG4xSSqmKiSSZ\nvwK8IiJLgXzg0uiEpJRSqrxiPgJUKaVU7MV0BKiIDBGRlSKySkTuiOW+3ExEmorIFyKyLDjA6o9O\nx+QkEfGKyCIR+dDpWJwmIpkiMik4+G65iPRyOianiMidwc/IUhF5Q0QqOR1TvIjIKyKyKVjpCD1W\nW0Q+E5HvReRTEck80jZilsxFxAuMAYYA7YARItI2VvtzuQLgVmNMe6AXcGMKtwXAzcByQL8Wwr+A\nqcaYtkBHYIXD8ThCRLKAq4EuxpgOgBe4wMmY4uxVbK4s6q/AZ8aY1sCM4HKpYnlk3gNYbYxZGxxY\n9BZwZgz351rGmI3GmNzg/TzsB7aRs1E5Q0SaAKcDLwMpfXJcRGoC/YwxrwAYYwqNMbscDsspu7EH\nPVVFxAdUBdY7G1L8GGO+BHYUe/gM4LXg/deA3x9pG7FM5o2BX4osrws+ltKCRyCdgXnORuKYp4G/\nABW4ymHSaQFsEZFXRWShiLwkIlWdDsoJxpjtwJPAz8CvwE5jzHRno3JcA2PMpuD9TUCDI60cy2Su\nX6GLEZFqwCTg5uARekoRkaHAZmPMIlL8qDzIB3QBnjXGdAH2cpSv0slKRI4FbgGysN9aq4nIRY4G\n5SLBecSPmFNjmczXA02LLDfFHp2nJBFJA94F/m2Med/peBzSBzhDRH4E3gQGisgEh2Ny0jpgnTFm\nfnB5Eja5p6JuwFfGmG3GmEJgMvb9kso2iUhDABE5Bth8pJVjmcwXAK1EJEtE0oHzgSkx3J9riYgA\n44DlxphRTsfjFGPMXcaYpsaYFtiTW58bY1J2fIIxZiPwS3DQHcAgYJmDITlpJdBLRKoEPy+DsCfJ\nU9kU4LLg/cuwAzVLFcmgoSMyxhSKyB+Aadgz0+OMMSl5ph7oC1wMLBGRRcHH7jTG/NfBmNxAS3Fw\nE/Cf4AHPD8DlDsfjCGPM4uC3tAXY8ykLgRedjSp+RORN4CSgroj8AtwLPApMFJErgbXA8CNuQwcN\nKaVU4tPLximlVBLQZK6UUklAk7lSSiUBTeZKKZUENJkrpVQS0GSulFJJQJO5UkolAU3mSimVBP4f\nbn7mmC/kIgwAAAAASUVORK5CYII=\n", 245 | "text": [ 246 | "
**\u00a1Importante!**: F\u00edjate que la funci\u00f3n del sistema recibe como primer argumento $\\mathbf{y}$ (un array) y como segundo argumento el instante $t$ (un escalar). Esta convenci\u00f3n va exactamente al rev\u00e9s que en MATLAB y nos equivocamos obtendremos errores o, lo que es peor, resultados incorrectos.
"
324 | ]
325 | },
326 | {
327 | "cell_type": "markdown",
328 | "metadata": {},
329 | "source": [
330 | "Condiciones iniciales:"
331 | ]
332 | },
333 | {
334 | "cell_type": "code",
335 | "collapsed": false,
336 | "input": [
337 | "y0 = 1"
338 | ],
339 | "language": "python",
340 | "metadata": {},
341 | "outputs": [],
342 | "prompt_number": 9
343 | },
344 | {
345 | "cell_type": "markdown",
346 | "metadata": {},
347 | "source": [
348 | "Vector de tiempos donde realizamos la integraci\u00f3n:"
349 | ]
350 | },
351 | {
352 | "cell_type": "code",
353 | "collapsed": false,
354 | "input": [
355 | "t = np.linspace(0, 3)"
356 | ],
357 | "language": "python",
358 | "metadata": {},
359 | "outputs": [],
360 | "prompt_number": 10
361 | },
362 | {
363 | "cell_type": "markdown",
364 | "metadata": {},
365 | "source": [
366 | "Integramos y representamos la soluci\u00f3n:"
367 | ]
368 | },
369 | {
370 | "cell_type": "code",
371 | "collapsed": false,
372 | "input": [
373 | "sol = odeint(f, y0, t)\n",
374 | "plt.plot(t, sol)"
375 | ],
376 | "language": "python",
377 | "metadata": {},
378 | "outputs": [
379 | {
380 | "metadata": {},
381 | "output_type": "pyout",
382 | "prompt_number": 11,
383 | "text": [
384 | "[El c\u00f3digo que tienes a continuaci\u00f3n est\u00e1 casi resuelto, busca `#completar#` y termina de resolver el ejercicio.
"
391 | ]
392 | },
393 | {
394 | "cell_type": "markdown",
395 | "metadata": {},
396 | "source": [
397 | "En primer lugar, cargamos los widgets interactivos:"
398 | ]
399 | },
400 | {
401 | "cell_type": "code",
402 | "collapsed": false,
403 | "input": [
404 | "from IPython.html.widgets import interact, interactive, fixed\n",
405 | "from IPython.html import widgets"
406 | ],
407 | "language": "python",
408 | "metadata": {},
409 | "outputs": [],
410 | "prompt_number": 5
411 | },
412 | {
413 | "cell_type": "code",
414 | "collapsed": false,
415 | "input": [
416 | "def #completar#(y, t, sigma, gamma, f, beta):\n",
417 | " \"\"\" Modelo simple epidemiol\u00f3gico para la\n",
418 | " propagacion del virus del \u00c9bola\"\"\"\n",
419 | " \n",
420 | " # Recuperamos las variables dependientes de entrada\n",
421 | " S, E, I, R, D = y\n",
422 | " \n",
423 | " Nn = S + E + I + R # Total de personas vivas\n",
424 | " \n",
425 | " # F\u00edjate que par\u00e1metros de las ecuaciones diferenciales \n",
426 | " # son son argumentos extra de la funci\u00f3n en esta ocasi\u00f3n\n",
427 | " \n",
428 | " # Definimos las ecuaciones diferenciales\n",
429 | " dSdt = -beta*(S/Nn)*I\n",
430 | " dEdt = +beta*(S/Nn)*I - sigma*E\n",
431 | " dIdt = sigma*E - gamma*I\n",
432 | " dRdt = (1.0 - f)*gamma*I\n",
433 | " dDdt = f*gamma*I\n",
434 | " \n",
435 | " return [#completar]\n",
436 | "\n",
437 | "def ebola(tiempo, recortes_sanidad=2):\n",
438 | " \"\"\" Solve the model and plot a basic graph of the results \"\"\"\n",
439 | " \n",
440 | " tiempo_simulacion = tiempo # d\u00edas\n",
441 | " t = np.linspace(0., tiempo_simulacion, 1000)\n",
442 | "\n",
443 | " #--------------------------\n",
444 | " # Par\u00e1metros del modelo\n",
445 | " #--------------------------\n",
446 | " \n",
447 | " sigma = 1./9 # Velocidad convertirse en contagioso una vez expuesto\n",
448 | " gamma = 1./8.5 # Velocidad de recuperaci\u00f3n o fallecimiento\n",
449 | " f = 0.4 # Fracci\u00f3n de gente que no supera la enfermedad (40%)\n",
450 | " \n",
451 | " # A mayor n\u00ba de recortes, mayor n\u00ba de personas contagiadas por infectado\n",
452 | " R_0 = #completar#\n",
453 | " \n",
454 | " beta = R_0*gamma # Velocidad con la que se propaga la infecci\u00f3n\n",
455 | " \n",
456 | " # Condiciones iniciales\n",
457 | "\n",
458 | " N = 350 # Poblaci\u00f3n total\n",
459 | " \n",
460 | " y0 = [N, 1, 0, 0, 0] # [S, E, I, R, D]\n",
461 | " \n",
462 | " # Llamamos al ode para resolver las eq. diferenciales \n",
463 | "\n",
464 | " solution = integrate.odeint(#completar#)\n",
465 | "\n",
466 | " S = solution[:,0] # S es el n\u00ba de personas sin infectar\n",
467 | " E = solution[:,1] # E es el n\u00ba de portadores del virus no contagiosos (a\u00fan)\n",
468 | " I = solution[:,2] # I es el n\u00ba de personas sintom\u00e1ticas e infecciosas\n",
469 | " R = solution[:,3] # R es el n\u00ba de personas recuperadas\n",
470 | " Ds = solution[:,4] # D es el n\u00ba de personas fallecidas\n",
471 | " \n",
472 | " # Creamos la figura para representar los resultados\n",
473 | " \n",
474 | " fig = plt.figure(figsize=(10, 5))\n",
475 | " \n",
476 | " # Representamos las curvas de evoluci\u00f3n con respecto al tiempo\n",
477 | " plt.plot(t, (E + I)/N*100 )\n",
478 | " plt.plot(t, I/N*100)\n",
479 | " plt.plot(t, R/N*100)\n",
480 | " plt.plot(t, Ds/N*100)\n",
481 | " plt.plot(t, S/N*100)\n",
482 | "\n",
483 | " # A\u00f1adimos leyenda a los datos\n",
484 | " plt.legend([\"% infectados\",\n",
485 | " \"% infectados contagiosos\", \n",
486 | " \"% recuperados\", \n",
487 | " \"% fallecidos\",\n",
488 | " u\"% poblaci\u00f3n sana \"], loc=1)\n",
489 | "\n",
490 | " #T\u00edtulo y anotaciones al gr\u00e1fico \n",
491 | " plt.suptitle(u'Variaci\u00f3n de la velocidad de propagaci\u00f3n del virus del \u00c9bola:',\n",
492 | " fontsize=14,\n",
493 | " fontweight='bold')\n",
494 | " \n",
495 | " plt.title(u\"Impacto de los recortes en el sistema sanitario.\",\n",
496 | " fontsize = 12)\n",
497 | " \n",
498 | " plt.text(10, 50, u\"A mayor n\u00ba de recortes, mayor \\n n\u00ba de contagios (accidentales) \\n y mayor velocidad de propagaci\u00f3n.\",\n",
499 | " style='italic',\n",
500 | " bbox={'facecolor':'red', 'alpha':0.5, 'pad':10})\n",
501 | " \n",
502 | " plt.xlabel(u\"tiempo / d\u00edas\")\n",
503 | " plt.ylabel(u\"porcentaje de poblaci\u00f3n\")\n",
504 | "\n",
505 | " plt.show()"
506 | ],
507 | "language": "python",
508 | "metadata": {},
509 | "outputs": [
510 | {
511 | "ename": "SyntaxError",
512 | "evalue": "invalid syntax (