Skip to content

Instantly share code, notes, and snippets.

@oroszl
Created August 6, 2017 14:37
Show Gist options
  • Save oroszl/9cb320b6bc2d92a473098c1337a3fbf6 to your computer and use it in GitHub Desktop.
Save oroszl/9cb320b6bc2d92a473098c1337a3fbf6 to your computer and use it in GitHub Desktop.
kpm finite size conductivity
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Populating the interactive namespace from numpy and matplotlib\n"
]
}
],
"source": [
"%pylab inline\n",
"import pybinding as pb\n",
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"pb.pltutils.use_style()"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"S0=array([[1,0],[0,1]])\n",
"S1=array([[0,1],[1,0]])\n",
"S2=array([[0,-1j],[1j,0]])\n",
"S3=array([[1,0],[0,-1]])\n",
"N=zeros_like(S0);"
]
},
{
"cell_type": "code",
"execution_count": 40,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"Tx=(S3+1j*S1)/2\n",
"Ty=(S3+1j*S2)/2\n",
"U=-1.5*S3"
]
},
{
"cell_type": "code",
"execution_count": 41,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"lat=pb.Lattice(a1=[1,0],a2=[0,1])\n",
"lat.add_sublattices(('A', [0,0],U ))\n",
"lat.add_hoppings( \n",
" ([1,0], 'A', 'A', Tx),\n",
" ([0,1], 'A', 'A', Ty),\n",
" )\n",
"model = pb.Model(lat,pb.primitive(a1=1000, a2=200))\n",
"kpm = pb.chebyshev.kpm(model, kernel=pb.jackson_kernel(),silent=True,energy_range=[-4, 4])"
]
},
{
"cell_type": "code",
"execution_count": 42,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"dos=kpm.calc_dos(energy=linspace(-5,5,500),broadening=0.01,num_random=1)"
]
},
{
"cell_type": "code",
"execution_count": 43,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYwAAAD+CAYAAAAppDI0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmYHFW5+PFv9Tr7ZA+ZLCYQOIRECAZCEBUSLhBFJShC\nUCFoLhIBAVkUkAE1oOAVcUG5PzYBrwJhk6CBiAQVFGRHCOGFhASykXX2tZf6/VGnZzqTSTJJuru6\ne97P8/QzNafrVL89mczbZ6lzHNd1UUoppXYl4HcASimlCoMmDKWUUn0SytSFjDGjgXuA4YAL3Coi\nvzDGDALuB8YCq4BTRaTO1rkCmAskgAtEZLEtnwLcBZQCi4ALRcQ1xkTta0wBtgCnicgqW2cOcJUN\n51oRuTtT700ppVRmWxhx4BIROQiYBpxnjDkIuBx4SkT2B56y32Ofmw1MBGYCvzHGBO21bgHOBva3\nj5m2fC5QJyLjgZuAG+y1BgHXAEcAU4FrjDEDM/jelFKq38tYwhCR9SLyij1uApYBI4GTgNSn/buB\nWfb4JOA+EekQkZXAcmCqMWYEUCUiz4uIi9eiSK+TutaDwLHGGAc4AXhSRLba1suTdCcZpZRSGZCx\nLql0xpixwKHAv4HhIrLePvUhXpcVeMnk+bRqa2xZzB73LE/VWQ0gInFjTAMwOL28lzpdLr30Ure0\ntBSA6dOnM2PGjD17g0op1X84qYOMJwxjTAXwEHCRiDQaY7qes+MQvs3jLS0tZf78+X69vFJKFbSM\nzpIyxoTxksXvReRhW7zBdjNhv2605WuB0WnVR9mytfa4Z/k2dYwxIaAab/B7R9dSSimVIRlLGHYs\n4Q5gmYj8LO2phcAcezwHeDStfLYxJmqMGYc3uP2C7b5qNMZMs9c8s0ed1LVOAZbYcY7FwPHGmIF2\nsPt4W6aUUipDMtkldRRwBvCGMeY1W3YlcD2wwBgzF3gfOBVARJYaYxYAb+HNsDpPRBK23rl0T6t9\n3D7AS0i/M8YsB7bizbJCRLYaY+YDL9rzfigiWzP43pRSqt9z+tPSILW1ta6OYSil1G7pGvTWO72V\nUkr1iSYMpZRSfaIJQymlVJ9owlBKKdUnmjCUUhkRW/oWnW+84XcYKouysjSIUqr/2Xj8CQCMXLt6\nF2eqQqUtDKWUUn2iCUMppVSfaMJQSinVJ5owlFJK9YkmDFWwOl9+hfZ/PON3GEr1GzpLShWsTZ8/\nCdBZOUrlirYwlFKZEY36HYHKMk0YSqmMcGzC6E8rYPc3mjCUUhnhRCLeQXu7v4GorMnYGIYx5k7g\ns8BGEZlky+4HUpt6DwDqRWSyMWYssAwQ+9zzIjLP1plC9+ZJi4AL7V7gUeAeYAretqynicgqW2cO\ncJW91rUicnem3pdSqm9SCSPZ2kqwtNTnaFQ2ZHLQ+y7gZrw/6gCIyGmpY2PMjUBD2vkrRGRyL9e5\nBTgb+DdewpiJt+PeXKBORMYbY2YDNwCnGWMGAdcAhwEu8LIxZqGI1GXwvSmldsUmDLe1FQYP9jkY\nlQ0Z65ISkX/gbZu6Hbs396nAvTu7hjFmBFAlIs/bvbrvAWbZp08CUi2HB4Fj7XVPAJ4Uka02STyJ\nl2RUEdN+8vzjlNgxjJYWnyNR2ZKrabWfBDaIyLtpZePs3t8NwFUi8gwwEliTds4aW4b9uhpAROLG\nmAZgcHp5L3W2sW7dOmprawGYPn06M2bM2Nv3pfyS6icv0Zk5+cIJBAFwW9t8jkRlS64Sxuls27pY\nD4wRkS12zOKPxpiJ2Q6ipqYG3dO7OCQbGwEIlJb5HIlKcRNxAJLawihaWZ8lZYwJAV8A7k+ViUiH\niGyxxy8DK4ADgLXAqLTqo2wZ9uvotGtW4w1+d5X3UkcVqVTCcMo0YeSNeAIAV2dJFa1cTKv9L+Bt\nEenqajLGDDXGBO3xvsD+wHsish5oNMZMs+MTZwKP2moLgTn2+BRgiR3nWAwcb4wZaIwZCBxvy1QR\nSzZ4CYOwLlaQL1ItDOIxfwNRWZOxhGGMuRd4zjs0a4wxc+1Ts9l+sPtTwH/sGMaDwDwRSQ2Ynwvc\nDizHa3k8bsvvAAYbY5YDFwOXA9h684EX7eOHaddSxaqz0/sai/sbh+qWamF0asIoVhn7eCYip++g\n/Kxeyh4CHtrB+S8Bk3opbwe+tIM6dwJ37ka4qsC5MS9huPppNm+48VQLQ5N4sdI7vVVB6voUqy2M\n/GETRSqZq+KjCUMVppiXMFz9NJs3uv4ttEuqaGnCUAXJ7ezwDjr102zeSNgxDE3iRUsThipIqS4p\n/eOUR1JdUprEi5YmDFWYbJcU8bguE5In3Hgcp7RUB72LmCYMVZDcWGf3siD6Byo/2IThxnQMo1hp\nwlCFqTNGoKwc0G6pfOAmk+C6XgtDu6SKliYMVZDcWKx7WRD9A+U/m7Sd0lJN4EVME4YqSG5nZ1fC\n0D9Q/nPtDCmnpEQHvYuYJgxVmGIxnLLSrmPlM/tv4JSW6r9HEdOEoQqSG4vhlGoLI1+4dh0pp7RE\n/z2KmCYMVZDczk4C5d6gt95ZnAcS3WMYOqZUvDRhqMKU1iWlCxDmgfRBb+2SKlqaMFRBSh/01gUI\n/edukzD036NYacJQBcntjHndH2gLIy+kJQx0tdqilbH9MIwxdwKfBTaKyCRb9n3gbGCTPe1KEVlk\nn7sCmAskgAtEZLEtnwLcBZQCi4ALRcQ1xkSBe4ApeFuzniYiq2ydOcBV9jWuFZG7M/W+VJ6KdeJE\nIhAK6SfaPNA1rVbvwyhqmWxh3AXM7KX8JhGZbB+pZHEQ3k58E22d36S2bAVuwUsy+9tH6ppzgToR\nGQ/cBNxgrzUIuAY4ApgKXGO3alVFzI3FIBzGCYf1E20+SLUwSkp0EkIRy1jCEJF/AH3dGvUk4D4R\n6RCRlXjbsU41xowAqkTkebtf9z3ArLQ6qZbDg8Cxdt/vE4AnRWSriNQBT9J74lJFxO2M4USjEA5r\nCyMPdE+rLdUNlIpYxrqkduJbxpgzgZeAS+wf9ZHA82nnrLFlMXvcsxz7dTWAiMSNMQ3A4PTyXups\nY926ddTW1gIwffp0ZsyYsXfvTPkn1okTDuOEQrr4YD6w40hONKoJvIhlO2HcAswHXPv1RuDrWX7N\nHaqpqWH+/Pl+vbzKoFSXFOGQDnrnATeegFDIdhHqv0exyuosKRHZICIJEUkCt+GNMQCsBUannTrK\nlq21xz3Lt6ljjAkB1XiD3zu6lipmnTGcSAQnWoLb3uF3NCoR91p7kbB2SRWxrCYMOyaRcjLwpj1e\nCMw2xkSNMePwBrdfEJH1QKMxZpodnzgTeDStzhx7fAqwxI5zLAaON8YMtIPdx9syVcTcVJdUaQlu\nW5vf4ahY3LYwInpfTBHL5LTae4FjgCHGmDV4M5eOMcZMxuuSWgWcAyAiS40xC4C3gDhwnogk7KXO\npXta7eP2AXAH8DtjzHK8wfXZ9lpbjTHzgRfteT8Ukb4OvqsC5cZiEAnjlJXhtrX6HU6/l2xuIlBZ\naac5a5dUscpYwhCR03spvmMn518HXNdL+UvApF7K24Ev7eBadwJ39jlYVfg6YzjhCE5JKW5bu9/R\n9HvJhkac6iqciI5hFDO901sVJLfTdkmVleG2agvDb25jI4GqKpxQWFsYRUwThipIbqwTImECpaU6\nhpEHko2NBKqq7aC3JoxipQlDFSS3pRWnrBynrFRbGHkg2dCIU1XpLdfSobPWipUmDFVw3GQSt7mZ\nQGWFHfTWFobf3MYGAlVVBKqqcdvbcTVpFCVNGKrguK2t4Lo4FZXeUhStmjD8lrRjGIGBA7zv6+t9\njkhlgyYMVXDcpiYAAlWVOGVlJLWF4btkozdLKjDAJoy6Op8jUtmgCUMVnKRNGE5lqoWhYxh+cxu8\nQW+ntBSnpEQTRpHShKEKjtvUDMEgTkmJlzDatYXht1SXFEBg4ECSddolVYw0YaiCk2xu8loXjqMt\njDyRbGzE2SZhaAujGGnCUAXHbbTLUACBsjK909tnbjKJ29REoNpLGM6AATroXaQ0YaiCk2xuxqmo\nAOyGPTro7Su3qQlct0eXlLYwipEmDFVw3KYmAlVeC0OXBvFfsrERAKe6GtCEUcw0YaiC47UwUgnD\nG8Nwk0mfo+q/kg1ewujqJqyswG1u8TMklSWaMFTBSTY2Eqj0uqQCVdXgurgt+gfKL25jI05Zmbfb\nHuCUl5Ns1X+PYqQJQxUct6GBgO3+cOxAa6pbROVesrGha4YU2G5CTeBFKZMbKN0JfBbYKCKTbNn/\nAJ8DOoEVwNdEpN4YMxZYBoit/ryIzLN1ptC9gdIi4EIRcY0xUeAeYAre1qynicgqW2cOcJW91rUi\ncnem3pfKP8mGBsL77APQNdDq1jfAyJF+htVvJRsau2ZIAQTKy3FbdFypGGWyhXEXMLNH2ZPAJBE5\nGHgHuCLtuRUiMtk+5qWV3wKcjbdt6/5p15wL1InIeOAm4AYAY8wgvN39jsDbM/wau1WrKlLJhgYc\nuwRFak+MZGODz1H1X25qaXPLKSsjqS2MopSxhCEi/8DbOjW97C8iktrg93lg1M6uYfcArxKR5+1+\n3fcAs+zTJwGplsODwLF23+8TgCdFZKuI1OElqZ6JSxWRZH09gQFpf6CqqrRLykfpN+0BOOW6bW6x\nyuUYxtfp3p8bYJwx5jVjzN+NMZ+0ZSOBNWnnrLFlqedWA9gk1AAMTi/vpY4qQsn6hq5F7gAC1VVd\nM3VU7vXsknLKtEuqWGVsDGNnjDHfA+LA723RemCMiGyxYxZ/NMZMzHYc69ato7a2FoDp06czY8aM\nbL+kygK3vr5r0Bu8mVKutjB8k9oLI8Ubw2jBdV0cx/ExMpVpWU8Yxpiz8AbDj7XdTIhIB9Bhj182\nxqwADgDWsm231Shbhv06GlhjjAkB1XiD32uBY3rU+VtvsdTU1DB//vxMvC3lE7e9Hbe9fZuEoV1S\n/ko2NhIeMaLre6e8DJJJaG+H0lIfI1OZltUuKWPMTOA7wOdFpDWtfKgxJmiP98Ub3H5PRNYDjcaY\naXZ84kzgUVttITDHHp8CLLEJaDFwvDFmoB3sPt6WqSKUbPAGt7fvktJBb7+k9sJIccrLvXK9A7/o\nZHJa7b14n/SHGGPW4M1cugKIAk8aY6B7+uyngB8aY2JAEpgnIqkB83Ppnlb7ON3jHncAvzPGLMcb\nXJ8NICJbjTHzgRfteT9Mu5YqMqnE4KQnjKoq7ZLyUWovjBSnzEsYbksLDB7sV1gqCzKWMETk9F6K\n79jBuQ8BD+3guZeASb2UtwNf2kGdO4E7+xysKljJrVshGsUpKekqC1RVEVu/3seo+rf0vTAAAuVl\nAHrzXhHSO71VQYmvWkVo7Ee2GUx1qqt1DMNHPafVEolAMEhSZ0oVHU0YqqDEl68gtN9+25QFqqtx\ndVqtL3ruhQF4G1uVl+PqelJFRxOGKijx5cu3Txg6S8o3PffCSAkOH05i9Zod1FKFShOGKiix5SsI\njx+/TZlOq/VPz70wUiIfO5TOV17xIySVRZowVMFwOztJfPABofE9u6SqcJuadE8MH/TcCyMlMmUK\nnS9rwig2mjBUwYi//z4kEr12SeG6XveIyim3sRGntLRrL4yUyJSPEX/3Xd3bu8howlAFI758OYHh\nw7b7NOvYewC0Wyr3klu2EBgyZLvy0AEH4FRW0vnKqz5EpbJFE4YqGN4MqfHblaf299a7vXMvsWlT\nrwnDCQS8cYyXX/YhKpUtmjBUweh8/XXCEyZsV+6EQjgVFTq11geJjRsJDhva63PhSZOIvf12jiNS\n2aQJQxUENxaj45lnKTnm6F6fD1RVaX+5D5KbNxMY0nvCCNbUkFi3LscRqWzShKEKQuw/b+B2dBA5\nclqvzwfHjPYGxVVOJTZu2mELIzhyJIl1umRLMdGEoQpC4sMPCQ4dSmAHy2WH9t2X+IoVOY5KJTdt\nJDC094QRqqkhuXkzbnt7jqNS2aIJQxWEZF0dgYE73qo9tN++xFe8l8OIFNgWxg4SRrDG2yMjoQtD\nFg1NGKog7DJh7LuftjByzE0kSG7eTHD48F6fdwYMwCkr026pIqIJQxUEL2EM2OHzkUmTSG7dSnzl\nyhxG1b8lPtwA8TjBUSN7fd5xHMITJtD5qt6LUSw0YaiCsKsWRrBmBKEDDe1P/y13QfVzibVrIBLZ\n4RgGQHT6MbQ//XQOo1LZlMkd9+7E27t7o4hMsmWDgPuBscAq4FQRqbPPXQHMBRLABSKy2JZPoXvH\nvUXAhSLiGmOiwD3AFLy9vE8TkVW2zhzgKhvKtSJyd6bel8oPybo6wgcdtNNzSo4+mo5nn6Xi61/L\nUVT9W2LtWoI1I3ACO/7cWXLsDJpu+jmJLVsI6u57BS+TLYy7gJk9yi4HnhKR/YGn7PcYYw7C22J1\noq3zm9Qe38AtwNl4+3zvn3bNuUCdiIwHbgJusNcahLcd7BHAVOAau7e3KiLJuvqdtjAAwh+dRGyZ\n3iiWK4k1awmNHLXTc8If/SjB0aNoe+xPOYpKZVPGEoaI/ANvr+10JwGpT/t3A7PSyu8TkQ4RWQks\nB6YaY0YAVSLyvIi4eC2KWb1c60HgWGOMA5wAPCkiW23r5Um2T1yqwO2qSwogfOCBJD74gGRzc46i\n6t/iq1cTHFmz03Mcx6H0+OPpePbZHEWlsinbYxjDRSQ1ReJDIDWdYiSwOu28NbZspD3uWb5NHRGJ\nAw3A4J1cazvr1q2jtraW2tpalixZsqfvSeWY67okNm8mMHjQTs8L7bcfhMPE3pYcRda/xV59jfDB\nB+/yvODo0d4AuSp4GRvD2BU7DuHm6vV6U1NTw/z58/0MQe2BxJo1uA0NhA88cKfnOZEIkYMPpuPp\np4keNiVH0fVPyYYGYsuWEZ06dZfnBocPJ7lBE0YxyHYLY4PtZsJ+3WjL1wKj084bZcvW2uOe5dvU\nMcaEgGq8we8dXUsVuLbFi2n5w73EXnudwPBhBEeM2GWdsi/PpuXee3UzpSxr+9OfcaqrCB1odnlu\nYPhwEhs36r9JEdhpwjDGhI0xhxpjhu3h9RcCc+zxHODRtPLZxpioMWYc3uD2C7b7qtEYM82OT5zZ\no07qWqcAS+w4x2LgeGPMQDvYfbwtUwWu6ee/oP6y79B89z1EDjmkT3VKPnU0yQ0b9RNtFrmuS+OP\nr6fqsktxgsFdnh/cZzjE4yS39hziVIVmm4RhjPlfY8xEe1wNvI438PyqMeb0nV3IGHMv8Jx3aNYY\nY+YC1wPHGWPeBf7Lfo+ILAUWAG8BTwDniUjCXupc4Ha8gfAVwOO2/A5gsDFmOXAxdsaViGwF5gMv\n2scPbZkqcE55BQCdzz1HyOz6kyxAYJ/hEI3qQoRZlPzwQ5J1dZQef0Kfzg8O8z5v6jhG4es5hvFJ\nEZlnj78GvCMis4wx++D94b53RxcSkR0llGN3cP51wHW9lL8ETOqlvB340g6udSdw545iU4UpuXVL\n13Fo3Lg+1XECAUJjxhB//wOi03pf2Vbtnfh7K3FKS73k3AdOJEJg8GBvqfNJE7Mcncqmnl1SnWnH\nxwF/BBCRD3MWkVJWcvMWQuO9HfZC+/YtYQAEx4whoS2MrIm/9x6hceN2esNeT+GDP0rniy9mMSqV\nCz3/xeuNMZ81xhwKHIXXXZQaZO59XWmlssBNJknW1RH52KFA31sYAKGxHyH+nq5cmy2xFSsI7bfv\nbtUpmTGDdp3KXvB6JoxzgPPx7tq+KK1lcSzw5xzGpfq5ZH09JJNEj/4UgX32IbAby0pEDjuMjn89\nh5tI7PpktdtiS9/q85hSSsn0Y4i/LcTX6gTGQrbNGIaIvEMvd0nbdZ505pHKidZHHqHj+RcAKDnu\nOEpPPBHHcfpcv2T6MdQ1NtL5yitEDz88W2H2S24iQez116k8/9zdqhcaN47QvvvSseRpQmd8NUvR\nqWzb7sY9Y8yn8WYgpUanlgI3iMiiXAam+q+2hY/R/pcnoSSKU1a2W8kCIFBZSfijH6XzlVc1YWRY\n/N13cVtaiEyevNt1ozOm0/7005RrwihY2yQMY8zZeN1S3wFessWHAdcbY0aJyK05jk/1Q4mNGyk/\n+7+JTJy428kiJWwOIP7uuxmOTHW+8iqh8eMJVFfvdt2SY2ewde7ZuB0dONFoFqJT2dazhfFt4BM9\n7mNYYlsdzwKaMFTWJT7cQHTqVEo/8+k9vkZo//1p/7M2ijOt89VXiRy6+60LgOgRR4Dr0v7kXyn9\n7IkZjkzlQs9Bb6e3m95EZEvPMqWywU0mSW7eTGDYni4u4AkfsD+xd9/Vge8M63z1VcKHHrpHdZ1o\nlMoLL2DrOfOo+/bFGY5M5ULPhNFojNluDQZb1pSbkFR/lty61dv2c/jeJYzIxz6GE43S9KubMxSZ\nSjY3E5d3uqY674mK88+j6ntX0vbEYk3mBahnwrgEWGiM+b4x5nP28QO89Zz0I4HKKrezk43TvYUB\ngnvZwghUV1P1nctove9+XNfXRZKLRuz1/0AkvMtVg3fGcRzKzzwDt7XVu54qKNskDBF5Fm/nugBw\nln0EgGn2OaWypv2vT3UtUJeJQdGSY44hsXo1iVWr9vpayo5fHHwwTji8V9cJVFQQmfIxWu5fQEx0\n75JCst20Wnuz3tXGmKH2+005j0r1S22LFlE2+zQqL/hWRq4XrBlBcNw4Ol54YbfuFFe963z9P4T7\nuGrwrkQ/+Umafnojrf/3f4QnTmTonx/b60Sksq/narWO7Y7aBAggxphNxpir/QlP9QfJtjbcjg7i\n77xL5JBDCH3kIxm7dni//Uis0nWlMiEuQnjChIxcq+SYo7uOY0uXktA7wAtCzzGMb+OtITVVRAaJ\nyCC8LqqjjDHfznl0ql+ov+hiGuZfS3zFCm+b1QwKjhlNfPXqXZ+odsptbye+ciXhPmyY1BeRQw9l\n6OPdqw21LXqcZEtLRq6tsqdnwjgDOF1EVqYKROQ94Kt4mxkplVHJhgba//Y3Wv5wL257O6HxmU0Y\nodGjSXygCWNvxZYvB9cldMABGbtmeOJEKPHGqhqv+xHNt92esWur7Og5hhEWkc09TxKRTcaYPepg\nNMYY4P60on2Bq4EBwNlAaozkytTyI8aYK4C5QAK4wK5lhTFmCt7CiKXAIuBCu1d4FG+jpyl427ae\nJiKr9iRelTvxtevYMPWIbcr29v6LnoIfGUP8gw8yes3+xk0mabzhJ0QOO4xAaeYWrXaCQcpPPZWW\ne34HQPsTi6m66MKMXV9lXs+E0dnrWbt+bodERIDJAMaYIN5+24/gbdB0k4j8NP18Y8xBwGy8taxq\ngL8aYw6wO/Ldgpdk/o2XMGbibew0F6gTkfHGmNnADcBpexKvyg23rY36K64kMGwYA274MfEV7xE5\n5JA9XgpkR0JjxpDctIlkUxOBysqMXru/aH3gATpfepnhT/0149ce8OMfkdi8hfZFi4i9+w5bv3Uh\nA3/584z/HqjM6JkwDjHGNPZyngOUZOD1jgVWiMj7ZsfLI58E3CciHcBKuyXrVGPMKqBKRJ4HMMbc\nA8zCSxgnAd+39R8EbjbGOHbPb5WHWh/7E7HXX2foIw/v1uZIuys0fjxOaSlNv7qZim+cTXDIkKy9\nVrFquetuKuedQ7BmRFauP+jmX+J2/pTEho1s+tznaVv4GGUnfT4rr6X2Ts/lzXe9o/vemc2227x+\nyxhzJt5Ch5eISB0wEng+7Zw1tixmj3uWY7+uBhCRuDGmARgMbNO9tm7dOmprawGYPn06M2bMyNDb\nUrsrtnQp0WnTsposAJxwmGBNDc2//g3JLVsYeONPd11JdUmsW0/sP28w8Fe/ytprONEoTjRKoLKS\n0k/PpONvf9OEkaf6vsfiXjLGRIDPAw/YolvwxjMmA+uBG7MdQ01NDfPnz2f+/PmaLHzkui7xZW8T\nmrDndwzvjshUb4nz1ocfIdncnJPXLBZNt95KaMIEwhmejLAj0SOPpHXBAzTfeltOXk/tnpwlDODT\nwCsisgFARDaISEJEksBtwFR73lpgdFq9UbZsrT3uWb5NHbudbDXe4LfKQ3XnnkfHP/9J+KDMzOnf\nlQHX/5gRb70JySRxeScnr1kMks3NtNz5W6qvvipnrxn91CcBaPjBD3VJlzyUy4RxOmndUcaY9A7R\nk4E37fFCYLYxJmqMGQfsD7wgIuvxFkecZoxx8Kb5PppWZ449PgVYouMX+SnZ2EjbwseIHDGV6JFH\n5uQ1nVCIQHU1oX331aUo+qjjn/9ivZmAE4kQ/fjHc/a6weHD2efVlwFIrFyVs9dVfbPd0iDZYIwp\nB47D25wp5SfGmMmAC6xKPSciS40xC4C3gDhwnp0hBXAu3dNqH7cPgDuA39kB8q14YyUqz9TXXk3s\nzTcJDB/GkAcW4ASzPWS2rfCBhtjbmjD6ovUBr+fYjcVwQjn5M9ElOGyYt6TLv/5FcNxYnTGVR3Ly\nmyAiLXiD0OllZ+zk/OuA63opfwmY1Et5O/ClvY9UZYsbj9P64EO4jY2Uf/1rOU8WAOHJk2n+zS2U\nf/XLhDN4A1oxii1fAUDl+ef58vqlM0+g/ruX03LffQz702O+xKC2l8suKdVPJZua6HjmGdz2dpzS\nUspmzfIljoq5Xyd8yCE0/uR/SDbp9i47kmxtJfbGGwx9bCFVl13qSwyl9nck9uprJD780JcY1PY0\nYaisSNbVkdi8ma0XXMT6Aw9iy1fPpOxLX2Kf114hMuVjvsTkhEJUnvMN2h9/gg+POJLON5f6Eke+\ni73yKk44TPij2zXmcyYyaSIjli0lOHIkH045nI4XX6LtL0/SMP9a2v/xjG9x9Xe57ZxU/ULHiy+x\n+eQvgJ3lEjrQUPGNsyk76SSckkzc/7nnokd9nKGLH6fxRz+m9cEHiUya6Gs8+ajjhRe8HQt9Xm48\nUFXF0EcfoeFHP2bLWWfh1jcA0LrgAfZ58d++/y71R9rCUBnXeP31lH3lK5TOOokBN/4Pw5/6K+Wn\nnZY3/8EjkyZRcswxdL7wgt+h5KWOZ58lcuQ0v8MAIDhiBJUXXYRb30D40EPZ56UXcEpKaL7zt36H\n1i9pwlAGLIdWAAAbsElEQVQZ5SaTxN54k9LPzGTQr2+mfHZ+TliLTDuC2Btvkqyr8zuUvJLYWkfn\niy9R8l/H+h1Kl/B++xIyB1A68wSCI0ZQcd43ab1/AeDdBNp86200zL+W2LJlPkda/DRhqIyJvfMO\n6/Y/ALelhfAk//q/+yI8aRLBMWNoWfDArk/uR1rvv5/gPvvk3b/fkAcWUHH2fwMQmTyZ+MqV3h4d\nK96j4Qc/pG3R47T87v98jrL4acJQGdNw1dXQ3gFAcPDgXZztLycQoPyMr9K6YIHfoeSNZHMzTT+9\nkaqrrsy7ex+Cgwd37fMeMgaSSWLLVxBbupRgTQ3lZ82h8z9v+Bxl8dNBb7XH4u+/j1NWRutDD5NY\nt46O555j4G9uJjAov5NFSumJn6Fx/rXEV67s93t+x1evpv7yKyAapfTEE/0OZ6cCpaUEP/IRYq+/\n7u0zPvEgIgd/lMaf/AQ3Hs/5jYb9if5k1R5x29vZOOO/cNvbtykvOe44AmVlPkW1e0KjRxOeOJG2\nxX+hct45u65QpFruX0Djj35McvNmQhMOLIg/uGWf/xz13/kuAJUXf7urC639icWUfvZE3GQSJ+B1\noMQ/+IDQmDG+xVpMtEtK7ZHWBQ/gVFd1fV9+5hkEhgwpmGSRUjLzBNoXL2bTyV+g9Y9/9DscX9Rf\nfAnJzZspOfFEBv78Jr/D6ZPKSy9h4G9+zdDFT1Ax7xwClZVUXXIJdZd9h87XXmPzaafTcO11JBsa\n2HDkUXrvRobk/0cJlTfan3mW4IgRdDz7LA21V1NdexXxdetoue12qn90HdU/2m41l7xX+umZNN34\nMwCckhLf7kL3U2DoUCKHHMKgX/48b6Y+74oTDG63Z0bFN+cRe/NNNp34OZyKCjr//W9C+3nLsifW\nr/cjzKKjLQzVpfPVV0m2tfX6nNvRwZbZp7Px6GNoqL2a0lmzKPvqV6i+upaa95bjOE7eDZT2RXjC\nBKp/8H0A3Lb2nZ9chNyODpKbNlF56cUFkyx2xHEcyr/yFQAqzvkGofH70fRLb+On2Guv6V4oGaAt\nDAVA2xNPsHXu2VR95zIqL7xgm+c2zfoCgSGDcSorGfzbOwgMHbbthjp29kqhqvjvuYTG78fWeefi\num5BJr49EX9vJcn6egCCI0ft4uzCEJl2BFVXfY/yr36FZH0DLXfcAUDLPb8j/t5Khtx/L+1/+xvh\nCRMIDh/uc7SFRxOGouX++2l77E+A7XYa+xGiRx5JcNgwElu20Pnii4A3TpGrPSxyLTxxIm5rK53P\nPU/048X5Hnvaeu55uLFOnLIyAgMH+B1ORjjBIJXfnAdA2amnkNyy2ZvJ94d76fjnP2l77E9snfdN\nAvsMZ8TLL/kcbeHRhNHPue3t1F/5PWjv8AaAn1hM53PPUTHvHNz2dlruuhuA4NixVJx3rs/RZk9w\n6FDKz5pDfW0tw5543Pd1lLLJ7eyk7uJLiL3h3bcQPuigomxVRSZNYtCvbyZZX0940iQarvweW+d9\nE4DkhxtIbNhA/N3lUFJC9LAp1F99DZXz5hGsGbGLK/dfudpAaRXQBCSAuIgcZowZBNwPjMXbQOlU\nEamz518BzLXnXyAii235FLo3UFoEXCgirjEmCtwDTMHbmvU0EVmVi/dWaNy2NlzXJVBWRuytZSQ2\nbui62a7085+j6tJLaH3oYVofepjk5s1UnH8ekSkfo+S444ryj0q6qksvYcPCx2j6+S98W9Y7Fzqe\neZa2R7pnhJV95cs+RpN9gQEDqJhzJoGKCpL19V7349yz2XjcCSS3eDs5D//3c7TccSftix6n8pKL\niRx+GKExY3AiEZ+jzy+5bGFMF5HNad9fDjwlItcbYy6333/XGHMQ3o55E4Ea4K/GmAPsrnu3AGcD\n/8ZLGDPxdt2bC9SJyHhjzGzgBuC0XL2xQrL1/G/R/sRi9nn1ZTYedzwEuuc9hMaMITxhAhVf/zrN\nt/wvAFWXf7foE0VKoKqKQb++mc2nnkbpF75AeL99/Q4pK9r+/GdKTjie6qtrSaxbT+SIqX6HlBNl\nX/xC9zfhcFeyAG/BRfBmU9VfehkAgX2GExw8hNC4cQz6f7fkNNZ85ecsqZOAu+3x3cCstPL7RKRD\nRFYCy4Gpdg/wKhF53u7XfU+POqlrPQgca/f97pcafnw98bVrtylzXZfWRxfS/sRiADZ/yebTZJLg\n2LEABO3NTcGaEQz548MMuv3WfpMsUqJHfZzQhAl0/P3vfoeSFW4sRtvixZR96RRCY8cS/fiRvux+\n6Lehjy1k2NNPMeSRh8BxqL/kMiJTD8cpLQWg5DOfJjhiBLGlS2n7059ovOEnxN5+m81nnoUbj/sc\nvX9y1cJw8VoKCeD/icitwHARSU2O/hBITVkYCTyfVneNLYvZ457lqTqrAUQkboxpwNsSNr1Fw7p1\n66itrQVg+vTpzJgxIzPvLo8ktmyh+eZfk3j/fQb97y203HcfsbfexgkFaf5/twJQVfs9GudfR3DE\nCBLr1zPgxz8iuWkTgUGDuq4TPfxwv96C70o+9UlaH/kjZafPJmD/gBSLjn/9CzpjlBxzjN+h+Cp9\nlt+wJxcTX76C6IzpOOEw6yd+lLJZs4ivWkXs1dcAaPrlr2j5v9+T3LqV9qf/Rulx/wV4/9/ctjZC\no4pjltmu5CphfEJE1hpjhgFPGmPeTn/SjkO42Q6ipqaG+fPnZ/tlfNX58ive11dfI7ZsGfWXXNb1\n3ICf/ZSSY44hMGgQ4QkTCI0Zw5azvk506uEFPwc/k8q/dhZtf15Ey223U3nBt/wOJ6Pa/ryIkmNn\ndH2SVt69OOEJE7q+H7bkrwRHjSKxbj2tDzxI/N13AUhu3Qp4K/qmEkbTT28kvmoVQ+79Q+4D90FO\nuqREZK39uhF4BJgKbLDdTNivG+3pa4HRadVH2bK19rhn+TZ1jDEhoBpv8LtfiC1fTn3t1bjJJJ0v\nv0zk8MNJNjay6eQvEv3UJxn4a+/mpdLPfIbg8OE44TAlRx9NaNw4hv/9aU0WPYRGj6Z01kk033ob\n9d//AW4i4XdIGeHGYrQ//kTeLy7ot9Do0TiOQ2hkDYPuuB2AoF2cMnLYYbQ//gTttsuy89XX6PzP\nf3Dd7s+7rusSe/vt7S9cBLKeMIwx5caYytQxcDzwJrAQmGNPmwM8ao8XArONMVFjzDhgf+AF233V\naIyZZscnzuxRJ3WtU4Aldpyj6HW8+CJb532Tljt/S9sjf6Rt4WOUzvo8Vd+5DLepifIvf5myWbOo\neW85gcpKv8MtGJFDJ5Osq6Plttu7phYXKjcWo+mXv6L5jjvBdYkeW3xdsdmSmmIbmXwIhEKUzT6N\nkuOPY8sZc2h74gliy5bh1jfQcvfdxNes8caIFi5k47HHkWxp8Tn6zMtFl9Rw4BFjTOr1/iAiTxhj\nXgQWGGPmAu8DpwKIyFJjzALgLSAOnGdnSAGcS/e02sftA+AO4HfGmOXAVrxZVkXLdV3q5p1L7N13\niK94j+iRRxKdexR1F1wIJVHKTj6ZQHU1ZSfPwqmuBujaS0D1TeTQQwGIHnUU7UuWUDH36z5HtOc6\n/vEMjTf8BIDq719TcAtE+ik1hhX95CeJfuIoSo47jvLTZ9Pwox9Td9HFEAxCPE7D92pp+F4t5V87\ni8TGTQC4zc1QXu5n+BnnpDelil1tba1bDGMYrY8+St255xOeNInqa64m+vEjcV2XlrvvIXLYFCJ5\ntltaoUo2NND5nzfYMvt0Kr51PtWXf9fvkHab67ps/cY5kExSdfl3Ce+/v98hFRy3vX27btv4ypVs\n+MSnKP3CFxhw7Q/58OOfwLXLrBAKQTzO8Gf+QWjfothnpWuqpC4+WICafnUzlRddyLDFj3ctY+E4\nDhVnzdFkkUGB6mqih00hcuSRNP/qZtqeeMLvkHZb67330fHPf1H57Ys0Weyh3sb4QuPGUX72f1N5\n7jz7e3JY17llJ3uz/ZMtxbfYoSaMAuPGYsSXvU3JzBP8DqVfcEpLGfrgAiq+Oa9rWnIhaX/6acpO\nOUU/SGTBgO9f0zW7qnTW56m87FJqVrzLwJ/fhFNSgluEYxiaMApMsrER8JY7ULlT/tWv0PnSy7Q+\n/IjfofSZG4/bxRSn+R1K0Ss7+WSqLrqw63unvBy3pdXHiLJDE0aBcRtswqiq2sWZKpNCY8dSfe18\n6q+4kmRr/v8hcONxtpxxJkQjRbvCcD5zysuLcv8NTRgFJtnYAI6Do1Nkc678K1/GqSin7bHH/A5l\nl1rvX0Bs2dsMW/RnAnamnModp7wMtwA+WOwuTRgFJtnQgFNV1bXBvcodJxSi/NRTaf39vX6Hsksd\nzz9P6cwTdJMgnwTKK3QMQ/kv2dConxh9VHb6bDpfeYXO11/3O5Sdir32OuHJh/gdRr/llJdpl5Ty\nn9vYqOMXPgqNGUPpF79I/RVXkti8edcVciz+/vs033ob8ffeI3KIJgy/OOXl2iWl/JfqklL+qb76\nKgAa5l/ncyTba77ztzT84IeU//dcQgce6Hc4/VagvFy7pJT/ko2NBAZol5SfgoMHU3X5d2l78EEa\nb/yZ3+Fso/OVV6mq/R4DfvD9freXST7xZklpwlA+cxsatEsqD0Q/8QkiUw+n6Ze/8juULlvPO5/Y\nK68QmTLF71D6Pa9LShOG8lmyqUmn1OYBJxBgwP/8BOJxkk1NfocDQNsfHyUwbBiRgw/2O5R+zykt\n1TEM5T+3s1NXns0ToZHeho+J1Wt2cWb2ubEYOA6D7/6t/n7kAScchljxbeWqCaPQxGI4oVxtlKh2\nxiktJTB0KPE1q/0OhcSHH4LrEhw5ctcnq+wLh3HjMb+jyDhNGAXGjcUgHPY7DGUFR4/OixZGYu1a\nKIlusy+78o8TDuN2Fl/CyPpHVWPMaOAevI2UXOBWEfmFMeb7wNnAJnvqlSKyyNa5ApgLJIALRGSx\nLZ9C9wZKi4AL7X7gUfsaU/C2Zj1NRFZl+735IhbHiUT8jkJZweHDSGzcuOsTsyyxdh3BETU6MypP\nOOEwxLVLak/EgUtE5CBgGnCeMeYg+9xNIjLZPlLJ4iC8HfMmAjOB3xhjgvb8W/CSzP72MdOWzwXq\nRGQ8cBNwQw7ely/cWKd2SeURp6QEt73d7zBIbNhAcB9dBiRvhEO4sU6/o8i4rCcMEVkvIq/Y4yZg\nGbCzjtaTgPtEpENEVgLLganGmBFAlYg8b/frvgeYlVYntfHyg8Cxdt/vouPG4hDRLql84USj0On/\nHwa3vR2nrLi2Ay1kTjiig957yxgzFjgU+Lct+pYx5j/GmDuNMQNt2UggfRRxjS0baY97lm9TR0Ti\nQAMwuOfrr1u3jtraWmpra1myZElm3lSuxWI4IU0Y+cKJRHA7OvwOA7ejQ2dH5ZNQyBtvLDI569sw\nxlQADwEXiUijMeYWYD7euMZ84Ebg69mMoaamhkLf09uNdeqgdz6JRnHtplZ+ctvacUo0YeQLJxIG\nTRh7xhgTxksWvxeRhwFEZEPa87cBf7LfrgVGp1UfZcvW2uOe5el11hhjQkA13uB38YnFvQE1lRec\naFRbGGp7oXBRtjCy3iVlxxLuAJaJyM/SykeknXYy8KY9XgjMNsZEjTHj8Aa3XxCR9UCjMWaaveaZ\nwKNpdebY41OAJXaco+i4sZgmjDzilJTgdvg/hkFHuyaMPOJEivM+jFy0MI4CzgDeMMa8ZsuuBE43\nxkzG65JaBZwDICJLjTELgLfwZlidJyIJW+9cuqfVPm4f4CWk3xljlgNb8WZZFSU3FtNB7zziRCKQ\nJy0MvQcjfzihMOh9GLtPRJ4FepuxtGgnda4Dtls7WkReAib1Ut4OfGkvwiwcOuidV7RLSvUqEsbV\n+zCU39x4XAe980g+JQw0YeQNJxSCzk5ct7h6xjVhFJrOTpyw3riXN6JR3M78SBhOSYnfYaiUsF2N\nIZHY+XkFRhNGgfHWktKlQfKFE82T+zDaO3RabR5x7Dhjsc2U0oRRQFzXhXhcWxh5xIlGcdv9Txh0\ndOBEtYWRL7rGGTVhKN/YXz5HWxh5w4lG82aWlA565xFtYSi/dc260Gm1ecOJluDmw1pSmjDyStcC\noZowlG/sHyZdrTaP6FpSqjd2CwJtYSjfdLUwtEsqbzglUYjHfZ9zr9Nq80t3C6O47sXQhFFIOlNj\nGNrCyBepT/V+d0t502o1YeSL1PI9xbYnhiaMApL65dO1pPJHVzeQj91SruvaWVKaMPJGV8LQFoby\nSdcvnyaMvJHaLtfXcQz72nrjXh4JBsFxQFsYyjfawsg7XV1SfrYw7Bax2sLIH47jQLj41pPShFFA\nuge9NWHkDfup3teEkWphaMLIK064+Fas1YRRSDpjEAziBPSfLV840SgEArjNLb7F4La1dcei8kc4\nVHR7YuhfngLixnXzpHzjBAIER40i/sH7vsUQf/99nIoKnAEDfItBbc8JR3CLrIVRNPMzjTEzgV8A\nQeB2Ebne55Ayzu2MaXdUHgqN/Qjxlat8e/348hWE9h/v9ZurvOGEQqAtjPxjjAkCvwY+DRyEt5vf\nQf5GlQW6PWteCo0b52/CePddQvuN9+311Q5EwtrCyFNTgeUi8h6AMeY+4CS8bV7zQrKtjeSGDd0F\n6RuruOAmk+AmIZl6uLhu0ltP334fe/tt0Jv28k5o3Dia77qL5tvvoPONN3HCIUJjx+JEIsQ/+ACn\ntBS3qQm3o4PIoYeS2LIFkkmckhKSDQ04JSXeOe3tJD7cQGBANW5zM05FBTgOwSFDcFtbcRMJnFI7\nyN7UTLKpGbe5ifa//52Kr33N55+C6skJhUl88AHx91b2vufoNifv+ITA0KEEysszG9weKpa/PiOB\n1WnfrwGO6HnSunXrqK2tBWD69OnMmDEjN9EBDVfV0nrf/XtW2XG8ed2BAJEpH8tsYGqvRY/+FO1/\nfYrWRxcS3m9f3I5OOp55lmRLC6GxH8HdsNFLCuEQzXfdRXDECHAc3PYOAgMH4HZ04ra14YRDBPfZ\nh9hrrxMYMpjk6jXgunQ2NOCUl0Ew2DWFNlBRgVNZSaCigrJTTqH05Fk+/xRUT6F9x9F4/Q00Xn/D\nXl2n5DOfZvBtt2Yoqr1TLAmjT2pqapg/f74vrz3ghuupvrp228L0TxWBgDcDynG849TDcbRvOs+F\nDziAIfff63cYKs8Mvuu3uLFY1yy2HdrFNq5OWVkGo9o7xZIw1gKj074fZcvyhhMK4VRX+x2GUiqH\nnHC4qMYdiyVhvAjsb4wZh5coZgNf9jckpZQqLkUxS0pE4sD5wGJgGbBARJb6G5VSShWXYmlhICKL\ngEV+x6GUUsWqKFoYSimlsk8ThlJKqT7RhKGUUqpPNGEopZTqE00YSiml+kQThlJKqT7RhKGUUqpP\nNGEopZTqE00YSiml+kQThlJKqT7RhKGUUqpPNGEopZTqE00YSiml+kQThlJKqT7RhKGUUqpPsrof\nhjHmK8B3AQdoAr4pIq/b51bZsgQQF5HDbPkg4H5gLLAKOFVE6uxzVwBzbZ0LRGSxLZ8C3AWU4u2J\ncaGI7HyjXKWUUrsl2y2MlcDRIvJRYD5wa4/np4vI5FSysC4HnhKR/YGn7PcYYw7C23p1IjAT+I0x\nJmjr3AKcDexvHzOz9H6UUqrfymrCEJF/pVoHwPPAqD5UOwm42x7fDcxKK79PRDpEZCWwHJhqjBkB\nVInI87ZVcU9anbyzZMkSv0PIG/qz6KY/i276s+iWbz+LXI5hzAUeT/veBf5qjHnZGPONtPLhIrLe\nHn8IDLfHI4HVaeetsWUj7XHP8rz09NNP+x1C3tCfRTf9WXTTn0W3fPtZ5GRPb2PMdLyE8Ym04k+I\nyFpjzDDgSWPM2yLyj/R6IuIaYzI2FrFgwYLXFyxYELXfvg98kKlr74YxCxYs8ON185H+LLrpz6Kb\n/iy65cPP4ioRuRaykDCMMefhjScAfAYYAtwOfFpEtqTOE5G19utGY8wjwFTgH8AGY8wIEVlvu5s2\n2iprgdFpLzXKlq1l266uVPl2RGTyXr49pZTqtzLeJSUiv7YD2ZPxEtLDwBki8k7qHGNMuTGmMnUM\nHA+8aZ9eCMyxx3OAR9PKZxtjosaYcXiD2y/Y7qtGY8w0Y4wDnJlWRymlVIZku0vqamAw3owm6J4+\nOxx4xJaFgD+IyBO2zvXAAmPMXLxuo1MBRGSpMWYB8BYQB84TkYStcy7d02ofZ9uxEqWUUhnguK7e\nrpBrxpjDgeeA2SLyoN/x+GFn9+j0J8aYmcAvgCBwu4hc73NIvjHGjMab5Tgcb1LMrSLyC3+j8o+9\nbeAlYK2IfNbveEDv9M45+0twA/AXv2Px2a7u0Sl69nfh18CngYOA0+39Rv1VHLhERA4CpgHn9fOf\nx4XAMr+DSKcJI/e+BTxE92B+v7SH9+gUm6nAchF5T0Q6gfvw7jfql0RkvYi8Yo+b8P5Y5u0U+Wwy\nxowCTsSbMJQ3NGHkkDFmJHAy3p3pqlvPe3T6ix3dW9TvGWPGAocC//Y5FL/8HPgOkPQ7kHSaMHLr\n58B3RSSvfgn8lHaPznf9jkXlB2NMBV4r/CIRafQ7nlwzxnwW2CgiL/sdS085uXGvP+txX0o1cJ+d\nHTYE+IwxJi4if/Qrvlzq6z06/ciO7i3qt4wxYbxk8XsRedjveHxyFPB5Y8xngBKgyhjzfyLyVZ/j\n0llSfjHG3AX8qR/PkhoDLAHOFJF/+R2PH4wxIeAd4Fi8RPEi8GURWeprYD6x91HdDWwVkYv8jicf\nGGOOAS7Nl1lS2sJQftnRPTr9hojEjTHnA4vxptXe2V+ThXUUcAbwhjHmNVt2pYgs8jEmlUZbGEop\npfpEB72VUkr1iSYMpZRSfaIJQymlVJ9owlBKKdUnmjCUUkr1iU6rVSqDjDEJ4I20ovt6W4HWGPNz\n4OGeu0ymPT8HmCkip6eVDcFbX2kU3v0KtSLybibjV2pnNGEolVltu9rZ0RgzGJi2i5vTHgFuNMaU\niUirLTsFeExEOowxt+CtNXT2Dq+gVIZpwlAq974IpDYMwxgzBfgZUAFsBs6yWxT/HfgccL89dTZw\nnT1+BrjLGBMSkXjOIlf9mo5hKJVZpcaY19Iep/VyzlHAy9C1dtKvgFNEZApwJ91J4V68JIExpgY4\nAG85FewClsuBQ7L5ZpRKpy0MpTJrl11SwAhgkz02wCTgSbtEShBYb5/7M97SKVV4WxU/lLYtMXh7\nqtRgk49S2aYJQ6nca8NbhRS8LWqXisiRPU8SkTZjzBN4e6jMBi7ucUqJvZZSOaFdUkrl3jJgvD0W\nYKgx5kjwuqiMMRPTzr0XL1EMx9sHPt0BwJtZjlWpLpowlMqsnmMY202pxetqOgbAbs16CnCDMeZ1\n4DXg42nnPonX7XS/iHStFGqMGY7X/fVhlt6HUtvR1WqV8oEx5lngsyJSv4f1vw00isgdmY1MqR3T\nFoZS/rgEGLMX9evxbt5TKme0haGUUqpPtIWhlFKqT/4/lyiagrj/ChYAAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f83cb625518>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"dos.plot()"
]
},
{
"cell_type": "code",
"execution_count": 47,
"metadata": {},
"outputs": [],
"source": [
"sigmaxy = kpm.calc_conductivity(chemical_potential=np.linspace(-4, 4, 500),\n",
" broadening=0.1, direction='xy', temperature=0,\n",
" volume=1000*200, num_random=1)\n",
"sigmaxx = kpm.calc_conductivity(chemical_potential=np.linspace(-4, 4, 500),\n",
" broadening=0.1, direction='xx', temperature=0,\n",
" volume=1000*200, num_random=1)"
]
},
{
"cell_type": "code",
"execution_count": 53,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"(-2, 2)"
]
},
"execution_count": 53,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYEAAAECCAYAAAAYfWtSAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXl8XFXduJ87d5YkTdO9TZp0bzmlTUopW4FaNlHABVBU\nFhXEDWQRXheQ+rrgW8FXfioqLiiLIC/iiwiovMriwg6F0jZpyynd96ZL0jRNMnNn5v7+mElMJ5Nk\nktztJOf5fPiQM3Pn5jnn3s4393zPYti2jUaj0WiGJiG/BTQajUbjHzoIaDQazRAm7NUvEkJMAh4A\nJgA2cLeU8s6cYwzgTuA8oAW4Qkq53CtHjUajGWp4+SSQBL4kpZwDLASuEULMyTnmXGBW9r/PAT/3\n0E+j0WiGHJ4FASnlrva/6qWUh4C1QGXOYecDD0gpbSnlq8BIIUSFV44ajUYz1PCsO6gzQoipwLHA\nazlvVQLbOpW3Z1/b1fmg66+/3h4+fDi2bbNo0SIWL15MLBYjlUqRTCYBiEaj2LaNZVkdZYBEIgFA\nJBLBMIyOcjgcxjRN4vF43rJpmoTDYQ58+rOYM2cw7CtfJhwOk0gksG2bUChEJBLpKBuG0fF6Op3G\nMAyi0SjJZJJUKgVALBbrUu5cB6/q1NTURCwW6yj3VKdoNIplWZ7Xad9pp1N0zRcoueiiguqUW/aq\nTvF4nPZ70+nr5FSdDl/1hcw9fNNXfb/3eqvToUOHiEajvt57hdQpHo9TUlLSY51av/JVzNGjid3y\nNd/+PZWWlhrk4HkQEEKUAn8AbpBSNvXnHMOHD2fp0qVdXo9EIl1ei8ViR5TbG6W7cu45upxz+zaK\nzz2H4uLiXo+3LKvX8/W1DM7Xqays7IjXnHZ2pE418witW09RUZEvjoWWY7HYET8fUYcBXienys1S\nEvvQhR33cKGfB+/rNHz4cFfvTXCmTp2ve3e/MzVlCkkpGTZsmKN16E+dOuPp6CAhRIRMAHhISvlY\nnkN2AJM6lauyrx1Be1T0Gtu2SW3ZSnjKlIKOb/+LIOio4GnMFlh1tX5r9ErQ2zJ1oIHUzp2kZ870\nW6Uggt6e7RTiaVZUkNq92wObvuFZEMiO/LkHWCul/EE3hz0JfFIIYQghFgIHpZS7ujnWc9J79mC3\ntRGeOtVvlSFHaO4crLrV2Om03ypKk1yzBqO4mNCUyX6rDDnMinJSu/f4rdEFL7uDTgU+AdQKIVZk\nX7sFmAwgpfwF8BSZ4aHryQwR/VS+E4VC/kxvSG7eDLEYofIJBR0fDvuScukzKnjGampoOXyY5KbN\nRGZM91unW4Lelom6OiJz5hDJ6S4IKkFvz3YK8TQrKkjv3YudSGD00D3jNZ61sJTyRaBLUiLnGBu4\nprdzGUaPp3GN5JYthCdPxigwCJmm6bKRM6jgGRk3DrOyEquuNtBBIOhtaa1eQ2TunMB7tjOYPM3y\ncrBtUvX1hKuqPLAqDCVnDPuVE0ht3kJ4amH5ABhc/Zl+E4/HidRUY9XW+a3SI0FvS2vNaiJz5wbe\ns53B5BkaMwYiEVK7gpUXUDII+EVy82bMApPCGueJVAc/CAQZu62N5DvriVTP9VtlSGKEQpjl5aR2\nBSbNCSgaBHztDupDUngwPcr6jWmaRGtqSNTVEuSVb4PclpaUYNtEhAi0Z2cGm6dZXk46YCOElAwC\nviWGt2wpeHgoDK6klt+Ew2EiNdXYjQdJbd/ut063BLktrdVrCM+ciVFcHGjPzgw2T7NCPwk4gh85\ngXRjI3bjwT51Bw2m/ky/icfjmBMmEBo/PtBdQkFuS2vVKiLV1UCwPTsz2Dx1d5DCJLdsgVCI8KTg\nZPWHIpm8QPAnjQWRxMqVROcf47fGkCYzYSxYcwWUDAJ+5ASSW7ZiTpzYp/G9g60/00/aHaM11STq\ngvskENS2tONxrLVvEzkmEwSC6pnLYPMM6ScBZ/AjJ5DqYz4ABl9/pp+0O0ZqqrFWBTc5HNS2tNau\nhXSayNzM6u1B9cxlsHmaFRWk9uwJ1Mx3JYOAHzmB5JYtmH2cat++cl/QUcGzY2XHmhrS+/aR3hOs\nR+p2gtqWiRUriQhBKLtoXFA9cxlsnubECrAs0vv3u2xUOEoGAT9I9mHhuHaC+tdqLip4tjualZUY\nI0eSCGhyOKhtaa1aRaRTPiConrkMNk9z/HgwjEB1CSkZBPzICfSnO8ivoax9RQXPdkfDMIjW1GAF\nNC8Q1LZMrFxJdN68jnJQPXMZbJ5GJEJo3LhArSaqRgvn4PWNYcfjpHbu7HN3UL51voOICp5HrClf\nE9wRQkFsy3RLC8l17xzxJBBEz3wMRk+zfAKpnfpJYEB4nRNIbtsOtt3nJ4HB1p/pJ50dg7x8RBDb\n0qqthUiEiBAdrwXRMx+D0TNo+wooGQS8JrV1K8bIkYRGjOjT5wZbf6afdHaM1tSQ2rmTVICSa+0E\nsS2tlauIzDn6iOHNQfTMx2D0zEwY00FAKTLLRfR9Ew6/1jjqKyp4dnY0p07BKC0NZF4giG2ZWLmS\n6DFHThILomc+BqOnWVFBWieGB4bXE0iSmzf3uSsIet7XM0io4NnZ0QiFiFTPDWSXUBDbMrFiZcck\nsXaC6JmPweipu4McIO3xRIvUli39WkLasiwXbJxHBc9cx6DmBYLWlunGRlKbNxM9Zt4RrwfNszsG\no2f7+kFB6epSMgh43XjJrX2fIwDeB6v+ooJnrmP7stJBI2htmVhVi1FSQjhnY/mgeXbHYPQMVVRg\nt7RgHzrkolHheDYnWwhxL/B+oF5KWZ3n/dOBJ4BN2Zcek1Le6pVfd9i2TWrrNsKTdU7AT3IdIzXV\npDZvIX3wYJ8T9m4StLa0Vq4kMq8GI6cLNWie3TEYPc2KcgBSu3YRKitzS6lgvFyY437gp8ADPRzz\ngpTy/b2dyMucQLqhAbu1FbOqss+fHYz9mX6R6xieMQOjqAhr9Rpip5zsk1VXgtaWiVWruiSFIXie\n3TEYPUMlJRgjRpDateuIYbt+4Vl3kJTyeeCAE+fy8hExtWMHGEZmk+g+kkwmXTByHhU8cx2NcJjw\nnDkkAjZpLGhtaa1YSSQnHwDB8+yOwepplk8ITHI4aEv0nSKEWAXsAL4spVyd76Dt27ezZMkSbNtm\n0aJFLF68mFgsRiqV6rgY0WgU27Y7EjbtkbpjIbJIBMMwOsrhcBjTNDs2h2gvH96wAWPsWNrSacKW\n1fG+aZqEw2ESiQS2bRMKhYhEIh1lwzCwLAvbtkmn05nlDqJRkslkx2S3WCzWpdy5Dm7VKbcOhw4d\nIpVKFVSnaDSKZVme16ndsXMdjNmCtpUrKcpzXfpynZysU0tLC6ZpunKd+lqn+K5dpHbuJHnUUViW\ndUQd2t39vvd6q1NzczOpVMrXe6+QOrW0tHS8Vsh1Ytx42rZsJWpZntaptLSUXIIUBJYDk6WUzUKI\n84DHgVn5DqysrGTp0qVdXs83dTsWix1Rzn1syy3nniOy/wBWVRXF2dUXu7zfQ7m5uZmSkpKCjy+k\nDA7UKadcUlJyxM0xUEc36pTP0ViwgOZf3k0kEvHcuadyu7vT16mv5aRchzFyBMNnz8YwjC73Zvs9\nXej5wPs6FRcXu3pvgnN16st3RKKqEvYfcOXehd7r1JnAjA6SUjZJKZuzPz8FRIQQY/Md62VOILVj\nB+HKvucDoOuFCCoqeOZzjNRUk1y/nnT2r7AgEKS2tFatIjpvXt6kZZA8e2KwepoVFYFZSTQwQUAI\nUS6EMLI/n0jGLe+6AJ7mBHbuxKyc2K/PDtb+TD/I5xg56igIh7FWr/HBKD9Bast8k8TaCZJnTwxW\nT7OigvRQywkIIR4GTgfGCiG2A98EIgBSyl8AFwFXCyGSQCtwsZQy74QAL+cJJHfsJHrCCf36rB+b\n3/QHFTzzORrRKJHZs7Hq6oidcLwPVl0JSlvato21ciXDLrsk7/tB8eyNweoZpA3nPQsCUsr8d+O/\n3/8pmSGkgSK1c0e/nwQ07hPkZaX9JLVzJ+n9+4nO0xvLBxGzoqJj+LmRk5vxmsB0B/UFr3ICdiJB\nek89ps4J+E53jkFbPiIobWmtWElo/HhCFfmHNgfFszcGq2eofAIAqfp6N3T6hJJBwKvuoNTu3WDb\n/Q4Cg/VR1g+6c4zW1GCtW4fd1uaxUX6C0paZlUPzJ4UhOJ69MVg9Q6NGQSwWiLkCSgYBrxLDqR07\nMIqKMhesHwzWpJYfdOcYOXo22DaWlB4b5ScobWmtXNVtUhiC49kbg9XTMAzMCcGYMKZkEPCK1J49\nhMrLlVm/ZChiFBcTnjUzUF1CfmOn090uF6EJDmb5hEBsLqNkEPAqJ5DaU485fly/Pz9Y+zP9oCfH\nSHUNiYAEgSC0ZXLTZuympiP2FM4lCJ6FMJg9zQkTSO/Z44JN31AyCHiVE0jX12OOn9Dvzw/W/kw/\n6MkxWlONFZBlpYPQltbKlZiTJmGOHt3tMUHwLITB7BkqL9fdQf3Fs5xA/V5CE8b3+/ODtT/TD3py\njNRUY619GzsAG5AEoS0TK1cSmdd10bjOBMGzEAazp6mDQPBJ19djjut/d5DGGyJz50I8TvKd9X6r\nBAJr5SqiPXQFaYKBWVFOarfuDuoXoZA32qm99YTG9/9JYDCuhe4XPS6AVVpKePr0QCwr7Xdb2skk\nVm1tr0lhvz0LZTB7mhMmkNqzx/dtJpUMAl4x0MSw3xe3UFTw7M0xUlONVed/ctjvtky+8w52WxuR\nmi6b9x2B356FMpg9zfJyiMdJNzS6YFQ4SgYBL3ICdjyO3dg4oMTwYNwk2y96c4zU1ARimKjfbZlY\nuZLwjBm9blvot2ehDGZPc0Lmu8XvheSUDAJekNq3D4DQAJ4ENN4Rqa7GWr0aW5HRJG5h9bByqCZY\nGMXFGCNHkNqjg0Cf8SInkN5TD6ZJqIdhdr0xmPszvaY3x2j1XOyWFpKbNnlk1I2Hz23ZvlxEb/jt\nWSiD3TMzQsjf5LCSQcALUnvrCY0dg+HhBjaa/hMaNQpz0qQhvaKoHY9jrX1bPwkoRBCGiSoZBLzI\nCaT3DGyiGPx7X8+go4JnIY6ZZaX9zQv42ZbW2rWQThOpntvrsSpccxj8npl9BXQQCCSpvXsHNDxU\n4z3R6urALB/hB4kVK4kIQcjn9ek1hROEpSOUDAJe5AQGOjwU8m8AHURU8CzEMVJTg1VX5+uwQj/b\n0lq1ikgB+QBQ45rD4PfU3UEBJl2/B3OATwKqrD6qgmchjpGaauymJlJbt3pglB8/2zKTFC4sH6DC\nNYfB7xmE9YO83GP4XuD9QL2UsstMluwm83cC5wEtwBVSyuX5zuVFTiC1dy+xAT4JJBIJJUY3qOBZ\niKM5bhyh8glYtXWEp0zxyOxI/GrLdEsLyXXv9LhyaGdUuOYw+D3NinLS+/ZhJxIYPtXTyyeB+4Fz\nenj/XGBW9r/PAT/3wKlbnEgMa7wnWl0TiOUjvMaqq4NwmIgQfqto+kD7hLHU3r2+OXgWBKSUzwMH\nejjkfOABKaUtpXwVGCmEqMh3oNs5ATudJrVvH6EBLh4XDnv2oDUgVPAs1NHv5SP8akurto7I0bML\n/mtShWsOg98zNHYshMO+jhAKUk6gEtjWqbw9+1oX3O4nTDc2gmVhDmAZafBu85uBooJnoY7tw0T9\nSg771ZZWbS2R6pqCj1fhmsPg9zRCIczx431dOkKNMJvD1q1bWbJkCbZts2jRIhYvXkwsFiOVSnWs\n6x2NRrFtu2NNj/b+uvbxvJFIBMMwOsrhcBjTNInH46Q2bwYgPWoUzc3NQOYih8Nh4vH4EeVEIoFt\n24RCISKRSEfZMAwsyyIajZJOpzEMg2g0SjKZ7NiAIhaLdSl3roOTdcpXbq9DQ0MDJSUlBdUpGo1i\nWZbndWpsbKSkpKTXOjFzFun9+zm0YQOh8nLP69TS0sLIkSNduU493XvxVauIXnYZzc3NBdUpmUx2\nePt57/X276mxsZHi4mJf771C6tTS0sLw4cP7/B0RjUZh/Hhatmwhffiw63UqLS0llyAFgR3ApE7l\nquxrXaisrGTp0qVdXs83TCt327fc5E1uORKJ0NbczOGyMqLDh5P7cJ37O3oqNzc3U1JS0u/P5ytD\n/+rUU7mkpOSIm2Ogjm7UqVBHe9gwDo0ZQ2TjRopnzvSsDp3L7e5OX6fuyunWVhreWU/p8ccRLfA6\nNjc3U5wzn8CPe6+3cnFxsav3JjhXp/b27LPTxArCjY0MGzasf5/vR506E6TuoCeBTwohDCHEQuCg\nlHJXvgPdzglkksIDnyg22PszvaRQR8MwfJ057EdbJteshVCoT0lhFa45DA1Pv+cKeDlE9GHgdGCs\nEGI78E0gAiCl/AXwFJnhoevJDBH9VHfncjsnkNq7d8BJYRj8/Zle0hfHSHU11ip/Rgj50ZaJ2loi\nQmAUFRX8GRWuOQwNT7O8HGvt2w7a9A3PgoCU8pJe3reBawo5l9ubT6f27BlwUhggHo8rMeNRBc++\nOEZramh59FGXjfLjR1tadXW9biKTiwrXHIaGp1le7uvSEUHqDgoMaYeeBDT+EKmpJr17D6n6er9V\nPMFaVdvnIKAJDuaECaR27/ZtRJuSQcD17qD6+o5JHANhKDzKekVfHM3JkzHKyrDqVrto1M3v9rgt\n7XgcS0qiNYWtGdSOCtcchoZnqLwcu6UF+9AhB4368Pt9+a0DxPXEcP1eTAeeBIZCUssr+uJoGEYm\nL+DDzGGv29KSEtJpwnOO7tPnVLjmMDQ8zfLsrGGfksNKBgHXcwL19Y4sI90xdj3gqODZV8doTTUJ\nH2YOe92W1qpawkfN6vPy0SpccxganqHSUozSUt92GFMyCLhJurUV+9ChAS8jrfGXIGww4wWJ2ro+\nzRTWBBM/h4kqGQTczAmks8nEkAOLxw2F/kyv6KtjpKaG1LZtpBsaXDLKj9dtadWuItqPpLAK1xyG\njqdZXu7b0hFKBgE3cwKp+nqIRAiNGjngcw2F/kyv6KtjeNo0jJISEh4nh71sS9uyMnsKz+v7k4AK\n1xyGjqef+wooGQTczAm0J4WdeNoY7PujeklfHQ3TJDJ3rucrinrZlsl174BlEZkzp8+fVeGaw9Dx\nNMsnkPJproCSQcBNUvX1hByYKAb4us1hX1DBsz+OmbyAtyOEvGzLRO0qwtOnE8qzKFhvqHDNYeh4\n6pxAH3EzJ5Das8eR4aHgzV7ITqCCZ38cM8NEvX0S8LItrdq6fnUFgRrXHIaOp1k+QQeBvuDmjZHe\nu9eRpDAM/k2yvaQ/jtGaGpIbN5L2cBKOl21p1dYRqe7fTGEVrjkMHU+zvJx0/V5sl4e/50PJIOBm\nTiBVX+/Y8NCh0p/pBf1xDM+aCbEY1po1Lhjlx6u2tJNJrNWridb070lAhWsOQ8fTnFAO6TRpH7aZ\nVDIIuEm6fq8jy0jD0OnP9IL+OBqRCJGjZ3vaJeRVWyY3bMBuayNSPbdfn1fhmsPQ8QyNHwehkC9d\nQjoI5JCZLezMk4Dbaxw5hQqe/XWMVNeQ8DAIeNWW1qpazKlTCI0Y0a/Pq3DNYeh4GuEwoXFjdRAo\nFLcmkNipFOl9+xx7EuhpN58goYJnfx2jNdVYdd6NEPKqLRO1df3uCgI1rjkMLc/MCCHvh4kqGQTS\n6bQ7592/H9JpxxLD7ft8Bh0VPPvrGKmpJrnuHdKtrQ4b5certrTqaokMIAiocM1haHm2LyntNUoG\nAbf6CdvXnzfHjXXkfG4FK6dRwbO/jhEhIBQi6dHOTV60pZ1OY9WtHtAeAipccxhann4tHaFkEHCL\ndP1eQqNGYTj0CDpU+jO9oL+ORlERkaOOIuHRpDEv2jK5cRP24cMDWjhOhWsOQ8vTr+4gTxfmEEKc\nA9wJmMCvpZS357x/OvAEsCn70mNSyltzz+NWTiBVv8expDAMrf5MtxmIY6Sm2rPlI7xoS6uuFrOq\nCnP0qH6fQ4VrDkPLM+TT0hFebjRvAncBZwPbgWVCiCellLmDuF+QUr6/p3O5lhOo34vpUD4AIJlM\nKjHZRQXPgThGampoeeT3Dhvlx4u2dGI7SRWuOQwtT7+WjvDySeBEYL2UciOAEOJ3wPlAn2fyuJkT\ncHJvYbc3v3EKFTwH4hiprsb69q3Y8ThGLFbQZ9Jpm0QyTTyZos1KEU/2/IdHOGQQMgza2loYngoR\nDhmY2dfM0L//c2Rhwto6ihadOqBzqHDNYWh5muXl2E1NpFtaCJWUOGBVGF4GgUpgW6fyduCkPMed\nIoRYBewAviyl9Gwt4FT9XsKTJ3n16zQeEa2eC+k01po1RI89tsv7m+qbWbZxP6u2NbLjQAt7D8XZ\n3xzHjb81DINMQMgGh1CnYDEsFmbmhFLmVo3k9KPHM2nMsC6ft20bq66O4Vdf5bycxlfM8nIA0rv3\nEJo+zbPfG7TFupcDk6WUzUKI84DHgVm5B+3cuZMlS5Zg2zaLFi1i8eLFxGIxUqkUyWQSyPTR2bbd\nMXSrvc+ufXp3JBLBMIyOcjgcJrVnN8yrprm5mXA4jGmaHdvGmaZJOBzuUk4kEti2TSgUIhKJdJQN\nwyAUCtHS0kI6ncYwDKLRKMlksuOvhlgs1qXcuQ5O1KlzHbqrUzKZpLm5uaA6RaNRLMvyvE6WZeW9\nLoVeJ1MIDr36GsVCEIlEaIvHeUHu43ev7eCdPc3MGD+MuZXDefecsVSMGsaokjDREMQiIYaXFBMi\nTTqdv07RaAwrmSSesEhYFtGiYlKpNPGERSptY4YjpNI28USCVNrGCJmkbTrexwiBYdBwqJUN9Yd5\ntnYXdz2zjuOmjuTyUydTM2lkR51SW7diNzVhzBYcPny439fJNE1aW1t9v/d6+/eUSqVobm729d4r\npE6WZdHa2trn74jOdcI0oaiIli1bOnKTTtepNM+Ks14GgR1A5z+zq7KvdSClbOr081NCiJ8JIcZK\nKfd1Pq6yspKlS5d2+QX5+uRiOY//uQmczuXGvXspraqipFND5Z6zL+XW1laKc/Z+Hcj52ulLnQo5\nZywWO8JzoI5u1CmVSg3IMbZgAfaatZSUlLCzoZXv/HENclcTH104hf++dAHlI/u2R2935Lvm/WHL\nvsM8+OImbny4lg8cW8kXzj6KstJSWjdsJFReTmziRHI7tty+N8H5e6+3cjQadfXeBGfq1Lk9B+LU\nXF5OuKFhQN9B0HudOuPlENFlwCwhxDQhRBS4GHiy8wFCiHIhhJH9+cSs3/7cE7mRE7BtO7uhjDOz\nhWFo9We6zUAdowvmY61Ywd9X7+bjP3+JkliY/73+XVx11izHAgA415ZTxg7j6xdUc/enT2TNjoN8\n8ucvs3p7I4na2n5tJ5mLCtcchp6nWeF9ctizICClTALXAn8D1gK/l1KuFkJcJYRo7+C8CKgTQqwE\nfgxcLKX0ZAUpu6kJu7WVULZfTjO4iM6fzxPF0/jPR1dx9VmzuOPSYxkzvLAksZ/MrRrJPZ9dyGlH\nj+eqe1/nia1xIvPm+a2lcQk/5gp4mhOQUj4FPJXz2i86/fxT4Ke9nceNeQLt0descC4I5D6SBRUV\nPAfiaNs2d2+2eeSED3Pr3DBnnTTFQbMjcaMtI+EQN557NNVVI/nOIwn2D4vwxbRNKNT/kUYqXHMY\nep7mhAkkt2135FyFouSMYTe6g1K7dmGMGOHo0Kyh9ijrJv11TKbSfOePdTzx1k7+a9vTnLDd3ZnD\nbrblGcPjfPvP3+dvDWGW/O9K2qz+/y4VrjkMPU8/lo5QMgi4MVkstXs3ZrlzE8WAjox90FHBsz+O\nLfEkX/6f5by5+QC//PSJHDN7ItZbK1yw+zdutmVixQqOjsb59edOZsOeQ1x7/zIOtvRvMxMVrjkM\nPU8/JowpGQTcILVrN2ZFhd8aGofY3xznC/cvY9+hOL/6zElMG1dKdP58Em+9pcxGJbkk3lpBZP58\nqkaX8KvPnIRhGFx93zL2NrX5raZxiFB5Oan6emwPF85TMgi4khPYtbtjsoZTDLX+TDfpi+PKrQ1c\n8ctXGBYL84srT2R8WREA0eOPI71/P6lNm12ydLctrRUriR47H4ARJVF+/MnjGDs8xlX3vs6OAy19\nOpcK1xyGnqdZPgEsi/SBA46crxCUDAKu5AR2O/8kMNT6M92kEEfbtnnopc1cc/8y3lNTwZ2fOI7S\non+PoTbHj8ecNo34a6/56tkfbMsiUVfbEQQAiqNh7rh0AUdVlPH5e19nY31zwedT4ZrD0PM0J2S6\npL3sElIyCLiRE0jv2uX4k8BQ6890k94cdxxo4brfvMH9z2/gux+dz3XvEYTNrrd3bOFJJFwMAm61\npbV2LSSsLhvJRMMhvnPRPE6eOZar7n2d1dsbCzqfCtcchp6nEY0SGjOGtIfDRJUMAm7gxpOAxn1S\naZtHXtnCZT97maKoyf9ccyqLZ3c/4S920knEX3vdQ0NniL/yKpHquYSGdV1PKGyGuOX8ubxv/kSu\n/c0bLNvYZX6lRiG83mGsoHkCQohhQJuUMhDPZqGQs7HLbmsjfeCA408CQ2ktdLfJ57hyawN3/vVt\ndjS0csv5czm7urzXVTqjC08itXUryR07CVdO9MTTCeIvvkRs0aJu3zcMg+vfKygrjvClh5bz7Q/P\n44w53Y92U+Gaw9D0DHk8QihvEBBChMgs63AZcAIQB2JCiH3AX4BfSinXe2bpMu0bOYQcnCgG7i15\n7TQqeLY7JlNplm9u4PevbeGVd/bxwQVV3HHpAkaXFpaYM6uqMCsrSbz8MuGPXOSap6PntCwSr71G\n6Weu7PE4wzD41GkzGFES5ev/u5Lr3nMUH1s4JW9gVOGaw9D09HrpiO6eBP4BPAt8DaiTUqYBhBCj\ngTOA7wkh/iil/K03mkfidE4gtWsXxGKERvV/p6Z8WJalxOiGIHomU2k27W1m/Z5mdjW0smVvE/XN\nFhv2HKLNSnN2dTkPXn0K08d3XRWxJwzDIHbGGbQ99xwlLgQBN9oysWIldiJB9IQTCjr+QydMYuzw\nGN9+rJblmxu48dzZVOSsjxTEa56PoehplpeTWL7ckXMVQndB4N1SSiv3RSnlAeAPwB+EEMHf7qdA\nUjt3YVZNg8zcAAAgAElEQVRUKLOf6WDmcFuS+1/YyJ+Wb+dgq0XlqGIqR5UwrjTMyTPHcsnJUzl+\n2mhKYv1f8aTo3WfRcN312JaFocCuVfEXXyR63II+zWZfPHs8939+Id//y1o+9pMXeffccs6YO4G5\nlSMYNUyNLpahSnjSJFr/+Lh3vy/fi+0BQAgRAz4MTO18rJTy1nxBwiuczgkkN20iPG2qo+eEodmf\nORDW7DjI1x5ZQWkszA3nzuZdR41nWFHmtkskEo55xhadim0lSCx7g9gpJztyznbcaMu2Z56h+Nxz\n+/y5SWOGcecnjuPNTQf44xvbufWxWg61JQmbBqWxMKVFYYbFIlRXjeADC6qYPbHMcfeBEpR7szec\n9DSnTSO5dSt2MokRdn95t95+wxPAQeBNMnmBQUkmCHi3k4+mK2t3HOT6B97gnHkVfPG9s4mE3Ru4\nFiouJnbKqbQ984zjQcBpktu2Ya1cxaif/KRfnzcMg+Onj+H46WNIpW32HGxl36E4Dc1tJFJwsNXi\n1fX7+PSvXuW8+RO5/j2C4cXBfzoazISnTYVkktSOHYSnuLfYYcfv6+X9KinlOa5b9BGncwLJTZso\nWbDA0XOCs3+9uonfnm/vbOL6B97gvGMmcuO5s/N2yzntWHzeuRz6fz+g7OtLMBycge60Z8ujfyBS\nU0NkxvQBn8sMGUwcVcLEUSU0Nzd37DJ10YmTWbvjILf9aTWX/exlvnFhNcdPHzPg3+cEft+bheKk\nZ2j0aIyyMpKbN3sSBHr7c+tlIURNL8cojW3bJDdtJjx1qt8qQ5J3dmcCwHvnVXQbANyg+P3vI93Q\nQPyFFzz5ff3BtixaHv4dwy69xPXfdXTlCO75zELOmz+RLz74Jj/669vEB7BKqab/GIZBeNpUkus3\nePL7uhsiWgvY2fc/JYTYSKY7yABsKaWvu1o4mRNI79qNffAgYXGUY+dsJ9+2b0HELc9U2qa5zaI4\nGiaap3vnzU0HuOl3b3F2dQVfOu/oHgOA046h4cMp/tCFNP/q1xSdfrpj53XSs+UPf8BubaX4og87\nds528nlGwiGuOmsWp8way7cfq+X1Dfv55odqEBX+5QoGy78h27bZ0dDK2NIYRdHenzwjc+ZgrV7t\nlF6PdNcd9HlgF+DdUnY+kVi1ktDYsZgTnZ84pMpoI6c9l23cz8Mvb+b1jftJpjLjp8tHFlFdNZLq\nqhGMKY2xYksDTyzfzuXvms5nTp/Rq4MbbVn6+c9Tf9a7ib/6KrGFCx05p1Oe6cZGmr73fYZfd62j\ne1y005PnvMmjePDqU7jzb5Ir787kCj53xkzGZRfi60wyleYvK3bydO0u9h2KM2N8KWfMncBZc8oH\ntOlNIZ5BojfP+/61kbv/sZ7p40v51WdOYlgvo9si1dW0/M/DTip2S3cmnwBOAtYBfwX+KqX0dpHr\nHnAyJ2CtXEVk3jxXbrah2J95//Mb+fU/1/PBBVXc+YnjGV8WoyWRYmN9M6u3H+Rvq3ZxsNVi2rhS\nfnDZcZw4o7C+ZzfaMjJzBsMuv5yGL32F8X/5E6GRIwd8Tic87bY2DnzhGsyJFQy78lMDdspHb54l\nsTBf++BcPrCgkp/8TfKRH7/IpadM4bJTpzEsFsa2bV5Zv4+fPr2OhsMJLjiuisrRJby98yC3P7mG\ne/+5gStPm8GZc8sxBxAMgv5vaOu+w7QkklSUhrr1XLeriXv+tYH/vuRYfvq05JFXtnDl6TN6PG+k\nuhpr3TrstjaMoq7B10mMnma6CSFmA+cC7wVGkJlE9lfgJT+XkLj55pvt22+/3ZFz7b3gQ8ROW0zZ\njTc4cr7OdE6+BRmnPJ9/u55bfr+CH378OE5wOLHoVlvabW3sveij2C2HGbFkCbFTTsYo7rrxvG3b\nkExCMomdSkEqBeFwZp5BJNLxR8ShQ4cYZhikDzZhNx/Cbm3FbmsDM0yobDih4WUYo0YSyvc7Egni\nr7xC023fI93czLj/fcS19az60p62bfP82/Xc9cw6DhxOUDNpJPVNbWzb38JHTpzMp06bfsRqrU2t\nFr97ZTO/f20rpUVhzjh6AovEOOZNGtXnUV+Feu4/FOfZut20JJIsEuOYVe5uF1Yyleb2P63hz2/t\nIGRA1ehiln70WGaVDz/iONu2ufq+ZZSPKOJbH57Hk29u5xd/f4fHbzwtbxdpx+fa2thVcwyjf303\nRaed5qR6l4jcYxDojBCimMxs4XOBk6WUx/f1twshzgHuBEzg11LK23PeN7Lvnwe0AFdIKbtMnbvl\nllvs7373u3399V1I7drF7uNPZPxzzxCZPXvA58ulra2NojxRfNnG/fzulS1MHlPCZ86Y2eujodt0\n9nxz0wH+tmonI0qifGzhFMYWuBn77sZWPvmLl7n8XdO57FTnh9t215ZOkD58mKbv/Bctf3gMu6UF\nimIY0RikUtjZL356Wyo4EsGIRLDj8S7HGkVF/z5PO0UxQiNHERo1CqO4GLupidTOndipFMM++hHK\nbr7JkSeT7uhPeyZTaV5+Zx9yZxNlJRHOmDOhY6+GfBxuS/JM3S5ekHt5Y9N+zJDBpSdP5eOLplEU\nKWxEVm+eiWSaR17dwn3/2kDFyGLKiiOs2tbIxSdP4QvvPmpATyE98cirW7j/+Y385PLjKR9RxA/+\nsoZ/yX3cfvH8I/4Aerp2F7c/uZrfX/8uxg6PEbdSXPDD57n2PUfxvvmVR5zz7Z1N/Pmt7RxVXsZ5\n8ydy8KqrMMeNY+R3lzqp3v8gMFCEECaZ7qWzge3AMuASKeWaTsecB1xHJgicBNwppTwp91xLliyx\nly4deMM0fv0/Sbz2OuOe/mu/uoPiVormtiRjuvmitCyrS8Lo9Q37+Y+H3uSDC6pYvvkAo0qi3PnJ\n43v8q8Bt2j0femkTP3/uHc6aW872Ay3sbGjlF1eeyJSxXVeu7Ewylebq+5YxvCizvr0TfcHdObqJ\nnUhkHsGbm7ETCQwzDGETIxw54v+EIxhmCDuZAiuBbVnYCQsSCVKmSWTMGEIjRxAqLYVYDMMwMk8S\nbW2km5pINzaSbmjI/teI3dpKaMQIzPJyosctyPsk4jRetGdn2qwUL8q9/OzZdRRHTL5/6bFMHNV7\nriOfp5VMs2X/YV5bv58/vrGNeDLFde8RHQsIrtzawM2/W8G8ySP5zkXHdPtv61CrRZuVypvr6InW\nRJIP/egFPnfGTC48YRKQ6bb67SvbuPdfG/jK++bwwQWVbN3fwlX3vs5lp2QCXzu/+vt6npf1PHDV\nyR3fO7XbGvnig28wb9Io5K4mKkYW8eWxBxn51esZ/8+/E66q6pNjD/QeBIQQZwMfBe6SUq4QQnxO\nSnn3QH+zEOJk4FtSyvdmy18DkFLe1umYXwL/lFI+nC1L4HQp5a7O53rvez5iH3XMuYRsGwM7WzMb\nI1OhzP+xwYYQNnR6rf194nHspibGvOcMzj5jHvMmH7lu0MGWBE+t2MmY4THOnDOhy9r0z9Tt4vYn\n13A4nuT46aP55oU1XW6m3EfZlniSS+56ifOOmcjnz5pFw+EEn/7Vq5w0Yww3fWBur22YTKX50/Id\nJJJpzj+uqqBRBoXQ3NzMq5ub+fZjq/jvS47l5FnjSKdtvv3HWlZubeDXn1nY4xPBXc+s46+rdvLg\nVacw0qUlCYZa15rb+OXZEk/y7cdqWbG1gds+Np8FU0f3ePyBxib2HLaRu5p4e2cTb+9qYv2eQyRT\nNtPHl3LuMRO56MRJFEePfJre2ZDZX6JqdAnfu/jYI/6ttCaS3PXMOv74xnZSaZvZE8v4/JkzWThz\nbEF/DP7m+Y08uXw7j1y3qON7ob09/7pqJ7c/uYay4ggNh+OcObecb1xYc8QTyf7mOBf84F98/9IF\nLJw5llVbG7jhwTf58ImT+cK7Z9HUavGD/3ubZ2t3M7VtP8MP1BMeVpz5o2QAfOKMWZz4vnd1qWC+\ns14JXA18Pbtg3Pw8x/SHSmBbp/J2Mn/t93ZMJZmRSh1ESRBp3AG2TdnIkZSNGIERMknbNmkbbAMM\nI4QNpG07Ex6MzMVKkznGiEQJVVWyMx3m6vuWcezkkVx52nTGDQvx5Fu7+cObOxk7PEZjS4LHXt/C\n0g/PpWxYEfF4nFc3HODbj7/NtWfP4piqUn7y7AY+/atX+X+XzmdiWQTbtjEMg2QySUtLC+l0GsMw\n+MU/txALG3zkuAmZmyYW41vnH821v13BjLFFnH/8ZFKpVMcGFbFYrKOcStvc+uQ61u5sIho2+Pvq\nnfzgsmOJmCESicxm45Fs33R7ORwOY5om8Xg8b9k0TcLhMBt3NfDdJyRXnzmD46eO5PDhw9i2zX+8\nZzpfe3QNNz38Jnd8rJpo2CQajWJZVkedlm89xP+8vJk7PjaXsJ3AsjL1bt9pqXMd8pWj0Si2bWNZ\nVkcZ6FKnlpaWPtUpt5xIJLBtm1AoRCQS6SgbhtGlTtFotEsdCq1TS0tLwXXq63Vysk7JZJLW1lbH\nr1NvdUpbcZa8fyYPvryd63/zBlefMY0PHFtOJBwmHA7TcKiFt7Y0snzLQVbvPMSmvc3YNkwZW4Io\nL+Pdc8Zy9elTmD5uGCOHl2TqlGijOXHkdSqLwE8/uYAbfvsW19z3Grd8QDB5XBmrtjbw3SfXYoYM\nfnjpMYwbHuOxZdv46sNvUV01gssXTWVuRQmGYeStU4uV5sEXN3LtWdNJWgnsdKaOLS0tmKbJWUeP\nY+6EBazYdpAJI4qZP2U0ba0tR1ynIiPJh46byHcfr+PC4yt54KUtfHB+OZefPJFkMknUSPGV907n\ngvnjWbf7MPvr3iZ1oAE7mSQUCmHbdseqpfnK2DbpzmUglUpTMjJ/0M/3JHC3lPJz2Z9vB86SUha2\nfGEPCCEuAs6RUn4mW/4EcJKU8tpOx/wZuF1K+WK2/Bxwk5Tyjc7ncion0M6OAy3c+68N/N/KnaRt\nmDpuGJ89YyZnHD2BhpYENz74Zuam+fhxbN53mC8++AZXnTmLS06ZCmT+Qv/uE6t5cV09d1y6oOOp\nonN/Zt32Rj5/z+vcdcUJzJ9y5FPHUyt2cPuf1nD9ewRnzJ3AmDzLIv/82XU8uXwH931uIZFwiCt+\n+QoXHjep21EG+w7FueefG1i9o5EZ44dz5WnTmTSma7eOlUzz2V+/yoQRxdx+8fwufwkdbElwxS9f\nZeHMrk8rm/Y28/l7XuOSk6fyqdN6Hu0wUNzMCTiJ9iycZ+t2c9uTqyktClM1uoSDLQk27T3MsJjJ\nCdPHcuzUUUwfU8ScyWMKziHk0ng4wbceq+WNTfupHFXC1v2H+ciJk/nC2Ucdcc7dja3c868N/HXl\nTqaMHcbNH5hL9aSuOZlfPPcOz79dz4NXn3LEX/d9bc9EMs2P//Y2K7c2cv5xVXz4hEleDYct6Eng\nL+0/SClvFkJc59Av3wFM6lSuyr7W12Mcb6zK0SX854U1XHP2URxqSzJpdElHv/aY0hh3XXECX/6f\n5Vz4o+dJJNNc/q5pHQEAMjs7/eeF1fz82Xe47jdv8JX3z+GceRWY2eUIEsk0tz2xmvOPq+oSAADO\nm19J2s580d/x1FpGDYtyVPlwPrpwCqfMGsvTtbt46OXN/PTyEyjPLgn81ffP4ZZHVnDGnAlMy1lO\n+e2dTdz42zeZPKaEDx5bxcvv7OVTd7/Kdz86v8uQzB/839s0HE5w5yePz9uuI0qifO/i+Xz2ntcY\nX1bEFYunYxgGW/cd5sbfvsmio8ZzxeKBL2nQG6aDSzu4ifYsnHdXl3PijDG8tG4v+w7FKYmazKkc\nwVEVZR1fsJmcQP9dRw6L8sOPL6Bu+0G27DvM3MoRXf69AJSPLGbJ+dV84d1Hcc8/N3D1fa/z3Y/N\n513i37vUHWiO88irW/jmh2q6JJz72p7RcIgvv29O/yrlML0NETWAjwPTpJS3CiEmA+VSyj7vzyeE\nCJNJDJ9F5ot9GXCplHJ1p2PeB1zLvxPDP5ZSnph7LieHiBZKMpXmjU0HGFkS7XG1xceWbePnz66j\ntCjCB+dPYO7ksTz00iZ2NLTym8+f3LEqZj7SaZtdja1sqG/mzU37efzN7ZREwzS1Wtz0/jl88Lgj\nk0PfeHQlOxta+eWnT+q4KVdsaeBLDy3n3GMquPHcozFDmaTkfc9v5N5/beCr7/v3eR55dQs/e3Yd\nP7qkhmNn9Lyhzmsb9nHLIyuZPKaE8pHFvLRuL6cfPZ5vXFiTdy9fp9F97c6iPXvmkVe38NOnZUe/\nPcAdf1nL6u2N3Pu5hV3+YFKlPSnwSaAzPyMza/hM4FbgEJn9BPrcPSSlTAohrgX+RmaI6L1SytVC\niKuy7/8CeIpMAFhPZoioOzNl+kHYDHXcDD3xoRMm8d55FTz55nb+tHwbv3lpG/OnjOKuK47vMQAA\nhEIGlaNLqBxdwuLZ47li8QxWbm1g2rjSvCN0bjz3aD7+s5f47z+v4YZzBP9Ys4fb/7SGj586lc+e\nMbPjRjUMgytPm0HV6BL+6/E6nlu9GzNk8MamA3znonnMmtDz6B+Ak2aM5eFrT+XPb+3gUKvF9y85\nlhNnjFFmRqdG0xc+tnAKVjLNTQ+/xR2XLcBKpfnjG9v42adOGHT3fG9PAsullAuEEG9JKY/NvrZS\nSnmMZ4Z5+NrXvmbfdtttvR/oM62trRS7PNxP7mriqw+/xZ6DbZRETb54zmzOP6774WSb9zbzxJvb\nSds25x83ienjSz3xHCgqOIL2dBq/Pe/71wZ++ffMTrrXnn3UEUM9O+O3Zx/o85OAlR3fbwMIIcYR\ngPWEnN5Uxi3CHmwIISrK+P11i9i87zBVo0p6fdqYOq6UL55z5MQ4LzwHigqOoD2dxm/PT502g/fU\nZHYdnDiq+y95vz0HQm/fpj8G/giMF0IsBV4EnBuW009Svc3eDAjtw+PcJhYxERVlvQaA7vDKcyCo\n4Aja02mC4Fk5uqTHAADB8OwvPX5rSCkfEkK8SSaZawAXSCnXemKm0Wg0Gtfpbj8BQ0ppA0gp3wbe\n7ukYr1ElMROEYXiFoIKnCo6gPZ1Ge7pPd08C/xBC/AF4Qkq5tf1FIUQUWARcTmZF0ftdN8yDzgk4\niwqeKjiC9nQa7ek+3ZmfQ2b5iIeFENOARqCYTA7haeBHUsq3vFHsiio5gUQiocTOSCp4quAI2tNp\ntKf75A0CUso2MnMEfiaEiABjgVYpZaOXcqrj1QqtA0UFTxUcQXs6jfZ0nx77VYQQLcBsMquKXiaE\n6H2pSw9QJSegSreVCp4qOIL2dBrt6T69mRcBDwCjgXLgESHENa5b9YIqDa7K46EKnio4gvZ0Gu3p\nPr1lM/aR2UWsDTpWFX0FuMttsZ7QOQFnUcFTBUfQnk6jPd2ntz+p93PkDOEUeaYda/KjSj+hCp4q\nOIL2dBrt6T69PQk8AbwmhHicTMC4gEzCWFMAquQuVPBUwRG0p9NoT/fp8UlASnkzcA1wENgLXCal\n/LkXYj2hysSM9t2Xgo4Knio4gvZ0Gu3pPr3OcJBSvgy87IFLwaTTvq9hVxBeb+bdX1TwVMERtKfT\naE/3UWOYTQ6q9L+pEqxU8FTBEbSn02hP91EyCKiCKv2EKniq4Aja02m0p/soGQR0TsBZVPBUwRG0\np9NoT/dRMgio8uiVTCb9VigIFTxVcATt6TTa0308Wfouu2H9nWT2D24BrpBSLs9z3P3AaWRGI5E9\nbkXucarkBFSZ1KaCpwqOoD2dRnu6j1frn54LzMr+dxLw8+z/8/EVKeWjHnlpNBrNkMar7qDzgQek\nlLaU8lVgpBCior8nUyUnEIvF/FYoCBU8VXAE7ek02tN9vHoSqAS2dSpvz762K8+xtwkhvgE8B9ws\npeyyeefWrVtZsmQJtm2zaNEiFi9eTCwWI5VKdfTNRaNRbNvGsqyOMmTW+IDMgk+GYXSUw+Ewpml2\n7BWaWzZNk3A43KWcSCSwbZtQKEQkEukoG4ZBOp3GsizS6TSGYRCNRkkmkx2PjrFYrEu5cx28qtPB\ngwcpKioqqE7RaNSXOjU1NVFUVOTKdXKyTm1tbZSVlfl+7/VWJ9u2A3Hv9VanpqYmYrGYr/deIXVq\na2tj2LBhvt57hdSptLSUXIK2Hc7XgN1AFLgbuAm4NfegyspKli5d2uXD+SZr5Ebo3Cx+bjn3HAMp\nNzc3U1JS4uj5wfk6FRUVHXFzONkGhZSh9zr57diXcru7n/deb+Xm5maKi4sLPr4dr+sUi8Vcve7t\nv6Mz/a1Te3v6fS9C73XqjGtBILvk9GezxWXApE5vVwE7cj8jpWx/MogLIe4DvuyWn0aj0WhczAlI\nKe+SUs6XUs4HHgc+KYQwhBALgYOdvvA7aM8TZEcTXQDU5Tu3zgk4iwqeKjiC9nQa7ek+XiWGnwI2\nAuuBXwFfaH9DCPGUEGJitviQEKIWqCWzpeV/5TuZHiLqLCp4quAI2tNptKf7eJITkFLaZFYjzffe\neZ1+PrOQ8+nJYs6igqcKjqA9nUZ7uo+SM4Y1Go1G4wxKBgGdE3AWFTxVcATt6TTa032UDAI6J+As\nKniq4Aja02m0p/soGQR0TsBZVPBUwRG0p9NoT/dRMghoNBqNxhmUDAKhkBraqqwxroKnCo6gPZ1G\ne7qPGt+miqJK7kIFTxUcQXs6jfZ0HyWDgCo5gfZFnIKOCp4qOIL2dBrt6T5KBgGNRqPROIOSQUDn\nBJxFBU8VHEF7Oo32dB81vk01Go1G4wpKBgFVcgLtGzkEHRU8VXAE7ek02tN9lAwCGo1Go3EGJYOA\nKjmBfDv+BBEVPFVwBO3pNNrTfdT4NlUUwzD8VigIFTxVcATt6TTa032UDAI6J+AsKniq4Aja02m0\np/soGQQ0Go1G4wxKBgFVcgLhsCcbtw0YFTxVcATt6TTa0308MRdCzAbuAxYAS6SUd3Rz3DTgd8AY\n4E3gE1LKLs9ZqvS/qbL5jQqeKjiC9nQa7ek+Xv1JfQC4Hsj75d+J7wE/lFLOBBqAT+c7SJUNHOLx\nuN8KBaGCpwqOoD2dRnu6jydBQEpZL6VcBnS7ypIQwgDOBB7NvvQb4AIP9DQajWbIEqSOrDFAo5Sy\nfYue7UBlvgN37NjBkiVLsG2bRYsWsXjxYmKxGKlUqmOHn2g0im3bHav7ta/t0Z7Fj0QiGIbRUQ6H\nw5im2RHRc8umaRIOh7uUE4kEtm0TCoWIRCId5fYuq5aWFtLpNIZhEI1GSSaTHU8ysVisS7lzHYJY\np2g0imVZntepuzoErU7xeFyJ6wTQ2tqqxL3X3Nzs671XSJ3i8bjv914hdSotLSWXIAWBgpk8eTJL\nly7t8nq+CRu5G0DnLvSUW849x0DKlmU5er52nK5TWVnZEa857exEnfx2LLQci8WO+LkzXt57vZX7\nc2+C93UaPny4q9cdnKlT5+vu970IvdepM651BwkhrhFCrMj+N7GAj+wHRgoh2gNTFbAj34E6J+As\nKniq4Aja02m0p/u49iQgpbwLuKsPx9tCiH8AF5EZIXQ58IRLehqNRqPBo8SwEKJcCLEd+A/g60KI\n7UKIsux7T3V6UrgJ+A8hxHoyOYJ78p1PDxF1FhU8VXAE7ek02tN9PMkJSCl3k+neyffeeZ1+3gic\n2Nv59GQxZ1HBUwVH0J5Ooz3dR41v0xx0TsBZVPBUwRG0p9NoT/dRMghoNBqNxhmUDAI6J+AsKniq\n4Aja02m0p/soGQR0TsBZVPBUwRG0p9NoT/dR49s0B1VyAqqsMa6CpwqOoD2dRnu6j5JBQBVs2/Zb\noSBU8FTBEbSn02hP91EyCKiSE1Cl20oFTxUcQXs6jfZ0HyXNVWlwVTafVsFTBUfQnk6jPd1HjW/T\nHHROwFlU8FTBEbSn02hP91EyCKiCKv2EKniq4Aja02m0p/voIOAiquQuVPBUwRG0p9NoT/dRMgio\nMjGjpzW8g4QKnio4gvZ0Gu3pPkoGgXQ67bdCQbTv7hN0VPBUwRG0p9NoT/dRMgio0v+mSrBSwVMF\nR9CeTqM93UfJIKAKqvQTquCpgiNoT6fRnu6jZBDQOQFnUcFTBUfQnk6jPd1HySCgyqNXMpn0W6Eg\nVPBUwRG0p9NoT/dRMgiokhNQZVKbCp4qOIL2dBrt6T6erH8qhJgN3AcsAJZIKe/o5rj7gdOAg9mX\nrpBSrvDCUaPRaIYiXi2CfQC4HriggGO/IqV8tKcDVMkJxGIxvxUKQgVPFRxBezqN9nQfT7qDpJT1\nUsplgCODaXVOwFlU8FTBEbSn02hP9wnidji3CSG+ATwH3Cyl7LKD8/bt21myZAm2bbNo0SIWL15M\nLBYjlUp1XIxoNIpt2x2TONqz9+0LPUUiEQzD6CiHw2FM0+zYMDq3bJom4XC4SzmRSGDbNqFQiEgk\n0lE2DAPLsrBtm3Q6jWEYRKNRkslkR/9hLBbrUu5cB6/qdOjQIVKpVEF1ikajWJbleZ3aHd24Tk7W\nqaWlBdM0fb/3eqtTu7vf915vdWpubiaVSvl67xVSp5aWlo7X/Lr3CqlTaWkpuRheJlmFEN8CmnvI\nCVQAu4EocDewQUp5a+5xN998s3377be7qeoIzc3NeRs9aKjgqYIjaE+n0Z6O02VCg2tPAkKIa4DP\nZovnSSl39vYZKeWu7I9xIcR9wJfzHadzAs6igqcKjqA9nUZ7uo9rQUBKeRdwV18+I4SokFLuEkIY\nZJLIdfmOU2mIqAqbTajgqYIjaE+n0Z7u49UQ0XLgDaAMSAshbgDmSCmbhBBPAZ/JPik8JIQYR+aR\nZQVwVb7z6cSws6jgqYIjaE+n0Z7u40kQkFLuBqq6ee+8Tj+f6YWPRqPRaDIoOWNY5wScRQVPFRxB\nezqN9nQfJYOASjkBFVDBUwVH0J5Ooz3dR8kgoHMCzqKCpwqOoD2dRnu6j5JBQKPRaDTOoGQQCIXU\n0FZljXEVPFVwBO3pNNrTfdT4NlUUVXIXKniq4Aja02m0p/soGQRUyQmosvm0Cp4qOIL2dBrt6T5K\nBqIysk0AAAe2SURBVAGNRqPROIOSQUDnBJxFBU8VHEF7Oo32dB81vk01Go1G4wpKBgFVcgLta3gH\nHRU8VXAE7ek02tN9lAwCGo1Go3EGJYOAKjkBVZaWVcFTBUfQnk6jPd1HjW9TRTGMLpv4BBIVPFVw\nBO3pNNrTfZQMAjon4CwqeKrgCNrTabSn+ygZBDQajUbjDEoGAVVyAuGwJ3v2DBgVPFVwBO3pNNrT\nfdT4Ns1Blf43VTa/UcFTBUfQnk6jPd3Hqz2GLwNuIrN38CHgainlyjzHTQN+B4wB3gQ+IaXs0tmm\nygYO8XhciVEDKniq4Aja02m0p/t49SSwCThNSlkDfAe4u5vjvgf8UEo5E2gAPu2Rn0aj0QxJPAkC\nUsqXpZQN2eKr5Nl0XghhAGcCj2Zf+g1wQb7z7d692w1Nx3nxxRf9VigIFTxVcATt6TTa01mEEB/I\nfc2PnMCngf/L8/oYoFFK2b5P23agMt8JVAkCL7zwgt8KBaGCpwqOoD2dRns6TpcgYHi5GYIQ4gzg\nZ8AiKeX+nPfGAq9mu4IQQkwC/k9KWZ3nPBuBeLa4Bdjqqnj/mUxw3TqjgqcKjqA9nUZ7OkuxlPIT\nnV9wLTEshLgG+Gy2eB4wFvg1cG5uAMiyHxgphAhnnwaqgB35zi2lnO6Cskaj0Qw5XOsOklLeJaWc\nL6WcTybYPEZmtM+6bo63gX8AF2Vfuhx4wi0/jUaj0XiXE/gGmT7/nwkhVggh3mh/QwjxlBBiYrZ4\nE/AfQoj12ePv8chPo9FohiSe5gScQghxAvAKcLGU8tHejvcaIcT5ZIbCpoEkcIOUMnDDBwqdv+E3\nQojZwH3AAmCJlPIOn5XyIoQ4B7gTMIFfSylv91mpC0KIe4H3A/X58m1BIJsPfACYANjA3VLKO/21\n6ooQogh4HoiR6e14VEr5TX+tukcIYQJvADuklO9vf125GcPZinwPeNpvlx54Djgm2xV2JZlcSBAp\ndP6G3xwArgcC+eUPHfflXcC5wBzgEiHEHH+t8nI/cI7fEr2QBL4kpZwDLASuCWhbxoEzpZTHAPOB\nc4QQC3126okvAmtzX1QuCADXAX8A6v0W6Q4pZXM2xwEwjMxfM4GjkPkbQUBKWS+lXAZYfrv0wInA\neinlxuws998B5/vs1AUp5fNkgmpgkVLuklIuz/58iMwXV97h4n4ipbSllM3ZYiT7XyD/rQshqoD3\nkecPUqVWPRJCVAIXAmcAJ/is0yNCiAuB24DxZBo/6HQ3f0NTGJXAtk7l7cBJPrkMGoQQU4Fjgdd8\nVslL9gnwTWAmcJeUMpCewI+ArwLDc99Q7UngR8BNUsrAbyggpfyjlHI2mVnP3/Hbpyey8zc+TSY/\noNEEAiFEKZmn/huklE1+++RDSpnKdvtWAScKIQKXZxFCtOeA3sz3fuCfBHLmG4wAfieEgMy8g/OE\nEEkp5eN++bWTOy9CSrkTMo/fQojpQoixUsp9/hlm6Mf8DV/orj0DzA5gUqdyt/NcNL0jhIiQCQAP\nSSkf89unN6SUjUKIf5DJt9T57ZPDqcAHhRDnAUVAmRDit1LKj4Oio4MAhBD3A38O6OigmcAGKaUt\nhFgA/Amo6pQnCARCiMnA34FPSilf9tunN4QQ3wKagzg6SAgRBtYBZ5H58l8GXCqlXO2rWB6yXSx/\nDvDoIIPM2mEHpJQ3+O3THUKIcYCVDQDFZAarfE9K+Wef1bpFCHE68OXOo4MC/ySgKB8GPimEsIBW\n4GNBCwBZOs/fAEhKKY/3V6krQohyMkPbyoC0EOIGYE6QugiklEkhxLXA38gMEb03oAHgYeB0YKwQ\nYjvwTSll0ObjnAp8AqgVQqzIvnaLlPIpH53yUQH8JpsXCAG/D3IA6A5lnwQ0Go1GM3BUSwxrNBqN\nxkF0ENBoNJohjA4CGo1GM4TRQUCj0WiGMDoIaDQazRBGBwGNRqMZwuggoNFoNEMYPVlMo+kn2Vmi\nfyWznHCqm2P+Adwupfxbp9duAASZpX2fzX4+6YGyRtMF/SSg0fSfK4HHugsAWR4GLs557WLg4eyS\n088BH3PJT6PpFf0koNHkkLu2jhDiy0CplPJbOYdeBlza6XMfJ7P5TZTM0sdfAB4F/ksIEZVSJrLn\nngi8kP3Y42SWHH/IrfpoND2hnwQ0mn4ghIgC06WUm7Plo8n8RX9qdmnhFHCZlPIA8DqZHccg8xTw\n+05rSdUR8L0xNIMbHQQ0mv4xFmjsVD4LOA5Yll307Cxgeva9zl1CF2fLQGY9eiAhhOiy2YdG4wW6\nO0ijyY/R6edInvdbyazN3vn430gpv5bn2CeAH2aXFS/Js7lHDGgbiKxG01/0k4BGk58pQohxQogQ\nsJjM8tAdZPdmNoUQ7YHgOeAiIcR4ACHEaCHElOyxzcA/gHvp9BSQPW4MsE9KGeT9kzWDGB0ENJr8\n7AceILOPQR2Z/SFm5BzzNLAIQEq5Bvg68LQQYhXwDJn15tt5GDiGnCBAZr/svzhur9EUiO4O0mjy\nc0hKeW6n8lfyHHMXcCOZsf5IKR8BHsl3suwWqEaety4Fbh6YqkbTf/STgEbTT6SUy4F/ZHeW6jPZ\nEUaPSynXOWum0RSO3llMo9FohjD6SUCj0WiGMDoIaDQazRBGBwGNRqMZwvx/Ve6z0KbLeVwAAAAA\nSUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f83c65532e8>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"sigmaxx.plot()\n",
"sigmaxy.plot()\n",
"grid()\n",
"ylim(-2,2)"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.5.2"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment