commit 966e68707b992c4ee03d1280a098938f7bdc1b39 Author: zara <912104620111@njust.edu.cn> Date: Mon Feb 21 09:46:08 2022 +0800 init with success diff --git a/.idea/.gitignore b/.idea/.gitignore new file mode 100644 index 0000000..13566b8 --- /dev/null +++ b/.idea/.gitignore @@ -0,0 +1,8 @@ +# Default ignored files +/shelf/ +/workspace.xml +# Editor-based HTTP Client requests +/httpRequests/ +# Datasource local storage ignored files +/dataSources/ +/dataSources.local.xml diff --git a/.idea/inspectionProfiles/profiles_settings.xml b/.idea/inspectionProfiles/profiles_settings.xml new file mode 100644 index 0000000..105ce2d --- /dev/null +++ b/.idea/inspectionProfiles/profiles_settings.xml @@ -0,0 +1,6 @@ + + + + \ No newline at end of file diff --git a/.idea/linearLoss.iml b/.idea/linearLoss.iml new file mode 100644 index 0000000..d0876a7 --- /dev/null +++ b/.idea/linearLoss.iml @@ -0,0 +1,8 @@ + + + + + + + + \ No newline at end of file diff --git a/.idea/misc.xml b/.idea/misc.xml new file mode 100644 index 0000000..d56657a --- /dev/null +++ b/.idea/misc.xml @@ -0,0 +1,4 @@ + + + + \ No newline at end of file diff --git a/.idea/modules.xml b/.idea/modules.xml new file mode 100644 index 0000000..fc6539e --- /dev/null +++ b/.idea/modules.xml @@ -0,0 +1,8 @@ + + + + + + + + \ No newline at end of file diff --git a/main.py b/main.py new file mode 100644 index 0000000..5596b44 --- /dev/null +++ b/main.py @@ -0,0 +1,16 @@ +# This is a sample Python script. + +# Press Shift+F10 to execute it or replace it with your code. +# Press Double Shift to search everywhere for classes, files, tool windows, actions, and settings. + + +def print_hi(name): + # Use a breakpoint in the code line below to debug your script. + print(f'Hi, {name}') # Press Ctrl+F8 to toggle the breakpoint. + + +# Press the green button in the gutter to run the script. +if __name__ == '__main__': + print_hi('PyCharm') + +# See PyCharm help at https://www.jetbrains.com/help/pycharm/ diff --git a/notebook.ipynb b/notebook.ipynb new file mode 100644 index 0000000..19e5e52 --- /dev/null +++ b/notebook.ipynb @@ -0,0 +1,367 @@ +{ + "cells": [ + { + "cell_type": "code", + "execution_count": 1, + "metadata": { + "collapsed": true + }, + "outputs": [], + "source": [ + "import numpy as np\n", + "from matplotlib import pyplot as plt" + ] + }, + { + "cell_type": "code", + "execution_count": 104, + "outputs": [], + "source": [ + "file = open(\"tempwithbright.txt\")\n", + "data = file.read()" + ], + "metadata": { + "collapsed": false, + "pycharm": { + "name": "#%%\n" + } + } + }, + { + "cell_type": "code", + "execution_count": 105, + "outputs": [], + "source": [ + "x = []\n", + "y = []\n", + "for line in data.split(\"\\n\"):\n", + " x.append(np.float(line.split(\",\")[0]))\n", + " y.append(np.float(line.split(\",\")[1]))" + ], + "metadata": { + "collapsed": false, + "pycharm": { + "name": "#%%\n" + } + } + }, + { + "cell_type": "code", + "execution_count": 106, + "outputs": [], + "source": [ + "x = np.array(x[0:150])\n", + "y = np.array(y[0:150])" + ], + "metadata": { + "collapsed": false, + "pycharm": { + "name": "#%%\n" + } + } + }, + { + "cell_type": "code", + "execution_count": 107, + "outputs": [], + "source": [ + "xy = np.concatenate(([x],[y]),0)" + ], + "metadata": { + "collapsed": false, + "pycharm": { + "name": "#%%\n" + } + } + }, + { + "cell_type": "code", + "execution_count": 108, + "outputs": [ + { + "data": { + "text/plain": "[]" + }, + "execution_count": 108, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "text/plain": "
", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD4CAYAAAAXUaZHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAkbElEQVR4nO3deXhV5b328e8vI4QAARLmGRkEmYNTnVBbsXWq9nBAa8VSkdaeHluPWq2v1fY9p7a2ttW2KirFATkiirOtaC1qW9Awz/OUkEAgkIHMye/9IxvfiImEZIe1s3N/ritX9n72kDsh62bl2Wuvx9wdERGJLjFBBxARkfBTuYuIRCGVu4hIFFK5i4hEIZW7iEgUigs6AEBqaqr3798/6BgiIi3KsmXLDrh7Wl23RUS59+/fn4yMjKBjiIi0KGa2q77bNC0jIhKFVO4iIlFI5S4iEoVU7iIiUUjlLiIShVTuIiJRSOUuIhKFIuI4dxGRaFJV7RwsKmN/YRlmkJwYR2yM4Q75JRXkFpWRW1jz0a9LEpeN6hn2DCp3EZF6VFc7ZmBmn17fcfAIazLz2ZNXTG5RGT1T2jIwtR0bsgv557YD7M4rZn9hGVXVDVsr4/LRPVXuIiJNVVZZRcbOQ2zLLeJAUTmD0trRKSmBxZtzydiZR0WVU1ldzYGicvKOlAOQEBtDXKxR7U5pRfWnz9W+TRyFpZUAmMHIXh05e1Aq3Tsm0r1DG7p2aANAYWkl1dUOBh3axJPWPpGu7RNJTU6kbUJss3yfKncRiRruzs6DxSzZfpAl2w+yr6CU2BijW/s2nDmwC7lFZcz5505yC8s+99iE2BjG9UshLTSFMqF/Z7okJ4I7FdVOZVU11Q5Du7VndJ8U+nVJok18LIeLy9mWW8SA1GQ6t0sI4Luu23HL3cz6AM8A3QAHZrn7783sQeByoBzYBtzo7ofNrD+wAdgUeool7j6zOcKLSOtWXF7JovX7WJuVz+68YlbtySenoBSAtPaJDOjSjvLKKhZvzuXlFVkAnDckjV98vR8je3ekU1IC2w8Usb+gjHH9OpGceOL7uylJCYzv1zms31c4NOQ7qQRuc/flZtYeWGZmi4BFwF3uXmlmvwTuAu4MPWabu49plsQi0irlHSnn3fX7eGttNoWllSTGxbBqz2GOlFeRGBdDn85JpPfvxFmDunDmwC4MTG336Vy5u7Mtt4jYmBgGpLb7zPMO696BYd2D+I6a13HL3d2zgezQ5UIz2wD0cvd3at1tCfCN5okoIq2Bu7Nkex4vZuxhf2EZZZVVlFZUU1pRRU5B6adz2307J9G3cxJHyiu5bFRPrh7Xiwn9OxMTY/U+t5lxStf2J+tbiQgn9DdIaMplLLD0mJu+DbxQ6/oAM1sBFAD3uPuHdTzXDGAGQN++fU8khohEkU05hby9Npu31mSzeV8RKUnxDExtR2JcLF2S42gTF8tZg7rQu1NbTh/QhdG9O366Ry71a3C5m1ky8BJwq7sX1Br/CTVTN3NDQ9lAX3c/aGbjgVfMbETtxwC4+yxgFkB6enrDjhkSkahxuLic/3lrA/MzMjGD8X078atrRnHFmJ60iW+eI0hakwaVu5nFU1Psc9395Vrj04DLgIvc3QHcvQwoC11eZmbbgCGAVuMQaeXKK6t5fdVePtiSy9835VJUVsnM8wfxnXMHkJqcGHS8qNKQo2UMeArY4O4P1RqfBNwBnO/uxbXG04A8d68ys4HAYGB72JOLSIvh7vx9cy4/f3092w8cITU5kYlD07j5/EGc2qND0PGiUkP23L8EXA+sMbOVobG7gYeBRGBRaP7r6CGP5wE/M7MKoBqY6e554Q4uIpGvqKySV1Zk8dySXWzMKWRAajtmT0tn4tCumjdvZg05WuYjoK5/hbfquf9L1EzhiEgrlXeknN+9u5mXl2dRVFbJ8B4d+J+vj+Sa8b1IjNN8+smgd6iKSFitycxn5nPL2F9YyuWjenLdmf0Y1zdFe+onmcpdRMKivLKa2f/YwUOLNpPaLoEFM89mdJ+UoGO1Wip3EWmSI2WVvL02h8cXb2PL/iK+PLwbD1w9sua8LBIYlbuINIq78/Q/d/LgXzdxpLyKgWnteOqGdC46tVvQ0QSVu4g0QkFpBXe9vIY3V2dzwdA0bpl4Cun9OmlePYKo3EWkQQpKK/jbhv28tSabxZtzqaiq5o5JQ5l53qAvPK+LBEPlLiJfqLKqmp+/sZ55H++hvKqabh0SmTKhD98Y34eRvTsGHU/qoXIXkXoVl1dyy9zlvL8plykT+vBv6b0Z26eT9tRbAJW7iNRpe24R339+BRtzCvjvr5/GdWf0CzqSnACVu4h8hrvz0vIs7n11LQlxMTx1wwQmDusadCw5QSp3EflU1uES7n55DYs353L6gM78fsoYenRsG3QsaQSVu4gA8ObqbH780mqq3Ln/ihFcf2Y/za23YCp3kVZu6/5C/vj+NhauyGJs3xQenjKWPp2Tgo4lTaRyF2mlyiurufOl1SxckUViXAy3TBzErRcPIT42JuhoEgYqd5FWqKraue3FVby+ai8zzx/ETecO0LlgoozKXaQVKa+sJmNnHnOX7ubNNdncOWkY371gUNCxpBk0ZJm9PsAzQDfAgVnu/nsz6wy8APQHdgKT3f1QaFm+3wNfBYqBae6+vHnii0hD7ckrZsqsJWQdLiEhNoZbLx6sYo9iDdlzrwRuc/flZtYeWGZmi4BpwHvu/oCZ/Rj4MXAncCk166YOBs4AHg19FpGAZOeXMPWJJRSVVfLodeM4b0ga7RL1h3s0O+4rJ+6efXTP290LgQ1AL+BK4OnQ3Z4GrgpdvhJ4xmssAVLMrEe4g4tIw3ywOZfJj/+L/OIKnp1+OpeO7KFibwVO6F/YzPoDY4GlQDd3zw7dlEPNtA3UFP+eWg/LDI1l1xrDzGYAMwD69u17orlF5DiW7z7En97fxrsb9tG/SxLPTD+dUb1Tgo4lJ0mDy93MkqlZ+PpWdy+ofd5md3cz8xP5wu4+C5gFkJ6efkKPFZH6bcgu4N5X1/LJzkO0bxPH7ZcM5TvnDtDC1K1Mg8rdzOKpKfa57v5yaHifmfVw9+zQtMv+0HgW0KfWw3uHxkSkGVVXO797bwt/en8rKUnx3Hf5cP4tvY+mYFqp4865h45+eQrY4O4P1brpNeCG0OUbgFdrjX/LapwJ5NeavhGRZuDu/N83N/Dwe1u4fHRPFv3wfKZ9aYCKvRVryL/8l4DrgTVmtjI0djfwADDfzKYDu4DJodveouYwyK3UHAp5YzgDi0gNd2dNVj4FJZX8a/sBZv9jB9PO7s9PLx+u5e7k+OXu7h8B9f2mXFTH/R24pYm5ROQLlFZU8X9eWcuLyzI/HbtqTE/uvUzFLjX0N5tIC5KTX8q7G/bx3JJdbMwp5PsTT+G8IWnExxqjeqfoLI7yKZW7SAtQWlHFQ4s28+SH26l2GJDajj9P0yIaUj+Vu0iEW7brELcvWMX23CNMPb0P088ZwKC0ZE2/yBdSuYtEqPLKah7860ae/GgHPTu25bnpZ3DO4NSgY0kLoXIXiUBV1c6P5q/kjdXZXHdGX+766qkk67BGOQH6bRGJMO7O/a+v443V2dx16TBuPl9nbpQTp3IXiRBb9hXy7JJdLFq/j+z8UmacN1DFLo2mchcJ2M4DR/jNos28sXovCbExnD8kjTsnDeOK0T2DjiYtmMpdJCDV1c5zS3fxi7c2YkZoubuBdG6XEHQ0iQIqd5GTrKraWbQ+hz/9fRurM/M5f0gav7xmFN07tgk6mkQRlbvISXS4uJybn13G0h159OuSxIPfGMU3xvfWMesSdip3kZNk6/4iZjybQWZeCQ9cPZJ/S+9DrE4XIM1E5S7SjEorqnh1ZRYvZmSSsesQHdvG8+z00zljYJego0mUU7mLNIOKqmqe/HAHT3y4nbwj5QzumsztlwzlmnG9NbcuJ4XKXSTM1mTmc/uCVWzMKWTi0DRuOm8gZw3sonl1OalU7iJhtGzXIa5/aint28TxxLfS+fLwbsd/kEgzOG65m9ls4DJgv7ufFhp7ARgauksKcNjdx5hZf2ADsCl02xJ3nxnu0CKRaG1WPtP+/DFd2ycy/+az6NpB0y8SnIbsuc8B/gA8c3TA3f/96GUz+w2QX+v+29x9TJjyiUS8kvIqHlu8jccWbyM1OZG5N52pYpfANWSZvQ9Ce+SfE1o8ezJwYZhziUS84vJK5n28hyc+2E5OQSmXjerBPV8brhdMJSI0dc79XGCfu2+pNTbAzFYABcA97v5hXQ80sxnADIC+ffs2MYbIybX7YDHXPrmEzEMlnDGgM7+bMoYzdXijRJCmlvtUYF6t69lAX3c/aGbjgVfMbIS7Fxz7QHefBcwCSE9P9ybmEDlpdhw4wrVPLKGkoornbzqDswdpAQ2JPI0udzOLA64Gxh8dc/cyoCx0eZmZbQOGABlNzCkSEfKLK7juiSWUVVbz/HfOZHjPDkFHEqlTTBMeezGw0d0zjw6YWZqZxYYuDwQGA9ubFlEkMrg7d7+yhv2FZcy5cYKKXSLaccvdzOYB/wKGmlmmmU0P3TSFz07JAJwHrDazlcACYKa754Uxr0hgXlmZxZurs/nhl4cwqndK0HFEvlBDjpaZWs/4tDrGXgJeanoskcix+2Axjy7eyoJlmUzo34mZWh1JWgC9Q1XkC8zP2MM9r6wFYMqEvtx68WCdyVFaBJW7SB3KKqv42evrmbt0N2cP6sJDk8fo+HVpUVTuIsfIyS/lu3OXsWL3YW4+fyC3f2UocbFNOfZA5ORTuYvU8u76fdz50mpKK6p49LpxXDqyR9CRRBpF5S6tXk5+Kct3H2LR+n0sXJHFqT068MjUMZzStX3Q0UQaTeUurdrHO/K47sklVFQ5CbExfO+CQdx68RAS4jQNIy2byl1arcLSCn40fyU9U9ry8JSxDOvRnsS42KBjiYSFyl1arZ+/sZ69h0t4ceZZjO6TEnQckbBSuUurU1ZZxX+/uYH5GZl874JBjO/XOehIImGncpdWJe9IOdP+/DGrM/O56dwB/OjLQ4KOJNIsVO7SapRVVnHzsxlsyinksW+OZ9Jp3YOOJNJsVO7SKrg7d720hk92HuIP145VsUvU0/Fe0ir86e/beHlFFj/68hAuG9Uz6DgizU7lLlHvzdXZPPjXTVw1pif/ceEpQccROSk0LSNRa3tuES8uy2T2RzsY368TD1wzipo13UWin8pdotJji7fxwNsbiY0xJg7tygPXjKRNvN6gJK2Hyl2izsPvbeGhRZv52qge/PSy4XTtoFP1SuvTkGX2ZpvZfjNbW2vsPjPLMrOVoY+v1rrtLjPbamabzOyS5goucqxDR8q5c8FqHlq0mWvG9ebhKWNV7NJqNWTPfQ7wB+CZY8Z/6+6/rj1gZsOpWVt1BNATeNfMhrh7VRiyitTrbxv3cdv8VRSUVvLdCwZx+1eGEqMVk6QVa8gaqh+YWf8GPt+VwP+6exmww8y2AqdTs8C2SLNYv7eAW+auYEBqO+bNGM2w7h2CjiQSuKYcCvl9M1sdmrbpFBrrBeypdZ/M0NjnmNkMM8sws4zc3NwmxJDWLO9IOTOezaBD2zjm3DhBxS4S0thyfxQYBIwBsoHfnOgTuPssd0939/S0tLRGxpDWrKC0gm/P+YT9hWU8fn265tdFamlUubv7Pnevcvdq4Alqpl4AsoA+te7aOzQmElb5xRVc/+RS1u3N55GpYxmjU/aKfEajyt3Mai8s+XXg6JE0rwFTzCzRzAYAg4GPmxZR5LOy80v491n/YkN2IY9eN55LRug8MSLHOu4LqmY2D7gASDWzTOCnwAVmNgZwYCdwM4C7rzOz+cB6oBK4RUfKSDhtyinkhtkfU1RWyZM3pHPeEE3pidTF3D3oDKSnp3tGRkbQMSTC7Tp4hGse/SexMcafp53O8J568VRaNzNb5u7pdd2md6hKi3CwqIwbZn9MZbXzws1nMSgtOehIIhFNZ4WUiJdfXMG0P39Cdn4pT90wQcUu0gAqd4lo+SUVXD976aerJ43v1+n4DxIRlbtErtKKKm7888dszC7ksevHMXFY16AjibQYmnOXiOTu3LFgNct3H+bR68Zx4bBuQUcSaVG05y4R6XfvbuG1VXu5Y9JQLh3Z4/gPEJHPULlLxJn90Q5+/94WvjG+N989f1DQcURaJJW7RJT5GXv42RvruWRENx64eqSWxRNpJJW7RIxVew7zk4VrOHdwKg9PHUtcrH49RRpLW49EhPySCr4/bzld27fhkaljSYzTeqciTaGjZSRw1dXOnQtWk324lPkzzyIlKSHoSCItnvbcJXC/fXczf1mXw48vHca4vnqTkkg4qNwlUC8vz+SRv21lyoQ+TD9nQNBxRKKGyl0C8+rKLG5fsJqzBnbhZ1eepiNjRMJI5S6BmP/JHm59YSXp/TrxxA3pJMTpV1EknPSCqpx0qzMPc9fCNZw7OI3Hvzmetgk6MkYk3I67u2Rms81sv5mtrTX2oJltNLPVZrbQzFJC4/3NrMTMVoY+HmvG7NIClVVWcdv8VaQlJ/LI1LEqdpFm0pC/hecAk44ZWwSc5u6jgM3AXbVu2+buY0IfM8MTU6LF797dwpb9RfzimpF0bBsfdByRqHXccnf3D4C8Y8becffK0NUlQO9myCZRZun2gzy+eBuT03szcahO3yvSnMLxKta3gbdrXR9gZivMbLGZnVvfg8xshpllmFlGbm5uGGJIJDtcXM6tL6ykb+ck7r18RNBxRKJek8rdzH4CVAJzQ0PZQF93Hwv8CHjezOpcxdjdZ7l7urunp6VpBftodvTc7AeKynhk6jiSE/U6vkhza3S5m9k04DLgOnd3AHcvc/eDocvLgG3AkDDklBbsuaW7eWf9Pu6cNIyRvTsGHUekVWhUuZvZJOAO4Ap3L641nmZmsaHLA4HBwPZwBJWWaWNOAT9/Yz3nD0nj21/SO1BFTpbj/n1sZvOAC4BUM8sEfkrN0TGJwKLQuwqXhI6MOQ/4mZlVANXATHfPq/OJJeodOlLO959fQYc28fxm8mhiYvQOVJGT5bjl7u5T6xh+qp77vgS81NRQ0vLtLyzlm08uZXdeMXNunEBqcmLQkURaFb2yJWFXVFbJlMeXkFNQypxpEzh7UGrQkURaHZW7hN2f3t/K9gNHeP47Z3D2KSp2kSDobE0SVrsPFvPkhzu4emwvFbtIgFTuElb//dZ64mKNOy8dFnQUkVZN5S5hMz9jD39dt49bJp5Ctw5tgo4j0qqp3CUslu3K456FaznnlFRuPm9g0HFEWj2VuzRZbmEZNz+7nB4pbfjDtWOJi9WvlUjQdLSMNNl9r62joLSCud85g5SkhKDjiAjac5cm+svabN5ck81/XjSYod3bBx1HREJU7tJoh46Uc88r6xjRswMzNM8uElE0LSON4u7cvXAN+SXlPP3tCcRrnl0komiLlEZ5MSOTt9fm8F9fGcqInjqNr0ikUbnLCduYU8B9r6/j7EFduOlcTceIRCKVu5yQPXnFfOupj2nfJk6n8RWJYJpzlwY7XFzO9U8tpayymhdnnkWPjm2DjiQi9VC5S4P95p3N7DlUwvybz2RINx32KBLJGjQtY2azzWy/ma2tNdbZzBaZ2ZbQ506hcTOzh81sq5mtNrNxzRVeTp71ewuYu3QX15/Zj/H9OgcdR0SOo6Fz7nOASceM/Rh4z90HA++FrgNcSs3aqYOBGcCjTY8pQXJ37n99HR3bxvPDi7XeuUhL0KByd/cPgGPXQr0SeDp0+Wngqlrjz3iNJUCKmfUIQ1YJyGur9rJ0Rx63fWUoHZPig44jIg3QlKNlurl7duhyDtAtdLkXsKfW/TJDY9IC5R0p5/7X1zO6TwpTT+8bdBwRaaCwHArp7g74iTzGzGaYWYaZZeTm5oYjhjSDn7+xnoKSCn55zUhiddijSIvRlHLfd3S6JfR5f2g8C+hT6369Q2Of4e6z3D3d3dPT0tKaEEOay3sb9rFwRRbfu2AQw7p3CDqOiJyAppT7a8ANocs3AK/WGv9W6KiZM4H8WtM30kKszjzMD+at4NQeHbjlwlOCjiMiJ6hBx7mb2TzgAiDVzDKBnwIPAPPNbDqwC5gcuvtbwFeBrUAxcGOYM0sz23ngCDf++RNSkhKYc+MEEuNig44kIieoQeXu7lPruemiOu7rwC1NCSXBqa52fjR/JVXuPDP9dK2FKtJC6dwy8hkLlmWyfPdhfvLVUxmUlhx0HBFpJJW7fOpwcTkP/GUj6f06cc243kHHEZEmULkLcPRdqOvJL6ng51edprM9irRwKncBYN7He1i4Iov/uPAUTu2hwx5FWjqVu7AmM5/7XlvHeUPS+MGFg4OOIyJhoHJv5YrLK/nB/66gS3ICv/v3MZqOEYkSOp97K/fA2xvZceAIz990Bp3bJQQdR0TCRHvurdgHm3N55l+7mH7OAM4elBp0HBEJI5V7K5VfXMEdC1ZzStdkbr9kaNBxRCTMVO6t1L2vreVAURm/nTyGNvE6vYBItFG5t0JvrN7Lqyv38oOLBjOyd8eg44hIM1C5tzKZh4q5++U1jO6TwvcuGBR0HBFpJir3VqSyqpr//N+VuMMjU8YSF6t/fpFopUMhW5EH39nEsl2HeHjqWPp2SQo6jog0I+26tRKPL97G44u3880z+3LF6J5BxxGRZqZybwWe/ddOfvH2Ri4f3ZP7rzgt6DgichKo3KPcgmWZ/J9X13HxqV15aPJoLXIt0ko0es7dzIYCL9QaGgjcC6QANwG5ofG73f2txn4daby312Rzx4JVnHNKKn+4dhzxegFVpNVodLm7+yZgDICZxQJZwEJq1kz9rbv/OhwBpXGKyyu5e+EaRvZOYda3xuuNSiKtTLh25S4Ctrn7rjA9nzTRixmZHCqu4J6vnUpSgg6KEmltwlXuU4B5ta5/38xWm9lsM+tU1wPMbIaZZZhZRm5ubl13kUaqrKrmyY+2M65vCun96vzxi0iUa3K5m1kCcAXwYmjoUWAQNVM22cBv6nqcu89y93R3T09LS2tqDKnl7bU57MkrYcZ5gzDTC6girVE49twvBZa7+z4Ad9/n7lXuXg08AZwehq8hDVRSXsXD721hYGo7vjy8W9BxRCQg4Sj3qdSakjGzHrVu+zqwNgxfQxrA3bnjpdVszS3i3suH67BHkVasSa+0mVk74MvAzbWGf2VmYwAHdh5zmzSjJz/cweur9nL7JUO5YGjXoOOISICaVO7ufgTocszY9U1KJI2yNiufX/5lI5NGdNfZHkVE71CNBmWVVfzXi6vo3C6BB64ZqRdRRURnhYwGv393CxtzCpk9LZ2UJC1yLSLac2/xVuw+xGOLtzE5vTcXDtPRMSJSQ+XegpVWVHHbi6vo3qEN91w2POg4IhJBNC3Tgv36r5vYnnuE56afQYc28UHHEZEIoj33Fur9Tft58qMdXH9mP84ZnBp0HBGJMCr3Fmjv4RJ+9MJKhnVvz0++dmrQcUQkAqncW5iqaucH81ZQXlnNH68bp1P5ikidNOfewsz+aAcZuw7x0OTRDEpLDjqOiEQo7bm3INtzi/j1O5u4+NRufH1sr6DjiEgEU7m3EFXVzh0LVpMYF8P/fP00vQtVRL6Qyr2FmPPPnWTsOsRPLx9B1w5tgo4jIhFO5d4C7DxwhAf/upGJQ9O4epymY0Tk+FTuEa66uuYc7fGxMfzi6lGajhGRBlG5R7hXVmbx8Y48fvLVU+neUdMxItIwKvcIdqSskgfe3sjo3h2ZnN4n6Dgi0oI0+Th3M9sJFAJVQKW7p5tZZ+AFoD81qzFNdvdDTf1arc0f39/K/sIyHrt+PDFaMk9ETkC49twnuvsYd08PXf8x8J67DwbeC12XE/DXdTk8+eEOrh7bi3F9OwUdR0RamOaalrkSeDp0+Wngqmb6OlHphU92893nljG8ZwfuvVyn8hWRExeOcnfgHTNbZmYzQmPd3D07dDkH+NwqEmY2w8wyzCwjNzc3DDGiw3sb9nHnS2s4Z3Aaz990hlZWEpFGCce5Zc5x9ywz6wosMrONtW90dzczP/ZB7j4LmAWQnp7+udtbo+z8Em57cRXDe3Rg1vXjdVIwEWm0Ju+5u3tW6PN+YCFwOrDPzHoAhD7vb+rXiXaVVdX8YN4KKiqr+cO1Y1XsItIkTSp3M2tnZu2PXga+AqwFXgNuCN3tBuDVpnyd1mDu0t18svMQP7/qNAbqbI8i0kRNnZbpBiwMvWsyDnje3f9iZp8A881sOrALmNzErxPVDhSV8Zt3NvGlU7robI8iEhZNKnd33w6MrmP8IHBRU567NfnVXzZSXF7F/VeM0OkFRCQs9A7VgC1Ylsn8jEymnzOAU7q2DzqOiEQJrcQUoDn/2MF9r6/n3MGp3HrxkKDjiEgUUbkH5Nklu7jv9fV8ZXg3Hrl2LIlxOjpGRMJH5R6At9Zkc++ra7loWFf+eN044mM1OyYi4aVyP4nWZuXz/Me7WZCRyfi+nfjDtSp2EWkeKvdmtmVfIW+tyeGtNdls2ldIYlwMV47pyT1fG07bBE3FiEjziMpyd3cOFVfQsW08sY08VW5VtWPwmVPtHi4uZ21WATkFpVRUVQOQGBfDwaJyVu45TElFFeP7dSItOZF1e/P5x7aDbN1fhBlM6NeZ+68YwVVjetExKT4c36aISL1afLm7O8t2HeKd9fvYe7iE/YVlbNlXyKHiCtrGxzK0e83hhQWlFXRKSqB3p7aktI2nbUIcSQmxtI2PZXdeMev25lPtkJqcyIGiMjbmFFBeWU3ndonExRjF5ZUUlFbWm6N3p7a0iY/lbxtrzrSQlBDLmD4pfOusfkwa0V2LWovISdWiy33VnsN8f95y9uSVkBAXQ6+UtqQlJzLptO4MTE1mb34JG7MLiYs1eqW0Je9IOct3H6KgpJKS8irKQ3vfyYlxDO/ZgcS4GDIPFdOhbTzXnt6PdomxHCgqo6raaRsfS/eObRnVuyN9OiWREBeD45RXVtMuMY7U5EQA8o6Uk19SQb/OSVpgQ0QC06LLvV+XJAamJnPrRUO45LTuJCee2LdTUVVNSUUVyQlxYSvizu0S6NxOp+kVkWC16HJPSUrg6W+f3ujHx8fG6GgVEYlKajYRkSikchcRiUIqdxGRKKRyFxGJQip3EZEopHIXEYlCKncRkSikchcRiULm7kFnwMxyqVlIu7FSgQNhitMcIj0fKGM4RHo+UMZwiKR8/dw9ra4bIqLcm8rMMtw9Pegc9Yn0fKCM4RDp+UAZwyHS8x2laRkRkSikchcRiULRUu6zgg5wHJGeD5QxHCI9HyhjOER6PiBK5txFROSzomXPXUREalG5i4hEoRZd7mb2QzNbZ2ZrzWyemQW+UKmZzTaz/Wa2ttZYZzNbZGZbQp87RWDGB81so5mtNrOFZpYSSflq3XabmbmZpQaRrVaOOjOa2X+Efo7rzOxXQeULZanr33mMmS0xs5VmlmFmjV/tpun5+pjZ+2a2PvTz+s/QeMRsL1+QMWK2l/q02HI3s17AD4B0dz8NiAWmBJsKgDnApGPGfgy85+6DgfdC14M0h89nXASc5u6jgM3AXSc7VC1z+Hw+zKwP8BVg98kOVIc5HJPRzCYCVwKj3X0E8OsActU2h8//HH8F3O/uY4B7Q9eDUgnc5u7DgTOBW8xsOJG1vdSXMZK2lzq12HIPiQPamlkckATsDTgP7v4BkHfM8JXA06HLTwNXncxMx6oro7u/4+6VoatLgN4nPdj/z1LXzxDgt8AdQOBHAdST8bvAA+5eFrrP/pMerJZ6MjrQIXS5IwFuM+6e7e7LQ5cLgQ1ALyJoe6kvYyRtL/VpseXu7lnU7BntBrKBfHd/J9hU9erm7tmhyzlAtyDDNMC3gbeDDlGbmV0JZLn7qqCzfIEhwLlmttTMFpvZhKAD1eFW4EEz20PN9hMRe5xm1h8YCywlQreXYzLWFnHbC7Tgcg/Nw10JDAB6Au3M7JvBpjo+rzn2NPA9z/qY2U+o+VN0btBZjjKzJOBuaqYRIlkc0JmaP99vB+abmQUb6XO+C/zQ3fsAPwSeCjgPZpYMvATc6u4FtW+LlO2lvoyRuL0c1WLLHbgY2OHuue5eAbwMnB1wpvrsM7MeAKHPgf65Xh8zmwZcBlznkfUGiEHU/Ce+ysx2UvMn8HIz6x5oqs/LBF72Gh8D1dScZCqS3EDNtgLwIhDYC6oAZhZPTWnOdfejuSJqe6knYyRvL0DLLvfdwJlmlhTaO7qImvmwSPQaNRsVoc+vBpilTmY2iZr57CvcvTjoPLW5+xp37+ru/d29PzUlOs7dcwKOdqxXgIkAZjYESCByzh541F7g/NDlC4EtQQUJbbdPARvc/aFaN0XM9lJfxkjeXj7l7i32A7gf2AisBZ4FEiMg0zxqXgOooKaEpgNdqHnVfwvwLtA5AjNuBfYAK0Mfj0VSvmNu3wmkRuDPMAF4LvT7uBy4MAIzngMsA1ZRM3c8PsB851Az5bK61u/dVyNpe/mCjBGzvdT3odMPiIhEoZY8LSMiIvVQuYuIRCGVu4hIFFK5i4hEIZW7iEgUUrmLiEQhlbuISBT6fxhrzLYd0/FkAAAAAElFTkSuQmCC\n" + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "plt.plot(xy[0],xy[1])" + ], + "metadata": { + "collapsed": false, + "pycharm": { + "name": "#%%\n" + } + } + }, + { + "cell_type": "code", + "execution_count": 109, + "outputs": [], + "source": [ + "def linear(x,y):\n", + " xybar = np.average( x*y)\n", + " xbar = np.average(x)\n", + " ybar = np.average(y)\n", + " x2bar = np.average(np.power(x,2))\n", + " k = (xybar -xbar*ybar)/(x2bar - xbar**2)\n", + " b = ybar - k*xbar\n", + " return k,b" + ], + "metadata": { + "collapsed": false, + "pycharm": { + "name": "#%%\n" + } + } + }, + { + "cell_type": "code", + "execution_count": 110, + "outputs": [], + "source": [ + "def loss(x,y):\n", + " k,b = linear(x,y)\n", + " yy = k*x+b\n", + " l = np.average(np.power(y-yy,2))\n", + " return l" + ], + "metadata": { + "collapsed": false, + "pycharm": { + "name": "#%%\n" + } + } + }, + { + "cell_type": "code", + "execution_count": 111, + "outputs": [], + "source": [ + "xy = np.sort(xy,1)" + ], + "metadata": { + "collapsed": false, + "pycharm": { + "name": "#%%\n" + } + } + }, + { + "cell_type": "code", + "execution_count": 112, + "outputs": [], + "source": [ + "startlinearloss = []\n", + "endlinearloss = []\n", + "for i in range(120):\n", + " startlinearloss.append([x[i+1],loss(x[0:i+2],y[0:i+2])])\n", + " endlinearloss.append([x[150-i-2], loss(x[150 - i -2:150],y[150 -i -2:150])])\n", + " # print(i,200-i-2,x[200-i-2],loss(x[200-i-2:200],y[200-i-2:200]))" + ], + "metadata": { + "collapsed": false, + "pycharm": { + "name": "#%%\n" + } + } + }, + { + "cell_type": "code", + "execution_count": 113, + "outputs": [ + { + "data": { + "text/plain": "150" + }, + "execution_count": 113, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [], + "metadata": { + "collapsed": false, + "pycharm": { + "name": "#%%\n" + } + } + }, + { + "cell_type": "code", + "execution_count": 114, + "outputs": [], + "source": [ + "startlinearloss = np.array(startlinearloss)\n", + "endlinearloss = np.array(endlinearloss)" + ], + "metadata": { + "collapsed": false, + "pycharm": { + "name": "#%%\n" + } + } + }, + { + "cell_type": "code", + "execution_count": 115, + "outputs": [ + { + "data": { + "text/plain": "[,\n ,\n ]" + }, + "execution_count": 115, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "text/plain": "
", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD4CAYAAAAXUaZHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAA4RElEQVR4nO3dd3xUVf7/8ddJMumNEBJKgFBC74TeRRBBKaKIioKriw2su19c/e6qu35/uruuru4qiqKIAiICglIkoNKUEnqH0FJIQhpJSJ+Z8/vjBozS0u/M5PN8PPLI5M7MvR9C7nvOnDn3HKW1RgghhGtxM7sAIYQQ1U/CXQghXJCEuxBCuCAJdyGEcEES7kII4YI8zC4AIDQ0VEdGRppdhhBCOJVdu3ala60bXO0+hwj3yMhIYmNjzS5DCCGcilLq7LXuk24ZIYRwQRLuQgjhgiTchRDCBUm4CyGEC5JwF0IIFyThLoQQLkjCXQghXJBDjHMXQghXYbdrsvKLSb9YTEZeEQAWdzcs7m54uClyC61k5ReTkVdMVl4x3ZsFMyjqqtchVYmEuxBCAFprcgqtaK0pttnJzCsmp8CKuxt4uLnh7qaw2jXHU3I5cT6XlJwisgtKaBLsTcNAH44k57A34QLncwuxV2CZjMeGtpJwF0KI8tBak5JTeDmcwwK9CfS2kF9sZc2BFH46mUFCZj55xVYi6/vh6eHG9lMZnMsuLNf+vS1uNAz0JtDHwsGkbDLzimkW4kvfliFE1PMl1N+T0AAvQnw9UUphtdspsdmx2jT+3h6E+HkS4utJsK8nnh410zsu4S6EcGp5RVYSsvJJvlDI4eQc9sRfYF/iBdJyi371uJYN/DifU8TFIiuh/l60DPUj1N+Lg+eyySuy0btFPaY1Db7cfRLi50Wgjwd2DTa7nRKbxk0pWof50zzEFzc3dXnfhSU2vC3utf1Pvy4JdyGEU7DbNQeSstkSl06x1Y6Hm2LHmUy2ncqgxPZLP0jLUD8GtQ6la9NgQv29sNrtxGfksy8xm+jm9bizZ1N6RdZDKXWdo1WMowU7SLgLIRyE3a7ZdiqDmCOpZBeUUGLTFFttFJbYScstIjErn5xC66+e0yLUj98NaEGnJkE0CvImKiyAIF+LSf8CxyLhLoQwTbHVzneHUth8Io1Nx9NJySnEx+JOiJ8nXh5ueJZ+NQzypluzYHpF1mNImzCCfSwU2+x4ebhVawvclUi4CyFq3cUiKzGHU3gz5jgJmQUEenvQv1UoY7o0YkSH8HJ1c3i7OV5XiCORcBdC1Ioiq40lsYks2B7P0ZQctIb2jQL5eFpHhrQJw91NWuDVScJdCFGj4s7nsnJfMktiE0jOLqRrRBBP3hRFr8gQ+req/6tRJ6L6SLgLIapddkEJi3bEs2LvOY4k56AU9G9Vn3/c2YWBrUOln7wWSLgLIarVD0fP8/yy/aTmFNGtaTB/ua0Dt3VpRFigt9ml1SkS7kKIanE0JYe3Yo7z3aFU2oT7M+f+aLo2DTa7rDpLwl0IUWnZ+SUs35PI2kMpbD+dib+nB8+NaMP0IS3x8pDRLGaScBdCVMqaA8n8ecUh0i8WERXmz9PD2zC1f3OCfT3NLk0g4S6EKCe7XXM4OYfNJ9L54eh5dpzJpGPjQD6eFk2XiGCzyxO/IeEuhLiuwhIb//zuGCv2JpF+sRiAdg0DeHF0ex4cEImHu6z544gk3IUQ13Qq7SJPLNzDkeQcxnRpxE1twxgUFSojX5yAhLsQ4golNjsfbj7F2+tP4OvpzicP9mJY2zCzyxIVIOEuhLhMa80Px87z2uqjnDh/kVs6hvPK2E40DJKWurORcBdCABB3/iIvrTzI1rgMIuv78uED0YzoEG52WaKSJNyFqOMSs/L5fFs8H285jbfFjVfGduTePs2wyAelTk3CXYg66kJ+MX9Ysp8NR1MBGNe1MS+O6UCDAC+TKxPVQcJdiDro3IUCpn68g7MZ+cwc1pq7opvSNMTX7LJENbphuCulmgLzgXBAA3O01m8rpUKAxUAkcAaYpLXOUsZ0b28Do4F8YJrWenfNlC+EKC+tNct2J/Hj8TQ2n0jDZtPM+10v+rcKNbs0UQPK03K3As9prXcrpQKAXUqpGGAasEFr/bpS6nngeWAWcCsQVfrVB5hd+l0IYZJiq51ZS/ezfE8S4YFeDI5qwGNDW9G+UaDZpYkacsNw11onA8mlt3OVUkeAJsA4YGjpwz4FfsQI93HAfK21BrYppYKVUo1K9yOEqGUnUnP584qDbDuVyR9GtuGJYa1lPvU6oEJ97kqpSKA7sB0ILxPYKRjdNmAEf0KZpyWWbvtVuCulpgPTAZo1a1bRuuscrTWHMw4TERBBkFeQ2eUIB2e3a7afzmTB9rOsOpCMj8WdNyd15Y4eEWaXJmpJucNdKeUPLAWe1lrnlH3l11prpZSuyIG11nOAOQDR0dEVem5ddOLCCSavmsxL/V7izjZ3ml2OcFA2u+brPUm8veEE8Zn5+Ht58OiQVvx+UEtC/GS2xrqkXOGulLJgBPsCrfWy0s2pl7pblFKNgPOl25OApmWeHlG6TVRBVHAUEf4RrI9fL+EurupMeh6PL9jN4eQcOjcJ4t93d+OWjg3x8ZR51euiG16lUDr6ZS5wRGv9Zpm7VgJTS29PBVaU2f6AMvQFsqW/veqUUtzc/Ga2J28npzjH7HKEgzl0Lps73/+J5OwC3p7cjRVPDGB89yYS7HVYeVruA4D7gQNKqb2l214AXge+VEo9BJwFJpXetxpjGGQcxlDIB6uz4LpseLPhzDs0j02Jm7it5W1mlyNMpLVmSWwiH2w6SbHNTlpuEfV8Pflieh9ah/mbXZ5wAOUZLbMFuNZH68Ov8ngNPFHFusRVdGnQhTCfMDac3SDhXkdl5RWz7VQGC7bHsyUuna5Ng2kV6oe3pztPDGtNk2Afs0sUDkKuUHUibsqNm5rdxNdxX1NgLcDHQ07kuiIlu5A31h1j6e5EtIYAbw/+Nr4T9/VuhpubDGsUV5JwdzI3N7+ZL459wU9JPzG8+RVvnISLyS+2MmfTKT7YeAqbXfO7AS0Y3bkRXSKCZGIvcV0S7k6mZ3hPgr2CiYmPkXB3cT/FpfPMl3tJzSliTOdGzBrVjmb1Zf4XUT4S7k7Gw82Dm5rdxHdnvqPYVoynu4xddkVrD6bw5KI9NKvvy7v39iA6MsTskoSTkXB3QiOaj2DZiWX8fO5nhjQdYnY5oppkXCziu0OpbI1LZ83BZLo1Debjab0I9pUXcFFxEu5OqE/DPgR4BrDu7DoJdxdwJDmHDzef4tt9yRTb7DQM9ObePs14YXR7fD3lFBWVI385TsjibmFY02H8kPADJbYSLO4Ws0sSlZCWW8SbMcf4YmcCvhZ3Jvduyr19mtE2PEAm9hJVJuHupEY0H8HKkyvZlryNQRGDzC5HVMDehAt8tPkU3x1KQWt4sH8LnhoeRZCvvEiL6iPh7qT6Ne6Hn8WP9fHrJdydhM2u+e/3cby94TgB3ham9G3O/X2b07KBXFEqqp+Eu5PycvdiSMQQvo//nj/3/TMebvJf6aiOpeTy3aEUVh9I5mhKLuO7NebVCZ3x95L/M1Fz5K/LiY1sPpLVp1cTmxpL30Z9zS5HlFFstfNlbAKLdsRz6FwOSkG3psG8OakrE7o3kT51UeMk3J1Y/yb98fHwYf3Z9RLuDmTdoRT+b/URzmbk07FxIC/f3oExXRrTIMDL7NJEHSLh7sR8PHwY1GQQ68+u50+9/4S7m0zvaiatjT71f8Ucp214AJ882IuhbRpIK12YQiancHIjIkeQUZjBnvN7zC6lTkvLLeKPX+3nXzHHuaN7E76ZOZBhbcMk2IVppOXu5AY3GYyXuxfrzq4jumG02eXUOSfTLvL5trMs2hFPsdXOjGGteW5kGwl1YToJdyfna/G93DXzfO/ncVPyZqw2nM8t5I9L9rPxeBoeboqx3RozY1hrGdYoHIaEuwsYGTmS9fHr2XN+Dz3De5pdjss7mJTN7+fHciG/hD/e0pa7oiMIC/A2uywhfkXC3QUMiRhidM2cWSfhXsMOJmUz6YOfCfax8NVj/ejYOMjskoS4KnkP7wJ8Lb4MbDKQmLMx2LXd7HJc1vncQn4/P5ZgHwvLnxggwS4cmoS7ixjZfCRpBWkyaqYGaK3ZE5/Fw58aXTEfTo0mPFC6YYRjk24ZFzGkqdE1892Z76RrphrtTbjAX1YcZH9iNn6e7rw9uZu02IVTkJa7i/Cz+DE4YjDrzqzDZreZXY7TKyyx8fqao9zx3lbScov42/hObH/xZkZ2bGh2aUKUi7TcXcgtkbcQczaG2NRY+jTqY3Y5TmtPfBZ/WLKPk2l5TO7VlBfGtCfQW6bjFc5Fwt2FDI4YjI+HD2vPrJVwr4TMvGLeWHeML3bE0zDQm/m/683gNg3MLkuISpFwdyE+Hj4MazqMmLMxvNDnBSxu0tq8keyCEt79IY6DSdnsT8ymoMTGA/0ieW5kGwKktS6cmPS5u5hRkaPILspme/J2s0txeFabnRkLdzN3y2nyim3c1qURa54axMtjO0qwC6cnLXcXM6DJAAIsAaw5vYaBTQaaXY5De23NUTafSOcfE7swqVdTs8sRolpJuLsYT3dPbmp2E9/Hf0+xrRhPd0+zS3I42QUlvLb6CF/sTODBAZES7MIlSbeMC7ol8hZyS3LZmrTV7FIczu74LEa8uZEvYxN4ZHBLXhzd3uyShKgREu4uqG/jvgR5BbH2zFqzS3Eou+OzeGDuDnw83Vk5YyB/Gt0eD3c5BYRrkr9sF2Rxs3Bzs5v5MeFHCq2FZpfjEC4Fe6i/J19M70unJnKVqXBtEu4u6pbIW8i35rMlaYvZpZiubLAvmt6XRkE+ZpckRI2TcHdRvRr2IsQ7hG9OfmN2KaY5k57H+xtPSrCLOklGy7goDzcPxrYay+eHPye9IJ1Qn1CzS6o1xVY7L608xKId8QBEN6/Hf+7tLsEu6pQbttyVUh8rpc4rpQ6W2fayUipJKbW39Gt0mfv+pJSKU0odU0rdUlOFixubEDUBq7ay8uRKs0upNekXi5jy0XYW7Yjn94NasGXWML56rL8Eu6hzytNynwf8F5j/m+1vaa3fKLtBKdUBmAx0BBoD65VSbbTWMk2hCVoGtaRHWA+WnVjGgx0fdOlFm3MLS/h4yxk+2nKKYqudtyd3Y1y3JmaXJYRpbthy11pvAjLLub9xwBda6yKt9WkgDuhdhfpEFU1sM5GzOWfZlbrL7FJqTMbFIia89xNvrT9Ov5b1+WbmQAl2UedV5QPVGUqp/aXdNvVKtzUBEso8JrF02xWUUtOVUrFKqdi0tLQqlCGuZ0TzEfhb/Fl2YpnZpdSI7PwS7p+7g8SsfD5/qA9zHoimTXiA2WUJYbrKhvtsoBXQDUgG/lXRHWit52ito7XW0Q0ayLSqNcXHw4dRLUYRczaGi8UXzS6nWiVnF3DvR9uIO3+RD+6PZmBU3fnQWIgbqVS4a61TtdY2rbUd+JBful6SgLITdUSUbhMmGt96PIW2QmLOxphdSrXZE5/F2P9u5WxGPh880JMhMu+6EL9SqXBXSjUq8+ME4NJImpXAZKWUl1KqBRAF7KhaiaKquoR2ITIwkq/jvja7lCrTWjN3y2kmffAz3hY3lj3en2Ftw8wuSwiHc8PRMkqpRcBQIFQplQi8BAxVSnUDNHAGeARAa31IKfUlcBiwAk/ISBnzKaUY13ocb+9+m/iceJoFNjO7pEopLLHxzOK9rDmYws3tw3njri4E+8qsl0JcjdJam10D0dHROjY21uwyXFpqXiojl47k4c4PM7P7TLPLqbC8IivTP4tla1wGL45uz8ODWrj00E4hykMptUtrHX21+2T6gToi3C+cvo36surUKuzabnY5FZJ+sYgpc7fz88kM3rirK78f3FKCXYgbkHCvQ25reRtJF5PYc36P2aWU24nUXMa/u5XD53J4774e3NkzwuyShHAKEu51yPBmw/Hx8HGaycQOJmVz5/s/U1hiZ/Ej/RjVqdGNnySEACTc6xRfiy83N7uZdWfWUWQrMruc6zqYlM19H23H38uD5Y/3p1vTYLNLEsKpSLjXMbe1uo3cklw2Jmw0u5Rr2nE6k3s/3Ia/lwdfTO9L0xBfs0sSwulIuNcxfRr2IcwnzGG7ZtYeTGbK3O2EBnhJsAtRBRLudYy7mztjWo5hS9IWMgvLOx9c7fjh2HmeWLiHTo0DWfpofwl2IapAwr0Our3V7Vi1lTWn15hdymUHErN5YsFu2jUMYP5DfajnJxcnCVEVEu51UFS9KNqHtHeYRTzizufy4Lwd1PP15JNpvfD3kgXChKgqCfc6alzrcRzOOExcVpypdcSdz2XynO2AYv5DvQkL9Da1HiFchYR7HTUqchQeysPU1nvsmUwmz9kGwBfT+9Kqgb9ptQjhaiTc66j6PvUZGDGQb059g9VurfXjL9wezz2lwx0XP9KX1mES7EJUJwn3OmxC6wmkF6SzJWlLrR535b5zvLD8AP1bhbLiiYHSYheiBki412GDIgZR37s+y08sr7VjHknOYdZX+4luXo8PH4gmyNdSa8cWoi6RcK/DLG4WxrYay6bETaQXpNf48bLyinn08134e3vw3n098PSQPz8haoqcXXXc+KjxWLWVb09+W6PHKbLaeOSzXSRfKOT9KT1kVIwQNUzCvY5rGdSSbg26sSxuGTW1cIvWmv/5aj87zmTyxqSu9GweUiPHEUL8QsJdcEfUHZzOPl1j87y/FXOcFXvP8cdb2jK2a+MaOYYQ4tck3AW3RN6Cn8WPpSeWVvu+v9qVyDvfxzEpOoLHh7aq9v0LIa5Owl3ga/FldIvRrDuzjpzinGrb7+Kd8Ty/dD8DWtfn/yZ0lqXxhKhFEu4CgIltJlJoK2T1qdVV3pfNrnlpxUFmLT1Av1b1mT2lJxZ3+VMTojbJGScA6Fi/I+1D2vPV8a+q/MHqnE2n+PTnszw8sAWfTOtFoLeMZReitkm4i8vubHMnx7KOcSjjUKX3cTQlh7dijnNrp4a8OKY9HtJiF8IUcuaJy0a3GI2Phw9fHf+qUs8vstp4dvE+An08eHV8J+ljF8JEEu7iMn9Pf25tcSurT6/mYvHFCj23xGZnxsI9HE7O4bU7ulDf36uGqhRClIeEu/iVO6PupMBawOrT5f9g1WbXPL14LzGHU3llbEdGdAivwQqFEOUh4S5+pVNoJ9rWa8uXx74s9werb6w7xqr9ybwwuh1T+0fWbIFCiHKR9czEryilmNxuMq/8/Aq7z++mZ3jP6z5+3aEUZv94knt6N2P64Gq+SKkwB87thszTkJsChdlgLQAPH/ALheDmENYOQtuAh3QDCVGWhLu4wpiWY3hr11ssPLLwuuF+Oj2P55bso3OTIF66vUP1HPzcXjj6LRz/DlIOAJfePSjwCjBCvKQAyn4moNyhfmuIiIZ2Y6DlMPD0rZ56hHBSEu7iCj4ePkxsM5H5h+aTkpdCQ7+GVzwmO7+Eh+btxMNN8d59PfC2uFf+gDYrHFkB22ZD4k5QbtC0Lwx9HiJ6QYN24B8G7mXGyxfnQ9YZSDsC549A6iHjRWHvArD4QuuboeMEaD8W3OXPXNQ98lcvrmpy28l8euhTFh9bzFM9nvrVfSU2O48v3EVCVj4LHu5L05BKtpJtVtgzH7b8Gy6chZBWcOs/oPNd4HuDmSM9fSG8g/F1eX8lcGYLHPnGCPojK42um4FPQ9d7wSLTDIu6Q9XUNK8VER0drWNjY80uQ/zGMz88w87UnaybuA5fixHghSU2Zi7aQ8zhVP55Zxfuim5auZ2f/AHWPg9pR6FJTxj4LLQdDW7V9Bm/3Q7H18Lmf0FSLPiFQZ9HoNdD4FOveo4hhMmUUru01tFXu09Gy4hrmtpxKtlF2Sw7sQyAi0VWpn2y4/KQx0oFu90G61+Bz8aDtQjuXgAPb4D2t1VfsIOxr3aj4eH1MPUbaNgZvv8bvNkR1syC7KTqO5aoPvmZsGIG5KaaXYnTk3AX19QtrBs9wnow//B8Suwl/OXrg+w8k8Xbk7tVbshjwQVYeDdseRN6TIXHtxmhXpNXsioFLQbD/cvg0S3QYSzs/Aje6QbfPA0ZJ2vu2KJiTv0Is/vDvi8gYbvZ1Ti9G4a7UupjpdR5pdTBMttClFIxSqkTpd/rlW5XSql3lFJxSqn9SqkeNVm8qHkPdX6I5Lxk/r7pC5btSWLGsNaM69ak4js6fxQ+vMk4gce8CWPfqf0+8IadYcL7MHM3dJ9ifPj6n56w6F5I3le7tYhfFOfB6j/C/HHGiKiH1xsvwqJKytNynweM+s2254ENWusoYEPpzwC3AlGlX9OB2dVTpjDLoCaDaBHYisUn5tM5IoAZN7Wu+E6Ofwcf3QxFuUYXSa+Hqr/QiqjXHG57C54+CIP/CGe3wgeDYcmD0pKvbSfWG631HXOgz6MwfSM07mZ2VS7hhuGutd4EZP5m8zjg09LbnwLjy2yfrw3bgGClVKNqqlWYwK6hJHMIeKYwZWhhxeZl1xp+fg8WTYaQFjD9R2jer8ZqrbCAcLjpRXh6vxHyx9fCu71h1R/gYprZ1bm2rLPwxX2wYKJxncK01XDr3+X6hGpU2T73cK11cuntFODSZCJNgIQyj0ss3XYFpdR0pVSsUio2LU1OJEf17/XHOXKiFYGWBqxNXFT+J9rt8N2L8N2fjAuLfrcWgirRnVMbvIPgpv+FJ/dCjwcg9mOjT/7Hv0NRxSZQEzdQUgCb/gnv9oGT38Pwl+DxnyFygNmVuZwqf6CqjbGUFR5PqbWeo7WO1lpHN2jQoKpliBrww9Hz/Of7OCZFN2d616nEpsZyIO3AjZ9ot8GKJ2Dbu9D7EbhrPnj61XzBVRUQbnTXPLEdWt0EP/4/eKc77PgQrMVmV+fcrMXG7/HtbvD9q9DmFpixEwY9K1NH1JDKhnvqpe6W0u/nS7cnAWXHx0WUbhNOJuZwKo98vosOjQJ5ZWwn7mxzJwGWAD459MmNn7zuz7BvIQx9wXirXZ1DHGtDaBTc/Rk8tN64vfoPRnfNga+MdySi/KzFsGse/KeH8XsMaWF0wUz6FIIizK7OpVX2rFsJTC29PRVYUWb7A6WjZvoC2WW6b4ST+Hb/OR79fBftGwaw4OE++Hi642fxY1LbSaw/u57T2aev/eTd840We59HYeismh3mWNOa9oJpq+DeL413Hksfgg8GwdFVxucJ4tryM2HLW/B2F/jmKfAPhylL4cE10gVTS254hapSahEwFAgFUoGXgK+BL4FmwFlgktY6UxlL7/wXY3RNPvCg1vqGl57KFaqO42TaRW57ZwsdGwcy73e98ff6ZYaKjIIMbl12K0ObDuUfg/9x5ZNPbYTPJ0KLQXDvEtea08Vuh4NL4cfXIPMkNOoGQ/9kdC848wtYddLaGJ++cy4cXgG2ImMSt/4zjW4u+T1Vu+tdoSrTD4jLiqw27njvJ85dKGDt04MJD7xyHPq/d/2bjw9+zNKxS4mqF/XLHcn74JMxxlvt360Fn+DaK7w22ayw/wvY+A9jPpzG3WHw/0DbW+tueBXnGy98O+ZAyn7wCoKud0PPaRDe0ezqXJpMPyDK5bXVRzl0Loe/T+xy1WAHmNZxGr4WX2bvK3MJQ9ZZ+PxOI9DvX+a6wQ7Gu5HuU2DmLhj7XyjIgi/ugdkDYP8SI/zrAmsxnIiBlU/Cm+1g5Qxj4rbb3oLnjsDof0qwm8yF3jeLqvhw0ynm/XSGBwdEMrLjlVP8XhLsHcz9He7n/X3vczTzKO0CmsPiKcY8MdO+hcDGtVi1idwt0ON+6HoPHFhi9C8vexi+/yv0mwnd73OOEUIVYSuB+G1GK/3Qcii8AJ7+xruW6N9Bs351992LA5JuGcHyPYk8s3gfYzo34p17uuPudv0TNKc4h1FfjSK6YTTvZJcYI2Pu/dLof66r7HY4vsaYvjhxB3gHG90SvR6G4ErOnGm24jxI2GHMqpm4y7iStyjHWAmr/W3Q6U5oOVSmUjbR9bplpOVex+06m8Wsrw7Qr2V93ry76w2DHSDQM5D7O97Pe3vf41BSMh2HPF+3gx1KZ6EcY3zFb4Nt78FP78DWtyFqBHS7F9qMAouP2ZVeW36msVjK2a1w9ic4twfspd1MoW2MxU+iRhiB7hVgaqnixqTlXoclZxdw+3+24uflzoonBhDs61nu5+ae+oFRP86gu0cA/53yk/ONZa8NFxKMoaF7PoPc5F+6MDqMh9bDzQ36kgJjLv3UQ0brPGG78TOAm8WYY795f2PYYkQv4ype4XCk5S6uUFhi45HPdlFQbGXh7/tUKNjJTSVg2aNM9ffiP+557M84SJcGXWquWGcV3NSYu2bo83Bms3ER1NFvjT56Dx+IHAithhnfwzuBWxWWKrye4jxjPdq0Y8ZXwnZI3vtLq9w7CCJ6Q+c7jeUNm/SUOV5cgLTc6yCtNc9+uY/le5KYc3/P636AegW7zVhoI2EnedO+YfTWP9AquBVzR85FyYdpN2YrMYL++DqIi4GMOGO7xddYK7ZBO2MOnqCmENLS+ApoeOPgL86Hi6nGV3aiMUzz/FFIPWi0yHXplbXuXsbwzeb9jO9hHY1jyDsvpyQtd/ErH24+xfI9STw3ok3Fgh2MSZ9Ob4Jx7+IX0YvpXabz+o7X2XpuKwObDKyZgl2Ju8W4oKfVTcDrxopQZ7YY/dvnDxvz3V9M+SWMAdw8wL8h+DcwlghUbsYFQ9pujFLKOm10+/xWYISxxmz720uDvL3xolFT7xCEQ5GWex2z9mAyjy3YzehOjfjvvd0r1to+sxXmjYGuk2H8bFCKElsJt399O/4Wf768/UvclLQAq8xmhZwkyDxlBHd2IuScg7w0YzWrS/P0KXfjxSK4OdRvZQxD9QszLiQLaiIfetYB0nIXAOyJz+KpL/bSrWkw/5rUtWLBXlII3zxpLHQx5l+XxzNb3C3M7D6T5zc/z8qTKxnfenzNFF+XuHsYv+d6zYFhZlcjnJQ0s+qInWcymfbJTsIDvfnwgWi8LRV8a77lLaN/eMybV1ycc2uLW+kS2oV/7/o3F4tl/nMhHIGEex2w5kAy9320nfp+nix4uA+h/hWcPzvtmLGodee7jCF8v+Gm3PhTnz+RUZjBnP1zqqlqIURVSLi7uHlbT/P4wt10ahzIV4/1p2lIBYe4lRQaU916+sMt/++aD+sU2okJrSfw2ZHPrj8lsBCiVki4u7B/rTvGy98cZkT7cBb+vi8hfhUYy35JzF+MMdLjZ4N/2HUf+mSPJ/Fx9+G17a/hCB/UC1GXSbi7qH0JF/jP93Hc2TOC2VN6VryPHYxFKXZ8AH0fh7ajbvjwUJ9QZnSfwc/JPxNzNqYSVQshqouEuwvSWvPamiPU9/Pk5bEdyzVfzBVyU2DFDGjUFW5+udxPm9R2Eu1C2vH3nX8nvyS/4scVQlQLCXcX9OOxNLadyuTJ4VG/Wkmp3LQ2FrguKYA7PqrQAsYebh682OdF0vLTeCP2jYofWwhRLSTcXczFIiv/t/oIzev7ck/vZpXbyc6PIG49jPwbNGhT4ad3C+vG1I5TWXJ8CRsTNlauBiFElUi4u5CCYhsPzdvJ6fQ8/jauE54elfjvzU2B9a8Yl8f3erjStczsPpM29drwl5/+QnpBeqX3I4SoHAl3F2G12Xl8wS52nMnkzUldGdymQeV2tO7PxsLGo9+o0qo6nu6evD7odS4WX+SvP/9VRs8IUcsk3F3EP747xg/H0vjbuE6M69akcjs5sxUOfAkDnjLmKqmiqHpRPNnjSX5I+IEVJ1dUeX9CiPKTcHcBK/YmMWfTKe7v25wpfZtXbifWYlj1nDFr4MBnq622Ke2n0DO8J6/veJ1zF89V236FENcn4e7kDp3LZtbS/fSKrMefb+tQ+R399DakHTEmBavGhRrc3dx5dcCrAMzaNIsSe0m17VsIcW0S7k4sK6+YRz7bRbCPJ+/e16NyH6ACpMfBxn8ay7/VwFqoEQERvNTvJfam7eW9ve9V+/6FEFeScHdSJTY7Mxbt5nxOEbOn9CAsoJIr0GsNq54FD2+49e/VW2QZt7a4lYlRE5l7YC6bEzfX2HGEEAYJdyekteaFZQfYGpfBqxM60b1Zvcrv7MhKOL0Rhv/ZWM6tBs3qPYs29dowa9MsmVxMiBom4e6E3oo5zpJdiTw5PIpJ0U0rv6PifPjuRWNx5p4PVl+B1+Dj4cM7N72Dxd3Ck98/SU5xTo0fU4i6SsLdyazcd453vo9jUnQEz9wcVbWdbX0bshPg1n8Yq//Ugsb+jXlz6JskXkxk5oaZFFgLauW4QtQ1Eu5O5HhqLrO+2k9083q8Or5zxZbJ+60LCbD139DxDogcUG01lkfP8J68Nug19qbt5ZkfnqHYVlyrxxeiLpBwdxI5hSU88tku/L09qjYy5pINrxjfR/y16sVVwqjIUbzc72W2ntvKrE2zsNqtptQhhKuScHcCdrvmuS/3kZCZz7v39iA8sJIjYy5J2AEHlkD/mRBchT77KpoQNYFZvWaxPn49L/30EnZtN60WIVxN7XS0iiqZvfEkMYdT+cttHejdIqRqO7PbYe2fwL8hDHi6WuqriikdppBbnMt7+97Dx8OHF/q8gJuSNocQVSXh7uB2nM7kX+uOMbZrYx4cEFn1HR5YAkmxMO498PKv+v6qwaNdH6XAVsAnBz9Ba82LfV+UgBeiiiTcHdjFIivPLdlLRD1fXrujih+gAhTnwfqXoXF36HpPtdRYHZRSPNPjGdxwY+7BuRTaCnm538tY3C1mlyaE06pSuCulzgC5gA2waq2jlVIhwGIgEjgDTNJaZ1WtzLrp/1YdJjGrgCWP9MOvMisq/daWf0PuObhrHrg5VstYKcVTPZ7C28Obd/e+S2p+Km8NfYsAzwCzSxPCKVXHGT5Ma91Nax1d+vPzwAatdRSwofRnUQFWm53XVh9h0Y4EHhnciujIKvazA2QnwU//gU4ToVmfqu+vBiileLTro7w64FV2pexiyuopJOQkmF2WEE6pJppv44BPS29/CoyvgWO4rOyCEu79aDsfbDrFlL7NeHZExZe5u6of/g+0HYa/VD37q0HjWo/jgxEfkFGYwT2r72Fb8jazSxLC6VQ13DWwTim1Syk1vXRbuNY6ufR2ChB+tScqpaYrpWKVUrFpaWlVLMM1FFltTJ8fy574LN66uyuvju9c9fHsAMn7Ye9C6PMI1KvkfO+1rHej3iwcvZD63vV5JOYR3t/3Pja7zeyyhHAaVU2OgVrrHsCtwBNKqcFl79TG2mpXXV9Naz1Hax2ttY5u0KCSS8K5ELtd88cl+9l+OpM37urKhO4R1bfzmL+ATzAMeq769lkLmgU2Y9GYRdza4lbe3fsuj8Q8QmpeqtllCeEUqhTuWuuk0u/ngeVAbyBVKdUIoPT7+aoWWRfM3XKalfvO8T+j2lZ+mbyrid8Op34wgt0nuPr2W0t8Lb68NvA1/tr/r+xP38/EbyYSczbG7LKEcHiVDnellJ9SKuDSbWAkcBBYCUwtfdhUQBbPvIHD53L453fHGNkhnMeGVH3t0l/Z/Ab4hED076p3v7VIKcWEqAksvm0xjf0a8+yPz/Lsj8+SXpBudmlCOKyqtNzDgS1KqX3ADmCV1not8DowQil1Ari59GdxDQXFNp5evIcgXwuvT+xS9bHsZSXvgxProN/j4OlXffs1SYugFiwYs4CnejzFxoSNjPt6HMtPLMfo/RNClKUc4cSIjo7WsbGxZpdR6xKz8nnks10cTs5h3oO9GdKmmj97WHw/nNoIzxwA76Dq3bfJTmWf4pWfXmH3+d30adiHF/q+QMuglmaXJUStUkrtKjMM/Vcc60qWOmTX2Uxu/88W4jPy+eiB6OoP9qyzcOQb6PWQywU7QMuglnwy6hP+t8//cjjjMBNXTuTN2DfJKpTr5YQACXdTbDyexn0fbSfY15MVMwYwvP1VR4tWza55oJRT97XfiJty4+52d7NywkpGtxjNvEPzGLV0FO/sfofsomyzyxPCVNItU0usNjs/n8rg233JLNuTSFRYAPMf6k2ov1cNHKwI3uwATfvAPQurf/8O6uSFk8zeN5t1Z9bha/Hlvvb38UCHBwjycr13LkLA9btlZOKwamK3Gy+Sbm6/fCCamJXPpuPpbDx+np/iMsgtsuLv5cH4bk348+0dCPSuoYmxDq+E/HSjS6YOaRXcijeGvMGJLieYvW82c/bPYeGRhTzQ4QGmdZqGj4eP2SUKUWtcpuVeZLVxIb+EtNwiTqfnkZxdQONgH1qE+mG1afKKrYT4edIw0Bs/Lw8s7r/0SOUWlnAqLQ8Pd0Wwryfncwo5mZYHQD1fI4ALS+wkZOVzPDWXi4XGqkFKgUKRnFPI8ZRclIJOjYPwsrhxIvUiKTmFADQJ9mFwm1CGtg1jSJsGeFvcq/RvvaG5IyEvDWbscrgJwmrTscxjzN43mw3xG2js15jnez/P0KZDq3dEkhAmul7L3anDPelCAd/uO8f6I6nsOpuFvQL/FA83hbfFHYu7Iiu/pNzPCw/0op6vJwBag0YT4udJ+0aB2O2a/UnZlNjsRIUF0KlJEEPahNKqgX/tBcqpjTB/LIx6Hfo+VjvHdHCxKbG8uu1VTmafpEdYD2Z0n0F0eLSEvHB6Lhvuaw4k89iC3XRsHMiQNg1oUs+HEF9Pmtf3o0mwD4kX8onPyMfTww0fT3cy84pJyS6koNhGQYmNwhI7RVYbjYN9aB3mj9aQlV9MqL8XrRr44e5mBL8CvC3uNAz0JsjXgecY1xrmjoCcczBzN1iquByfCymxlbD0xFLm7J9DWkEa7UPac0+7exjdcjRe7jXwuYcQtcBlw72g2EZWfjGNg6UvFYBja2DRZLj9beg5zexqHFKhtZAVcSv44tgXxF2II9grmAlRExjebDgd63fEw00+hhLOw2XDXZRht8MHg6AkH57YAbKK0XVprdmZspOFRxfyQ8IP2LWdQM9ARkWO4o6oO+hQv4N02wiHJ6Nl6oKj30DqQbjjQwn2clBK0btRb3o36s2FwgtsS9nGjwk/suLkCr48/iVt67VlYpuJjG4xWoZSCqckLXdXcKnVbi2CJ7aDWw2PxnFhucW5rD61mqUnlnIk8wiebp4Mbzac8a3H06dRH9zldysciLTcXd3Rb41W+4Q5EuxVFOAZwN3t7ubudndzOOMwX8d9zapTq1hzZg1hvmHcGXUnd7W9i1CfULNLFeK6pOXu7Ox2+GBwmb52eb2ubsW2Yn5M+JFlccvYmrQVDzcPeoX3YkCTAfQI60G7kHZYpCtMmEBa7q7syApIPQATPpBgryGe7p6MjBzJyMiRnM05y1fHv2Jz4mbeiH3DuN/Nk06hnYw+/Ia96dagm4S9MJ203J2ZzQqz+4Fyg8d+ki6ZWpaal8q+tH3sT9vPrtRdHM48jF3b8fHwoVNoJ9qHtKddSDvah7SnRVAL6a8X1U5a7q5q/2JIPw6TPpNgN0G4Xzgj/YwWPRgfxu5M2cm25G0cTD/I4mOLKbIVAeBn8aNrg670btibwRGDaR3cWoZaiholLXdnZS2C/0YbS+hN/9GY6EY4FKvdyuns0xzJPHK5dR93IQ6AMN8wejXsRd9GfRnYZKB8QCsqRVrurmjnXLgQb1yNKsHukDzcPIiqF0VUvSjGthoLGF05m5M2sz15O9vObWPVqVUAtA9pT59GfejXqB+9GvaSPntRZdJyd0YFF+CdbtC4O9y/3OxqRCVprTmWdYxNiZv4+dzP7EvbR4m9hABLAEObDuWOqDvoGd5Tum/ENUnL3dVsedMI+JtfMbsSUQVKKdqFtKNdSDumd5lOgbWAHck72BC/gfVn1/PNqW+IDIzk5uY3M7TpUDrV7yQfyopyk5a7s7mQAP/pCZ3ugAnvm12NqCEF1gLWnl7Lt6e+ZVfqLmzaRoh3CAObDGRQxCAGNB5AgGeA2WUKk8nEYa5k6cPGwtczYiG4qdnViFqQXZTNlqQtbErcxNZzW8kuysbDzYMhEUMY12ocA5sMlD76Okq6ZVxF4i44sAQGPSfBXocEeQUxpuUYxrQcg9VuZX/aftbHr2fVqVVsiN9AoGcgI5qPYGrHqbQIamF2ucJBSMvdWWgNn9wKGXHw5B7wkrfkdV2JvYSfz/3M6tOr+T7+e4ptxYxtNZZpHafRMril2eWJWiAtd1dwYAnE/wy3vSXBLgCwuFkYHDGYwRGDSS9IZ+6BuSw+tpjlccvp3bA3d7e9m2HNhmFxky6bukha7s4gPxP+2wvqRcJD6+RqVHFNGQUZLI9bzpJjSziXd44wnzDubX8vd7e9G39Pf7PLE9VMPlB1ditnwp4F8MhGaNjZ7GqEE7DZbWxO2syCIwvYlryNAEsA97S/hwc6PCCLj7gQCXdnFrcePp8I/Z+EkX8zuxrhhA5lHOKj/R+xPn49vh6+3N7q9stLCQrnJuHurHJTYPYA8A+DhzeAp6/ZFQkndjzrOB8f/JiYMzEU24vpHtadhzs/zKAmg+QqWCcl4e6M7Db4bDwk7DQmBgtrZ3ZFwkVkF2Xzzclv+PTwp6TkpRDhH8GYlmMY12ocTQNliK0zkXB3NrYSWPZ7OLQcxv4XetxvdkXCBZXYSlh7Zi0rT65ke/J2NJreDXszvvV4hjcbjq9F3ik6Ogl3Z1JSCEumwfE1MPJV6D/T7IpEHZCal8qKkytYdmIZSReT8Hb3ZmjTodwSeQsDmgzAx8PH7BLFVUi4O4viPFh0D5zeCGP+Bb0eNrsiUcfYtZ295/ey6tQq1sevJ7MwEy93L3o17EXvhr3pFNqJjvU7SqveQUi4O4O8dPjiPkjcAePehW73ml2RqOOsdiuxqbFsTNjIlqQtnMk5A4C7cqddSDt6hPegV3gveoT3kOGVJjEl3JVSo4C3AXfgI63169d6bJ0P9xMx8PXjUHgB7vgQOo43uyIhrpBZmMnB9IPsS9vH7tTd7E/bT7G9GIAI/wg6hnakc2hnOoV2ok29NjJrZS2o9XBXSrkDx4ERQCKwE7hHa334ao+vU+Fus0J+BqQegIQdcORbOH8IwjrCHXOgYSezKxSiXIpsRRxIO8DetL0czjjMofRDnMs7d/n+cN9wGvk1ItwvnAY+DQjzDaO+T31CvEPwt/gbX57++Fp8sbhZ8HDzwF2546bcTPxXORcz5pbpDcRprU+VFvAFMA64arhX2enNsOo50Hag9MVKa7CXGGGKNn5GG4/R9tKfy1AKUL98/+32ilLKOIa2g91qfFkLwVZc9kHQrB+Meh16PggW74ofRwiTeLl7Ed0wmuiGv2RLekE6hzMOczzrOKezT5Oal8rRzKNsyt9EgbWg3PtWpeecUgqFQil1OfjdlTtKqcu33ZQbbrgZp2/pY1WZc1ajf7VfN+V21ReQssesTKO3bK3X/TeVfp/cbjKT202u8HHKq6bCvQmQUObnRKBP2QcopaYD0wGaNWtWtaN5+UNYe1Buv15P1M0C7h7Gdij9rso87tJjy4Z/2f/Usj9rrh7yV9teZh/K3ZgLxt3T+PL0A+8gCOtgtNJ96lX2Xy2Ewwn1Cb08mVlZWmvySvLIKMwgqzCLiyUXuVhykbziPPJK8iixl2C1W7FrOzZtM56Dvhyydm2/fN+l23ZtR6Ox2q2XH1v2OWVdCtyyzy0bwpee89sXgorQaOzafvX7rrL/et41e+6bNiuk1noOMAeMbpkq7axxd5j0aXWUJYSoAUop/D2Nbpjmgc3NLqdOqKnOrSSg7KVuEaXbhBBC1IKaCvedQJRSqoVSyhOYDKysoWMJIYT4jRrpltFaW5VSM4DvMIZCfqy1PlQTxxJCCHGlGutz11qvBlbX1P6FEEJcmwwoFUIIFyThLoQQLkjCXQghXJCEuxBCuCCHmBVSKZUGnK3k00OB9GospyZIjVXn6PWB1FhdHL1GR6qvuda6wdXucIhwrwqlVOy1Js5xFFJj1Tl6fSA1VhdHr9HR67tEumWEEMIFSbgLIYQLcoVwn2N2AeUgNVado9cHUmN1cfQaHb0+wAX63IUQQlzJFVruQgghfkPCXQghXJBTh7tS6hml1CGl1EGl1CKllOnr1CmlPlZKnVdKHSyzLUQpFaOUOlH63dTll65R4z+VUkeVUvuVUsuVUsGOVF+Z+55TSmmlVKgZtZWp46o1KqVmlv4eDyml/mFWfaW1XO3/uZtSaptSaq9SKlYp1dvE+poqpX5QSh0u/X09VbrdYc6X69ToMOfLtThtuCulmgBPAtFa604YUwvX3IKE5TcPGPWbbc8DG7TWUcCG0p/NNI8ra4wBOmmtu2Asbv6n2i6qjHlcWR9KqabASCC+tgu6inn8pkal1DCMtYK7aq07Am+YUFdZ87jy9/gP4BWtdTfgL6U/m8UKPKe17gD0BZ5QSnXAsc6Xa9XoSOfLVTltuJfyAHyUUh6AL3DuBo+vcVrrTUDmbzaPAy6tA/gpML42a/qtq9WotV6ntbaW/rgNY/UsU1zjdwjwFvA//GqRWnNco8bHgNe11kWljzlf64WVcY0aNRBYejsIE88ZrXWy1np36e1c4AjG+ssOc75cq0ZHOl+uxWnDXWudhNEyigeSgWyt9Tpzq7qmcK11cuntFCDczGLK4XfAGrOLKEspNQ5I0lrvM7uW62gDDFJKbVdKbVRK9TK7oKt4GvinUioB4/xxiBanUioS6A5sx0HPl9/UWJbDnS/gxOFe2g83DmgBNAb8lFJTzK3qxrQx9tT0lue1KKVexHgrusDsWi5RSvkCL2B0IzgyDyAE4+37H4EvlVLK3JKu8BjwjNa6KfAMMNfkelBK+QNLgae11jll73OU8+VaNTri+XKJ04Y7cDNwWmudprUuAZYB/U2u6VpSlVKNAEq/m/p2/VqUUtOA24D7tGNdANEK40V8n1LqDMZb4N1KqYamVnWlRGCZNuwA7BiTTDmSqRjnCsASwLQPVAGUUhaM0Fygtb5Ul0OdL9eo0ZHPF8C5wz0e6KuU8i1tHQ3H6A9zRCsxTipKv68wsZarUkqNwujPHqu1zje7nrK01ge01mFa60itdSRGiPbQWqeYXNpvfQ0MA1BKtQE8cZzZAy85BwwpvX0TcMKsQkrP27nAEa31m2Xucpjz5Vo1OvL5cpnW2mm/gFeAo8BB4DPAywFqWoTxGUAJRgg9BNTH+NT/BLAeCHHAGuOABGBv6df7jlTfb+4/A4Q64O/QE/i89O9xN3CTA9Y4ENgF7MPoO+5pYn0DMbpc9pf5uxvtSOfLdWp0mPPlWl8y/YAQQrggZ+6WEUIIcQ0S7kII4YIk3IUQwgVJuAshhAuScBdCCBck4S6EEC5Iwl0IIVzQ/wdda/g0GKXBOgAAAABJRU5ErkJggg==\n" + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "plt.plot(x[0:150],y[0:150],startlinearloss[:,0],startlinearloss[:,1],endlinearloss[:,0],endlinearloss[:,1])" + ], + "metadata": { + "collapsed": false, + "pycharm": { + "name": "#%%\n" + } + } + }, + { + "cell_type": "code", + "execution_count": 137, + "outputs": [], + "source": [ + "startlinearlossdiff = startlinearloss[0:50] +startlinearloss[1:51]+startlinearloss[2:52]+startlinearloss[3:53]+startlinearloss[4:54]-10*startlinearloss[5:55]+startlinearloss[6:56]+startlinearloss[7:57]+startlinearloss[8:58]+startlinearloss[9:59]" + ], + "metadata": { + "collapsed": false, + "pycharm": { + "name": "#%%\n" + } + } + }, + { + "cell_type": "code", + "execution_count": 139, + "outputs": [ + { + "data": { + "text/plain": "[]" + }, + "execution_count": 139, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "text/plain": "
", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX8AAAD7CAYAAACCEpQdAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAiqUlEQVR4nO3deXxU9b3/8ddnshASwhIStiyEJWwCIgRwQdwFlwq2amm1bq1Uq1brrVZrr/f6s6v2drG1VluXWluX1o1WBcWqKIoQFDXsYSdsYRHCFkjm8/sjo0YFFJIzJ5l5Px+PeczMOZMz7/OAx3vOfOcs5u6IiEhyiYQdQERE4k/lLyKShFT+IiJJSOUvIpKEVP4iIklI5S8ikoQCLX8zKzSzl81snpnNNbNrYtNzzOxFM1scu+8QZA4REfkkC3I/fzPrCnR197fNLBuYDYwHLgY2u/vPzexGoIO7/yCwICIi8gmBlv9n3szsGeD3sdvx7r429gHxirv3PdDf5ubmenFxcRxSiogkhtmzZ29097x9zUuNVwgzKwaOAN4COrv72tisdUDn/fzNRGAiQFFREWVlZXFIKiKSGMxsxf7mxeUHXzNrAzwBXOvu2xrO8/qvHvv8+uHu97p7qbuX5uXt88NLREQOQeDlb2Zp1Bf/39z9ydjk9bHhng9/F9gQdA4REflY0Hv7GHAfMN/df9Vg1iTgotjji4BngswhIiKfFPSY/zHAN4D3zWxObNoPgZ8Dj5vZN4EVwHkB5xARkQYCLX93fx2w/cw+Kcj3FhGR/dMRviIiSUjlLyKShFT+0mzV1NbxyMyVVGyoDjuKSMKJ20FeIl+Uu/PS/A3c9uw8VmzaSXarVP5wwVCOLdGxHiJNRVv+0qwsqdrOxQ/M4lsPlZEaMX47YQj5HVpzyQOzeGzWyrDjiSQMbflLs1C9ey+//08F909fRkZqCj86oz8XHV1MWkqEE/t14jt/e5sfPPE+Kzfv5L9O6Usksr+dyETki1D5S6jcnUnvruHHz86nqrqG80oLuH5MP/KyW330muyMNO6/eDi3PFPOXS8vYeXmXdxxzmAy0lJCTC7Ssqn8JTRLqrZzyzPlTK/YxOCCdvzpwlKGFLbf52vTUiL89OxBFOVk8YvJC1i3dRf3fKOUnKz0+IYWSRAqf4m73Xvr+P1/Krhn2hIy0lK4bfxAvj6iiJTPGcoxM644vhdFOZl87/E5XPLATP55xdGkpeinK5GDpfKXuHp5wQZumVTOqs27OPuIfH54ev9PDPF8EWcM7krE4Iq/vc2vX1zEDWP7BZRWJHGp/CUu6qLObf+ex4NvLKdXXhZ/v2wkR/fKPeTlnTaoK18bUcjdry5hVEluo5Ylkoz0fVkCt72mlsseKuPBN5Zz6TE9eP6a0U1S1v995gB65GZx3WPv8sHOPU2QVCR5qPwlUGu37uLcP77JKws3cNv4gdzypQGkpzbNf7vM9FTunHAEm3bUcOMT7xPPS5KKtHQqfwlMeeVWxt81nZWbdnDfxcP5xpHdm/w9Bua34/oxfZk8dx2PzlrV5MsXSVQqfwnE1HnrOfePb5Jixj+vOJoT+nYK7L2+Naono3rn8v/+NY+KDdsDex+RRKLylybl7tz3+jIu+2sZJZ3b8PSVx9C/a9tA3zMSMf7vvMPJSItwzaPvUFNbF+j7iSQClb80mfdXb+W8e97ktn/P49QBnXl04pF0apsRl/fu3DaD2885nLlrtnHL03PZumtvXN5XpKXSrp7SaOu37eaOKQt54u3V5GSm89OzBzFheGHcz79zyoDOfGtUD/78+jL+9d4avjK0gIuPKaZXXpu45hBpCayl7CFRWlrqZWVlYceQBnbvreO+15dx18sV1NY5l4wq5soTetM2Iy3UXOWVW3nwjeVMmrOGPXVRjuuTx8XHFHNcSZ5OCCdJxcxmu3vpPuep/OVg7KmN8s7KLby2eCNPvVNJ5Qe7GHNYZ354en+6d8wKO94nbNxew9/fWslfZ6ygqrqG/l3b8vi3jyQ75A8nkXhpluVvZmOB3wIpwJ/d/ecHer3KPxzuzrKNO3ht8UZeW1zFm0s2sWNPHSkRY1j3Dlx7ckmzP7p2T22Up9+p5IYn3uOqE3rz/TF9w44kEhcHKv9QxvzNLAW4CzgFWA3MMrNJ7j4vjDzJzN3ZvTfKjj21rNu6m2Ubd7C0agfLNm5n6cYdLKvaQXVNLQDdO2Zy9tB8ji3J46heHUMf3vmi0lMjnDe8kNcqNvLn15dywZHd6dIuPj9EizRXYf3gOwKocPelAGb2KDAOaPLyf2HuOqLuQP1Yr8WGfA1wwL2+AB2IuhN1qItG2bUnyu69dezaW1d/v6eOmtooe+ui1Ead2rooe2P3dVGnLlr/t/XLcKLRjx/XRZ06h+hHr/OP3ivqjjf4O4DUSISI1d+nRIyUiBGJGBbLX3//8fP6zI5/lP/jZX3I7OOx7praOnbW1LGjppYde2qJ7uPLX3771vTMy+Lsofn07ZLNqN65zW5Y52DdMKYvU8rX8asXF3L7OYeHHUckVGGVfz7Q8HDM1cDIT7/IzCYCEwGKiooO6Y2ueXQOu/Y2br9vM2idlkKr1AhpKRFSI0ZqSoTUFCMtVtCpKYaZETGImJESa+nUSIRWqfXlnWLUF7nZR/cWe/2Hf+fw0YdJXdSpjTp10Sh1sQ8piH1g4bEPrvplWoNlR+zjDweo/5D7kDu0SovQJj2VzFYpZKWnktUqlaxWKXTMakXPvCyKO2bROj3xLpRSmJPJhUd1577py7h0VA/6dQn2+AOR5qxZ7+rp7vcC90L9mP+hLOOpK48mGo0tj4/L80OfLuD6Eo3QOi2F1mkpZKRHSE+JfGLLWVquq07szeNlq/jZcwv4y6Ujwo4jEpqwyr8SKGzwvCA2rclp604aap+ZzlUn9uanzy3g9cUbGVXSvH+sFglKWEf4zgJKzKyHmaUDE4BJIWWRJHPhUcXkt2/NT5+bT3RfP3iIJIFQyt/da4GrgCnAfOBxd58bRhZJPhlpKdwwti/z1m7j6TmBfOEUafZCO7ePuz/n7n3cvZe7/ySsHJKcvjS4G4Py2/HLKQvZ3cgdAkRaIp3YTZJSJGLcdHo/1mzdzQPTl4cdRyTuVP6StI7ulcuJ/Trxh5cr2LJDl4GU5KLyl6T2g7H9qK6p5eEZK8KOIhJXKn9Jan27ZHNcnzwemrGCPbXRsOOIxI3KX5LepaN6UFVdw7Pvrwk7ikjcqPwl6Y0uyaV3pzbc9/oyWsopzkUaS+UvSc/MuOSYYsort1G2YkvYcUTiQuUvAnz5iALatU7j/teXhR1FJC5U/iJA6/QUvj6yiClz17Fq886w44gETuUvEnPhUd2JmPHQm8vDjiISOJW/SEzXdq05fVBXHp21iu2xq5eJJCqVv0gDl47qQfXuWp6YvTrsKCKBUvmLNDCksD1Di9rzwPRlOt2zJDSVv8inXDqqB8s37eQ/CzaEHUUkMCp/kU8Ze1gXurXL4P7p2u1TEpfKX+RTUlMiXHh0MW8s2cT8tdvCjiMSCJW/yD58bXgRGWkRne1TEpbKX2Qf2mWmcfqgrkyas4ade7TbpyQelb/IfkwYXkR1TS3Pvb8u7CgiTU7lL7Ifw4s70DMvi8dmrQw7ikiTC6z8zewOM1tgZu+Z2VNm1r7BvJvMrMLMFprZmKAyiDSGmTFheCGzlm+hYkN12HFEmlSQW/4vAgPdfTCwCLgJwMwGABOAw4CxwB/MLCXAHCKH7MtDC0iNGI/NWhV2FJEmFVj5u/sL7v7hL2UzgILY43HAo+5e4+7LgApgRFA5RBojt00rThnQmSfertRlHiWhxGvM/1Lg+djjfKDhZtTq2LTPMLOJZlZmZmVVVVUBRxTZt68OL2Tzjj1Mnb8+7CgiTaZR5W9mU82sfB+3cQ1eczNQC/ztYJfv7ve6e6m7l+bl5TUmqsghO7Ykj/z2rXlkpn74lcSR2pg/dveTDzTfzC4GzgRO8o8vjloJFDZ4WUFsmkizlBIxzi0t4LcvLWbV5p0U5mSGHUmk0YLc22cscANwlrs3vDTSJGCCmbUysx5ACTAzqBwiTeHc0vrtlX/oVM+SIIIc8/89kA28aGZzzOyPAO4+F3gcmAdMBq5097oAc4g0Wn771owuyeMfZauo06meJQEEubdPb3cvdPchsdvlDeb9xN17uXtfd3/+QMsRaS4mDC9k7dbdTFusnQ+k5dMRviJf0En9O9MxK51H9cOvJACVv8gXlJ4a4ZxhBbw0fwMbqneHHUekUVT+IgfhvOGF1Eadf+qHX2nhVP4iB6FXXhtG9c7lz68to3r33rDjiBwylb/IQbp+TF8279jDn6YtDTuKyCFT+YscpMML23PG4K786bVlbNimsX9pmVT+Iofg+lP7srcuym9eWhx2FJFDovIXOQTFuVmcP7KIx2atYknV9rDjiBw0lb/IIbr6pBIyUiPcPnlB2FFEDprKX+QQ5bZpxbeP68WUueuZvWJz2HFEDorKX6QRvjmqB7ltWvGz5xbw8YlrRZo/lb9II2S1SuXak0soW7GFqfM3hB1H5AtT+Ys00leHF9IzN4tfTF5AbZ0u9Sgtg8pfpJHSUiLcMLYvFRu267QP0mKo/EWawJjDunBEUXt+9eIidtTUhh1H5HOp/EWagJnxozP6s6G6hnt02gdpAVT+Ik1kWPcczhjclXunLWHt1l1hxxE5IJW/SBO6cWw/olG4Y8rCsKOIHJDKX6QJFeZkcsmoYp58u5L3Vn8QdhyR/Qq8/M3sv8zMzSw39tzM7E4zqzCz98xsaNAZROLpyhN60zErnR8/O18HfkmzFWj5m1khcCrQ8KKnpwElsdtE4O4gM4jEW9uMNL53Sh9mLtvMlLnrw44jsk9Bb/n/GrgBaLj5Mw54yOvNANqbWdeAc4jE1YThhZR0asPPnp/Pnlod+CXNT2Dlb2bjgEp3f/dTs/KBVQ2er45N29cyJppZmZmVVVVVBZRUpOmlpkS4+Yz+rNi0k4feXB52HJHPaFT5m9lUMyvfx20c8EPglsYs393vdfdSdy/Ny8trzKJE4u74vp0Y3SePO19azJYde8KOI/IJjSp/dz/Z3Qd++gYsBXoA75rZcqAAeNvMugCVQGGDxRTEpokknJtP78/2mlp+qyt+STMTyLCPu7/v7p3cvdjdi6kf2hnq7uuAScCFsb1+jgS2uvvaIHKIhK1vl2wmjCji4RkrWL5xR9hxRD4Sxn7+z1H/zaAC+BPwnRAyiMTNtSeVkJYS4VcvLgo7ishH4lL+sW8AG2OP3d2vdPde7j7I3cvikUEkLJ3aZnDpqGImvbuGuWu2hh1HBNARviJxMXF0L9q1TuP2yTrtgzQPKn+ROGjXOo0rT+jFq4uqeHPJprDjiKj8ReLlwqOK6dI2g9un6Hq/Ej6Vv0icZKSlcO3JJbyz8gNenKfTPki4VP4icXTOsAJ65mVxx5SF1EW19S/hUfmLxFFqSoTrT+3L4g3befJtXe9XwqPyF4mzsQO7MLigHb+Zupjde+vCjiNJSuUvEmdmxg/G9qPyg108PGNF2HEkSan8RUJwTO9cji3J5a6XK9i0vSbsOJKEVP4iIbnxtH7s3FPHxQ/Monr33rDjSJJR+YuE5LBu7bj7gqHMX7uNb/6lTOP/Elcqf5EQndivM/933uHMWr6ZKx6erat+Sdyo/EVCNm5IPj8eP5CXF1Zx3eNztP+/xEVq2AFEBM4f2Z3q3bX8/PkFZGek8dOzB2JmYceSBKbyF2kmLj+uF9t27eUPryyhbUYqN57WTx8AEhiVv0gzcv2YvmzbvZd7pi1lWPcOnHpYl7AjSYLSmL9IM2Jm3HrWQLIzUnl5YVXYcSSBqfxFmpmUiDG8OIe3luq8/xIclb9IMzSyRw5LN+5gw7bdYUeRBKXyF2mGjuzZEYC3lm0OOYkkqkDL38yuNrMFZjbXzG5vMP0mM6sws4VmNibIDCIt0WHd2tKmVSpvLdPQjwQjsL19zOwEYBxwuLvXmFmn2PQBwATgMKAbMNXM+ri7jm0XiUlNiTCsewdmLNWWvwQjyC3/K4Cfu3sNgLtviE0fBzzq7jXuvgyoAEYEmEOkRRrZM4eKDdvZqLN+SgCCLP8+wLFm9paZvWpmw2PT84FVDV63OjbtM8xsopmVmVlZVZV2e5Pk8uG4/0yN+0sAGlX+ZjbVzMr3cRtH/ZBSDnAkcD3wuB3k4Yrufq+7l7p7aV5eXmOiirQ4g/LbkZmeol0+JRCNGvN395P3N8/MrgCedHcHZppZFMgFKoHCBi8tiE0TkQbSNO4vAQpy2Odp4AQAM+sDpAMbgUnABDNrZWY9gBJgZoA5RFqsI3t2ZOH6ajbv2BN2FEkwQZb//UBPMysHHgUu8npzgceBecBk4Ert6SOybyN75AAa95emF9iunu6+B7hgP/N+AvwkqPcWSRSDC9qTkRbhrWWbGDtQJ3mTpqMjfEWasfTUCEOLOvCWxv2lian8RZq5I3t2ZP66bWzdqYu8S9NR+Ys0cyN75OAOM5dr61+ajspfpJk7vLA96akR7e8vTUrlL9LMZaSlMLSovc7wKU1K5S/SAozs0ZG5a7aybbfG/aVpqPxFWoCRPXOIOpRp3F+aiMpfpAUYWtSB9JSIdvmUJqPyF2kBMtJSGFLYnhka95cmovIXaSFG9syhvHIr22tqw44iCUDlL9JCjOzRkbqoa9xfmoTKX6SFGNq9PWkpxisLdWEjaTyVv0gLkZmeypmDu/HYrFW6tKM0mspfpAW56sTe1NTW8adpS8OOIi2cyl+kBemV14ZxQ/J56M0V2vqXRlH5i7QwV8e2/u/V1r80gspfpIXpmdeG8UPyeejN5dr6l0Om8hdpga46sTd7aqPa+pdDpvIXaYEabv1XVWvrXw6eyl+khfp4639J2FGkBQqs/M1siJnNMLM5ZlZmZiNi083M7jSzCjN7z8yGBpVBJJH1zGvD+CPy+euMFWyo3h12HGlhgtzyvx241d2HALfEngOcBpTEbhOBuwPMIJLQrj6xpH7r/1WN/cvBCbL8HWgbe9wOWBN7PA54yOvNANqbWdcAc4gkrB65WYw/Ip+H39LWvxycIMv/WuAOM1sF/BK4KTY9H1jV4HWrY9M+w8wmxoaMyqqqdD4TkX357okl7K1zbp00jy079oQdR1qIRpW/mU01s/J93MYBVwDfc/dC4HvAfQe7fHe/191L3b00Ly+vMVFFElZxbhbfOb4Xz5Wv5djbX+ZXLy5i6y5d7lEOzNw9mAWbbQXau7ubmQFb3b2tmd0DvOLuj8RetxA43t3XHmh5paWlXlZWFkhWkUSwcF01v5m6iOfL15Gdkcplx/bkkmOKyc5ICzuahMTMZrt76b7mBTnsswY4Lvb4RGBx7PEk4MLYXj9HUv+hcMDiF5HP17dLNndfMIxnvzuKI3t25FcvLvrom8B7qz8gGg1mQ09apiC3/EcBvwVSgd3Ad9x9duxbwO+BscBO4BJ3/9xNem35ixyc91Z/wK9fXMTLsfP/d8hMY1RJHseW5HJsSS5d27UOOaEE7UBb/oGVf1NT+Yscmo3ba5hesZFpizby2uIqNsSOCB7QtS23jjuM4cU5ISeUoKj8RQQAd2fh+mpeW7SRh2YsZ/WWXUwc3ZPrTulDq9SUsONJEwtrzF9Emhkzo1+Xtlw2uifPXzOaCcOLuOfVpZz1u+nMXbM17HgSRyp/kSTVplUqP/vyIB64eDhbdu5h/F3TuevlCmrromFHkzhQ+YskuRP6dWLKtaM59bAu3DFlIefe8ybz124LO5YETOUvInTISueurw/lzq8dwbKNOzjjzte46cn3dbGYBKbyF5GPnHV4N179/glcfHQP/lG2iuPveIV7Xl1CTW1d2NGkian8ReQT2mWmccuXBjDle6MZ2SOHnz2/gFN+NY3J5WtpKXsHyufTrp4ickCvLa7itn/PY9H67RTmtObk/p05pX9nhvfIIS1F24/NmfbzF5FGqa2L8sycNTz3/lper9hITW2UthmpnNivEycP6MxJ/TrTOl3HCTQ3Kn8RaTI799QybdFGps5fz38WbGDzjj0MLmjH498+iow0fQA0Jwcq/9R4hxGRli0zPZWxA7swdmAX6qLOpHcrue7xd/nBE+/xm68Oof70XdLcqfxF5JClRIyzjyigcssufvnCIvp3bcvlx/UKO5Z8Afq1RkQa7coTenPG4K78YvICXl6wIew48gWo/EWk0cyMO84ZTP8ubfnuo++wpGp72JHkc6j8RaRJZKancu+Fw0hPiXDZX8p0KclmTuUvIk2moEMmd18wjJWbd3LNo+9Qp6uHNVv6wVdEmtSIHjncOu4wbn6qnMsfnk3/rm1pm5FK24w0sjNSyc5Io6BDa4pzs8KOmtRU/iLS5M4f2Z1Vm3fx8IwVvDhv/Wfmm8E3j+nB98f01bEBIdFBXiISqLqos313Ldt276V6dy3Vu/fy7/fW8tcZK+jdqQ2/Pm8IgwrahR0zIQV2JS8zO9fM5ppZ1MxKPzXvJjOrMLOFZjamwfSxsWkVZnZjY95fRJq/lIjRLjONwpxMBnRry8ieHblt/EAeunQE23fXcvYfpvObqYvYq4vIxFVjf/AtB74MTGs40cwGABOAw4CxwB/MLMXMUoC7gNOAAcDXYq8VkSQzuk8eU64dzZcO78Zvpi7mK3e/QcWG6rBjJY1Glb+7z3f3hfuYNQ541N1r3H0ZUAGMiN0q3H2pu+8BHo29VkSSULvMNH791SHcff5QVm3eyRl3vs4jM1eGHSspBLWrZz6wqsHz1bFp+5u+T2Y20czKzKysqqoqkKAiEr7TBnVlyvdGM6JHDj96upxF6/UNIGifW/5mNtXMyvdxC3yL3d3vdfdSdy/Ny8sL+u1EJESdsjO4c8IRtGmVyv9OmqsLxwTsc3f1dPeTD2G5lUBhg+cFsWkcYLqIJLkOWel8/9Q+/Pczc5lcvo7TBnUNO1LCCmrYZxIwwcxamVkPoASYCcwCSsysh5mlU/+j8KSAMohIC/S1EUX065LNj5+dz649unZwUBq7q+fZZrYaOAp41symALj7XOBxYB4wGbjS3evcvRa4CpgCzAcej71WRASA1JQIt551GJUf7OKPry4JO07C0kFeItIsXf3IO7wwdx1TrzuOwpzMsOO0SIEd5CUiEpQfnt6PiBk/eXZ+2FESkspfRJqlru1ac9WJvZk8dx2vL94YdpyEo/IXkWbrm6N6UJSTyf/+a65O/9DEVP4i0mxlpKVwy5kDqNiwnb+8sTzsOAlFp3QWkWbtpP6dOK5PHndMWcjSjTv4+ogiBubrLKCNpS1/EWnWPrw+8JmDu/Hk26s583ev86Xf1Z8DaHtNbdjxWizt6ikiLcbWXXt5+p1K/v7WShauryYrPYWzhuRzbmkBRxS2x8zCjtisHGhXT5W/iLQ47s7bKz/gkZkr+fd7a9i9N0rP3Cy+MqyAs4/Ip1v71mFHbBZU/iKSsKp37+W599fyxOxKZi7fjBkc0yuXrwzL54xB3UhPTd7RbZW/iCSFFZt28MTblTz59mpWb9nF8OIO3H3BMHLbtAo7Wih0hK+IJIXuHbO47pQ+TLv+BH791cN5v3IrZ/3udcort4YdrdlR+YtIwolEjLOPKOCflx+NA+f88Q0mvbsm7FjNispfRBLWwPx2TLpqFAO7teO7j7zDLyYvoC7aMoa6g6byF5GElpfdir9fdiRfG1HI3a8s4bKHyqj8YFfSXylMR/iKSMJLT43w07MHMaBrW2791zyO+fl/aNc6jb5dsunbOZu+XbLp1yWbgfntyEhLCTtuXKj8RSQpmBnfOKqYkT078uaSTSxcX83CddU8/U4l1bEjhXvmZfHE5UfTISs95LTBU/mLSFLp0zmbPp2zP3ru7lR+sIuy5Vu44Yn3+PZfZ/PXb42gVWpifwPQmL+IJDUzo6BDJuOPyOf/zj2cmcs3c/0/3iOa4D8Ma8tfRCTmS4d3Y9WWndw+eSFFOZl8f0zfsCMFRuUvItLAFcf1YuWmnfz+5QqKcjI5b3hh2JEC0ahhHzM718zmmlnUzEobTD/FzGab2fux+xMbzBsWm15hZneaTsMnIs2ImXHb+IEcW5LLD596P2EvIdnYMf9y4MvAtE9N3wh8yd0HARcBf20w727gMqAkdhvbyAwiIk0qLSXCH84fSu9Obbji4dksXFcddqQm16hhH3efD3zmHNru/k6Dp3OB1mbWCsgB2rr7jNjfPQSMB55vTA4RkaaWnZHG/RcPZ/xd0zntt9NITdnHtrJ/eFf/4MPjxiIRo1+XbEq75zC8uAPDijvQKTsjTsm/mHiM+X8FeNvda8wsH1jdYN5qIH9/f2hmE4GJAEVFRYGGFBH5tG7tW/PIxCN58u3V7O/68R9u+1qD5zV7o7xfuZW/vbWC+6cvA6C4YyalxTlcdFQxgwrCvwzl55a/mU0Fuuxj1s3u/szn/O1hwC+AUw8lnLvfC9wL9ad0PpRliIg0Rq+8Nlw/pt8h/e2e2ijla7ZStnwzs5Zv4YW563hmTiU/OmMAFx7VPdQrj31u+bv7yYeyYDMrAJ4CLnT3JbHJlUBBg5cVxKaJiCSc9NQIQ4s6MLSoAxNHw5Yde/ivf7zL/0yay8zlm/n5lweRnZEWSrZADvIys/bAs8CN7j79w+nuvhbYZmZHxvbyuRA44LcHEZFE0SErnT9fWMoPxvZjcvk6zvr9dOav3RZKlkZdycvMzgZ+B+QBHwBz3H2Mmf0IuAlY3ODlp7r7htguoQ8Cran/ofdq/wIhdCUvEUkkby3dxNWPvMPWXXu5bdxAzhteyPaaWlZs2sHKTTtZsXknKzbtpC4a5fZzDj+k99BlHEVEmqGq6hqufewdpldson1mGh/s3PuJ+TlZ6ZR0asNj3z7qkJZ/oPLXEb4iIiHJy27FQ5eO5IHpy1hStZ2inCy6d8ykKCeToo6ZtA3w9wCVv4hIiFIixreO7Rn399VZPUVEkpDKX0QkCan8RUSSkMpfRCQJqfxFRJKQyl9EJAmp/EVEkpDKX0QkCbWY0zuYWRWwIuwcB5BL/RXMkkkyrjNovZNNS17v7u6et68ZLab8mzszK9vfOTQSVTKuM2i9w84Rb4m63hr2ERFJQip/EZEkpPJvOveGHSAEybjOoPVONgm53hrzFxFJQtryFxFJQip/EZEkpPJvJDO7xszKzWyumV0bdp6gmNn9ZrbBzMobTMsxsxfNbHHsvkOYGYOwn/U+N/bvHY1dkzrh7Ge97zCzBWb2npk9ZWbtQ4wYiP2s922xdZ5jZi+YWbcwMzYVlX8jmNlA4DJgBHA4cKaZ9Q43VWAeBMZ+atqNwEvuXgK8FHueaB7ks+tdDnwZmBb3NPHzIJ9d7xeBge4+GFgE3BTvUHHwIJ9d7zvcfbC7DwH+DdwS71BBUPk3Tn/gLXff6e61wKvUl0LCcfdpwOZPTR4H/CX2+C/A+Hhmiod9rbe7z3f3hSFFiov9rPcLsf/nADOAgrgHC9h+1ntbg6dZQELsJaPyb5xy4Fgz62hmmcDpQGHImeKps7uvjT1eB3QOM4zE1aXA82GHiBcz+4mZrQLOR1v+4u7zgV8ALwCTgTlAXZiZwuL1+wwnxBaRHJiZ3QzUAn8LO0u8uPvN7l5I/TpfFXaepqDybyR3v8/dh7n7aGAL9WOhyWK9mXUFiN1vCDmPBMzMLgbOBM735DxI6G/AV8IO0RRU/o1kZp1i90XUj/f/PdxEcTUJuCj2+CLgmRCzSMDMbCxwA3CWu+8MO0+8mFlJg6fjgAVhZWlKOsK3kczsNaAjsBe4zt1fCjlSIMzsEeB46k9vux74H+Bp4HGgiPrTbZ/n7p/+UbhF2896bwZ+B+QBHwBz3H1MSBEDsZ/1vgloBWyKvWyGu18eSsCA7Ge9Twf6AlHq/59f7u6VYWVsKip/EZEkpGEfEZEkpPIXEUlCKn8RkSSk8hcRSUIqfxGRJKTyFxFJQip/EZEk9P8B72/m9ZqdzQUAAAAASUVORK5CYII=\n" + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "plt.plot(startlinearloss[5:55,0],startlinearlossdiff[:,1])" + ], + "metadata": { + "collapsed": false, + "pycharm": { + "name": "#%%\n" + } + } + }, + { + "cell_type": "code", + "execution_count": 116, + "outputs": [], + "source": [ + "endlinearlossdiff = endlinearloss[0:50] +endlinearloss[1:51]+endlinearloss[2:52]+endlinearloss[3:53]+endlinearloss[4:54]-10*endlinearloss[5:55]+endlinearloss[6:56]+endlinearloss[7:57]+endlinearloss[8:58]+endlinearloss[9:59]" + ], + "metadata": { + "collapsed": false, + "pycharm": { + "name": "#%%\n" + } + } + }, + { + "cell_type": "code", + "execution_count": 134, + "outputs": [ + { + "data": { + "text/plain": "[]" + }, + "execution_count": 134, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "text/plain": "
", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAD4CAYAAAAJmJb0AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAfwklEQVR4nO3deXyU5b338c8vM5M9AbKwBkhYBVwAIyAotWrVatXW1lY9StVWa2t7Tk9Pn57HLk999XV6nu6L7emptPpYPdpdqscdobUVQQiC7EtYkxCyEMhC1pm5nj8y0BxNgDCZ3JOZ7/v1yiszc8/M/bsI+c6d677u6zLnHCIikphSvC5ARERiRyEvIpLAFPIiIglMIS8iksAU8iIiCczvdQE9FRQUuOLiYq/LEBEZUtavX1/vnCvsbVtchXxxcTFlZWVelyEiMqSY2YG+tqm7RkQkgSnkRUQSmEJeRCSBKeRFRBKYQl5EJIEp5EVEEphCXkQkgSnkRXpR09TO46v3s7my0etSRKISVxdDiXhtU+UxHn19H89tqiYY7l5r4aLiEdy9qISrZo3Gl2IeVyjSPwp5SQr76o+z/sBR8rNTGZmTxsicdPKzUklJMYKhMK9sq+HR1/dRduAo2Wl+llxczM2lRawqr+exN/bz6SffomhEBncuLOajF40nNz3gdZNEzojF08pQpaWlTtMayEDacPAoD7+2l5e3Head/9V9KUZBdiphB3XNHYzPy+CuhSXcXFpETo8QD4Udy7fV8Oiqfazd10BWqo8PzhnHrfMmcO64YYPcIpF3M7P1zrnSXrcp5CXRhMOOP++s5eHX9rJ2fwO56X7uuHgiN84eR3N7kLrmdmqbO6hpaqe2qYO2rhDXXzCWK2eMOm13zJaqRh5dtY/nN1XTEQwza2wut8ybwI2zx+roXjzjacib2TXAjwEf8Evn3Lf6eq5CXnoKhR0/+3M5Oel+5k/KZ/qoHFL6COGuUJith5pYu+8IvyurpLy2hXHDM7j7khI+dtF4stMGtmeysa2LZzZW8eu1FWyvbiI9kMJ1543lzoXFnFeko3sZXJ6FvJn5gF3A+4BKYB1wq3NuW2/PV8hLTw+t2M0Plu86eX94ZoB5xXnMn5TP/JI8WjqCrNvXwNr9Daw/cJTWzhAAs8bmcs+lk7ju/DEEfLEdQOacY3NVI79ZV8GzGw/R0hHkuvPH8KWrpzMxPyum+z4T9S0dVDS0kpsRICfdT256gPSAz+uyZIB5GfIXAw86566O3H8AwDn3f3t7vkJeTli7r4Fblq7mhgvG8sWrp/Pm3gbe3HeENXsbONjQevJ5ZjB9VA7zS/K4qCSPecV5jMxN96Tmlo4gS/+6l1/8dS9doTC3L5jI5y6fQn522qDWUdfcwUtbD/PCpmre3HeE8Dt+xVP9KeSmd4d+mj8l8uUjLZBCqi+F9FQfBVmpjMxNpzAn7eSJ6pG5aeRnpWJ2diOMnHPsrGnmhU3VvLy1hub2LnLSA2Sn+8lJ95Od5icnPUBmqo+AL4VUnxHwpRDwpxDwpZBi0BEM09EVpiMY6r4dDBEMOdL83XVnBCJfqT7SAz5Sfd2vDZx4r8jt1B5tPtH+VH/KyX0PhlDYsfNwM+sPHuWtA0eZPjqH+94z+azey8uQ/whwjXPuk5H7dwDznXOf7fGce4F7ASZMmHDhgQN9TossSeLo8U6ufehvpPlTeO4fL31XV0t1Y1vkBKif0uIRDM9M9ajS3tU2tfOjFbv57boKMgI+7nvPJO6+pITM1Oi6jA4da+Otg0fp6AqTHvB1B1vg7+G8qfIYz2+uZu2+BsIOJhVmcd15Y7igaDjHO4M0tQdpauuiqb2LprYgLR1BOk+EZVeYzlB3aLZ2hqhv7qCpPfiuGtL8KRSNyGB8XibjR2QyPi+DohGZjBmW3v2hkJ1Gqv/vIemcY1t1Ey9srubFzYfZW3+cFIN5JXmMG55JS0cXze3dtTS3d3+1dQbpCjs6g+E+/y0CPusOaX8KvhSjMxSmtTN0ytecqfRACtlpgR4fPH5S/Sl0BsN0BsN0hcJ0BLv/vZyD7DQ/uRl+ctIC5Gb4Ix+ggZM/l9TIB2lqpNZdkWDfePAYxyN/fRbmpHHzhUV86ZpzzqrmuA75nnQkL8457nm8jNd21bHsM4uG9OiV8toWvvPSDl7ZVkNOup+Fk/O5ZEoBi6YUUFKQdcojYucc++qPszbSHbV2XwOVR9tOu88pI7O59rwxXHfeGKaNyj7ro26A9q4Qdc0d1DZ3n6CuaWqn6lgblUfbqDjaSkVDG41tXe963YjMQOQvgHQqj7ay/0grvhTj4kn5vP+80Vw1czSFOaf/68Y5Ryjs6Ao5OkNhwmHXfXQeCcvehMKO9q4QbV0h2jpDdIXCdIVc5HuYYOTDozP0zr8IwnREXtMc+cBp6QjS0t5FS0eQjmD4ZGCn+rvDO+BPIcWMlvYumtqDNEc+PJvau052HfbGl2LMGJPD3AkjuHDiCOZOGEHRiIyoflanCvlYj5OvAsb3uF8UeUykV4+9sZ9Xt9fyfz4wc0gHPHQH7tIlpaw/0MDvyyr52+56Xt5aA8DYYeksmlLA+eOHc7wjSH1zB0eOd1Lf0kF9SyeHG9s42todoAXZqVxUnMcnLinhouI8ctL9dATDtHeFTn5v7wozMT+TaaNyBqz+9ICv+4g9L7PP5zS1d1HR0HpypFJtjw+F2uYOJuZncd97JvO+maP63W1lZvh9ht8HGZzZeQRfipGV5idrgE+091cw1P1BcuLo/8SRf1cozPgRmYNaX6yP5P10n3i9gu5wXwfc5pzb2tvzdSSf3LZUNXLTz95g8bQCfrGkNKojm3jknONgQyuvl9ezqryeVeVHTh4JpwdSKMhOIz87jcLsVApz0ji/aDjzSvKYdJqjfhHPjuSdc0Ez+yzwMt1DKB/tK+AlubV0BPnsU2+Rl5XKdz9yQUKGmpkxMT+LiflZ/MP8iYTCjsNN7QzPCHh+5CmJK+b/s5xzLwAvxHo/MnQ55/jqss0cbGjl1/csYERWfJ1IjRVfijFueIbXZUiC0yyU4rm/7KzjTxsP8Y9XTGX+pHyvyxFJKAp58ZRzjh+v2M244Rnc/94pXpcjknAU8uKp18vr2VhxjM+8d/KgXYQikkz0WyWecc7x0IrdjBmWzkcuLPK6HJGEpJAXz6zZ28C6/Uf51OJJpPk1n4pILCjkxTM/Wbmbguw0bpk3wetSRBKWQl48sf5AA2/sOcJ975mkWRFFYkghL554aEU5eVmp3DZfR/EisaSQl0H3dsUxXttVxycvjX5mRhE5NYW8DLqfrNzNsIwASy4u9roUkYSnkJdBtfVQI69ur+UTl5QM+JJ8IvJuCnkZVD9dWU5Omp+PLyz2uhSRpKCQl0Gz9VAjL245zJ2LihmWEfC6HJGkoJCXQbGx4hh3PLKW/KxU7l5U4nU5IklDIS8xt2J7DbcsXU12mp8/fHph0kwlLBIPdOZLYuqpNw/y1T9t5txxw3jk4xed0dqeIjJwFPISE845frh8Fw+tLOe90wv56W1ztfqRiAf0WycDrisU5oGnN/OH9ZV8rHQ83/zQufg1jbCIJxTyMuD+/YXt/GF9JZ+/cir/dMXUhFyvVWSoUMjLgGpq7+K36yr48NwiPn/lNK/LEUl6+htaBtQf11fS2hnirkXFXpciIijkZQCFw44nVh9gzoThnDtumNfliAgKeRlAq/bUs7f+OB/XxGMicUMhLwPmV28cID8rlfefN9rrUkQkQiEvA6KioZWVO2q4dd4ErdcqEkcU8jIgnnzzIIBWehKJMwp5iVp7V4jfrjvIVTNHM3Z4htfliEgPCnmJ2nObqjna2sWSiyd6XYqIvINCXqL2xOr9TBmZzcWT870uRUTeQSEvUdlYcYy3KxtZcvFETV8gEodiFvJm9qCZVZnZxsjXtbHal0SnMxjmy8s28/1Xdvb7tY+v3k92mp+b5hbFoDIRiVas5675oXPuezHeh0ShvSvE/U++xYodtQCMyk3n9gVn1rd+pKWD596u5pZ547Uot0icUndNEmvrDHHP42Ws2FHLN26cxeXnjOTBZ7eyZu+RM3r9b8sq6AyFdcJVJI7FOuQ/a2abzOxRMxvR2xPM7F4zKzOzsrq6uhiXIycc7why12Nreb28nu985HyWXFzMj26ZzYT8TD7z5FtUHm095etDYceTaw6ycHI+U0bmDFLVItJfUYW8mb1qZlt6+boR+E9gMjAbqAa+39t7OOeWOudKnXOlhYWF0ZQjZ6ipvYs7HnmTdfuP8qOPzeajpeMByE0P8MslpXSFwtzz+HpaO4O9vr69K8T3XtlJ1bE2lmieGpG4FlVHqnPuyjN5npn9Angumn3JwDjW2smSR9eyvbqJ/7htDtecO+Z/bJ9UmM1Pbp3D3Y+t43/9fhM/vW3OyVEzzjle2VbDvz2/jYqGNq6/YCxXzhjpRTNE5AzFcnRNz/T4ELAlVvuSM9MZDHP7I2+y43AzD99x4bsC/oTLpo/kX685h+c3V/Ozv+wBYHdNM3c8spZPPbGejICPpz45n5/cOkfL+onEuVgOifiOmc0GHLAf+FQM9yVn4LE39rGlqomf3z6Xy88Zdcrn3rt4Eturm/jeKzvZXt3Ei1sOk5Xq48HrZ3L7gokKd5EhImYh75y7I1bvLf1X29TOj1/dzRXnjOzzCL4nM+NbHz6fvfXHeX5zNbfNm8C/XDWdvKzUQahWRAaKBjcniW+9tIOukONrH5h5xq9JD/h46p4FNLR0MiE/M4bViUisKOSTwPoDDTz9VhWfuWwyxQVZ/XptdppfFzqJDGHqWE1wobDj689uZXRuOve/d4rX5YjIIFPIJ7jfrqtgS1UTX75uBlk6IhdJOgr5BHastZPvvryD+SV5XH/+6U+2ikjiUcgnsB8s30VjWxcP3jBL0wCLJCmFfILadqiJ/1pzgDsWTGTGmFyvyxERjyjkE5Bzjgef3crwzFS+8L7pXpcjIh5SyCegFzYfZu3+Br541XSGZQa8LkdEPKSQTzAdwRDffmkH54zO4WMXjfe6HBHxmEI+wTyx+gAHG1r58rUz8KXoZKtIslPIJ5BjrZ08tGI3i6cVsnia5uYXEYV8QnloRTktHUG+cu0Mr0sRkTihkE8Q++uP88Sa/Xy0dDzTR2s5PhHpppBPEN9+aQcBXwpfeN80r0sRkTiikE8AZfsbeHHLYT61eDIjc9O9LkdE4ohCfohzzvFvz29nVG4a9ywu8bocEYkzCvkh7rlN1WysOMa/XDWdzFTNMiki/5NCfgjreeHTh+cWeV2OiMQhhfwQ9vuySiqPtunCJxHpk0J+iHLO8fjq/Zw7LpdLpxZ4XY6IxCmF/BC1Zm8Du2paWLKgWHPFi0ifFPJD1BNr9jM8M8ANs8d6XYqIxDGF/BBU3djGy1tr+FjpeNIDPq/LEZE4ppAfgn795kHCznH7golelyIicU4hP8R0BsM8tbaCy6ePZHxeptfliEicU8gPMS9uqaa+pYMlC4u9LkVEhgCF/BDz+OoDFOdncukUDZsUkdNTyA8hW6oaWX/gKHdcXEyKLn4SkTOgkB9Cnlh9gIyAj49cqCkMROTMRBXyZnazmW01s7CZlb5j2wNmVm5mO83s6ujKlGOtnfxpYxUfnDOOYRkBr8sRkSEi2mkLtwA3AQ/3fNDMZgK3ALOAscCrZjbNOReKcn9J6/dllXQEwyy5WMMmReTMRXUk75zb7pzb2cumG4HfOOc6nHP7gHJgXjT7SmbhsOOJNQeYV5zHjDG5XpcjIkNIrPrkxwEVPe5XRh57FzO718zKzKysrq4uRuUMbU9vqOJgQytLFuooXkT657TdNWb2KjC6l01fcc49E20BzrmlwFKA0tJSF+37JZI1e4/ww+W7eHNfA1NHZnP1rN5+DCIifTttyDvnrjyL960Cxve4XxR5TM7Auv0N/HD5Lt7Yc4TCnDS+fv1Mbp03gYBPg6FEpH9itV7cs8BTZvYDuk+8TgXWxmhfCWN3TTPfeG4bf9tdT0F2Kl+9bga3L5ioSchE5KxFFfJm9iHgJ0Ah8LyZbXTOXe2c22pmvwO2AUHgfo2sOb0v/XETe+uO8+Vrz+H2BRO1ZquIRC2qFHHOLQOW9bHtm8A3o3n/ZNLY1sXbFcf47OVTuXfxZK/LEZEEoU7eOLFm7xHCDhZNzve6FBFJIAr5OPFGeT0ZAR9zJozwuhQRSSAK+Tjxenk980rySPXrRyIiA0eJEgcON7azp+44l2j6YBEZYAr5OLCqvB6AhVPUHy8iA0shHwdW7aknLyuVGaM1L42IDCyFvMecc6wqr2fh5HwtBCIiA04h77E9dcepaepgkfrjRSQGFPIeO9Efr5OuIhILCnmPrSqvZ3xeBuPzMr0uRUQSkELeQ8FQmNV7j+goXkRiRiHvoc1VjTS3B1k4WSEvIrGhkPfQG3uOALBQ89WISIwo5D30+u56ZozJJT87zetSRCRBKeQ90tYZYv2Bo1yiq1xFJIYU8h4pO9BAZyjMQp10FZEYUsh7ZFX5EQI+Y15xnteliEgCU8h7ZFV5PXMmjCArTUv8iUjsKOQ9cKy1ky2HGlmkoZMiEmMKeQ+s3nME5+CSqTrpKiKxpZD3wKo99WSl+ji/aLjXpYhIglOH8CBq6Qjy+7IKnttUzfxJ+QR8+owVkdhSyA+Cg0daeeyN/fyurIKWjiBzJwzni1dN97osEUkCCvkYWn+ggZ+/tpdXt9fgM+O688dw16ISZo8f7nVpIpIkFPIxUnm0lY89vIacdD/3XzaF2xdMZPSwdK/LEpEko5CPkRXbawmGHX/49EImF2Z7XY6IJCmd+YuRFTtqmVSQpYAXEU8p5GPgeEeQNXuOcPk5I70uRUSSnEI+Bl4vr6czFObyGQp5EfGWQj4GVm6vJSfNz0WafExEPBZVyJvZzWa21czCZlba4/FiM2szs42Rr59HX+rQEA47Vu6sZfH0Ql3sJCKei3Z0zRbgJuDhXrbtcc7NjvL9h5zNVY3UNXdwhfrjRSQORBXyzrntAGY2MNUkgBU7akkxuGy6Ql5EvBfL/oQSM9tgZq+Z2aV9PcnM7jWzMjMrq6uri2E5g2PljhrmThhBXlaq16WIiJw+5M3sVTPb0svXjad4WTUwwTk3B/gC8JSZ5fb2ROfcUudcqXOutLCw8OxaEScON7azpapJo2pEJG6ctrvGOXdlf9/UOdcBdERurzezPcA0oKzfFQ4hf95ZC8AV54zyuBIRkW4x6a4xs0Iz80VuTwKmAntjsa94smJ7LeOGZzBtlK5yFZH4EO0Qyg+ZWSVwMfC8mb0c2bQY2GRmG4E/APc55xqiqjTOtXeFWFVezxUzRupEtIjEjWhH1ywDlvXy+B+BP0bz3kPN6r1HaOsKaSoDEYkrulpngKzcXktGwMeCSVq3VUTih0J+ADjnWLmjlkumFpAe8HldjojISQr5AbDjcDNVx9p0lauIxB2F/ABYuaN76KT640Uk3ijkB8CK7TWcXzSMkbla3k9E4otCPkpHWjrYUHFMR/EiEpcU8lF6eWsNzsGVM3SVq4jEH4V8lJZtqGRyYRazxvY6NY+IiKcU8lGoaGhl3f6j3DS3SFe5ikhcUshHYdmGKgBunD3W40pERHqnkD9LzjmWbahifkkeRSMyvS5HRKRXCvmztLHiGPvqj3PT3HFelyIi0ieF/FlatqGKNH8K7z9vjNeliIj0SSF/FjqDYf777UNcOXMUuekBr8sREemTQv4svLarjqOtXdw0R101IhLfFPJn4em3KsnPSmXxtKG9Jq2IJD6FfD81tnaxYnst118wloBP/3wiEt+UUv30/OZqOkNhjaoRkSFBId9PJ6YxOG/cMK9LERE5LYV8P2gaAxEZahTy/aBpDERkqFHInyFNYyAiQ5FC/gxt0DQGIjIEKeTPwN66Fj731AZy0v2axkBEhhS/1wXEu02Vx7jz/63DgKc+uUDTGIjIkKKQP4VV5fXc+3gZwzNTeeIT85hUmO11SSIi/aKQ78MLm6v5/G82UlKQxeOfmMeo3HSvSxIR6TeFfC/+a80BvvbMFi6cMIJHPn4RwzLVRSMiQ5NC/h3+8y97+PZLO7j8nJH8x21zyUj1eV2SiMhZU8j38Ojr+/j2Szu4cfZYvnfzBZqATESGvKhSzMy+a2Y7zGyTmS0zs+E9tj1gZuVmttPMro660hj7zdqDfOO5bVwzazTfV8CLSIKINsmWA+c6584HdgEPAJjZTOAWYBZwDfAzM4vbfo9nNlbxwLLNXDa9kIdunYNfAS8iCSKqNHPOveKcC0burgGKIrdvBH7jnOtwzu0DyoF50ewrVl7eepgv/O5t5pfk8fPbLyTVr4AXkcQxkIl2N/Bi5PY4oKLHtsrIY+9iZveaWZmZldXV1Q1gOaf32q46PvfUBs4bN4xffvwi0gNx+8eGiMhZOe2JVzN7FRjdy6avOOeeiTznK0AQeLK/BTjnlgJLAUpLS11/X3+23tx7hE89Ucbkkdn86q55ZKfpHLSIJJ7TJptz7spTbTezO4EPAFc4506EdBUwvsfTiiKPxY0vL9vMmGEZPPGJeRoHLyIJK9rRNdcAXwJucM619tj0LHCLmaWZWQkwFVgbzb4GUnN7F3vqjvPhueMoyE7zuhwRkZiJto/ip0AasDyyUtIa59x9zrmtZvY7YBvd3Tj3O+dCUe5rwGyvbgZg1lgt4SciiS2qkHfOTTnFtm8C34zm/WNl66FGAGaOzfW4EhGR2ErK8YJbDzVRkJ3KyBx11YhIYkvKkN92qImZY4dpMW4RSXhJF/KdwTC7a5uZpa4aEUkCSRfyu2qa6Qo5hbyIJIWkC/lth5oAjawRkeSQdCG/9VAjWak+JuZlel2KiEjMJWHINzFjTC4pKTrpKiKJL6lCPhx2bK9uUn+8iCSNpAr5Aw2tHO8MqT9eRJJGUoW8rnQVkWSTZCHfhD/FmDoq2+tSREQGRdKF/NRROaT5tTiIiCSHpAl55xzbDjXqpKuIJJWkCfm65g7qWzoV8iKSVJIm5LfqSlcRSUJJFPLdI2tmjMnxuBIRkcGTRCHfxMT8THLStZ6riCSPpAp59ceLSLJJipBvau/iYEOr+uNFJOkkRchvj5x01ZWuIpJskiLkT46sGaOQF5HkkjQhX5CdxsjcdK9LEREZVEkS8rrSVUSSU8KHfEcwRHlti0JeRJJSwof8rsMtBMNOI2tEJCklRMg759hT19Lrtm3VmkNeRJKX3+sCBsLmqkZu+Okqzh2Xy01zirhh9lgKstOA7pOu2Wl+LdwtIkkpIY7kJ+Rl8vXrZ2IY33huG/P/fQV3P7aO5zYd4u3KRmaMydHC3SKSlBLiSH54Zip3LSrhrkUl7Kpp5um3qvjThipW7qgF4M6Fxd4WKCLiEXPOeV3DSaWlpa6srGxA3isUdqzZe4Tl22q4dd4Epo/W7JMikpjMbL1zrrS3bVEdyZvZd4HrgU5gD3CXc+6YmRUD24Gdkaeucc7dF82++suXYiyaUsCiKQWDuVsRkbgSbZ/8cuBc59z5wC7ggR7b9jjnZke+BjXgRUSkW1Qh75x7xTkXjNxdAxRFX5KIiAyUgRxdczfwYo/7JWa2wcxeM7NL+3qRmd1rZmVmVlZXVzeA5YiIyGn75M3sVWB0L5u+4px7JvKcrwBB4MnItmpggnPuiJldCPzJzGY555re+SbOuaXAUug+8Xp2zRARkd6cNuSdc1eearuZ3Ql8ALjCRYbqOOc6gI7I7fVmtgeYBgzM0BkRETkjUXXXmNk1wJeAG5xzrT0eLzQzX+T2JGAqsDeafYmISP9FezHUT4E0YLmZwd+HSi4GvmFmXUAYuM851xDlvkREpJ+iCnnn3JQ+Hv8j8Mdo3ltERKIXV1e8mlkdcMDrOgZAAVDvdREeULuTi9odPyY65wp72xBXIZ8ozKysr0uME5nanVzU7qEhIWahFBGR3inkRUQSmEI+NpZ6XYBH1O7konYPAeqTFxFJYDqSFxFJYAp5EZEEppCPkpk9ama1Zralx2OzzWyNmW2MzLA5z8saY6GPdl9gZqvNbLOZ/beZ5XpZ40Azs/Fm9mcz22ZmW83snyKP55nZcjPbHfk+wutaB9Ip2n1z5H7YzIbMkMIzdYp2f9fMdpjZJjNbZmbDPS71lNQnHyUzWwy0AI87586NPPYK8EPn3Itmdi3wJefcZR6WOeD6aPc64IvOudfM7G6gxDn3NS/rHEhmNgYY45x7y8xygPXAB4E7gQbn3LfM7H8DI5xz/+pdpQPrFO12dE9b8jDdP/eEmoDwFO0uAlY654Jm9m2AeP5560g+Ss65vwLvnJfHASeOYocBhwa1qEHQR7unAX+N3F4OfHhQi4ox51y1c+6tyO1mupe4HAfcCPwq8rRf0R0ECaOvdjvntjvndp761UPXKdo9pBZLinaCMund54GXzex7dH+QLvS2nEGzle7A+xNwMzDe02piKLKO8RzgTWCUc646sukwMMqrumLtHe1OGqdo993Abwe9oH7QkXxsfBr4Z+fceOCfgUc8rmew3A18xszWAzl0L/CecMwsm+4J+D7/zoVwImsqJGQf6Knancj6ancviyXFJYV8bHwceDpy+/dAwp147Y1zbodz7irn3IXAr4E9Xtc00MwsQPcv/JPOuRM/45pI/+2Jftxar+qLlT7anfD6anePxZL+4cRiSfFKIR8bh4D3RG5fDuz2sJZBY2YjI99TgK8CP/e2ooFl3YsmPAJsd879oMemZ+n+YCfy/ZnBri2WTtHuhNZXu/taLCleaXRNlMzs18BldE8/WgN8HdgJ/Jjucx7twGecc+u9qjEW+mh3NnB/5ClPAw/E+1FOf5jZJcDfgM10jyoB+DLd/bS/AybQPVX2RxNpkZxTtDsN+AlQCBwDNjrnrvaixlg4RbsforvtRyKPnVgsKS4p5EVEEpi6a0REEphCXkQkgSnkRUQSmEJeRCSBKeRFRBKYQl5EJIEp5EVEEtj/Bwt0DrmfLtTQAAAAAElFTkSuQmCC\n" + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "plt.plot(endlinearloss[5:55,0],endlinearlossdiff[:,1])" + ], + "metadata": { + "collapsed": false, + "pycharm": { + "name": "#%%\n" + } + } + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python 3", + "language": "python", + "name": "python3" + }, + "language_info": { + "codemirror_mode": { + "name": "ipython", + "version": 2 + }, + "file_extension": ".py", + "mimetype": "text/x-python", + "name": "python", + "nbconvert_exporter": "python", + "pygments_lexer": "ipython2", + "version": "2.7.6" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} \ No newline at end of file diff --git a/tempwithbright.txt b/tempwithbright.txt new file mode 100644 index 0000000..28437ce --- /dev/null +++ b/tempwithbright.txt @@ -0,0 +1,209 @@ +8,35.6181 +8.1,36.3403 +8.2,36.5208 +8.3,36.8819 +8.4,36.4931 +8.5,36.6319 +8.6,36.375 +8.7,36.2917 +8.8,36.2014 +8.9,36.6042 +9,36.4167 +9.1,36.8056 +9.2,36.6319 +9.3,36.7153 +9.4,37.1181 +9.5,37.1667 +9.6,37.1042 +9.7,37.5347 +9.8,37.0556 +9.9,37.4514 +10,38.0833 +10.1,38.7708 +10.21,39.6597 +10.31,41.8681 +10.41,43.7569 +10.51,47.1181 +10.61,49.8194 +10.71,51.6458 +10.81,55.0903 +10.91,58.1319 +11.01,61.1458 +11.11,63.8125 +11.21,66.7222 +11.31,69.7361 +11.41,72.2431 +11.51,76.6111 +11.61,78.6667 +11.71,80.8194 +11.81,83.5556 +11.91,86.5625 +12.01,89.375 +12.11,91.9861 +12.21,95.2917 +12.31,97.5069 +12.42,100.931 +12.52,104.368 +12.62,106.556 +12.72,109.382 +12.82,112.451 +12.92,114.896 +13.02,118.361 +13.12,120.278 +13.22,122.944 +13.32,125.778 +13.42,128.778 +13.52,130.681 +13.62,133.056 +13.72,135.743 +13.82,137.861 +13.92,140.236 +14.02,142.25 +14.12,144.5 +14.22,146.257 +14.32,148.229 +14.42,150.708 +14.52,151.924 +14.63,154.069 +14.73,154.826 +14.83,157.222 +14.93,159.569 +15.03,161.09 +15.13,163.986 +15.23,164.389 +15.33,165.771 +15.43,167.417 +15.53,168.924 +15.63,169.903 +15.73,172.757 +15.83,174.729 +15.93,175.688 +16.03,176.861 +16.13,178.965 +16.23,180.257 +16.33,181.139 +16.43,181.854 +16.53,183.604 +16.63,184.542 +16.73,185.91 +16.84,187.014 +16.94,187.639 +17.04,189.833 +17.14,191.229 +17.24,192.5 +17.34,193.111 +17.44,194.66 +17.54,194.646 +17.64,196.833 +17.74,198.132 +17.84,199.021 +17.94,200.215 +18.04,202.299 +18.14,202.882 +18.24,205.028 +18.34,205.938 +18.44,207.41 +18.54,208.708 +18.64,210.188 +18.74,211.972 +18.84,212.646 +18.94,213.278 +19.05,214.972 +19.15,215.729 +19.25,216.951 +19.35,218.535 +19.45,219.354 +19.55,220.181 +19.65,221.361 +19.75,221.889 +19.85,223.153 +19.95,223.778 +20.05,224.639 +20.15,224.694 +20.25,224.771 +20.35,225.028 +20.45,225.722 +20.55,225.632 +20.65,226.375 +20.75,226.639 +20.85,227.09 +20.95,227.646 +21.05,227.618 +21.15,227.118 +21.26,228.111 +21.36,228.424 +21.46,228.063 +21.56,228.681 +21.66,228.688 +21.76,228.889 +21.86,228.861 +21.96,228.806 +22.06,229.243 +22.16,229.472 +22.26,229.104 +22.36,229.66 +22.46,229.833 +22.56,230.083 +22.66,229.972 +22.76,229.986 +22.86,230.215 +22.96,230.361 +23.06,230.549 +23.16,230.833 +22.5,229.701 +22.4,229.521 +22.3,229.236 +22.2,229.09 +22.1,229.174 +22,228.556 +21.9,229.042 +21.8,228.569 +21.7,227.972 +21.6,228.618 +21.5,228.236 +21.4,227.715 +21.3,227.854 +21.2,227.813 +21.1,227.667 +21,227.639 +20.9,227.333 +20.8,226.819 +20.7,226.618 +20.6,226.167 +20.5,225.785 +20.39,224.854 +20.3,225.118 +20.2,224.375 +20.1,224.972 +20,224.389 +19.89,223.368 +19.8,222.917 +19.7,222.271 +19.6,221.958 +19.5,220.319 +19.39,218.785 +19.3,217.951 +19.2,216.118 +19.1,215.944 +19,214.083 +18.89,212.313 +18.8,212.764 +18.7,211.41 +18.6,209.042 +18.5,208.785 +18.39,207.861 +18.3,205.951 +18.2,204.722 +18.1,204.014 +18,201.861 +18,201.611 +17.89,200.313 +17.8,199.042 +17.7,198.514 +17.6,198.194 +17.6,197.361 +17.6,211.405 +17.6,211.41 +17.6,202.88 +17.6,203.222 +17.6,203.667 \ No newline at end of file