Skip to content

Instantly share code, notes, and snippets.

@LuxXx
Last active January 29, 2020 11:02
Show Gist options
  • Save LuxXx/cee38a4ddceb8fcef1d8c4e6f75a1c07 to your computer and use it in GitHub Desktop.
Save LuxXx/cee38a4ddceb8fcef1d8c4e6f75a1c07 to your computer and use it in GitHub Desktop.
Diskrete Brownsche Brücke
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [],
"source": [
"import numpy as np\n",
"import pandas as pd\n",
"import matplotlib.pyplot as plt"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [],
"source": [
"N = 10000\n",
"Xn = np.random.normal(0.2, 3, N)\n",
"Sn = Xn.cumsum()"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD8CAYAAAB0IB+mAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzt3Xd4FFX7N/DvSUJC76GXgPQiCJEiitKbiv0FURALNuwNVHwUbD/7o4+FYkMRKUoREAREKQISekIooYdQQi+BkHLeP3Y2md2d2Sk7bWfuz3VxkZ2dnTmzM3vPmVMZ5xyEEELcL8buBBBCCLEGBXxCCPEICviEEOIRFPAJIcQjKOATQohHUMAnhBCPoIBPCCEeQQGfEEI8ggI+IYR4RJzdCRCrWrUqT0pKsjsZhBASVdavX3+cc56otJ6jAn5SUhJSUlLsTgYhhEQVxth+NetRkQ4hhHgEBXxCCPEICviEEOIRFPAJIcQjKOATQohHUMAnhBCPoIBPCCEeQQFfpa2ZZ7D54Gm7k0EIIbpRwJcwfd1BjJgc2AHspv+txMDPV9mUIkIIiZyjeto6xYu/bLE7CYQQYjjK4RNCiEdQwCeEEI+ggE8IIR5BAZ8QQjyCAj4hhHgEBXxCXObYuUu4lFdgdzKIA1HANxjnHAtTD+NyfqHdSSEe1eGtpXj4h/V2J4M4kCEBnzH2DWPsGGMsVbSsMmNsMWNsl/B/JSP25XTLdx3HIz9uwMdLdtqdFOJhf+/MtnR/F3LzsTv7vKX7JNoZlcP/DkDfoGWjACzlnDcGsFR47XqnLlwGAGSdvmhzSkg0WZh6BDuPnrM7Gbrd/9069Pjwb7uTQRQYEvA558sBnAxaPBDA98Lf3wO4xYh9OR0HtzsJJAo98uN69P54ud3J0G3t3uCfP3EiM8vwq3PODwOA8H81qZUYYyMYYymMsZTsbGsfQ83E7E6ARhdy87H3+AW7k0EIMZHtlbac8wmc82TOeXJiYqLdyfGs4d+tQ7cP/rI7GboljZqP56ZvtjsZhDiamQH/KGOsJgAI/x8zcV+OwaO0ROdfFzyS/7Ih0+4kEOJoZgb8uQCGCX8PAzDHxH05DmPRVqhDCHE7o5plTgWwGkBTxlgmY+wBAO8C6MUY2wWgl/CaENv9sj4TczYdsjsZhFjOkPHwOeeDZd7qYcT2neSdBekY3b+53ckgEXhuhq+sf2Db2janhBBr2V5p62T/7D4esmz88j1hPxOtZfjEeu8sSMeHf+ywOxnERv5+O1ZxVcA/fj4XP67Zb9j27p64FofP6OtARSX4RMn45Xvw2Z8ZdieDWGjDgVPYccTXwW5h6mFcNW4x1u2zrsGEqwL+41M24NXZqdhjYBfvnMuhg1C9NX+bYduXs/f4BXz/zz7T9+N1zcb8jiembrQ7GcQjbvviH/T5xNfBbs0eX6DfmnnGsv27KuCfFB6P8grMLVeZuGKv7HtG7fn2L//Bf+amIa9A+yBsa/acwLFzlwxKiTXSss7gkA3DUVzKK8Rvm7Ms3y/RrrCQY/uRs3YnA8/P2IwBn66wOxm6uCrgm9ESUneZfIRpOXsxT/dnB01Yg1s//0fVujmX81FYaH/Fw4BPV6LLu3/anQziYBNW7EHfT1Zg08HTtqZj5vpMpGXZf+PRw1UB389N49noveGoyS1fvFyAFq8twju/p+vbCSEW2pLpC/SHTrlrYEIru+y4KuAzF1WVWnERnM/NBwDM2kht0s02a2MmMk/l2J0MW6zKOI58HUWTxHiuCvh+xjaNtPdpwU1PK17FOccz0zbjti/UFbO5yT8ZxzFk0lrbWyMdPJkD7rA203akx1UB358rtvO8GnUS5Z5WJq3YY+q4N2ZfhJfyCrD/RPhROXPzC5BzOd/UdNjh2LlcW/efdfqi5TOx+Y/ZzpFYNx08jeveW4Ypaw+o/kx+QaGp1+BPaw9gu9A808pyCVcFfDPojX+/bjCnmOTN+em4a/xqQ7dp5Q1y5E8bcP37f4V9xO/98XK0eG2RdYmKUn/vzEbSqPlF7brDybmcj2ve/ROjftliQcp8jp69hH0KN3c9lJ56x83bhqRR84te+5tpr99/SvU+Rv600dRr8OVZW22ZQ8CVAd/OYhCj9xxpMD5xPhdJo+bj539Dczd2jO/21w7fnAfhDmv/Cemy7gu5+fhlfabjHs3tsjD1CAAgZb9y4LiU57vBLtth3aC1Hd9eik+W7AIQeL53Z5/Hg9+nIDdf20Trauvovl4p32xarYVpRyLehhO5MuBHs+kpB/HNyr2GPecdOOkLnlPXHVS1vpmxNDe/APkRNAH9z9w0PDdjc9icGuccny/LsLzLulZOnGzmcn4hTueY/72NmZ2KJelHkbJPfY4bKM7IFRRyLN521LAb//Hzubgo0cFSyfSUg5ieou53FY6VI+tSwFdgRvz7bXMWkkbNx8GToTnZF2duwdh55vfkNcvZS3my7fq/WLY7om0fPevrTHYhzI9z9Z4TeH/RDoz+dWtE+zLbfd/+a9i2tMQ9DmD2xkOSueuHJqeg7djFhqUrYL8G5iQmrtiDhyanYJHKXLjSrpPfXIJbv1ilYjscHywqHvvoxZlb8OJM64rIjOCqgO+vBHH6E/9soRmkv/I1N78gpEet0j0/r6AQP6zZjwIDOk1xmb+1Op1zGVe+/gc+XrJT8v0zEXQmU8vfy/pCBBVuaVlnkDRqPlZlhA6eF+xcmBtcOPkm9wYP5g+4p3Py8PS0Tfjoj9Bz9PdO5SlGL14uwF1frVasNwgOxkYerb8dfvZ5bU8j4X5T21XUg5zOycP/lkX32EeuCvhKXpq5BSMmp9idjKJijedmbMa4edvQ9NWF6PHh35Lryt28JizfgzGzUzFNZVFNsIxj57Bmzwldn/Wli+Ot+duw8UDxY/lx4Qc4f8th3dt1grXCGCeLtx0Nu96pC5fR+vU/8InoBrdk21GcyTH/xiamVCKwMPUI2r+5JGCZUouhWRszkSXRee/ffSfx776TeFNhPKmHf1gfPlERcHh+TjMr69IMGQ8/WkwTytvO5+ajbIJJh67iaiwURXF/BdOBoOIdpYvAn1s+d0lbcLnps5XYd/wCzuVG1uSskPvGFJq4Yi/2vTsg8E339H8Lyz/K4fyth/Fs76bIPpeLByenoHPDKpg6opNl6VB6oh2nsYgwN78Az0zbjLqVS2HFi90jSJmIAVG6UHgIPm3xDdVNPJXD9xOXw7ld8O9s66EzYYO92rJWtTH9z+1H8fFi6SIeM1h5rxkRlIv1l4sH37yDZRw7r7mFSiSUzmlBIcfrc9NE6/v+P3Y28CngdM5l5OYZl26tRa9yxXRHzlxC0qj5+HO79BNZJPcat7UIc2UOX+kcXdJ40WqZDk9Nk1C9tfKfR1B+GG6Pai7q/IJCxDCGmJiwWwpZcv93viK0Z3o1Cdqn4i4BAP3+uwLZ5y4h5dVe6j5ggXOX8nD7l8W9ZndnX5D9DtOyzmDlruK6gDMX89DzI+niO738l9PbC9Lxy3p1E7mL0/vv3pP4TsVQ3EZV6BpdhOEfY+entQfRvVn1MDvWvu1CDsS66InVnQHfwFI+zoGnft5k2PYi8b7Mk0lhIQdj0jeSGUIxVva5XFzQUYyTczkfnd/5E2cu5qFjg8qY9nBnxc/4U2FEhXL6Yd+ohGN/24YVu5QrUf3MzJgtSjuKnUfPhyxrVbt8yLoDPl0Z8FpP879g5y7lISEuNmDZ7I2HMEFhNjY5cr+XQgO/xGgdIsR3Y3RPxPdkkY6Z1HQO0XL5qPmhNHx5AUbKTOIx9V9fwD90+iL6qxjDO3hvczdlFdUXaO0ZOFvHoGxyLUW+WRXYmWbzwdPIuZyPdftOhn0CO5OTh+dnbNZ1s5PzvDAnrtgRnTOjqcE5x+9bDxcNi9D69T9wz9drRe8DT08zPlNi5LwSC7YewY2f+a4/qfvIb5uzLB219eylPDz4/TpkGzDcxQszNuP3reEbKgyZtAbv/r494n1FyjUB3yllbUblZNT2KvTvTU3LGLkerOGKmMK1eQ9JS9Chh2saKbfLYd8ot08/ezEPAz9fhSenbsSdX60OeAIbGvT5z/7chZnrMzFlrbqpL6es3Y+v/vb1Fzh7KQ87jyo31wPUlxO/Mkt7/4AVu47j0SkbAua/1TKeklTawp1zqf4hRkg9FDiG/OEzF9Fg9HxsOngaT0zdiPF/h39CCU6z2oyTVGyYvu4glqQfUxxqYvG2o4qd0Wasz8SjUzaEXWdVxomi6yoYjaWjw27RtIZGxn67H0WtuI+JfxBm7s+oslt/HcxmDVPDSd1ApWYTe2VWalGTxV83HELvj5frTKV0oFm6PfzQBiclegifEgJOukxbcaPLxHupOGZ/8D1+Plf3vM8vzNwCzoH/yvTbsMIlUeV5YSEPOWePTtlQVDl/VmOLOCdyTRm+WZM2/ZEWvi12MKMCZvCPODgQ+Cc4sat00RnPU5GZvNq4Ce+lzvsPa7RvP3hQuV1Hz+EXYSC+5TLFXXquuUifiP2fTw5q36/k4Mkc/LNbf/8PybQovK+2kUTDlxfgjvZ18MGdbQKW+594ZqSoqxB3Mtfk8MWULwD12/rIhCaFWnNkJy9cRrtxgS0k/EU4UjnCSO0/caGo96iWwOBfU+4H9u2qfRGmTD01RTHnL2kr199/4kLRpDFqzDOgA1qvj5fLBnq1Dp9RmN9YR+zX2xa+6/vLQpYFXy/vLEjHt0F1NuG8MitV2I6uJAWYKdHKSfH7ixSNpWOuwkLg+3/2aW6eqdUbv6UFDNPqJ3d6pcqvOYCTF+QrloxoSXFK9OPdefQcrn//L3wpU96ohtTx/RHc1d7kR4TgnqQcHJNX78NtKsZMkXP9+3/h/4UZmlp8TLM2Zuqbt0DHb9/MeLFEprfx1kPqi9PEpM57cPLHL9+DN35T31ns+PnwFa9Oqd+Ts9/CgfRcGfDFJzi/oBDHzgXeoX9PPYz/zE2Tbeao1tL0o2EnSZDL0crlgMUtVOwqqvEXFUUyycquY+dDlgV3UrLC4TMXMUk0VO5rc9Kw4UDxBNha6mf815DayaufmRbakscsZsaz7Ufkj/eOL42ZwcuqDK5UPY6WaVGVxtHxNwzYe/xCwETrShnLSSv3YptFk6K7M+CL/h4zJw0d3loa0CwvR2h5EslgXjuPnsMD36fgZWFUxnlbsrBMoUJOD865LYPBcfjawL85X7mp3JS1+5E0aj5y84ydTek7DY/1Ukb+VNxU9ePFu0Le94/VDgBJo+aHnQ1q0Pg1YffFYc/8AobQme4UiWGqL+cXhrSUMoqeZK7bdxKXLJrl65VZqcgvKES3D/7CLZ8XP0mq6Xvx7HRr+vq4MuCLLd7mK0rIEX3pRsTPc0L5736hQmfkTxsx/Lt1sttemHoY09b5JiFRc+FqaQ7pN/zbfw0tpvo9Vd3ws/7AefqisfUJrys81ivdCMXDF1xU8b2Eqw/ZY9Fjd4e3lkq2HgpH7kaTV1CIN35Lk37TJOmHz+qoc9B3x0nZdzJg8L5/gkY33X/iAu78ajXGzE4tWrZ+/8mQ9dQYoKIPC2BtPZUergz4XyzLCDuFnpR/Mo4jLUtfuaQaj/y4AS/9or0NNgfw1gL5XLY46C3bkY1mYxbqSJ309uQiKuccP6pogaL0M865nI9PluzUfK7MmjAikia4RpYTrw7TikVqP3K7nrspS1UAsruIO9zpfGdBOjq9vbRoLgSx2ZuycKtoYvgLlwuQISpOvOmzlSGfuf3L1bh7UnGnNal9Z0gUSaotyhP/VvMLCiXTbSfXNMsUW5J+DL9uOIS7rq4r+b7U9eW/CEJGfpTl+5VsFJUJA8b/ePLyC4umBZQyyYDp3ORIHUrGsXPo+ZF0O22tx/7x4p2YuGIvalUopS1dqgd403ZjOHvRGROn71fZ8Unpvmfk0AhqGb3H8cJwET0+/Bu1KypfJ09PKy7GOxvUCkttPsGop+Sx87Zh8ur9+PuFGwzZnhFcmcMHijtU+Mdol6rJN+P38LKOnpThBI9jbiWpbufDvlmn+vPhvl7GiovZcjXm8MXbMJLcEA1qZkPKPHWx6FqzKtDanTM3iprTeD43HztUNLUN7s1rp6Xpvjq9cxqb/5rJlTl8KeLmZVFbuWYhzjl+lphc5ZDEpBhFnzEzQRKUxkHR2nRQbn7S4Kc4Kd/9s69oxMmI221zjs0HTyNWYmRScZCXGw44mkmNUxTt7vtWOZNk1by2rg34anI/4dq3E/WkvuszOXkKQzKblhzdjmucMs8sB07mYMwc5crWo2fDX79WTo4dCXEypTo+Wc3or02pn4CVXFOkExxAjp69hCteXhD2M8t2ZOORH9YHTOUmHpMnnNu/lO+Ao8TZv0PfF6knWGj5xKmcy5iy9oDmfXjBxBXG1MsYUZG86aC2p6QHvlNf5Oenta5FL6t/dk7s8OXaHP7yXdmqxmNfmHYECSWK73tyc8tGIrS9v6MjPgB9F6uWTwQUkzjwh2GlWRvV52q1fFNqb9rh1lqSfhQPfp+CO5PrqNrWCR1Dfcgl0+h+LVLfnZY5FrTKMntIBh1cG/DDxRBfjsK4IDP13/A51Tu/CuyR6OwcvjE4uOrBw7TmaKOlqEKt56arL7c+pTBUry4KX+eS9KNYkm5efYFcXcxwHU8L4azS0f7ebVxTpBMcA4LbzZqZhxz9a/iWOcGzI4UbjsFuuTp6Jfq/e/FTwaqMEyHHLUdpDlhSrO8n6joAAeGf0jh8xZ7Nxyy0rFu/HKkeu2YwfRC0CFiVhTE9h88Y2wfgHIACAPmc82Sz9yll0grR5ApB3+5vm7MsS8fGA6ewKsPY4WGN5J9MRG1HE6D4aUocXpRualp7k4q5sSWHWlIVgIvSpHtEKzUH/HP7MVzMKzB0mOhoEa5D2koTi3nsZlWRTjfOua3fYnAnDDGzxtKXIu4ZaLa9EQwHoGXYZX8Q0lLuH27cGq9hjEVUjyFXDj12nroRJ734hBWuefE7DpiK0CyuKdLRwivBptsHf1m6v/u/Syn6+4iDH58J8SorAj4H8AdjbD1jbETwm4yxEYyxFMZYSnZ2ZBM9EOf4UGHimOc8XCzjBAzR0FaMGM2KgN+Fc94OQD8AjzPGuorf5JxP4Jwnc86TExMTLUgOIc5iR+D1dkNY59l22CXj4XPOs4T/jwGYBaCD2fskhBASytSAzxgrwxgr5/8bQG8AqeE/pY/Z0xUSYpYCj3c8I9Yxu5VOdQCzhI4ycQB+4pxHNmC7jIdtmEKPECNQvCdWMTXgc873AGhj5j78nNypghCnOXAyB3M2Wdf/hDiDJ5tlEuJ1aoZ8Ju5DAZ8QQjyCAj4hhHiEKwL+jiPKU58RQojXuSLgp2Vpm6SBEEK8yBUBnxBCiDJXBHxqx0wIIcpcEfAJIYQoc0XApww+ISTaWTHpuSsCPiGEEGUU8AkhxAGsqIt0RcC34lGIEELMZMVk7q4I+IQQEu0yjp03fR8U8AkhxAG4Bc1PKOATQogDFFIZPiGEeAM1y1SJqmwJIdGu0IIsvisCPiGERLux87aZvg8K+IQQ4gBUhq8WlekQQogidwR8QgghilwR8K1ov0oIIdHOFQGfEEKIMgr4hBDiERTwCSHEIyjgE0KIR7gi4NPoyIQQoswVAZ8QQogyCviEEOIRrgj4VKJDCCHKXBHwCSGEKHNFwKdKW0IIUeaKgE8IIUSZKwI+jaVDCCHKXBHwCSGEKHNFwKcyfEIIUeaOgG93AgghJAqYHvAZY30ZYzsYYxmMsVGm7ISy+IQQosjUgM8YiwXwOYB+AFoAGMwYa2HmPgkhhEgzO4ffAUAG53wP5/wygJ8BDDR6J5S/J4QQZWYH/NoADopeZwrLijDGRjDGUhhjKdnZ2bp2QiU6hBCizOyAzySWBYRnzvkEznky5zw5MTHR5OQQQoh3mR3wMwHUFb2uAyDL6J1wyuITQogiswP+OgCNGWMNGGPxAAYBmGv0TijcE0KIsjgzN845z2eMjQSwCEAsgG8452lG7+fAyRyjN0kIIa5jasAHAM75AgALzNzH5oOnzdw8IYS4AvW0JYQQj3BFwC+kiE8IIYpcEfCpIT4hhChzRcCncE8IIcpcEfAJIYQoc0XApxIdQghR5oqAX7NCSbuTQAghjueKgP/oDVfYnQRCCHE8VwT82BipMdoIIYSIuSLgM8lBOQkhhIi5IuBzaphJCCGKXBHwCSGEKHNFwKdmmYSQaBcfa344dkXAJ4SQaPd4t0am78MVAZ8y+ISQaGdFY0NXBHxCCCHKXBHwaU5bQghR5oqATwghRBkFfEII8QgK+IQQ4hEU8AkhxAEYtdJRh6psCSHRzoq2J+4I+BTxCSFRrnalUqbvwxUBn/L4hJBo16ZuRdP34YqAX0jxnhCiQo9m1exOgiwrBnl3RcCnIh1CiBrDuzSwOwm2cknAp4hPCFHm9bkzXBHwqUiHEKJG+ZIl7E6CrVwR8AkhRA0rKkadzBUB3+uPaYQQooY7Aj7Fe0IIUUQBnxBCPMIdAZ+KdAghBnrtxhaW75NZMJiOOwK+QfF++sOdjdkQISRqtalbEb1bVrc7GaZwRcAv1BDxy5eMk32vQ4PKRiSHEMd4qW8zu5MQdf43+Cq7k2AaVwR8LRn865okom/LGqalhZCr6pnf9K+dyn1UKOW+dudjB7a0OwmmqF+5tOn7MC3gM8ZeZ4wdYoxtEv71N2tfVIRPnGTWY11M30c9lcHBjfVb93SsH/D62kZVbUqJsWJior8M/2POeVvh3wKzdqKlSAewZqIBt1AbWMx2dVIlu5PgKFZU8DlR1bIJIYGxaxP7A/7ckebf5I3gjiIdDfG+fb1Ktgf8hDhnf+1v3tLK7iSEYAaPJfjYDVcYuj1ivu3j+mLVqG4hy42+NhgDYjQEiTdubokr61TEvncHGJoOM5gdeUYyxrYwxr5hjJmWRbu6QWXUKF9Scb2Hr2+I4V2SDL9AtPrxwY627l/JdY2Lc0x23RwfuNa8UQ2rlUvAi1Fememl/P01V1QB4MsoJcTFhrxv9DVasXQ8alZQjicA0K1pIoZdk6R7X9XKJej+rB4RBXzG2BLGWKrEv4EAvgRwBYC2AA4D+FBmGyMYYymMsZTs7Gxd6ahQqgTWvNwDpeNDLwax+pXL+B6FLf61tKxVPuB1UpUy1iYgCt3dsZ7dSSAOMXFoMhY/0zWgGOvpno1l1+/aJDGi/ZVNiNNdZBZrQTl8JCIK+JzznpzzVhL/5nDOj3LOCzjnhQAmAuggs40JnPNkznlyYmJkJ0otZ58S+4mLyNR+V6P7mZxjppOmS4ua5ZVXcrgyCXFoXL1cwLKnezaRXb9UCfOLTL8eliy5PFbjjcLqKnUzW+nUFL28FUCqWfvy8weqezvVx08PdkTGW/0wcaj0iSHRpXZF+fk+Zz7iwQ5zKuLKqlHdcVW9SvhiSDu8fWtr89MkQW2z0J7N9c9EVT2oONfKIttoqzw381b4HmNsK2NsC4BuAJ4xcV8BRvdvhmsaVUVcbAwc/oQVwmmPhNdH+HisV3BF/O3t6siuq3XyZzc0VFQT1Pw3yf6ta9pWRKbmcq5aNh7/HaS/s1OVsvG6P+sErw5ojmGd6yuvaAD5bqcR4pzfa9a29fK3SVZbA/9Cn6Z4f9EOM5MUonbFUjhwMsfSfYqVCaoHGXNjC3S+ogpKxcdhb/Z5vP7bNlvSFRcrf85qVtAW8KNZXAxDvooZfwZ3qGtBaoxRsXQ8yiQYF4pszXSr2Hf7+pWwfv+potcPXtfQxAQFcnb7QJOovSAe79bIsP09qKLVyf/d3lqx4tlIvz52DTokBQ4nUSuo6CQuNgZ9W9XE9U0ScV+XBuguMwm0kT+ypkHltYDvR0KAzkKLlXJhhggZc2MLvHPblZakp2HVMqgaYQ470ilKS5YI/M3YOXqu05tcOzt1BisRa9/hjhJVaiaIKpViYxie7dUEW1/vjf93dT2M7t887HYaVA1s4VO/iv6OUe3qVcKEoe0DlukN3Eb9yJ7s3ghzJDqxaGkX7Wb+7znc12HVNzV2YEv8+fwNSHm1l6HbHTewJX4e0Un1+lfVrRjQMc/KSyX4ZnVlnQrW7VwHTwX8W6+qDcCcH8S+dwfgiyHtAABrRvcIeT9OdLMpX7IEZj/eBZtf6430sX3xZI/GKCfMtXl9k0S8d4d87mzy/YGNnayulDb7t1S+VImQHJtXNawa2nzXSUMllFVRDKOmUjN4nXs7J6FTwyqq08EYw+AOxXUUcpkPIztGyR2W1gpjq59GPBXw/Tl88QWWpJBDrlmhJOJUVqT2b10T+94dgBphOm20qu1rJte2bkVUKF0C8TY/AgZfcAwsbEgJfk9tB5Vg1zaqioeuC1fMZd4vIVwbbjWe7N4o7E1ZidpWRYnlEkKe6CJRuUx0V25Gg6uDiki7NQ1t9CDu2Gg1TwV8P/HYG0pFBSte7Ib0cX0N2e+2sX3w66PKY26ES1Fwcu0u6JB7hE15tWfYz7WrXwmvDAidZMLsZm7jbmmFIR0jaxExvEsD3JWsv1I0OUndMNyMAU2qlw1YVi5B/+iX4YYGN4uW8vl5T1yLz+9uV/R6ybNdMcPkJrdGjyb6RPdGWPb8DWHXebJ7ZBmOSHgy4DdK9FUKtqxVXjFixsXGRFz273/MKx0fF3GOXk9ADNe0MuQGonHzcoOrVS0baZfx4oSoqQi7T2X3djuCnpR/Xyku9guX4xN/v5/f3Q7JQlm1mW3NV74UOl6NlCsSi29GE+5tj/duD33qCW4EIMV/JK1qV8CAK4u77zSqVi4kx2w0tceqVkwMU3wqCxz8zdoyHVcGfHGmQiqDYXX9n9Zy13Brl9ZRvl0yTM/DSMsQH72hEfq1qoHBBrXzljo1as5XKZWtm25uU0tbggz26gBfpXy1cspFYQwMz/VuWvRaHAzllCoRiz6tpOd7kLspPtUjMMdZp1Jp3NymluSNVlws1KZu8Zj8vVvWwF1Xhz719G+tnGajafm9+evOvMJVAd+tDTnEnbEq6SiHFecG31dR9hy2SCnodeUy8fjynvYor/WHI3OnKW/AI/aXQ9rJvmd3z0itba61Pl2mj+sr2ytAgQpsAAAQHUlEQVT5vi4NJCsupZp4fjr4Kux4s1/I8iHCjV3uqWvi0GTMe+JaLUkmFnJVwA+XW+0h037c6V67sUXYyjY18Uuc41EKeErNQs307m2tcZvQkkpMzVNIddGog/1MyFVWKi1/Ixqi4ummZ/NqmiuLbRjnL8TCp6/DxKHJRbl5pfT0alEdrWqrb5rYpk4Fy4fj3vxa74i3ofbJeGjnpLDvv3mLtUNeuCrg+1mViRNXMJnl/msbWPqjt2sohfi4GAzqUE9UvqmtrMns2YKubZxY9N2UCMrdvnVra8XhnP9zU8uwA35JCVd3YdU13qxGefRqUR2d/c0kDd7xnJHXoqOGJph6vSLKyFQIc/PWSi4D5e9AqTRPdl+Z4jezuDLgqxVp+bWaMlXAusGc9DSRlPoKLKlG0hg47Mzpjr+3Pd69rTU+H9IOsx/voqr9uV6PiiZmef/ONqbtx23El5P/dy2+YT7U1brhC9ToqHAjMIvrA/4V1XwtCXq2qB7RduyY+PyO9r4Bw8LHxuI3V4/uYXrOz59j+XpYMv54pqumz1YpEx921EstpNo3m6VPyxookxCHsglxaFtX2wTlWo9XPJFP5C2djOe/vKwcAkQNceVwJL+Bhona+j2oaXY65/EuWPR04G/lu+EdsGpUd037MoLrA36DqmWQ9kYfDJJoQQCE5hzl5qb86t7iIQi+uS8Z3953tVFJDNGmji+o+IeMjbQpp5G9+UZ0bYjlL3RDj+bV0URizJtw1o/phTuT5Ue91OqJ7o0wXjgvdStZP/fuqlHdDa2gvKdTPdzbKXwfAfG5fPOWVgGT60zROZPa3R3r4S6V5yWhRAxe6d8cMx+9RnFdpZud2n2qIZ4JKz4uFhVKlcC4gdrqBn4e0QkzHlbX7l/tTYXD15qpaY1yaFStLPq09GU8S8XHGpb50cIZjZJNpmUkvivrKOfgujeL7GlBSdMa5bDrrX5FLTQm398R3T74y5BtB+dItA5cxRhDvQjG75Hdro7PiJssdrOhUr52xVKKP1otuc1eLWqorotgAO7pVB/3dKqPpFHzAQBdGqnrwfnZ4Kuw8+g5fPZnBgBf/5D37miD6SmZsp8RV/yrKR5Z9vwNqFw6Hr9tyZJ8f/fb/Q0furxn82pYkn4MHRtUxmeDfcMtv/jLFtWf1zKcgxKpQ1vy7PWGbV8vTwR8OXI/xhrlS6J6eeMep/U8Yoqb42npXs9g33jvH97ZBs/N2GzItupWLo1aFUoi68wlQ7ZntVmPXYMqZRJw96Q1hm63mnBd1owgd3iT0BfBH/C1UFsf5b9mB3eoh9gYhtG/bg1434x5HyYNuxqHTl9ELZ3DfXiB64t09Fjzcg/MGRkdbYnV3EyC15n6UCcsf8HYHoYAcHt76Uf02BhWNORxu3qVhP/DP0klxMVi5UvdUSY+Fq/f3NLYhJog+EHpqnqVDHkSCj53N7ephUlDkzFc1IlKPHCY08TGMEvTV7tiKdV9LSId1jkc5wxxF8jTOXw7x82OxH8HtcVTP29SvX7wcfrHVAdCxxI3w643+xUFrq5NErFhTC9VA3nFxDCkjfWNY1SoYtKPaKcmTDHGQhogvHNba7xzmz1TGEar6Q93Vhw40Qh296MIRjl8jUqEmXnJKgPb1i7qCKN10uRgZRLiTMnti8XEsIBclxWjNto5IiGxxtyRXfCGzqe/Dg0qo1p57xX9uCrgP9/HV4kXr6E7+gPXqm+fu/Dp67DqJe1NqaoYEOCGda4f0AV+wr3tMebGFkjSOHyu1ONucNFDpDMQ6fHrY8qtPrT44QF9LVb0krvvRtNT5O9PXYepD8lMPOLA47iyTkUMUzlonhEedlhbfj1cVaTzwLUNFHs8Bru7Yz3VTQWb1SivvFKQmhVK4qO72mr+XLA3BrbCG6JmZtXLl9R8rIC2YG7kWOxKWtZy9kxBXtC8pvL17dbxqtQY3b85xi/fY3cyIuKqHL5eJQwYAlnOoKvr6RrwTIs3bm6JupXlW204MHPmeEYVc3k5QEoZf297PNndmLmirdKpYWivWLX5Jqf99lyVw1crLtY/Pr2zegvqNeyapIgfbf96/gbk5hcakyATGBE4p6voVNO6dgVsPXRGcwubaCq6sVOfljXQx4Ze65H49r4OOJVzWfI9uRZBTr3PezLgN61eDi/0aVo0dIEbTRyajIcmpwBQd/FprQuIRkoDWQHAy/2bF7Vi6tm8miWTp4e7Wdg9nLOfl+9npeJjUSpe+gnajvquSHgy4DPG8Hi36Hqs1MqIsWai7WI2gji+ThqmffgM/wQnbuWM24/9HHIf1syTAd9IDRPLYE/2hZDlfVvWwMK0I5ZfGIyxkCyj1rBtVK5y+QvdcPxCriHbcrNoDR4k+lClbYTmP3EdNozpFbK8cdDk006hJvgblbOvV6V0Uc9aO02+vwNmPx46KF7/1qFlyUbH3hFCU74qJvbqtMKtwsQ0vaOs/F2Nv56/QXKcGyf3YNaLcvgR8pXvuaPyN5hTyo/16tKoCga2rY2uMpO6fDHEN9Kmf/CxSPg75AWPETO0c5LirEfRoHnN8pLTI7qBXP2VG3swU8AnrjXlQZlORGHovck90aMxLhcUasoV+ov9iPv4K/udVg9GAd8DHHbNOVopnWMLlU2Iw39u0tbN3z/Hwr1frw1Y/sWQdo6c/ISoN/2Rzpi7OcvU2dH0cFZqCLHRB3e2Qes61vf4Db4h9w+ahP3b4VcHTHRC1CtVIhZPaZw83gjNa5ZX1XPZahTwPUBcSuGWzmZmsLtfhlxpUrem5kzuMrJbI9SrbP1MYVZKH9fX7iQ4CgV8kzixGOWaK6qgfyt1E68T9/MPNki8g5plmswJ7Vz8N59h1ySpmkLPgfcqAL5WN4M7SM9NTIiV/GNvaZk+1Qko4LuUVNDWevOx4mZVNiFO1dj1jDFMebATerUwdz5hQtTockVVvNCnKd68RdtE6XaLrtsTUeSEJwotUt/oY3cSsOjprth3IrS3NDHOUz0a4+ylPLuTYZiYmOgcnoUCPokqTYS5cQe2rW3YNpvWKIemNcoZtj2t1BSzRbtnejWxOwkEFPA9oU4l30h/5UuVsDklkatTqbTreny+d/uV+Orv3ejcsIryyoREIKIyfMbYnYyxNMZYIWMsOei90YyxDMbYDsaY/c/tHjaqXzN8OaQdOlkQUGJjGLo3M6cZoVvVqFASr9/cEnEmTcJDiF+kOfxUALcBGC9eyBhrAWAQgJYAagFYwhhrwjkviHB/UYPb1Nblrqvr4qe1BwLGcS9ZIhb9WlvTHHP32/0t2Q8hRLuIAj7nPB2QHH9kIICfOee5APYyxjIAdACwOpL9RSOrxx8bN7AVXh3QHLExDI93uwJZpy9p3ob/ZhEfF305TjWTnBDiVWaV4dcGsEb0OlNYRkwWG8NQOt53Wl/o00zXNpKqlMZTPRrb3vNUq/SxfYumrySEhFIM+IyxJQCkBsF+hXM+R+5jEsskyzgYYyMAjACAevXcN/50NGKMRWWrCrcOU02IURQDPue8p47tZgIQd4msAyBLZvsTAEwAgOTkZKd28tSsfElfixinjZZHCPEus6LRXAA/McY+gq/StjGAf03alyPdf20DxMfF4J5O9e1OCiHEID+P6IRDpy7anQzdIgr4jLFbAXwGIBHAfMbYJs55H855GmNsOoBtAPIBPO6lFjqAb6yN4V0a2J0MQoiBrGjabKZIW+nMAjBL5r23ALwVyfYJIYQYJ/ra3RFCCNGFAj4hhHgEBXxCCPEICviEEOIRFPAJIcQjKOATQohHUMAnhBCPYJw7ZzQDxlg2gP06P14VwHEDkxMN6Ji9gY7ZGyI55vqc80SllRwV8CPBGEvhnCcrr+kedMzeQMfsDVYcMxXpEEKIR1DAJ4QQj3BTwJ9gdwJsQMfsDXTM3mD6MbumDJ8QQkh4bsrhE0IICcMVAZ8x1pcxtoMxlsEYG2V3evRijNVljC1jjKUzxtIYY08JyyszxhYzxnYJ/1cSljPG2KfCcW9hjLUTbWuYsP4uxtgwu45JLcZYLGNsI2NsnvC6AWNsrZD+aYyxeGF5gvA6Q3g/SbSN0cLyHYyxPvYciTqMsYqMsZmMse3C+e7s9vPMGHtGuK5TGWNTGWMl3XaeGWPfMMaOMcZSRcsMO6+MsfaMsa3CZz5ljGmbxJlzHtX/AMQC2A2gIYB4AJsBtLA7XTqPpSaAdsLf5QDsBNACwHsARgnLRwH4P+Hv/gB+h28O4U4A1grLKwPYI/xfSfi7kt3Hp3DszwL4CcA84fV0AIOEv78C8Kjw92MAvhL+HgRgmvB3C+HcJwBoIFwTsXYfV5jj/R7Ag8Lf8QAquvk8A6gNYC+AUqLze5/bzjOArgDaAUgVLTPsvMI3c2Bn4TO/A+inKX12f0EGfMGdASwSvR4NYLTd6TLo2OYA6AVgB4CawrKaAHYIf48HMFi0/g7h/cEAxouWB6zntH/wzXm8FEB3APOEi/k4gLjgcwxgEYDOwt9xwnos+LyL13PaPwDlheDHgpa79jwLAf+gEMTihPPcx43nGUBSUMA35LwK720XLQ9YT80/NxTp+C8kv0xhWVQTHmGvArAWQHXO+WEAEP6vJqwmd+zR9p18AuBFAIXC6yoATnPO84XX4vQXHZvw/hlh/Wg65oYAsgF8KxRjTWKMlYGLzzPn/BCADwAcAHAYvvO2Hu4+z35Gndfawt/By1VzQ8CXKsOK6qZHjLGyAH4B8DTn/Gy4VSWW8TDLHYcxdiOAY5zz9eLFEqtyhfei5pjhy7G2A/Al5/wqABfge9SXE/XHLJRbD4SvGKYWgDIA+kms6qbzrETrMUZ87G4I+JkA6ope1wGQZVNaIsYYKwFfsJ/COf9VWHyUMVZTeL8mgGPCcrljj6bvpAuAmxlj+wD8DF+xzicAKjLG/HMui9NfdGzC+xUAnER0HXMmgEzO+Vrh9Uz4bgBuPs89AezlnGdzzvMA/ArgGrj7PPsZdV4zhb+Dl6vmhoC/DkBjobY/Hr4Knrk2p0kXocb9awDpnPOPRG/NBeCvqR8GX9m+f/lQoba/E4AzwiPjIgC9GWOVhJxVb2GZ43DOR3PO63DOk+A7d39yzocAWAbgDmG14GP2fxd3COtzYfkgoXVHAwCN4avgchzO+REABxljTYVFPQBsg4vPM3xFOZ0YY6WF69x/zK49zyKGnFfhvXOMsU7CdzhUtC117K7gMKiSpD98LVp2A3jF7vREcBzXwveItgXAJuFff/jKLpcC2CX8X1lYnwH4XDjurQCSRdu6H0CG8G+43cem8vhvQHErnYbw/ZAzAMwAkCAsLym8zhDebyj6/CvCd7EDGlsv2HCsbQGkCOd6NnytMVx9ngG8AWA7gFQAP8DX0sZV5xnAVPjqKPLgy5E/YOR5BZAsfH+7AfwPQRX/Sv+opy0hhHiEG4p0CCGEqEABnxBCPIICPiGEeAQFfEII8QgK+IQQ4hEU8AkhxCMo4BNCiEdQwCeEEI/4/zbDSD8dnnkhAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.plot(Xn);"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAD8CAYAAAB+UHOxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzt3Xd8FVX+//HXh9AxSIdI2SAGBAUpAQsWkA7uV9eyX7Fh2XV/il913dUFe0OxrLu6VlZxddeyuuqKCiIiCDYgoPQWIEKoQbpAIMn5/XEnl3vTSW7Jvff9fDzyyMxnzsw9k4HzuXPmzIw55xARkcRTI9oVEBGR6FACEBFJUEoAIiIJSglARCRBKQGIiCQoJQARkQSlBCAikqCUAEREEpQSgIhIgqoZ7QqUpVmzZi41NTXa1RARiSnz58/f7pxrXl65ap0AUlNTycjIiHY1RERiipn9WJFy6gISEUlQSgAiIglKCUBEJEEpAYiIJCglABGRBKUEICKSoMpNAGbW1sxmmNlyM1tqZrd48fvNbKOZ/eD9DA9YZ6yZZZrZSjMbEhAf6sUyzWxMeHZJREQqoiL3AeQBf3DOLTCzZGC+mU3zlv3FOfdkYGEz6wJcCpwEHAd8bmYdvcXPAYOAbGCemU1yzi0LxY6IiMSDb9f8xObdBxhyUisa1AnvrVrlbt05txnY7E3vNbPlQOsyVjkfeNs5lwusM7NMoI+3LNM5txbAzN72yioBiIh4Rv79O29qIVnjR4T1s44qvZhZKtADmAP0BW4ys6uADHxnCTvxJYfvAlbL5kjC2FAkfmqlai0iEmdSx3xSLOacw8zC9pkVvghsZscA7wG3Ouf2AC8AHYDu+M4Q/lxYtITVXRnxop9zvZllmFlGTk5ORasnIhKz3pq7vlgsc9ywsDb+UMEEYGa18DX+bzjn3gdwzm11zuU75wqAv3OkmycbaBuwehtgUxnxIM65Cc65dOdcevPm5T7LSEQk5o19f3HQ/KqHh1EzKfyDNMvtAjJfCnoFWO6ceyognuJdHwD4FbDEm54EvGlmT+G7CJwGzMV3BpBmZu2BjfguFF8Wqh0REYl1X/zhHI5vfkzEPq8i1wD6AlcCi83sBy92JzDSzLrj68bJAn4H4Jxbambv4Lu4mweMds7lA5jZTcBUIAmY6JxbGsJ9ERGJKT/+9DO3v7vIPx/Jxh8qNgroK0ruv59cxjrjgHElxCeXtZ6ISCI554mZ/umLeraJ+OfrTmARkSgoOurniYu7RbwOSgAiIhG2ZOPuoPlbBqRRo0Z4R/yUpFq/EUxEJB5d//qRNx3Ou2sgzZPrRKUeSgAiIhF081vfs2n3QQDWPjI8Kt/8C6kLSEQkTJxzHDiU75/P3rmfSQuP3P4UzcYfdAYgIhI246es4KVZa0tcturhYRGuTXE6AxARCZPSGv83f3MqtWtGv/mNfg1EROLQRwuLPekGgMFdWnLGCc0iXJuSKQGIiIRB5rZ9ADx9aXfm3jnAH3/+8p7RqlIxugYgIhIG2/flAnBet+NIqmGsengYSTWMpChf+A2kBCAiEkIFBY4Pvt/IG3PWU6dmDX+DXx36/ItSAhARCYGDh/PZsvsg/Z6c6Y/l5hVEr0IVoAQgIlIJ2Tv3893aHZyd1ow+j0wvsczLV6VHuFZHRwlAROQolfT6xqIuO7UdA7u0jEBtKk8JQETkKJT0+sZAtw/pxOj+J0SoNlWjBCAiUgG3vfMDrRrW5fmZa0pc3iWlIbv2H4qZxh+UAEREypVf4Hh/wcag2KtX96Z5ch2aNKjNcY3qRalmVaMEICJSjoysHUHzc+8cQIuGdaNUm9CpfgNTRUSqmf0BT/RM/0XjuGj8QWcAIiIl+mHDLi547uug2B8Hd+TXvdtGqUahpwQgIlKCoo0/wE3npkWhJuGjLiARkSKytv9cLPbKqOp9U1dl6AxARKSIKyfO8U+f0aEpr4zqTb3aSVGsUXgoAYhIwsrLL+CEu6YAkDV+BAD7D+WxYceBoFi8UgIQkYQV+H7e1DGfMLBzCz5fvi2KNYosJQARSThbdh8kqYZx2zsLg+KBjf/dIzpHuloRpwQgInFty+6DnPbodF64vCdfr9nOv74r+1k+ha47s32YaxZ9SgAiEpeccwx86kvW5PhG9NzwxoISy90yII0b+nVgz4HD/sc6x3vffyElABGJS09+ttLf+JekQe0kHrrgZC7s2QaAurWSEqbhL6QEICJxoaDAsf3nXJ7+fDWbdh1gxsoc/7Kz0poxe/V2AJY9OIT6tdX0QQVuBDOztmY2w8yWm9lSM7vFizcxs2lmttr73diLm5k9Y2aZZrbIzHoGbGuUV361mY0K326JSKJ5atoq+oybzhtz1gc1/q9f24cJVx65iUuN/xEV+UvkAX9wzi0ws2RgvplNA64GpjvnxpvZGGAM8CdgGJDm/ZwKvACcamZNgPuAdMB525nknNsZ6p0SkcTy075cnp2RWSz+zu9Op0/7JkDi9OsfjXITgHNuM7DZm95rZsuB1sD5QD+v2GvATHwJ4HzgdeecA74zs0ZmluKVneac2wHgJZGhwFsh3B8RSTD7cvPo9fDn/vlYeiNXtB3Vs4DMLBXoAcwBWnrJoTBJtPCKtQY2BKyW7cVKi4uIVFr2zv3+6Rv6dVDjfxQqnADM7BjgPeBW59yesoqWEHNlxIt+zvVmlmFmGTk5OSWsIiJyxMotewG4oPtx/GnoiVGuTWypUAIws1r4Gv83nHPve+GtXtcO3u/CW+iygcAHZrcBNpURD+Kcm+CcS3fOpTdv3vxo9kVEEsz+Q3nc8vYPANx9Xpco1yb2VGQUkAGvAMudc08FLJoEFI7kGQV8GBC/yhsNdBqw2+simgoMNrPG3oihwV5MROSofZO5nS73HmlCmjaoHcXaxKaKjALqC1wJLDazH7zYncB44B0zuw5YD1ziLZsMDAcygf3ANQDOuR1m9hAwzyv3YOEFYRGRitp78DBd7/8sKPbBjWfg+64qR6Mio4C+ouT+e4ABJZR3wOhStjURmHg0FRQRCVS08Qfo0a5xFGoS+/RGMBGJCTt+PsQDHy31z796TW9ObJXMBzeeEcVaxTbdEiciMaHnQ9P8078f2JH+nVrQv1OLMtaQ8ugMQESqvWemrw6aH9EtJUo1iS86AxCRamX3gcOc8oCvn3/oSa34dOkW/7KmDWrz2rV9OKHFMdGqXlxRAhCRqEsd8wkAyx8cysMfL/PHAxt/gIy7B2q0TwgpAYhIVG3efcA//dmyLbw7P7vEcnqYW+jpGoCIRM3eg4c5/9mv/fOFd/UC/HDvIP/07Dv6R7ReiUJnACISUf/67kfu/u8SzujQlG/W/FRimdl39KdR/dr61h9mOgMQkYgpKHDc/d8lAMUa/xev8L076pJebWjbpH7E65aIdAYgIhEztchFXYC3rz+N7m0bJeQ7eaNNCUBEImLDjv3c8MaCoFjDujU57fimUaqRqAtIRMJu2rKtnPX4DP9860b1APjwpjOjVSVBZwAiEgG/fT3DP31t3/b89uz2fLV6O+2bNYhirUQJQETCqqDgyIv/bujXwf/WrkvS25a2ikSIEoCIhNXLX60F4PGLuvHr3mr0qxNdAxCRsJm1KodHJq8A4My0ZlGujRSlMwARCbmDh/O54LmvWeG9sL1fp+Yc5134lepDCUBEQu7Eez4Nmu/YMjlKNZGyqAtIRELqkcnLi8Wa6IXt1ZLOAEQkpCbM8l30vfncE6hTK4nZq3P43dnHR7lWUhIlAJEEU1DgeOHLNVx1+i9IrlsrpNv+5d++8k/fNrgTAKP7nxDSz5DQUReQSIJ5evpqnpi6kq73fxbS7a7YsofFG3cDcNmp7UK6bQkPnQGIJJinA96vW/gmLoDLT23HNX3bV+p1i4HbARh3wcmVr6BEjM4ARASAN+asZ+BTXx7VOvkFjnczNgTF1j06XK9tjBFKACIJZFH2LgBGdE0ptcyyTXvI2ZvL1j0Hy93eb16bx+3/WeSfb5FcR41/DFEXkEgCKeyjH9EthZF92pGe2pjaSTV4fmYmq7buY9LCTQx/ZnbQOpnjhlEzqQa79h/ipje/59VrelMryffdccbKHH+5Wbf3p9WxdSO3M1JlSgAiCeQ/3gvXB3RuQZ2aSf74TeemcfBwPpMWbiq2zgl3TeGGfh14YeYaAJ6atorzuqVQt9aR9W8+9wTaNdVbvGKNEoBIAtm+LxcgqPEvVLdWEv+8rg9XvjK32LLCxr9wOnAejgz5lNiiBCCSAHYfOMwpD/iGffZp36TUcmelNefVq3vTsmFdOqck037s5HK3ve7R4SGrp0SWEoBInNuy+yCnPTrdP59WzjDP/ie28E9njR/B6q17GfSXWaWW10Xf2KUEIBLH9uXmBTX+AA8f5Rj9tJbJfDv2XFok1+X79Tu5+MVv+eDGM3jpy7Vc0OO4UFZXIqzcYaBmNtHMtpnZkoDY/Wa20cx+8H6GBywba2aZZrbSzIYExId6sUwzGxP6XRGRor5YsS1oPnPcsEp9Y085th5JNYz01CZkjR9Bj3aNefHKXgw9ufThpFL9VeQM4B/As8DrReJ/cc49GRgwsy7ApcBJwHHA52bW0Vv8HDAIyAbmmdkk59yyKtRdRMqwcdcBbn7rewA+uPEMerRrHOUaSXVTbgJwzs0ys9QKbu984G3nXC6wzswygT7eskzn3FoAM3vbK6sEIBImkxdt9k+r8ZeSVOVO4JvMbJHXRVT4r6s1EHhfeLYXKy1ejJldb2YZZpaRk5NTUhERqYBx3nP5n760e5RrItVVZRPAC0AHoDuwGfizFy+pc9GVES8edG6Ccy7dOZfevHnzSlZPRAqd373E71oilRsF5JzbWjhtZn8HPvZms4G2AUXbAIW3FpYWF5EQW7f9ZwBuH6IbtKR0lUoAZpbinCvsYPwVUDhCaBLwppk9he8icBowF98ZQJqZtQc24rtQfFlVKi4iR3y/fie/ev4bHr2wK2PfX+yPN0+uE8VaSXVXbgIws7eAfkAzM8sG7gP6mVl3fN04WcDvAJxzS83sHXwXd/OA0c65fG87NwFTgSRgonNuacj3RiTObdtzkBYNgx+4lpuXz5j3fI1+YOMP8MtuGqcvpTPnSuyKrxbS09NdRkZGtKshElWvf5vFvR8Gf1/KGj8CgAOH8ul876clrvev607lzLRm4a6eVENmNt85l15eOd0JLFKN5Re4Yo0/wDeZ27ns5Tmlrrd63DD/I5tFSqMEIFKNzVpV8lDokhr/+37Zhe37crllQEc1/lIhSgAi1diLX/oeu/zRTWeS1tL3ELcT7wnu8nn2sh6cp75+qQQlAJFqbM66HQB0bXNssWUDO7dgwpXp1Kihp3FK5eg8UaSa+n79zhLjC+8bzHndUrjvlyep8Zcq0RmASDW0LzePXz3/DQDjL+watOzYerV49rKe0aiWxBmdAYhUQ58v899sz4U920SxJhLPdAYgEiUrt+ylcf1aNDumDocLCrjkxW9ZlL2bV0alc+u/fwDg+ct7UrumvqdJeCgBiETQN2u2k9YimUXZu7juNd9NjvVqJXHgcL6/TGEcYOhJrSJeR0kcSgAiYXAor4Atuw/Srml9wNenf+2r85ibtaNY2cDGvyhd5JVwUgIQCbFzn5zJWu9pnI9f1I012/fx0pdrj2obQ05qyV3Du4SjeiJ+SgAilVRQ4Ph8+VYGdm7p/6bunPM3/gB3vLeo1PXXPjKc0W8uYMqSLQB8eutZHMoroFXDusUe+CYSDkoAIkU8Onk5L81ay+WntuPCnm147NMVLN+0h725efzzuj6ceUIzzIzj75zsX6dfp+YkmdGvU+kvMfp1ehsevbAbA5/6kofOP5kaNYwXrujFRws3MbBzS+rVTorE7on46WmgIkWkjvmkzOUj+7Tj/51zPOc8MbPUMncN78y8rB18tmwrw05uxdOX9tBoHomYij4NVAlAxDN+ygr/s3eqat2jw3EOvl6z3X/GIBIpFU0A+koiCWPB+p1s2X2wxGXOuaDG/7GLgu++bd2oHs9e1qPYer87+/gSt2dm1KhhnJXWXI2/VFu6BiBxzznHZX+fw7drfwJ8384DG+Wf9uXS6+HPg9a5uFdb/rd3u2Lb6tQymUF/meWfv2PoibRtUp+Le7WhTs0a/OXz1Vx3Zvsw7YlIaKkLSOLexK/W8eDHy4Ji53RszmvX9uGcJ2bw40/7g5bd0K8Dfxp6Yqnbc87x/oKNjOiWQt1aunAr1Y/eCCYJbfSbC/hk0WaWPjCEd+dnF1v+5aoc7vjPwqDGf8KVvUht1oDUpg3K3LaZcVEvPZ9HYp8SgMSV/AJHh4DhmSfdN9U//emtZzH0r7P98+9kHEkM7Zs1YLAeuyAJRglA4srLs0u/4/bEVg3JGj8C5xztxx5JEs+M7MGQk1pGonoi1YoSgMSN/YfyaHZMnRKXBY7qMTOyxo/glre/Z3T/E+jYMjlSVRSpVpQAJCbl5Rfw2rc/csVp7di2J5ezHp8RtHzeXQPpPc43smfunQNKfLTC05cWH9YpkkiUACTm5Obl0+lu34vRHyoyuqdQ8+Q6/HDvIL7K3K7n6oiUQjeCScx5efa6Mpf/zynHAdCofm3O63ZcJKokEpN0BiAx5a2563li6spi8Vev6U3/Ti2YunQLZ3RoGoWaicQeJQCJGYfzCxj7/mL//JpHhvPZ0i2c27kFdWr6bsgaoqGcIhWmBCAxYf+hPLrce2RMf+/UxiTVMIZ1TYlirURimxKAxITAxn/qrWfTqZWGbopUlS4CS7X2ylfrgp7P36B2khp/kRApNwGY2UQz22ZmSwJiTcxsmpmt9n439uJmZs+YWaaZLTKzngHrjPLKrzazUeHZHYknzrmgYZ4juqWw9MGhUayRSHypyBnAP4Ci/+vGANOdc2nAdG8eYBiQ5v1cD7wAvoQB3AecCvQB7itMGiKlmTAr+LEOz47UjVsioVTuNQDn3CwzSy0SPh/o502/BswE/uTFX3e+Z0x/Z2aNzCzFKzvNObcDwMym4Usqb1V5DyRuPTcjE4A/DT2RerVq6MUqIiFW2YvALZ1zmwGcc5vNrIUXbw1sCCiX7cVKixdjZtfjO3ugXbviL+SQxLHnYB7gez6/iIReqC8Cl/QVzZURLx50boJzLt05l968efOQVk5ERI6obALY6nXt4P3e5sWzgbYB5doAm8qIi5Ro4YZd0a6CSNyrbAKYBBSO5BkFfBgQv8obDXQasNvrKpoKDDazxt7F38FeTKSY7J37ufnt7wH4+1XlvtVORCqp3GsAZvYWvou4zcwsG99onvHAO2Z2HbAeuMQrPhkYDmQC+4FrAJxzO8zsIWCeV+7BwgvCIoGenLqSZ72LvwCDuuhFLSLhUpFRQCNLWTSghLIOGF3KdiYCE4+qdpIwcvPyeW/+xqDGX0TCS4+CkGrhsSkrmfh18GOep956dpRqI5IYlACkWghs/Gff0Z+2TepHsTYiiUEJQKJu58+H/NOL7h9Mw7q1olgbkcShBCBRUfiAt7l3DaDPuOkA3D2isxp/kQjS00Al4nxjBXwKG39Az/YXiTAlAIm4Bz4q/iL3Ti2Tad2oXhRqI5K4lAAk4v7xTVbQ/DV9U5n6e434EYk0JQCJiIUbdvHq1+s4eDi/2LK7R3SJQo1ERBeBJew+W7qF6/85Hwju/nn5qnS27c0lqYYe8ywSDUoAEnb/99b3xWL/Hd2X7m0bRaE2IlJIXUASVtv2HiQ3r6BYXI2/SPQpAUhYnfrIkWGef77kFABSm+ouX5HqQF1AElaFQ/7n3jmAZsfUoWG9WvTvpBf9iFQHSgASFoOe+pLV2/b551s0rOuL6/HOItWGEoCEVEGB48VZa4Ia/9sGdYxijUSkNLoGICE1eclmHv90ZVBMozxFqiclAAmpm948MuTz89vOplXDulzdt30UayQipVEXkFSZc472YycHxb6/ZxCNG9TmuzuLvThORKoJnQFIlb23YGPQ/G/Pak/jBrWjVBsRqSglAKmSQ3kF/PHdhUGxO4d3jlJtRORoKAFIlUxZsrlYzExXfUViga4BSKUdPJzPLW//AEDj+rWYe9dAamrIj0jMUAKQSjvxnk/90xl3D9JTPUVijLqApFLezdjgn/7v6L5q/EVikBKAVMoni319/7cP6aQne4rEKCUAOWo5e3OZuTIHgNH9T4hybUSkspQA5Kjk5RfQe9zn0a6GiISAEoAclRPumuKffuyirlGsiYhUlRKAVNjeg4f907cN6sj/9m4XxdqISFUpAUiFXfdahn/65gFpUayJiIRClRKAmWWZ2WIz+8HMMrxYEzObZmarvd+NvbiZ2TNmlmlmi8ysZyh2QCJj/o87mLtuBwArHhoa5dqISCiE4gygv3Ouu3Mu3ZsfA0x3zqUB0715gGFAmvdzPfBCCD5bImD3gcNc9MK3/vm6tZKiWBsRCZVwdAGdD7zmTb8GXBAQf935fAc0MrOUMHy+hNAf3lnIKQ985p9//8YzolgbEQmlqiYAB3xmZvPN7Hov1tI5txnA+93Ci7cGNgSsm+3FgpjZ9WaWYWYZOTk5VayeVEVG1g7eW5Dtn591e396tmscxRqJSChV9VlAfZ1zm8ysBTDNzFaUUbakZwW4YgHnJgATANLT04stl8i5+MUj3T7JdWrSrmn9KNZGREKtSgnAObfJ+73NzD4A+gBbzSzFObfZ6+LZ5hXPBtoGrN4G2FSVz5fweWb6av/0+Au7cm7nFmWUFpFYVOkuIDNrYGbJhdPAYGAJMAkY5RUbBXzoTU8CrvJGA50G7C7sKpLqJWdvLk9NWwXAPed14dI+7WiRXDfKtRKRUKvKGUBL4APv5R81gTedc5+a2TzgHTO7DlgPXOKVnwwMBzKB/cA1VfhsCZP9h/KCHvUwsk/bMkqLSCyrdAJwzq0FTikh/hNQ7E3gzjkHjK7s50lkdLl3qn96+YNDqVdbQz5F4pXuBBa/X7/0bdC8Gn+R+KY3gsWBg4fzj+rmrA079nPBc18z4apejJwwh0P5BcXKZI0fEcoqikg1pAQQo3L25rL7wCEO5TmGPzObXr9ozHs3VOwmrbMenwEQdHdvoD9fUqxnT0TikBJADPnxp58Z9vRsCpzj4GHft/axw04EYP6PO3HO4V2UL9X6n/aXuXzR/YNpWLdWaCosItWargHEiHczNnDOEzPZfyjf3/gDPDrlyL13m3cf5MtVxe+e3nPwML3Hfc5bc9dz9hMzgpadldbMP505bpgaf5EEojOAGLB59wFu/8+icsud97ev2PHzITq2PIbTj2/KA+efDMCvX/yWnL25jH1/sb/sOR2b06NdI24d2JG8/AIKHNRM0vcBkUSi//HV1L/nrWfAn2eSm5fP6Y9+4Y9PueUsOrY8hhevKP407R0/HwJg1dZ9vPbtjxzOL+BwfgErtuwNKjf0pFa8dm0fbh3YEfA1/LVr6p+CSKIx3/D86ik9Pd1lZGSUX7Aacs5x5mMzeOyibpwZ0M1Slgc/WsaFPVsz8et1vL9gY7HlH47uyyltG/nnt+45yKqte+nethFd7/+sWPkv/nAOQ5+ezaG8I11GV5+Ryh1DO1G/tk7+ROKVmc0PeER/qdQKhMm8rJ1s3HWAK16Zw4w/9uNwfgEdWyaXWNY5R/uxkwGY+PW6UrcZ2PgDtGxYl5YNgx/R8OIVvXjhyzUs3LCLc//8pT9eO6kGq8YNq+zuiEgcUgIIk8Cbqvo/OROAN39zKmecUPxsYN32n8vd3tdjzi1zedb4Efy0L5emx9Sha5tj6Tv+SLdR60b1mHVH/wrWXEQShTp+w+Dl2WtLjF/28hz/9NY9B7ni5Tnsy80rdgcuwMf/dyaZ3jf2RvVr0bpRvXI/t+kxdQBfg39jvw4AXND9OL68vR9JNcoeHioiiUfXAELsUF4BHe+eAsB9v+zCAx8tC1rePLkOT15yCqMmzi227qqHh9Hx7ilcd2Z77jmvS5Xrcji/gFoa2SOScCp6DUCtQwh8tnQLz36xmgOH8v2Nf+P6tbimb3ueu8w3WufqM1IB3x28JTX+g7u0pHbNGmSNHxGSxh9Q4y8iZdI1gBC4/p/zAXjys1X+2OAurQAY0S2FEd1GkF/g+Mc3WaVuY8JV5SZrEZGQUgKogp9z87jhjQUlLhvjPaKhUFINo0VyHQZ1ackbc9bTqH4t5t89iO37ckmuq8MgIpGnawCVFDh0M1DKsXX5dmyx1yGIiESM7gMIg5y9ufQe9zl9UpvwiyIvSF/18DDdTSsiMUUJ4CgUvipxbtYO5mbt8MdfurKXGn8RiTlqtSqotK6yF6/oyeAuLSNcGxGRqlMCqKAF63eWGB9yUqtyn8EvIlIdqQuoHP+et54vVmxj6tKtALxweU+6HNeQzG37ODOtmRp/EYlZSgBFrMnZx8otexneNQWAP723OGj54JNakVTD+EXTBtGonohIyCR8AtiXm8cDk5by+MXdgoZ1zvxjPy6d8F1Q2R7tGumZOiISNxI+AZx831QA3p2fHRTv5z3BM9CdwztHokoiIhGRkBeBnXNs3HWA/IKjuwmuS0rDMNVIRCTyEu4MIL/A0eHO4nfwFlrzyPCg5cO7tuL5y3tFomoiIhGVMAlgTc4+DhzK57FPVxRbtuCeQTRpUNs/v/DewZzyoO8Vi2OGqttHROJTQiSAFVv2MPSvs0tdHtj4AxxbvxZZ40eEu1oiIlEV9wng7v8u5l/frS8W/2bMufzti9XcMeTEEtYSEYl/cZ0Adu0/FNT4v35tH+Zl7eDmAWnUSqrBoxd2i2LtRESiK+IJwMyGAk8DScDLzrnx4fqsjxdt9k9PvfVsOrVK5uyOzcP1cSIiMSWiw0DNLAl4DhgGdAFGmllo3n9Ygrfn+b79/21kDzq1Sg7Xx4iIxKRInwH0ATKdc2sBzOxt4HxgWZlrVULqmE/800NPbhXqzYuIxLxI3wjWGtgQMJ/txUIqa/vPQfN6ObqISHGRbhlLepBO0O24Zna9mWWYWUZOTk6lPiS1WQNSjq1Lo/q1WHz/4EptQ0Qk3kW6CygbaBsw3wbYFFjAOTcBmAC+dwJX9oP0Xl4RkbJF+gxgHpBtAq1HAAAFG0lEQVRmZu3NrDZwKTApwnUQEREifAbgnMszs5uAqfiGgU50zi2NZB1ERMQn4vcBOOcmA6U/jU1ERCJCw2NERBKUEoCISIJSAhARSVBKACIiCUoJQEQkQZlzlb7XKuzMLAf4sQqbaAZsD1F1YkWi7XOi7S9onxNFVfb5F865ch99XK0TQFWZWYZzLj3a9YikRNvnRNtf0D4nikjss7qAREQSlBKAiEiCivcEMCHaFYiCRNvnRNtf0D4nirDvc1xfAxARkdLF+xmAiIiUIi4TgJkNNbOVZpZpZmOiXZ+qMLO2ZjbDzJab2VIzu8WLNzGzaWa22vvd2IubmT3j7fsiM+sZsK1RXvnVZjYqWvtUEWaWZGbfm9nH3nx7M5vj1f3f3uPEMbM63nymtzw1YBtjvfhKMxsSnT2pGDNrZGb/MbMV3rE+PQGO8e+9f9NLzOwtM6sbb8fZzCaa2TYzWxIQC9lxNbNeZrbYW+cZMyvppVulc87F1Q++x0yvAY4HagMLgS7RrlcV9icF6OlNJwOrgC7A48AYLz4GeMybHg5Mwff2tdOAOV68CbDW+93Ym24c7f0rY79vA94EPvbm3wEu9aZfBG7wpm8EXvSmLwX+7U138Y59HaC9928iKdr7Vcb+vgb8xpuuDTSK52OM71Ww64B6Acf36ng7zsDZQE9gSUAsZMcVmAuc7q0zBRh2VPWL9h8oDH/w04GpAfNjgbHRrlcI9+9DYBCwEkjxYinASm/6JWBkQPmV3vKRwEsB8aBy1ekH35vipgPnAh97/7i3AzWLHmN875Y43Zuu6ZWzosc9sFx1+wEaeo2hFYnH8zEufD94E++4fQwMicfjDKQWSQAhOa7eshUB8aByFfmJxy6giLx4Phq8094ewBygpXNuM4D3u4VXrLT9j6W/y1+BO4ACb74psMs5l+fNB9bdv1/e8t1e+Vja3+OBHOBVr9vrZTNrQBwfY+fcRuBJYD2wGd9xm098H+dCoTqurb3povEKi8cEUO6L52ORmR0DvAfc6pzbU1bREmKujHi1YmbnAducc/MDwyUUdeUsi4n99dTE103wgnOuB/Azvq6B0sT8Pnv93ufj67Y5DmgADCuhaDwd5/Ic7T5Wed/jMQGU++L5WGNmtfA1/m845973wlvNLMVbngJs8+Kl7X+s/F36Av9jZlnA2/i6gf4KNDKzwjfYBdbdv1/e8mOBHcTO/oKvrtnOuTne/H/wJYR4PcYAA4F1zrkc59xh4H3gDOL7OBcK1XHN9qaLxissHhNAXL143ruq/wqw3Dn3VMCiSUDhaIBR+K4NFMav8kYUnAbs9k4zpwKDzayx9+1rsBerVpxzY51zbZxzqfiO3RfOucuBGcDFXrGi+1v4d7jYK++8+KXe6JH2QBq+C2bVjnNuC7DBzDp5oQHAMuL0GHvWA6eZWX3v33jhPsftcQ4QkuPqLdtrZqd5f8OrArZVMdG+QBKmiy7D8Y2WWQPcFe36VHFfzsR3WrcI+MH7GY6v/3M6sNr73cQrb8Bz3r4vBtIDtnUtkOn9XBPtfavAvvfjyCig4/H9x84E3gXqePG63nymt/z4gPXv8v4OKznK0RFR2NfuQIZ3nP+Lb7RHXB9j4AFgBbAE+Ce+kTxxdZyBt/Bd4ziM7xv7daE8rkC69/dbAzxLkYEE5f3oTmARkQQVj11AIiJSAUoAIiIJSglARCRBKQGIiCQoJQARkQSlBCAikqCUAEREEpQSgIhIgvr/4HxaLFgRrm4AAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.plot(Sn);"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [],
"source": [
"BB = (lambda t: Sn[t] - (t/(N-1) * Sn[-1]))(np.arange(len(Sn)))"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYMAAAD8CAYAAACVZ8iyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJztnXd8FHX+/1/v3fSQAIEAgRBCR3qJiAVFQYoFu4K9fVFPz/bz7vDUs3Knd+rZ5bCcngXFCkoTFEUUiAHpvRNqqOnJls/vj53ZnZmdbdkyu7Pv5+ORR2Y+85mZz+zszns+70pCCDAMwzDJjcXoATAMwzDGw8KAYRiGYWHAMAzDsDBgGIZhwMKAYRiGAQsDhmEYBiwMGIZhGLAwYBiGYcDCgGEYhgGQYvQAgqV169aiuLjY6GEwDMMkDCtWrDgihMgPpm/CCIPi4mKUlZUZPQyGYZiEgYh2B9uX1UQMwzAMCwOGYRiGhQHDMAwDFgYMwzAMWBgwDMMwYGHAMAzDgIUBwzAMAxYGTIIjhMDnK8pRb3MYPRSGSWhYGDAJzY+bK/DQZ6vxr/mbjR4KwyQ0LAyYhOZPn68BAGyvqDZ4JAyT2LAwYBKW9ftP4kh1AwDXDEGPBrsDS7cfjeWwGCYhYWHAJCxHqhv9bl+06TBuercUE99ahq2HqmI0KoZJTBImUR3DaHE4nX633/Leb+7lynpbtIfDMAkNzwyYhKVlVlrQfYWI4kAYxgSwMGASkpoGO+r8uJNWVDWo1t/4cTsAoKrehvs++R1f/74vquNjmESD1URMQtLn8fl+tx+qrFetbzpQCQC444MV+HX7UcxctR+XDuoQtfExTKLBMwPGlFQ32FXrdTYHNh6oxK/sWcQwurAwYExJ62Zqe4KFCDNX7TdoNAwT/7AwYEyJQ+NodOtZnXFa5zxjBsMwCQALA8Y0fFK6B3PXHgAA2DTSIMVCaLD7d0VlmGQmIsKAiN4losNEtE7RlkdEC4hoq/S/pdRORPQKEW0jojVENDgSY2CYyV+uxV0frcSXK8u9hIHN4USDnZPZMYwvIjUzeA/AWE3bZADfCyG6A/heWgeAcQC6S3+TALwZoTEwSYLT6T9o4MEZq/HAp6tUbTsqajizKcP4ISLCQAixGMAxTfMlAN6Xlt8HcKmi/X/CxTIALYioIBLjYJIDmyby+NnL+3n12XW0VrX+5e/7UG9jNRHD+CKaNoO2QogDACD9byO1dwCwV9GvXGpjmICcrLVh1Z4TqrZeBblB7audGTgCzDAYJpkwIuiMdNp0f5VENAkuVRKKioqiOSYmQRj/+hLs1rz1W/S+UTpoZwbVDXY0z0yN1NAYJqGJpjA4REQFQogDkhrosNReDqCjol8hAF0HcCHENADTAKCkpIRf45KY79YfRL3d6SUIAFcMQTD8e+EW1brghEUM4yaawmAWgJsAPCv9n6lov4eIPgFwGoCTsjqJYXwx6YMVuu2vThzkVxhYLeRTHcRqIobxECnX0ukAlgLoSUTlRHQbXELgfCLaCuB8aR0A5gDYAWAbgLcA/CESY2CSj8FFLXDxgPaw+tETlf51JAYVtVC13TuyOwDg2zX8DsIwMhGZGQghJvrYNFKnrwBwdyTOyyQ3lw0uBADk56T77JORakWKRli0b54BAHh81nrcdEZx1MbHMIkEZy1lEpYbhnUCAORlp+H3x87Hit3H0b+wOYb+/Xt3H5cwUE+A/c0kGCZZ4XQUjClomZ2GUb3bok1uhqrdaiFUVHtqG6RaCSlWjzDYo2OQZphkhIUBE/ccq/Ff6zgQ2w5Xu5fTrBYcq/GUwDxUVa+3C8MkHSwMmLinSqd+8ZonRvvs78/TtKbRgfLjntmA1p7AMMkKCwMm7tHzAM3N8B0s5m8bAJzdI9+9rLUnMEyywr8EJu7ZcqhKtd6zbY7f/q9f6z8RbprV87VftoMrnzEMwMKASQCm/rRdtf7nsT399u/dXp2r6HJNrWNleuspczYGzILKMMkACwMm7rE71A/rHgFmBqlWtR3AqUk7YdMcr1FbFo1hkhAWBkzcoy1UEyhOINWq/lorn/0ts1Jxfu+2eHJ8H5/Hl3l85jo8OGOV7jaGMRssDJi4JzPNqloP5AGUphEGshrouSv64dfJrqD4sX3bubdrZwoy7y/djS9X7uMcRkxSwMKAiXsmnNpRtR5oZmCRtmdJQmTPMZcraYrF4hYsymMcq2mAP/Ye48A0xvywMGDiHtKUwEixBv7avnTNQMy+dzgAYO2+kwCAJduOKI7pYXqpstaSN1sVQWsMY1ZYGDBxj0NjAA4mTuzSQR3QuXU2AOD+Ua4spfec1829PS87zb1c2DJTte/eY7XYqnBnDUL2MEzCw4nqmLjHLunshxbnoXTXMS8DcSDuG9kd/ze8C7LTPV93UoQpF7bMUvUf/s9FqvVb3yvDrmcvDHXYDJNQ8DsPE/c4JG+fVyYOwqKHRiAj1RpgDzVEpBIEMs9fNQAA0GB3eG1jmGSDZwZM3CM7+2SlW9GueYb/ziEwtDgPANBgCxxnYHM4Q56RMEwiwd9uJu55/9ddAABrkLWOgyU91fX1rw9iZnDxq0siem6GiTeiPjMgol0AqgA4ANiFECVElAfgUwDFAHYBuFoIcTzaY2ESE9k1NNJFaTJSXOqmusbAwmDTwaqAfRgmkYnVzOBcIcRAIUSJtD4ZwPdCiO4AvpfWTcsHy3bjt13HjB5GwhNpNY08M3hm9kbc+G4piifPxpryE7p9Wym8jxjGjBilJroEwPvS8vsALjVoHDHhsa/X4aqpS40eRsIi6/YjPTNQRiov3lIBABj/2i+6fS/oVxDRczNMvBELYSAAfEdEK4hoktTWVghxAACk/21iMA4mQWmRlYpe7fwnp2sKlhCES05GCjYeqIz4GBgmXoiFMDhTCDEYwDgAdxPR2cHuSESTiKiMiMoqKiqiN8IocqTaf6oDJjBGevLcdlZnAMC0xTsw7uWfsWqvvhqJYRKdqP/ChBD7pf+HAXwFYCiAQ0RUAADS/8M+9p0mhCgRQpTk5+frdYl76m3swx4uNofwSksdK6wWQnqKxR34tu94nSHjYJhoE1VhQETZRJQjLwMYDWAdgFkAbpK63QRgZjTHoWTJ1iP4cNnuWJ2OiQA2hzOofERNITvNfwBbXnYaGuxc74AxP9GeGbQFsISIVgMoBTBbCDEPwLMAzieirQDOl9ZjwvXvLMejX6/zaj9cVY87PihDY4R/+PwgCZ8th6pQ22iPyrFrAriVymoiGQFOZ82Yk6gKAyHEDiHEAOmvjxBiitR+VAgxUgjRXfpviN/l4cp6XDX1VzTanbjxnVLMX38IN71bGtFzbD3EGS/D5XitDev2GWO81doqtBlUGcYsJE0E8o3vluKpbzao2ob+/Xv8tus4znj2B3fR9aURLpB+54crInq8ZMOI+sR/GNHV57YXF2yO2iyFYYwkaYTB4i0VePeXne71XUdq3MsnahvRq52riHrfDrle+wbL1kNVePvnHU0fJONFbZQN8N894O3cdsfZvoXB9ooa9P7b/GgOiWEMIWkT1Y14/kf38gPn90BdowMbDlRidO92vncKwPn/XgwAuH14l3CHx0hU10f3Lbx5ZqpXW2oKq4KY5CNpZgb+SLNa3IbBXxTVsJqK7E5q91FonQme6gYbAOCOc6IjYC06ye+CiWkQgg3JjLlgYQBg7/Fadxrj5TvDt2VXSW+zJ+tsYR8r2ZE/y2FdWkXl+K2beeccStFEJndokenVJ9K2JYYxmqQQBoHe4qob7Ph5a/Azgj1Ha1E8eTb+89N23e02aUZwtKZR1b7/BAcsBYvd4cTdH63ESwu3AgBydIrTRALSmRlo23760wivPidqWdAz5iIphEEgh5TCFpnYLHkTdWiRicVbKrBXSpusxx2Sh9A/5m7S3W6XqrEcrVYLg+O1jXrdGR2GPLMQs9cewE9SArlQq5uFSnErV+lLeRZw3WlF+O8tpwLQT5A31ceLAMMkKklhQLY7/evu7QppUdtox41SrMH3/+8cdM1v5tV/RM98v0nLZPVQmSZt9eaDVejTvnnQ405GhBDYerjaS8XWNjdyFc60lD4yEtlpKVix+zh6SgnxplzWz72diPDwuF4q4b+m/GTUxsMwRpAUMwP5TR0Aerb1zn65QfFgP66Y/o984SdU1XurA7QZNIUQKlXUxa+5qmK9sGALAOC8Xq6krA/OWN2U4ScVb/28A6Mlrywl+TnpUTtnm5wMZKen4Owe+T6Fzh3n+HY3ZRgzkBQzA+VbvN4s4cfNvjOi2hzeOqb7PlmlWh/5wk/YoYhb0NKjbQ5+2KSbi4/RMG1xYsRpTDi1o9FDYJiIkhQzg30Kw22zEA2RgSqU7T1W61MQDO3sKsoyoJBVQ8FypDox7Cqf/LbX6CEwTERJCmGgTBb3xvVDQtpX62aoZe9xfUNzVb0NpZKbaptcj4qjePLskM5vdq7+z1Jc/oZ+dTGGYWJHUggDZSbSDi0yMbBjC91+n04a5tUWqNSiXuKyi/oXqPTeQzrlBTvUpKN05zGs3JMYBWNm33uWar2mIfFzFM1bdwAn2MstLhn70uKYvjwmhTCQi5mfWtwSADD9/1wP/acv6aPupxOAFAi9oivfrjmAAyfrAQBtpVlBuyh6wzCxoU/75rhvZHf3+h8+WmngaMLnUGU97vxwJe76MLGvw6xsOlgV0/MlhTBoKQmDB0b1AABkplmx69kLccPpxap+hS2zvPa16xiQlXxWVu53e6bkH59iUKUuJrI8cH4P93JpBKLVjUSeMe/xE1PDxJ4Vu49hs0IQOGKUuTcpvInkiOC0FP+yL00nJ02gGIVAWTXlYCm9YzOJTaDvRqJgluswC1e8uVS1vvdYLYpbZ0f9vEnxhJKFgV4CsrF9PFlKLTr2Adm19N8LtuCOD8q8tgcqmyinNlAKoh82HQpi1MmFXo0AnUwRcYXdgFoL0eBQZYPRQ2D8EMhuGSkMEwZENJaINhPRNiKaHK3zHDhZh1vfcz3E9YTBFUMK/e4vvzW9/P1WzF/veoinWS0Y0sllfwjkYijfxiyF0Pjz52uDGnsysfGAt35UrjERb9wlFb9J9MSlsVI/MOHhjNEXzRBhQERWAK8DGAegN4CJRNQ7Gud65fut7uU0nTz1WqGrrXl78KT6relwZT0aHU6kB1A5ycjRzcrI5vooF2xJRLQuvGd1a42Zd59p0GiSg0e+5peSeENPa6BNeBktjJoZDAWwTaqR3AjgEwCXRONEi7d4spFaLd6Xu0RTv+Cxi3pjx98vQEFzl/fPc/PUyeiml7pmApk6idMu7F+AlY+drzuOnYrANKWrK+OiUVP74cPbT0NaigVrnhiNskdHGTQqfWL1phZtftnGabjjjZcXbvVqu/yNX2NybjKiSAcRXQlgrBDidmn9BgCnCSHu8bVPTuccMeTx0ALGAGCZIu98/8IWKnUN4Howr9xzHO1yM1RGmkaHEyt3HwfgyqUvH8dqITicAs3SU1Ct8TNvlZ2G7m1zVOdskZWGXu3UbfIxGc/9OaUgF80zU93r8fz5lO065rYXxPM4A6H8TibydZgJ7XNCZmjnPN1CTIH46ZafVgghSoLpa9TMQO+qvKQSEU0iojIiKrPZws8fr+fRk5ZiweCill7Wel8fu6xn1ctZpCdWY2T7SUiUb9jKdxI9F994IkXne7TveJ1pDMpM/NEUQRAqRrmWlgNQZvoqBLBf20kIMQ3ANAAoKSkRP978Y8gnUkbwLbntwqD3O1bTiMFPLwAA/HjzhV6RgMU5Wdh1VO2fPa5HO7x5/RBV3wt7FeD1awfjH3M24j+KJGw/3hz8WMyEEMLtYXWyzoYBT34HAHj27CE4r1cb9Hh0Lu4f0AP3nNfd32EMZdvhKox6cTG65Gfjh5tH4Kvfy/HAp6shGjLx4+TzjB5e0Ci/p8n6fYw3fEUcN/X+0C3BCxGjZga/AehORJ2JKA3ABACzDBqLLkL11ur9xqdXSF2W3soyieP6ulxX1+/3Xf8gWTha3YDOD8/B9NI9AOBWwwEutZw/F+B4olubHORlp2FHhcsO9MCnrtTksp2JYZrKAJ1UOaNOaROTcxvyqxNC2AHcA2A+gI0AZggh1kfjXE0N9lJW1io/7l2usrdOkRp5JjeoyHVD5943HBf1bw+Aq5wBnkhXWRjc/j9P3MbHy3e701cHCg6MB47peHjkZCRWDGenVvGtjktG9F4828QolY1hvzohxBwhRA8hRFchxJRonaf0kZFN2i87PQU3DOsEABj+z0Ve268q8Y5PkAuwPHtFf0y7YQhOKfD4yVfqFMlJNmT10Jryk6hrdKj83JftOIaXJTfgRFK9v/2zR/W3yE9djHhEWchn5Z7jfnoysULPn0fr7h4t4v8VLExaZKVh7n3D8fOfzw1536I8329Og4taqtZvObMYfx7TC4CrZsJoRWQzAIzoEZupXjyj1F4+Mcv3RHDO2gPRH0yEeGb2RqOH0GQaFPEutQ0c+xIPCI0bSufW2bqld6OB6YUB4HJb7Ojnwe6LgOmrFZsnnFqETD+pKZ4c3wf/vLK/e31vEiYHU37Ndx31XRmO8zjFhnqbJ7YjM40/83hAmybq8YujEourC38D/CCnp544VL/EobKecqBAJIuFcHWJ5zjaIKtkwK64Zpuf69fGgjDRod7uQOtmLtVmEn4d4xI5Y8HbN5bgzesGY0TP2GkUEsviFWNkf/LaRgdSreQVW6BMSRFqnhdnIinGI4RWALbJScfhqgYQqXWlZ3RrHeORJSf1NgdyM1NwpLqBo+LjjFG928b8nDwz8IOcQ2jmqv3ISPF+W33t2sHu5VADufWC1syO9poHdGyBXu1yvNRC1yRAsflFD40weghh02B3IifD5SLtb6bGRB8hhOGV81gY+GGfwqVUz92xY14W+kvF7h0hSoNkVBPZFG+fRAS7w4m0FIuqRjUAZOnkfYo3Oscgv3y0qbc5kCu5wx6qrDd4NMnNGz9uR5/H5xs6BhYGfkhP9Xw8SmFwQT+Pp9C/rhyAsX3aoXdBcOmWn760L4DkTFb3sRRfAACbD1bB5hBe2UrP69VGt64EE1mEEKi3Od2xEZO/XKuy6TCx5ZvVXgkYYg7bDPyQYtEXBiN7efR5PdvlYOoNwSfQ69XOZXRusCefK98Pmw67l6sb7Gh0OL2ijd++MaicWkyYVEkqiS2Hqt1tNoeAjjaUSRJ4ZuCHWQppvftordu7KJzKQ7LRORlnBlpKdx7Dck0d4USdFcgV8xLl7Xr2Glcsx7bDHmHw2Qr/hZoYc8PCwA+tstNU65Ew+sozDK2e3OwcNrlOupmkblHagsqP16Js1zFfuxhKvw4uW5dctQ0A/jZzPeavP2jUkJKWFbuPY9NB70p/sYaFgR9evHqgbns4qYrTpXm48o0sGThcZb46u89d0c+9LOveaxocbs+cs55bhCunLtXd12jkMXZvo45u5YSKseeKN2NTvCYQLAz8UKRI5KXUXoQTIGuVwpZfXLCl6QdJQNT1CjJ1+7yVYPaCZumezLU/bXHlJTp1ykL0emyeUUMKGnmWq1fsiUlOWBgEyZ3neKbT/IMJnVRF/ekzuupX1erTPjiPrHjhpy0eg/i5ikhRbQDilNkb4s6PX/4OZ6WpfUjMUtKTCR0WBgF4aHQPr7avft/X5OO1be4K/+9iAj/1UFA+IH1VbYr3OgZarhjsyVwrx5vo8dbPOzFvXXzp4mXhpJ0ZTFMUYGIih83hxLx1B3RTVGs5s5sxJUgT69dnAHLaZeUtDOehlZ5iRU56Cg5V1qO20diIw1iiNL5/8pvHa+Xe87q5lxOhjoESpRoxXeOTqf3R19viy5VYdmDwl1yRiRy3v1+GOz9c6fUiqed99sZ1odd6jwSJ9eszgK+lm7d0u6dQdUaYEbJVDXbUNDrw1DcbwjpOIrG2/IRu+3mneGI2EjlbqTJAEfBWFcVb+hF5ZpCeYsHmZ8YaPBrzI9uU5qxVzxDrdVTOelUUY0Hi/vpixFbJ62d7RTUyJSHgS+cdKvtPmtvdUsljM131C27XFOpQ2gkSbWYg69uJ4P5uyGg9zuJtFijbDFKtFq9ZTTIGRMaKXE01vLrG+Pmso/brI6IniGgfEa2S/i5QbHuYiLYR0WYiGhOtMUSCe0e6CrNnplqx/skx+M8NQ3DzGcURObY1MeOrwuKC/gWqdaXKLZxgPiNonpmK3x4ZhW1TLnAHJMpoDcbxVgRHHp+eAJ5RVh7r4SQt8fSSEO1XsX8LIQZKf3MAgIh6A5gAoA+AsQDeIKK4VVw+MKo7rinpiPn3nw2LhTCmTzu3HSFcEq1MYiRIsRBevHqAqu3pS/uiR9vYVHOKNPk56bBayMuOZI8ztZCWGumNVM/+1RBn9g0z8aXGZvDx8j0+esYeI+bllwD4RAjRIITYCWAbgKEGjCMoiAjPXdkfLTXRyNFk15EanKw1R83kcS//jDs/WOFetxC5C6rI3DCsE7574JxYDy2iKPNYAUBVvfqN7+IB7WM5nIA8/a3LXiULg7+M7eXelmgztETj121H3MvtW+jH3BhBtIXBPUS0hojeJSK5aHAHAMokKOVSW1JS02BH8eTZWL//JADgtR+2YsTzP2LAU98ZPLLIsPFAJeYpUhwcPFmPlCTQj/3lizWq9fYtMnz0NJYMyfB95zld8IyUUffJBHdsOFlnwwfLdgflxhkrlO6iyihvXwGYRhCWMCCihUS0TufvEgBvAugKYCCAAwBekHfTOZTuXSOiSURURkRlFRXmVKnI6QoufGUJAOD578wRmVxZb0Px5Nle7SXFLb3eos2AtpD50h0u77NubZohI9USd2qjG4Z1QvPMVLfxmIjQt4PvWIlE4slZ6/HY1+u8kiAaSWaqx3Cs/K58prDP/GlMTyx7eGRMx6UkrF+lEGKUEKKvzt9MIcQhIYRDCOEE8BY8qqByAMpSVoUAdJN5CyGmCSFKhBAl+fn54Qw1btl4wJy5YA5X6uciapGVZsqZga9Auk55WchMtcZlBHKGxh1WWVvi7Z8TN/isUlLRbT1cHTflZX/e6nmZ3XmkBoAreaM8a75kYHvcdU5XtGtu3Awymt5ESreRywCsk5ZnAZhAROlE1BlAdwCl0RoHYwx6aue+HVxupKkmnBm0zHLZlFpkqX3E01MtSLFaUNMQvlH22zX7MfalxRFRfzRKVeaUKI3J8eb9FApyzMdjX6/Da4u2AQDunf47vv59H1btPRHzdDJ2h1OVpXh6qUtLrnQtf2h0T8PTt0fzV/lPIlpLRGsAnAvgAQAQQqwHMAPABgDzANwthGD3BSROLvxg0Hshkx82EXLGiivaNc/AvPuHY9Qp6kLmVosFFVUN+GJledhuhH+c/js2SRXiwmHZjqM4WtPoFeRnFsOx8ioWbXblj5q1ej/u/3QVLn39F/R4dG7EzvXukp0onjwbVfW+HT7kgDMlM1ftw+3vl7nX4yEVS9QqnQkhbvCzbQqAKdE6dyLSulm6bl3kmav2Ic1qwbh+BTp7xS96b1/yFz6O7HoRpVe7XGw9pM5Lv3L3cfdydb3dKzFcKMifWzjJ5E7W2TBh2jIAwCmaUq3aEqSJSqXCk6uyzhZVQ/JTkldWvye+w2+PjEJ+TrpXn2qdQvdvLNqOI9UeVWo8BFwaP4Ik5PVrB3u1dW6dhQab+gH6zer9uO+TVbjro5WxGlrYOJ0CNodTV7CZ5WHjjwdH91StTx7ncdmM1CNJm+oiFJS2C+0DSO+eJRpHqxtw8GSde/1knQ1Hqhsjeo7iybN1nSNkj0AtORneLwBaG5I2aNEIWBgYwIX9C7zC0ru1yUGNRo3wx+m/x3JYEeHRmevQ/ZG5uioROUVDxzyXO93fLuod07HFgmJF8rpUK6niC8qP10bkHLIBsil8UuoJckrXqCaKW6kz6caTa2awDHlmoaqu85HqRrwu2Q2ija+PSy/X0Hm92qjWw813FglYGBjE/w3volq3OZy4OsiqWDaHE6v36id+MxIhhDuictL/Vnhtl20iLbLSsP3vF+BWTZ4iM9BJ8UDVPhyueHMpahvtOF7TiBveWY7n5m1q0jkuenUJ9h4LXbBUVDWoXJe1MwPtejgV/eKJD5btjspx1+1TzwR8zdj0bDxzNSnN48FmYPwIkpRrTytSrdscTpV3gValIj9Idx6pQfdH5uKS13/BL4pIxnhge4XnjUypJ5VdGJVpPMxirPSHnm6/99/mY9DTC/Dz1iN488ftTT72yj3HA3fScOqUhap1bWwEALw6cZB7Od5iI5pKOGo1f3y34ZD6PD6mBnpuxftO1On0NBYWBgaRqngLI/L+wmjfym6VPA/eUEx5N8RZvVpfKuee7VyGSqsZ3Yj8EM1HaahGZD2Vzy/bjnq1KdVaZrAhRJMNGhuBr5gG+bd997lddbd3UqgWjYSFgUHkpCsiEoV3nnMtiyX3tM9WeCIWp8yJL19wXy/7tdIswYThBbo8cbHLFhJNlfsDn67Gf3/ZGXT/g5Whp0s3k6tzNFi48bBq3dfMoNHuar+gXwE2POWdpHnq9cYUs9GSJD/P+EMv8+mYPm3Rs21O0Mdoo+PGZiR615STkYKJQ10qMV9Rumbjt11qFc77t4afh3GLxmUVAD4KIeNlKIJp4lBXgoB4K8gTDxyt1o+sB4B7PtZ3+JCr3KVZLbqpWPTcUY2AhUEcMX/9Iew+VoMpl/UNqv/4OMuEqdXNtm6WjrVPjHGH2KfEgZEsFozoqU6dkhoB+4iejnnb4Wqdnvro6c0LfKQ+GFzkyilpRAqN2ka7rl9+MGhVYbPvPSsSQ1JR24RiNB9LHlyNDqeXLbBdboZXFl+jSI5fZ5zyysRBeGXiIGQr6tDW25zYc1TfU2TzQfXb4a/bj+K+T+LH/dTuVD885KAaeUaQDHEGAHBOD7UwiIQQbBVmCnW96mX3SYWbtMieLUYIg6FTvkffx+c3aV/tTKZVtu+H7InaRhzSqM6e+mZDQKeMQMFh7yzZ6WU7OKtbawBAYcssr5QTbXPjQxAALAwMZfyA9hg/oL270IjMCR+1DK6c+qtqfcOBSsxctR8VVb4GwxLwAAAgAElEQVSnrrHk5636PyT5B6asI21mcjV+5YFkwecryjGjbK/fPvW28B7MDToR4edqfN1lPMIgumqiZTuOoutf56hSOTR1VgCoDd5/GtNT90H7pzGuoMCBTy3AaX//XrXt3V924rq3l+see/GWCtTbHAE9k57+dgNO+ds81XXIn6deje/V5fqBakbAwiAOKW6drduuLZgis2K38al6nU6BZ+fq+83LLqd1SVJBKyPVij+P7Ykv7joDgG8vK5mHPluNP3++xm+fcD87rTC4pqQj2ubqq4nkt99oJ3SbMG0ZHE6BWau9kxY3JY7CphjvH0Z09bJhlXRq6VWrWv7t+Auw23KoCje+W4onv1kflJtqg92JkS/86F6XDfF62XpP7xKZeuqRgIVBHHJ1SaFqvUfbZl56aCV3frjS8MLa+09667Tlt7BJZ7sC7AYVtYjpmIzkDyO6YUgnl+5dqz5TskDjq+6Leh1hEIraTWtfyEr3HfHqFgYxUhM98tU6r7bDTZjtKtVasiBQpoL4/K4zvCJ9r3hzqbSv74d8ZZ1r5jK9dG/QLr2HFCncbZIA0btf0ycNC+p4sYCFQZxx3WlFqhgEAMhKS9F9GCgxOohFT5d6pqQrlaNyz+5uzpoUgfAXvPV//yvzuU2JfP/n3Dvc3ebLlVFLVb3Na+Zx5ZBCH7096oxYp3pW8sr3W0PeR5793HxGsbttxh2nAwCmXu/KB6at4eDaz6FKn6LV+VcpVD6hxK+t3nsCe4/Vwi4ZjiNVOz1aRC1rKRM8X9x1uvsN5alL+qrecM7tmQ+7U3ipiFo3S1Ml4Kr0k0I3Fui9/MpueJ1bZ+OnP41Ax5bxEVwTa0qKW2JgxxZ45tK+OF7biOJW2Rj+z0Ve/YQQPh8YstokN1MdnxIMSg+Yj28/DRYLoU9731XNZFXIlkNVOL2rMWqMJqmJpN/NwI6eGegpBbnY9eyF7nW9HEA9H52HAYWez6PO5kC2Ig7oVYVg0qqJvv3jWbjo1SW647nk9V8AuGbGiVDQiWcGccCQTnnuZauFVHlK/nvLUGSkWlUzgzO7tcIFmpTWesapWHJAR02Uk+ExpHZqlW148Q6jyEpLwdd3n4m+HZpjePd8dMzTF4rKBGta5JxCGalWnN0jtBmW6gFGwLAAeuoCqV7zxyHEMUQC5Xe8Kb73sqrHn8fPYR/Bd5kKjz6t++jKPZ48YFo1Ud8Ozd1xGUqU7qLTFu8I2wEgFrAwiEO0eXsyUq0qA2D/whZeAVw2hxMfLd/t5S4XK6aXej84hnbO0+nJAMDo3m292oLJaZ+RasVL1wxEWool6AemcqYZjG1JNrJu1gl0ixRaN2kAeP/XXe7lSwZ2CPmYslrLX9K37RX6GV/lwEjA+zMaoJhp6NkM/nF5f6+2I36C02TuOLtLwD6xhIVBApCRYlG9NdkdTnRt00zVp/x4HR75ap2qelIsiURZx2RCT13h8GNolklPsSAvOw3XDi0KWqevNI72CCLCPRaFVk7UqmsM3PDOcpXQCtYeokQ2ePurDXCTwp6gRPkZ1drUKtkR0kysc+vsiCa9i4e01UrCuutEdBURrSciJxGVaLY9TETbiGgzEY1RtI+V2rYR0eRwzm8m5t0/XGUcVKJVE9mdApcMVEcfywawtfuM8Vvu3rZZ4E6Mm+5tvD+vWasPBNxP9khpsDtxsi44O5HSm8mXO6mSWAgDbSDez1uPqLL22h1OCCFQPHk23lkSXA4mObDO3/iV6iAlyjxMlXVqYfChlAJ7/4k6rN4bud/XLWcWR+xYkSDcu74OwOUAFisbiag3gAkA+gAYC+ANIrISkRXA6wDGAegNYKLUN+np1S4Xvdvn6m47cLIOxxWBaBNOLUKW5q1CL6golrRvkWno+RONO0d4Z7AMxoNGNjDLarndRwMXulHOIIKpqBUL+5Oem6XSRmF3CPfb+t+DTMj41DeuEpT+JhW+rn/yl2vdyy98t1m17WiNaxbTYHfir1+tRbhsfmYstjwzDi2ywosqjzRh3XUhxEYhxGadTZcA+EQI0SCE2AlgG4Ch0t82IcQOIUQjgE+kvoyGGXec7p4paLMj9myX4/VmpS2ZGWsCBU0xasItZiLblcqPB3Ypfvtnz5t1MO6NsRAGgYZhczrdaiN/qpkGu8NdQW6TZIfQS70ho7w2PcMvACzf6QniDDYi+uPbTwuqHwCkp1jjouaxlmiNqAMAZXx9udTmq53RMLRzns+Zgh5GRvfuj8NCHWajVXaaKvDw/FNcBuhg3lQ3HAit7oXFQhh1SltdVVakCKR7d80MPC84367xjlKutznQ89F5OOu5Rdh22GOQ9md6UQrh0b3bBRxnsA4ZZ3RrrXJhTUQCCgMiWkhE63T+/L3R68l94afd17knEVEZEZVVVFQEGqppuX5YUcA+TcmmGCnu+mile7n0ryMNG0eicu953dzLxZNn41iNdwH3ozWNqlq6Xdu4AvmC8ccPJbupTOnOo9h6uFq3lnW4nKhtxG4fyRhl7A4n9h7zvGT8e8EWrz5KtdKoFz2a6mY6BehllMIgmNmZv2DPcX0DC5NEIuCnIYQYJYToq/M3089u5QCUc7BCAPv9tPs69zQhRIkQoiQ/PzmjVwHgwfN76rY/d0U/d8qHaPxog0VZj7lNbkZSlLSMJEWaQvTaB7w885q5yvNTaSnpm4NxbpEjcn+ZfF7QY6qUghw/Whb5WIMRz/+I+z9d5bePzSnwvEJ3r+cS+tS3G3T3Pc2PS7PSZpBqJZUgViInVXxpob4dZ80To/HatYN9nicRiZaaaBaACUSUTkSdAXQHUArgNwDdiagzEaXBZWSeFaUxmAZf+sVrTi3CVVIeo0BvWrGk9K8jQ3rwJDuj+6hjDgTU3i16MwU5z1PrZmoj5K4jNSiePFuVBl3+/jRLCz3hQKRnnHPWHvCZlVeJ3eFscgZTf3YR5bbCvCxVmU8lE99aBsD3zCs3I1X3pWfzM2NDGWpcEa5r6WVEVA7gdACziWg+AAgh1gOYAWADgHkA7hZCOIQQdgD3AJgPYCOAGVJfxg/pCmHwwW3qqllyXeGfthijRtObRrdqlo4O7F0UFG1z071qQ3+xohzdHpnrVwUkR60rU5IAwNtLdgAAnp3n8cCZttjV1pSUCP9e6K2eCYdgcw699fNOVCtSsPQu8NjPGu1O/LApuAR/evQuyMXAji3QoUVmQF//7iFUHgRcxmE9mqXHf+afsEYohPgKwFc+tk0BMEWnfQ6AOeGcN9lQuuFpo3r1yuj5y3ETaYL1AWe8WfrwechKS/F6w/xA8mvffLAKHfOyAjoH7D9Rh3a5GbBYyJ3/qU2OdzxBU9V3vR6bi01Pj2vSvlo2aSKPR53SxstbDgCy0qwY3KmlOwpa6UL7wneb8R9JwGnp4iP9u5I593nieVo18+/eOaxLHr5ZvR9WC4UVcPb13Wc0ed9YEX/+TYwXyge79s2jeZbLqNgyy2NcbErx86aiTKC35C/nxuy8ZqCgeSaaZ+qrGwBXRK3TKXDV1KV+j3PGsz/gzZ+2AwBOSIFo7ylSO8iE4jKao3iTjW5eHf1rr210qFKcKAtAle0+rrcLAOCHh0aEdHZtfQMt8rU3xbNqkWIsvmYM8QQLAxPQLD1FFZR2uDJ2lc+UAU2FSZqVNFxSrRbd1MpV9TZVTYG3b1QF+WOwoj7Ez1srcLymEXk6gUxndmuFTq28Sy7649JBsfH4bhNC2UdZJbnCjzAIFeWL1va/X+DznJPH9Qr6mF3yXbMTh1OgSEpKaI9gGotowcLABGijKv251kUaf4VbmOC59czOXm31NqfqIaK9r8qssHuP1WHQ0wswRSda1+EUaKujNvLHX0J4+EWSG0/v5HObXpr20b3b4sHze4R1zpcnDMTc+4brztD+Nd/l0XRG19ZBF66X7WVEwHNX9Ef/wuYJYUOLf6sGAwDo2yEX6/bpBw8d13hnxLKQufzm8+iFp8TsnGZEG1EOuGw/NlUqCXUfpS3JX3GjeptTVfErGDI0HmwHT9ajXfPQBEowKB+/ctDW/5bu9tHX+2HdoWUmuoUZHKfMkNosPUXXiyktxYKyR0fhn/M2YURP/drRMi9PGIT56w+ia34zdM1vhln3nBXW+GIFzwwShBl3nI7SR4IL6LLZYzclzZLcFS/qr++ixwRHqs5bqUO40jLIpGse0MEahOttjpAzZGqF063v/RbS/r6QS4HKhOLnoNc3M9Xq9bmEg/JFSi8N9Z/H9gqYmj0vO02VEjtRYGGQIGSlpeh6iOhhi6HqpkZ6iwomARrjGz19fpqVVKmVta6h/lxFe7XzuEQ22J1hp0sONaWFL5SuyJcP7hB0tTbAVYu4ePJsVVuq1RJR46yyHOhlb/wSseMmAiwMTIgthhlMZR21nosrEzyyblpJq2bpqnupraWsLXCkRPm2XG9zeKl9jELpJntOj3zfuWgUyGkf9Cqv/bbrGNJ1jO9N5cnxffDQaJcNQpkOIxmIj28IE1EM8VzgiUHE+cNHK7Fggye4Svt2/+0a3/UPGhWCoylqomhRr3ARLWyZ5Z4ZjPcRCXzvyO4YKwmDt3ViWnYdqYmomijFavEyFF/Qz1w5iHzBwsCENMbQgDy8e2sAQG4MPZjMyPu3DtVtl4OrxvVtF5KhdOOBSrdKpd7m1HVdNYI6mwPXDyvCdw+cjSGdWrqzgvZQFEe6S1Hr4Y/ndfMbH0FEXsVowmWjRiU2Z+3BiB4/XomPbwgTUWKpJurRNgfZadaYRTybFV+Rs7IRU/baCpUlW4+g3h4/M4M6mwNZaSnu8psr97hiBpR1GW5WlKZMtVr85v7Pz0lH+xaR9XIyMgOwkbAwMAFKYyEQWzVRnc3BWUojQCBVR1OF7fXvLIcQgesH6LH4T+qI8rXl4ZV8dDqFNEvxCKYzu7VW/Qe8S3P6Ky4/rEsrdMlvhn9c3i9iyREjaYNIJJLzqk3Gi1cPVK3HKs7g9UXb8PHyPbo+8kxoRPszXLc/dG+golZZ+PIPnpw6ry0KLsmcLyqkh7oyBYRTElL+jOF6LzdDi13unXLBn4lDiyIW2DWub0FEjpNosKLXBPRun4uFD56NI9WNmDBtmcodMZrIHjB6KZaZ0MhO96/GuXxweOkhFjcxq+3gIk9cwJZDoRfJUXLd28sBALNW73fbBZxCFga+98vSKWL/6R3DYHeKsMuH6qGNI1AWFTIz/EpnErq1yUGnVi69ciwjkJnIkJ5ixSsTByEvO023sl17nbfeWKc4CPd7JVdcUxpo5ejfPu2bq/peMbgQHfNc16cX0EhEUREEgDqyGwCev2pAVM4Tb/DMwETIP45YCIMDJ5PLBzsWjB/QHuMHtIcQAh9qKozpBfVNHtcLf5z+OwDgzesGo7Lehr98oV8TWemtE2uml+5R1Wy451xPdbEL+hXo1g5+4WrPA1j70P8ohOLzTUFrn0kWLSgLAxPhEQbRVxPt95MLhwkPPWNxqk5Q31hFDd5x/Vx6bl/CoKTYfwoFf1w6sD2+XrU/pGhhmU0HK/Hwl+oxhZvjSGlsjgV52cFnVk1kkkTmJQfy22MsZgbsShpb9NJV6KlJpl4/BGd2a+XV/tiFvZt87kwp/1SDPXSXywadWgjXJlDenqnXD8HAji0CdzQB4Za9vIqI1hORk4hKFO3FRFRHRKukv6mKbUOIaC0RbSOiV4ifKhFDfjjYYyAM/vvLrqifgwmOAYqH1di+7XR17Jk6RthgkeWQtsRmMDz02WrVOpG+YAvE2D7GRAGf37tt4E4mIdyZwToAlwNYrLNtuxBioPR3p6L9TQCTAHSX/hK3gnScIRu+GmOgJvpm9f6on4MJzK5nL8TMu89UtYVS0SwYrjm1Y5P33XpY7YHUFFUTAEy9YQgAxPwtPZliaMKtgbwRCF5lQEQFAHKFEEul9f8BuBTA3HDGwbhweVgQjlY34MZ3S/HuTSUcA5CEaCN2wzUe9y+MDzXJ7HvPQscmRmKHyn9vPhVlu4/F5FzxQjQNyJ2J6HcAlQAeFUL8DKADgHJFn3KpjYkQqVYLPpKyOz4+az2mXNYv6uf84q74L/adaBTlZSHFQthxpCZwZx9c2K8A3ds2C+vNXktdoyNolZPQmQb069Bcp2dwaN1Po8m5vdrg3F7+i9iYjYCvjUS0kIjW6fxd4me3AwCKhBCDADwI4GMiyoV+bkufE0cimkREZURUVlHRtKCZZEPrIx1tpv/fMK+CJUz4LP7zuVj44DnhHYSA+0f1QEHzyMUjPD5rXdB91+tEPSdLAFciEnBmIIQYFepBhRANABqk5RVEtB1AD7hmAoWKroUAfCqfhRDTAEwDgJKSkvivKB0HKFUEZ3SNvgtefo53AXYmMjTF0AoAbXJcrpBdfSS/C4fNB6uC7vvYTG/B4a8gD2MsUVETEVE+gGNCCAcRdYHLULxDCHGMiKqIaBiA5QBuBPBqNMaQrCjdDaPlYlrb6EkZzM5g0WdkiOqK07q0wse3nxawPGNTWB1CsrotOoIjWlHDTPiEJQyI6DK4Hub5AGYT0SohxBgAZwN4iojsABwA7hRCyNaYuwC8ByATLsMxG48jiPLNqzFKqaynl+51L+vphZnIsfpvo5vkFnpGjAOz9KjRSQUdaU8nJnKE6030FYCvdNq/APCFj33KAPQN57yMb5RvXk0JEgqGOsXMoLhV5FURjIfmWfGnYz9Za2vyuFhNFL+wmDYZyrQFDVGaGRyvtQEAstOs7LqahAx46juU7QrsdnlWt9ZIs1pwYT9PSmh/qaoZY+Ffssk4VuuJEo2WMOgpFdN5cHTPqByfiT8+u/N01fot//0tYOryFCuhV0EORFBl7xmjYWFgMiqqPFWhomUzkN1X5frHjPk5tTjP7aUEAFUNdjz1zXq/+zTanUizWuBUfA2Xbj8arSEyYcLCwMS8/H14lan0WLX3BB6c4co3c5yL2iQVh6vU5Sf9GbYdToFftx9F2e7juGKIx5v8YGV91MbHhAcLAyYkLn39F/cyGwOTm91Ha31u+2DpLvfy+b3bYvwA7+R5THzBwsDkRDOd9YA4yVnDGIO/SmtPfLNBtT5LSmzYIg69oxgXLAxMxssTBqrXF0ZeVSTDnkTJTdc23gnwft12BO8s2elzHw5LiV/412xyXlu0LSLHWVN+Auv2BR99ypiPSweqVT31Nu84lmvfXo6nv93g1d4s3RXSlJvJxRXjFRYGScCKCKTiHf/aL7jo1SURGA2TqNx2VhfVeijeaq9OHAQAuG9kj4iOiYkcLKZNhl6uoJZZnEyOCZ+MVPW7Yygan3N7tcHCB89B13yOWI9XWBiYjFSdTJfJVK2JiR7pKWpX0sOVDT56evjh/3nScHfTsTEw8QOriUzGqN5tcdtZnVVtDmd4Vrto5ThiEgvtzOCLlZ46VYer6jHpf2Ve+3Ti3FUJAwsDk5FqteCxi3ojJ8Mz6QtXGPR8dJ5XWzgVq5jERDkzIAIuURiU7/pwJb7bcMhrH56VJg4sDEzKr5PPcy87IuzPd0bXVphxx+mBOzKmIl0xM2jfPBMpiqSInVp51yZ+6ZqBXm1M/MI2A5OSk+EJ7rE7IisM3rtlqFfRdcb8pCvuudVCcCpeMmobvFWJ/Qt59phI8C86CXBGeGaQymkokhKlp1qKhWBXqB9TdV4OivK8ZwtM/MLCwMRMHtcLAFQ/WiXfbzyE4smzdYOH/MGlLpObO87uAquF4FCkIy1s6Z2agiPUE4uw7hYR/YuINhHRGiL6iohaKLY9TETbiGgzEY1RtI+V2rYR0eRwzs/4p0/7XACAUyEMhBDuUpW3ve/y/rjjgxWxHxyTkOx69kI8fMEpsFpIpX5s7ydPEZMYhCu6FwDoK4ToD2ALgIcBgIh6A5gAoA+AsQDeICIrEVkBvA5gHIDeACZKfZkoIHtyKGcGnR+egydmqfPQ/7Slwucx6nTq2DJMVb0dmxQF7x/7ep2Bo2EiQVjCQAjxnRBCLoi7DICcuPwSAJ8IIRqEEDsBbAMwVPrbJoTYIYRoBPCJ1JeJAlZJnSO7lu6RUg6/v3Q3ACAv2xWZfN1pRT6P8dBnq6M5RCZB2XeiDnuO1eLqqUtRq6iJzSQukVTq3QpgrrTcAcBexbZyqc1XOxMF5HoDsjB4cMYq1fbBRS0B+I8Z+GHT4SiNjjEDpbuOofff5hs9DCYCBHQtJaKFANrpbHpECDFT6vMIADuAj+TddPoL6Asfn64uRDQJwCQAKCry/fbK6GPRzAxqNSof2XZg8xOUVheicZlhmMQkoDAQQozyt52IbgJwEYCRQrh9GMsBdFR0KwSwX1r21a537mkApgFASUkJZ0IPETkoSBYGB07WqbbLH+jmg5WwOZxI1Xh/KLNSTjq7C6Yt3oFzeuRHb8BMwjB+QHt3wRotlw/uoHJaYBKDsILOiGgsgL8AOEcIoayBNwvAx0T0IoD2ALoDKIVrxtCdiDoD2AeXkfnacMbA+MZjQHY91I/X2lTbZSHx4bI9AIBnLu2n2l5R7UlEZiHCpqfHIoXTCzAAuvtJOvePy/t5JbVj4p9wbQavAcgBsICIVhHRVAAQQqwHMAPABgDzANwthHBIxuZ7AMwHsBHADKkvEwU8NgPvbZsPVqm8iD5ctgf7TtR5d5RwCoGMVCv7jjMAgIxU3w97ZZoKJnEIa2YghOjmZ9sUAFN02ucAmBPOeZngkGcGlfWuGcGF/Qowe+0BAMDEt5Z59T/z2R8w654z0V+qbfzlCk9WynCT3THmoqreptu+/skxnJwuQWERbmJkA/LDX64FAOTnpLu3Hatp1N1ne0W1+/8LC7a42+160wsmaZHdk7Vkp3O6s0SFhYGJ0eYksjuDf6DXNKh9x3u0y4nImBhz8MfzfCoFmASFhYGJERphYLMHVvXIu1gU+YeeHN8HE09l117Gw8COLQJ3YhIKFgYmprClJ2vkriM1+LRsr5/eapS56Hq2y4GF9cCMAr3XCmWJSybxYGFgYpQeHyOe/zGofWQhoExed/BkfSSHxZgAm44NqW1uhgEjYSIFCwNGhawmKj/ucTNt15x/5IwavYJJbDxObPjuMSr+OW8zVu89oWrjQDNGC7samw+eGZic3IzQ5P3Bynovt0E57oBhZFpKGW9bZqUG6MkkCiwMTM59o3qEfQyud8xoGdixBT687TQs/6vf1GVMAsG/cpOT6SdtQHoQD/m7z+0ayeEwJuKs7q35RcFE8J00OVlp/nLIENY9OQaLHhrhs8+fxvSKwqgYhok3WBiYHL3a9Z/feToAwGIhNEtPQefW2TEeFcMw8QYLgyTj39cMQJd8V/phi56kYBgmKWHXUpOj9Afv2TYHlw0qhM3hROtm6fjbxb0NHBnDMPEEzwxMTl6zNPeykJIIpFotKHt0FMYPaO/e9u0fz4r52BiGiR9YGJicET3y8d9bTkXn1tl4/OI+Pvv17dA8hqNiGCbeYDWRySEinNuzDc7t2cbooTAm5myujZ3whFsD+V8ALgbQCGA7gFuEECeIqBiuspabpa7LhBB3SvsMAfAegEy4Kp7dJ7S5lhnDuX5YEcb0aWf0MJgEYNezFxo9BCYChKsmWgCgrxCiP4AtAB5WbNsuhBgo/d2paH8TwCQA3aW/sWGOgYkwfTvk4plL+2F4d37bY5hkISxhIIT4TipyDwDLABT6609EBQByhRBLpdnA/wBcGs4YmMjzzk2nGj0EhmFiTCQNyLcCmKtY70xEvxPRT0Q0XGrrAKBc0adcamPiCM5LzzDJR0CbAREtBKCnPH5ECDFT6vMIADuAj6RtBwAUCSGOSjaCr4moDwC9KCef9gIimgSXSglFRVx2MdrMvW84lm4/avQwGIYxgIDCQAjhNy0hEd0E4CIAI2VDsBCiAUCDtLyCiLYD6AHXTECpSioEsN/PuacBmAYAJSUlbGSOMqcU5OKUglyjh8EwjAGEpSYiorEA/gJgvBCiVtGeT0RWabkLXIbiHUKIAwCqiGgYERGAGwHMDGcMDMMwTPiEG2fwGoB0AAtcz3a3C+nZAJ4iIjsAB4A7hRDHpH3ugse1dC7UdgaGYRjGAMISBkKIbj7avwDwhY9tZQD6hnNehmEYJrJwOgqGYRiGhQHDMAzDwoBhGIYBCwOGYRgGLAwYhmEYAJQoCUOJqALA7ibu3hrAkQgOJxHgazY/yXa9AF9zqHQSQgSVcTJhhEE4EFGZEKLE6HHEEr5m85Ns1wvwNUcTVhMxDMMwLAwYhmGY5BEG04wegAHwNZufZLtegK85aiSFzYBhGIbxT7LMDBiGYRg/mFoYENFYItpMRNuIaLLR4wkHIupIRIuIaCMRrSei+6T2PCJaQERbpf8tpXYioleka19DRIMVx7pJ6r9VqkcRtxCRVaqY96203pmIlktj/5SI0qT2dGl9m7S9WHGMh6X2zUQ0xpgrCR4iakFEnxPRJul+n27m+0xED0jf6XVENJ2IMsx4n4noXSI6TETrFG0Ru69ENISI1kr7vCKVCQgeIYQp/wBYAWwH0AVAGoDVAHobPa4wrqcAwGBpOQfAFgC9AfwTwGSpfTKA56TlC+BKD04AhgFYLrXnAdgh/W8pLbc0+vr8XPeDAD4G8K20PgPABGl5KoC7pOU/AJgqLU8A8Km03Fu69+kAOkvfCavR1xXgmt8HcLu0nAaghVnvM1xlb3cCyFTc35vNeJ/hSu0/GMA6RVvE7iuAUgCnS/vMBTAupPEZ/QFF8YM/HcB8xfrDAB42elwRvL6ZAM4HsBlAgdRWAGCztPwfABMV/TdL2ycC+I+iXdUvnv7gqoT3PYDzAHwrfcmPAEjR3mMA8wGcLi2nSP1Ie9+V/eLxD0Cu9HAkTbsp77MkDPZKD7cU6T6PMet9BlCsEQYRua/Stk2KdlW/YP7MrCaSv2Qy5VJbwiNNjQcBWA6grXBVkIP0v43Uzdf1JwAtFc8AAAKXSURBVNLn8hKAPwNwSuutAJwQQtildeXY3dclbT8p9U+k6wVcM9kKAP+V1GNvE1E2THqfhRD7ADwPYA9ctdNPAlgB899nmUjd1w7SsrY9aMwsDPT0ZQnvOkVEzeAqHHS/EKLSX1edNuGnPa4goosAHBZCrFA263QVAbYlxPUqSIFLlfCmEGIQgBq41Ae+SOjrlnTkl8Cl2mkPIBvAOJ2uZrvPgQj1OsO+fjMLg3IAHRXrhQD2GzSWiEBEqXAJgo+EEF9KzYeIqEDaXgDgsNTu6/oT5XM5E8B4ItoF4BO4VEUvAWhBRHKFPuXY3dclbW8O4BgS53plygGUCyGWS+ufwyUczHqfRwHYKYSoEELYAHwJ4AyY/z7LROq+lkvL2vagMbMw+A1Ad8krIQ0uY9Msg8fUZCTPgHcAbBRCvKjYNAuA7FFwE1y2BLn9RskrYRiAk9I0dD6A0UTUUnorGy21xRVCiIeFEIVCiGK47t0PQojrACwCcKXUTXu98udwpdRfSO0TJC+UzgC6w2Voi0uEEAcB7CWinlLTSAAbYNL7DJd6aBgRZUnfcfl6TX2fFUTkvkrbqohomPQ53qg4VnAYbVCJsrHmAri8brYDeMTo8YR5LWfBNe1bA2CV9HcBXPrS7wFslf7nSf0JwOvSta8FUKI41q0Atkl/txh9bUFc+wh4vIm6wPUj3wbgMwDpUnuGtL5N2t5Fsf8j0uewGSF6WBh0vQMBlEn3+mu4vEZMe58BPAlgE4B1AD6AyyPIdPcZwHS47CI2uN7kb4vkfQVQIn2G2wG8Bo0TQqA/jkBmGIZhTK0mYhiGYYKEhQHDMAzDwoBhGIZhYcAwDMOAhQHDMAwDFgYMwzAMWBgwDMMwYGHAMAzDAPj/RoPLHpqdL/gAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.plot(np.arange(N), BB)\n",
"plt.axhline(0, color='green');"
]
}
],
"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.7.3"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment