Skip to content

Instantly share code, notes, and snippets.

@kozo496
Last active June 1, 2022 14:59
Show Gist options
  • Save kozo496/968cc804554b90b45f1e3510a5268ec7 to your computer and use it in GitHub Desktop.
Save kozo496/968cc804554b90b45f1e3510a5268ec7 to your computer and use it in GitHub Desktop.
最初に閾値を超える値のindexを求める
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"id": "187eb089",
"metadata": {},
"outputs": [],
"source": [
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"plt.style.use('seaborn-whitegrid')"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "74dbc0cb",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAD0CAYAAACVbe2MAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAABFBUlEQVR4nO29eWBTZd72f2Vt0yRNuu87LbS0pbQFRCgIiorLDKNs7QguDD4zI+ooLy7zjOiPUdCZFx9fnUcdl0GnjlpBZxQddQQRZKfpQltStrZpm+5b2qTNfn5/pAmWli5pknNycn/+anJyzrly7vS67/v7vRcORVEUCAQCgcBauHQLIBAIBIJ7IUZPIBAILIcYPYFAILAcYvQEAoHAcojREwgEAsshRk8gEAgsh0+3gKtRKBR0SyAQCASvJC8vb8z3GWf0wLXF0olSqUR6ejrdMkbBVF0A0eYsTNXGVF0A0QaM30gmoRsCgUBgOcToCQQCgeUQoycQCASWQ4yeQCAQWA4xegKBQGA50zL6yspKbNiwYdT733//Pe6++26sW7cOn3zyCQDAarVi+/btWLduHTZs2ACVSjWdWxMIBAJhkjg9vPLtt9/GF198AZFINOJ9k8mEXbt2Yd++fRCJRCgsLMSyZctQXl4Oo9GIkpISVFRU4MUXX8Qbb7wx7S9AIBAIhPFx2ujj4+Px2muv4Yknnhjx/uXLlxEfHw+ZTAbANia+tLQUFRUVKCgoAADk5OSgurp6GrIJdGA0W/HBSRU+KW3C5U4tpP4CLE0Lw0PLUjAjXEq3PIKL6NIa8OYPl/F1dRva+vWIlvvj9qxoLI+y0C2N4CROG/0tt9yC5ubmUe9rtVpIpVf+6cViMbRaLbRaLSQSieN9Ho8Hs9kMPn+0BKVS6awst6HX631aV5fOjB2H2nGx24D0MD/8fFYg+vQWfFPVgi8rW/DQdSG4JTWQFm3OQLSNTVXbEF443I4BgxULYgOwKFaGhj4j/nr4MvaKuHhmwISZof60aBsPUp7j4/KZsRKJBDqdzvFap9NBKpWOet9qtY5p8gAYOcONqTPvPKGro1+P3751Ep1aC968Jw+3ZkY6jnVpDXispAKvHO9CaHgk7rkuwaPanIVoG83Jum5s//40YuQi7P1NHtIirjTYqtUa/GrPSTxzoAMf/GoB5sTJPa5vPEh5enhmbEpKClQqFfr6+mA0GlFaWoq5c+ciNzcXR44cAQBUVFQgLS3N1bcmuAGj2Yr/+kCB9n493n9g3giTB4BQiR/euTcfN84Kx/bPq3H8chdNSgnTobl3EL/+QIHYoACU/NfCESYPAJkxMvzp1mjIxQL86u+l6Bww0KSU4AwuM/r9+/ejpKQEAoEATz31FDZt2oT169fj7rvvRkREBFasWAGhUIj169dj165dePrpp111a4Ib2f3deZQ39uFPq7ORlxA85mf8+Dy8WjgXyWESPPJROXp0Rg+rJEwHi5XCwx+Vw2Kh8M7GfIRK/Mb8XJiYj7c35qN/yITHSipAtpv2HqYVuomNjXUMn7zzzjsd7y9fvhzLly8f8Vkul4sdO3ZM53YED3OupR/v/FiPwvlxuCM7etzPiv34+EvRXNzx6lHs+rcSf14zx0MqCdPlg5MqlDf24f+tz0FiqHjcz86KDMQzd2TgD/+qxj/L1bgrN9ZDKgnTgUyYIowJRVH4w7+qIBcJ8NStk4svzooMxOYlydiraIZC1eNmhQRX0KU14M/fnkdBaih+Nmf8ytxO0fx4zI2X44WvlOjXm9yskOAKiNETxuSgsgNljX3YdstMyAIEkz7v4eUzECrxw5++OU+69l7Amz9cxqDRjGfvnA0OhzOpc7hcDnb8LBPdOiP+drTezQoJroAYPWEUFEXh5e8uICEkAHfnTa1rHiDkY8uyFJyq70FFm95NCgmuoE2jR/FJFe7KjcWMcMnEJ/yErFgZbpkdgXd/rIdmkLTqmQ4xesIo/nOuHeda+/HI8lQIeFP/iRQuiEe0zB8fVJDwDZN58/BlWKwUHr0x1anzf3dTGgYMZrx7tM7Fygiuhhg9YRTv/liPuGARfp4zuZjt1fjxefhVQTLOdRhQ2dTnWnEEl6AZMuGT0ib8LCcaccEBTl0jPSoQN2dEoPikCnoTmTXLZIjRE0ZQrdbgdEMP7l2YCL4TrXk7a/JjIRJw8N7xBteJI7iMvaVNGDRa8MCipGld5/5FSegdNOGLihYXKSO4A2L0hBG8f7wBIgEPa/LjpnUdqb8AN8+Q4suzLejoJ7F6JmGxUnj/RAPmJQYhM0Y2rWtdlxyMmRFS7DneQJLvDIYYPcGBZtCEzytb8IvcGMhEkx9pcy3unCWDyUJhr2L0mkgE+jh2qQtNPUO49/rEaV+Lw+HgvkWJULb2Q6Hqnb44glsgRk9w8FVVK4xmK9bPm15r3k5MoADzk4KxT9FMWnsMYp+iGTKRACsyIlxyvTvnREMk4OHTMlKhMxVi9AQHn5U1IzVcgqxpdud/yuq8WNR36VDWSFp7TKBfb8K3NW342Zxo+PF5LrmmxI+PlVmR+LKylSRlGQoxegIAoKFLh1JVL+7KjZ30xJnJcHtWFAKEPOwtJa09JvDV2VYYzNYpz4+YiNW5sRgwmPFtTZtLr0twDcToCQCAz8rV4HCAX8yNcel1xX583JoZia+qWmEwk9Ye3XyqaMaMcAnmxLqu1wYA1yWHIEYuwmdlapdel+AaiNETQFEUPq9QY1FKKCJlrt9U4o7sKAzozTh2iSxhTCctfUMoVfXiF3NjXNprA2zLItyRHYVjl7rITFkGQoyegNq2Aai6B3F7dpRbrr94Rhik/nx8dZZ06+nkm2rb81951Z4CrmJlVhTMVgrfKdvdcn2C8xCjJ+Dr6jZwOHDZKIyrEfK5WJERge/OtcFotrrlHoSJ+aa6DbMipUgOm9q6NpNlTqwM0TJ/fF3V6pbrE5yHGD0B31a3YV5i8DU3nHAFt2VGoV9vxjGyAxUtdAzocUbVM2qHMFfC4XCwMisKP17sIssXMwynNh6xWq147rnncP78eQiFQjz//PNISLDtFdrZ2YnHH3/c8VmlUomtW7eisLAQq1atcmwcHhsbi127drngKxCmQ12nFufbB/DsnRluvU9BWiikfnz8+2wrls0Md+u9CKP5T007KApYmeme8Jyd27Ii8e7RehxUtuMXc8mmJEzBKaM/cOAAjEYjSkpKUFFRgRdffBFvvPEGACAsLAzFxcUAgPLycvzP//wP1q5dC4PBtsek/RiBGXwzPBzultnua+kBtoXObpgVjkPnO2C1UuByXZsMJIzPtzVtSA4VIy3CPWEbO3PjghAu9cOBcx3E6BmEU6EbhUKBgoICAEBOTg6qq6tHfYaiKPzxj3/Ec889Bx6Ph9raWgwNDeGBBx7Axo0bUVFRMS3hBNdwUNmBrBgZouUit99r+awwdGmNqFJr3H4vwhV0BjNO1fXgpowIl4+2uRoul4NlM8Nx5EInTBaSj2EKTrXotVotJJIrLQMejwez2Qw+/8rlvv/+e6SmpiI5ORkA4O/vj02bNmHNmjVoaGjA5s2b8c0334w4x45SqXRGllvR6/Ws0zVgsKC8sRfrsuRu+W5Xa4uCBRwAJUdrIMwZe6NxT8HU8gRcr+1kow5GixXJ/kPTuu5kdaVKjBgwmPHZkQpkR7q/AQH4Vnk6g1NGL5FIoNPpHK+tVusow/7iiy+wceNGx+ukpCQkJCSAw+EgKSkJcrkcnZ2diIoaHTNMT5/cHqWeRKlUsk7X/soWWCkVVi/KQHpCkIuVja0t92Q/qrqstD9LppYn4HptxbVVEAt5uGtJDoR858dfTFZXXLIZL/34H1weEmGdh56xL5XntVAoFNc85lSp5+bm4siRIwCAiooKpKWljfpMTU0NcnNzHa/37duHF198EQDQ3t4OrVaLsLAwZ25PcBGHzndAHiBATpzcY/dcPiscVWoNWbrYQ1AUhcPnO7E4NXRaJj8VJH58LEgKwfe1HR65H2FinCr5FStWQCgUYv369di1axeefvpp7N+/HyUlJQCAnp4eiMXiEfHA1atXY2BgAIWFhXjsscewc+fOMcM2BM9gtVI4cqETBalh4HkwMbp8lm3EzaHzxAQ8wcUOLdR9Q7jBwyOdbpgZhksdWjT1DHr0voSxccppuVwuduzYMeK9lJQUx9/BwcH4/PPPRxwXCoXYvXu3M7cjuIGaln50aY24Ic2zvapZkVJEBPrhyMUurJsX79F7+yI/DFeoN8z0bDkvnxWO579S4ocLndhwXYJH700YDZkw5aPYDWCJh42ew+Fg0YxQHL/UBauVrFHvbg7VdmJWpBRRMs8kRe0khYoRIxfh2EUyQY4JEKP3UQ5f6ERWjAxhUvfNhr0Wi2eEonfQhHOt/R6/ty8xZLSgVNWDpR6uzAF7hR6CE3XdsJAKnXaI0fsgOoMZFU19WJwaSsv9F82w3ZesZuleSlU9MFkoXD+DvnLWDJlQ00LmTdANMXofpFTVC7OVwsLkEFruHxHoj9RwCY4So3crxy93g8/lIN8NQ2cnw/UptgqGlDP9EKP3QY5f7oKAx0F+Ij0GANhae2caesjWc27kxOVu5MTJIfajZ3RbmNQPsyKlpOfGAIjR+yAnhw0gQEjf8NbFM0KhN1nJXrJuol9vwtnmPlyfQk+vzc71KaE409BLKnSaIUbvY/TrTahSa2gL29hZkBwMHpeDo2RUhls4U98DKwUsTKEnPm9ncWoIjGYrFCpSodMJMXof43QdMwxA6i9AdqwMp+p7aNXBVo5f7oaQz8XceDmtOuYnhYDP5ZDwDc0Qo/cxTtQxwwAAYH5SMM4295FuvRs4frkb+QlB8BfwaNUh8eNjdowMZxpIhU4nxOh9jBMMMQAAWJAUDJOFQnljH91SWEWPzghlaz/t8Xk7C5KCUdmkIRU6jRCj9yF6dUaca+2nPT5vJy8hGBwOcJqEb1zKqbpuAMBChhj9vMRgGC1WVDb10S3FZyFG70PYu8/XMcQAZCIB0iMDcbqhm24prOJMQy/8+FxkxcjplgIAmDc8jJdU6PRBjN6HKFX1QsjjIitGRrcUB/OTgqFQ9cJoJrsRuQqFqgc5cXKPLUs8EfIAIWZGSHGaxOlpgxm/BIJHKG3oQVasjBHxeTvzk4KhN1lRTabJu4RBoxnVLf2Yl0jvDl5XMz8pGGWqXpjJ9oK0QIzeR9CbLKhSa2ibDn8t7IZ0hnTrXUJFUx8sVgp5NM56Hot5ScHQGS1kITuaIEbvI5xt1sBkoZDPsJZemNQPyWFiEr91EaUNveBwgNx4Zhn9/OHfHSlneiBG7yOUqmz/YHkMa9EDtuF3pxt6yHK2LqBU1YuZEVLIRAK6pYwgUuaP+OAAYvQ04ZTRW61WbN++HevWrcOGDRugUqlGHN+zZw9uv/12bNiwARs2bEBdXd2E5xDci6KhF8lhYgSLhXRLGUV+QjAG9GZc6tDSLcWrsVgplKl6aV2sbjzmJwWjVNULiiIVuqdxyugPHDgAo9GIkpISbN261bHpt52amhq89NJLKC4uRnFxMZKTkyc8h+A+rFYKpapezEtgVtjGTu5wL4MscDY9zrcNQGswI5+h5ZyXEIQenREN3WQfWU/jlNErFAoUFBQAAHJyclBdXT3ieE1NDd566y0UFhbir3/966TOIbiPy51aaIZMjEvQ2UkMCUCwWIgysvDVtLCH55jaorfnDUg5ex6n1qnVarWQSCSO1zweD2azGXy+7XK33347ioqKIJFIsGXLFhw6dGjCc36KUql0RpZb0ev1Xqvr6wu2kQ5B5h4olZ4Lj0zlmc0I4uPkpXaPPWOmlifgvLbvz7YjJICH/tYGKNs4jNFlx0pRCBBwcLCyDhkBAy5Uxs7ydCVOGb1EIoFOp3O8tlqtDsOmKAr33nsvpFIpAGDp0qU4d+7cuOdcTXp6ujOy3IpSqfRaXW9XVSBELMRNC7LB4bjeAK7FVJ7Z0jYB/vzteUQlpEAe4P48AlPLE3Be24V/tWDhjHBkZGS4QZVrnlluwgAaBowuf/ZsLM+polAornnMqdBNbm4ujhw5AgCoqKhAWlqa45hWq8Udd9wBnU4HiqJw6tQpZGZmjnsOwb0oVL3ISwjyqMlPFftqmuVkPRSnaNUMoUWjZ9w8iavJjZejtq0fOoOZbik+hVMt+hUrVuDYsWNYv349KIrCzp07sX//fgwODmLdunV47LHHsHHjRgiFQixcuBBLly6F1WoddQ7B/XRrDVB1D6JofjzdUsZlTqwcXA5QrurFspnhdMvxOiqGVwCdy7Dx81czNyEIVgqobO5z7ClLcD9OGT2Xy8WOHTtGvJeSkuL4e9WqVVi1atWE5xDcT2VzHwAgJ05Oq46JEPvxMSsyEGVkyWKnqGjqg5DHRXpUIN1SxiU3zlYRlTcSo/ckZMIUy6lo0oDLATIZtJDZtchNkDum8BOmRnlTHzKiAxmzkNm1kAUIkBImRjkZSutRmP2rIEybiqY+pEVIIfajbyPwyZIbHwStwYyLHa4dkcF2zBYrqpo1jO+12ZkbH4Syxj4yccqDEKNnMRRFobKpz2sM4Mo46z56hXgZF9q1GDJZGLE95GTIjbdNnFKRiVMegxg9i2noHoRmyOQ1Rp9gnzhFuvVTomJ4pJK3lHNughwAUN5EytlTEKNnMRXD/0g5XtLS43A4yI2Xk/jtFKlo6kVQgADxwQF0S5kUqeFSSPz4pOfmQYjRs5jKJg0ChDykhkvpljJpsmPlqOvSYUBvoluK11DZpMGcODmj50n8FB6Xg+xYmaMnQnA/xOhZTHlTH7JiZOBxvcMAACA7VgaKAqrUZMepyaA1mHGhY8BrwjZ2smNtE6cMZgvdUnwCYvQsxWC2QNnS7zVhGzvZsXIAto1SCBNztrkPFOU98Xk7c2JlMFkoKFvJCCtPQIyepShbB2C0WJEzbJzeQrBYiPjgAJwdnuhFGB9vS8TayR7WS8rZMxCjZykVjd6ViP0p2bEyVDaRFv1kqGjsQ2JIgEcWgnMl0TJ/hEqEpJw9BDF6llLR1IdwqR8iA/3pljJl5sTKoe4bQpfWQLcURkNRFCq8aJ7ET+FwOMiOlZMWvYcgRs9SKodnSnrLSIyfkh1rW66BmMD4tPXr0TFg8EqjB2zlfKlTCy1ZydLtEKNnIX2DRtR36bwybAPY1uXhckC69RNgX7Eyh+ErVl6LObFyUBRQTUZYuR1i9CykcnjEirclYu2I/fiYES4hLfoJqGzWQMDjID3Ke+ZJ/BTSc/McxOhZSNXwP05mLPNXrLwWtvithix8NQ7Vag1mRkrhx+fRLcUpQiR+iJGLHA0TgvsgRs9CqtQaJIWKEegvoFuK08yJlaFbZ4S6b4huKYyEoihUqTXI8oLlp8djTpyMtOg9ADF6FlKt7veK9efHY45jnDVp7Y1FU88QNEMmry/n7Fg5mnqG0KMz0i2F1Ti1SLnVasVzzz2H8+fPQygU4vnnn0dCQoLj+Jdffon3338fPB4PaWlpeO6558DlcrFq1SrHpuGxsbHYtWuXa74FAfjTn4B589AzfxHUfUO49/oE4NAh4MwZ4Ikn6FY3ZWZFBkLI46KyuQ+3ZUXRLYdx2JeIyI6R0ytkmvw0Tn8D2ULSbThl9AcOHIDRaERJSQkqKirw4osv4o033gAA6PV6vPLKK9i/fz9EIhEef/xxHDp0CIsXLwYAFBcXu0494Qrz5gFr16Jp91sAhLi+sQp4bDPwySd0K3MKIZ+L9CgpzpKRN2NSpbYlYtMiJXRLmRZZMTJwOLaeGzF69+FU6EahUKCgoAAAkJOTg+rqascxoVCIjz/+GCKRCABgNpvh5+eH2tpaDA0N4YEHHsDGjRtRUVExffWEKyxbBnzyCdIe2YTHfvwAGXaTX7aMbmVOkx0rR5VaAyvZWnAU3p6ItSP1FyA5VEzi9G7GqRa9VquFRHKlJcHj8WA2m8Hn88HlchEaatv0t7i4GIODg1i0aBEuXLiATZs2Yc2aNWhoaMDmzZvxzTffgM8fLUGpVDr5ddyHXq9nvq7ISCivvwOPfl2Mzt/8Bl2RkQCNmqf7zMJ4g9AazPju1FnEy107xZ+p5QlMrI2iKFQ09mBxgtij38FdzyxRyoGioRvnzp1zeoKfN5enJ3DK6CUSCXQ6neO11WodYdhWqxV//vOfUV9fj9deew0cDgdJSUlISEhw/C2Xy9HZ2YmoqNHx1/T0dGdkuRWlUsl8XYcOIfLIF/j6F7/Cyr17EbZmDa0t+uk+M05QP14+1okh/1Ckp8e4UBlzyxOYWFtj9yC0xnoUZCYgPT3hmp/ztC5nKejxx8G6cwiOSUakzLklO7y5PF2FQqG45jGnQje5ubk4cuQIAKCiogJpaWkjjm/fvh0GgwGvv/66I4Szb98+vPjiiwCA9vZ2aLVahIWFOXN7wlgcOgTrmrX47Z1PovHRp2xhm7VrbQlZL2VGmAT+Ai5Zm/4q2JKItWMfOURmyLoPp1r0K1aswLFjx7B+/XpQFIWdO3di//79GBwcRGZmJvbt24f8/Hzce++9AICNGzdi9erVePrpp1FYWAgOh4OdO3eOGbYhOMmZM6h6+S2cOCfEwzEyYIYtZo8zZ7w2Ts/ncZEeFUgM4CrYkoi1kxEdCA7H9r1uyoigWw4rccppuVwuduzYMeK9lJQUx9+1tbVjnrd7925nbkeYDE88gaOHLgHnzmO2fWz1smVea/J2MqNl+Ge5GlYrBa4X7ZTlTtiSiLUTIOQjJUyCmhZSobsLMmGKRVSrNUgICYBM5L0zYq8mK0YGrcEMVc8g3VIYAVtmxF5NVoyMhOjcCDF6FlGl1nj9TMmrsX8fYgI27DNis1gSn7eTGSNDe78BHQN6uqWwEmL0LKFXZ0Rz7xDrWnqpERII+VwSpx/GXuGxrZwzowMBADXqfpqVsBNi9CyhuoWdBiDgcZEeKSVGPwzbErF2ZpORN26FGD1LsLf0MqPZZfSAzQSq1WTJYoB9iVg7Ej8+kkPFJETnJojRs4RqtQbxwQGQBbAnEWsnK0aGfr0ZjT6ekGVrItZO5nCFTnA9xOhZApsNIMvRrfft+C1bE7F2smJkaNHo0U02hXc5xOhZQN+gEU09Q6wbcWMnNUICAY/j8916tiZi7cyOsSVkq1t8u0J3B8ToWYC9pctWA/Dj8zAzUurzE2rYmoi1MzuaJGTdBTF6FuBIxA63iNiIfUKNLydk2ZqItSMTCZAQEkCM3g0Qo2cB1WoN4oJFkAe4dilfJjE7Woa+QROae31zD1m2J2LtZJIZsm6BGD0L8AUDsH8/Xw3fsD0RayczWobm3iH0DZI9ZF0JMXovZ8BgQWPPIGsTsXZmRkrB5/puQpbtiVg7ZISVeyBG7+Vc6rG1fNhuAP4CHlIjpKjyUQNgeyLWTqZj5I1vVujughi9l3Op2zbmmI0zYq8mKyYQNT6akGV7ItaOPECI2CCRz/bc3AUxei/nUrcBsUEiBInZm4i1kxkjQ7fOiFaNb61wSFEUqlvYn4exkxUjQw0xepfilNFbrVZs374d69atw4YNG6BSqUYc//7773H33Xdj3bp1+OSTTyZ1DsE5LnYbfMYAfHXLOVty0sT6PIydzBgZGroH0a830S2FNThl9AcOHIDRaERJSQm2bt3q2AsWAEwmE3bt2oW//e1vKC4uRklJCTo7O8c9h+AcmiETWgfMPmMAGVGB4HE5Pmf01T6SiLXjqxW6O3HK6BUKBQoKCgAAOTk5qK6udhy7fPky4uPjIZPJIBQKkZeXh9LS0nHPcRVag9nl12QyNT5mAP4CHmaESXwuflul1oDP5WBmpJRuKR7BV9em1xnMbss/OWX0Wq0WEsmV7D+Px4PZbHYck0qv/CDFYjG0Wu2457gChaoXuTu+Q2O376xw6CtD7n5KZozM59ZCqVJrkBbB/kSsnRCJH6Jl/j5VoVusFJbv/gF/P+GekLZTm4NLJBLodDrHa6vVCj6fP+YxnU4HqVQ67jlXo1Qqp6yps9cIo8WKL05U48YU17d89Hq9U7rcyTFlO8ICeGhrvIw2usWMgTueWRh/CJ0DBhxVVCEkwKmfLwBmlqedn2qjKAqVjT1YGC+mXa8nn1lCIBdl9Z2Tvp+3lOe1UPUa0d5vwEBPB5RK18/+duo/JTc3F4cOHcJtt92GiooKpKWlOY6lpKRApVKhr68PAQEBKC0txaZNm8DhcK55ztWkp6dPWVOalYLo61Z0U2Knzp8IpVLplutOh8Yv25AW6sc4XXbc8cx0oh68efoE9KJwpKdHOH0dJpannZ9qU/cNod9Qj8UZ8UhPT2SMLnezsIWPkwcuIC45FRK/iW3KW8rzWpxTNANoxi3z0pEa4VxDVaFQXPOYU0a/YsUKHDt2DOvXrwdFUdi5cyf279+PwcFBrFu3Dk899RQ2bdoEiqJw9913IyIiYsxzXAmPy8Hs6ECfSeD0601o6B7E0vgguqV4lIzoQHA5tgk1N2U4b/TeQlWzfcE63wnPAbZwJEUB51r6MT8pmG45bqe6RQORgIfkMPdMiHPK6LlcLnbs2DHivZSUFMffy5cvx/Llyyc8x9VkxshQcqYJFisFHpfj1nvRjT1RlRriR7MSzxIg5CMlTOIzFXpNiwY8LgfpUexdmXQs7BVblVrjG0av1iAjOtBtvsWqCVNZMTIMmSyo69TSLcXt2I0uxceMHriyZLEvUKXWIDVcAn+BbyRi7YRJ/RAR6OcTFbrFSqGmpd+tgyrYZfSxV1oBbKdKrUG0zB9yf98yAMDW2mvvN6Cjn90zZCmKQrVa43NhGzu+UqHXd+kwaLS4tZxZZfQpYRKIBDyf+HH4tAH4SIXe3m9Al9boGFfua2TGyHC5Uwsdy+fH2Hst7tw4iFVGz+NykOEDCdkBvQl1XTqfGj//UzKiAsHhsN/oHfMkYn2znB0J2VZ2z5uoUmvgL+BihpsSsQDLjB4YXhCppR8WK3tXOKwZnjCU6aMGIPbzjYRstVoDLgc+l4i1Y2/I2EcesZVqtQbpUYHg89xnx6wz+swYGQaNFtR3sTch6+jq+cDSxNfCF+K31WoNUsIkCBA6PzHMmwkP9Gd9QtY6nIh19/8y64w+K4b98dtqtQaRgf4Ik/reiBs7WT6QkPWlpYmvRVaMDGdZ/L+s6hmE1mB2ezmzzuhTwsTwF3BR1czeuF6VDydi7bA9IdsxoEd7vwGzfbyc2Z6QrVJ7ZkIc64yez+MiI4q9CVmtwezTiVg7bE/I+trSxNciO5bdCdlqtQZCPhepEe7dIpJ1Rg/YE7IaWFmYkD3X0g+KArJifTNBZ4ftCdlqdT84HNuSD76MvaV7lqUJ2Wq1BumRUgjcmIgFWGr0mTEy6IwW1HXpJv6wl1Htoa6eN8DmhGyVWoOkUPGkFvRiM+FS9iZkPTkhjpVGb4/fsvHHUa3WIFzqh3CpP91SaIfNM2Sr1SQRa4etFXpjzyD69Z7ZIY6VRj8jTGJLyLLwx1FFDMBBNksTsn1DFrRq9D49fPanZMXIcblTy7od5KqHFyb0xP8zK42ez+MiPSqQdQYwaDTjcqeWhG2GYWtC9lKPAQAJz9nJig10LFnMJqrUGgh4HLcnYgGWGj1gqyXPtfSzKiF7rqUfVoqMxLDD1oTspW6b0c9249on3kQmS+fGVKs1mBnpmS0iWWv0mTEyaA1m1HezJyHr62ufjAUb47eXug1IDAlAoL+AbimMIFzqj8hAf1Q199EtxWVQFOXRCXGsNXr7A2RTa69a3Y9QiR/CfXhG7NWwMSF7sdtIwjZXkcmyCr25dwh9gybM9lAexqmxW3q9Htu2bUN3dzfEYjFeeuklBAeP3AXmvffew1dffQUAWLp0KbZs2QKKorBkyRIkJiYCAHJycrB169bpfYNrkBougR+fi6pmDX6eE+OWe3ga20iMQHA47N49ayr8dMmLGwO9fyRSr86IDp37p8R7G1kxMhysbYfWYGbFkFNPT4hz6ol99NFHSEtLw8MPP4yvvvoKr7/+Ov7whz84jjc1NeGLL77A3r17weFwUFRUhJtuugkikQizZ8/Gm2++6bIvcC3YlpAdMlpwsWMAt8xm/z6pU2F29JWE7I3T2CycKVS3kHkSY2FPyNaoNViQHEK3nGlTpdaAz+VgZqRzG4FPFadCNwqFAgUFBQCAJUuW4MSJEyOOR0ZG4p133gGPxwOXy4XZbIafnx9qamrQ3t6ODRs2YPPmzairq5v+NxgH+5LFbEjInmu1JWKJAYyEbQlZ+wxQMrRyJGxLyFYNJ2I9tUXkhC36vXv34v333x/xXkhICKRSW00kFosxMDAw4rhAIEBwcDAoisKf/vQnZGRkICkpCV1dXXjwwQexcuVKlJaWYtu2bfj0009H3VOpVE7nO13RydVBazDj4OmziJUJp3UtvV7vMl3OcEBp+4H7D3VCqex1vE+3rvHwlLZ4CVDe0D2lezH1uR1TtiFKwkOL6hJa6BZzFXQ/s5AAHo4pm7Ao1DDqGN3axuNqbRRFoVzVg4JEscc0T2j0a9aswZo1a0a8t2XLFuh0ttEsOp0OgYGjh4EZDAb8/ve/h1gsxrPPPgsAyMzMBI9nq8Hy8/PR3t4OiqJGxZzT09Od+zZXQcn68crxHzHoH4r09OnF6ZVKpct0OcO71ZUIEQtRkJc54nnRrWs8PKVtUbc/vq87h5CYJIRPMk7P1OdW/68WzAzzZ6Q2up/Z3AQd6ru0Y2qgW9t4XK2toUsHrbEeSzMTkZ4e77L7KBSKax5zKnSTm5uLw4cPAwCOHDmCvLy8EccpisJvf/tbzJw5Ezt27HCY+1/+8hdH76C2thbR0dFuTSymRkgg5HNZ0a0/29yHOXFykogdA7bsQdAxoEerRo+0EDKqaiyyY2Wo69J5/QzZyuFhotmxco/d06lkbGFhIZ588kkUFhZCIBBg9+7dAIA9e/YgPj4eVqsVp0+fhtFoxI8//ggAePzxx/Hggw9i27ZtOHz4MHg8Hnbt2uW6bzIGApYkZLUGMy52aHFbVhTdUhgJWxKyZ5tsv9PUUO8fPeQO7HvIentC9myzbY/YNA/MiLXjlNGLRCK8+uqro96///77HX9XVVWNee5bb73lzC2dJismEJ+Xt8BqpcDlemdruFqtAUUBc+LkdEthJGxJyJ5t7gOXA8wInl4+ia38NCHr3Ubfh9nRMrfuEXs1rJ0wZScrRoYBgxmqnkG6pThNZVMfAGCOB7t63gYbZshWNmuQFiGFv4D1/5ZOESb1s82Q9eJyNlusqFb3Oxbk8xSs/0VlxcgB2GpRb+VsswZxwSIEi0lL71p4+wxZiqJwtrnP4wbgbWTFeneFfqlTiyGThRi9q0mLkEAk4KG8sY9uKU5T0dTn0cSNNzJn+B+nYrj342009w6hd9BEynkCsmNkqO/SoV9voluKU9jzMJ4uZ9YbPZ/HRVaMzJHp9ja6tAao+4aQQwxgXDJjZOBzOV5bznbdJDw3PnPi5KAooMpLtxasbO6D1I+PpBCxR+/LeqMHgJx4OWpa+mE0W+mWMmXOOoZikS79ePgLeJgVJfXaFv3ZZg2EPK7HpsR7K/YBCd5czlmxMo8PDPEJo58TK4fRbEVtm/dtXFDZpAGXQ5Y+mAw5cXKcbfLOTeErm/qQER0IId8n/iWdRiYSIDlM7JWhWIPZgtq2flrCcz7xq8qJlwPwzlZAZXMfUsOlELNgxT53MydWjgGDGXVdWrqlTAmL1bZJ9BzSa5sUOXFyVDT1gaK8q0JXtg7AZKFoKWefMPpomT9CJX5eZ/S2kRgazIkjBjAZ5g5X6N7W2qvr1EJntJBE7CTJiZM7clfehCMMS8N8GJ8weg6H42gFeBPNvUPo0RmJAUyS5FAJpP58ryvnyuHEIqnQJ0eOl8bpK5s0CJUIES3z/MxnnzB6AMiJk6GuUwfNoPcMy7KPxMghM2InBZfLwZxY76vQK5v6IPHjIznUc1PivZlZkbZcRoWX9dxs8yToWa/Kh4w+CABwVt1Hr5ApUNnUByGfjMSYCnPiZKhtG4DeZKFbyqQ529yHzJhAr12iw9MI+VxkRgd61VDafr0Jlzq1tDXafMbo7Rtqe1MroLJZg9nRgRB4cE0MbycnLsiR3PQG9CYLalr6MTc+iG4pXkVOXBCq1BqYLN4xZLqyqQ8UBeTSVM4+4yAykQApYWKvaQXY1sTQkAk0U8Tb4rdVag3MVoo2A/BWcuLl0JusON82MPGHGUCZqg8cDn15GJ8xesA22cJbhmXVtg1g0GhxjCQhTI4wqR9i5CKvMfoylW23MFLOU2Oul1XoZY29SAuXQuovoOX+PmX0c+Pk6NIavWJYVlmjzQDyEkhLb6p40wirssZeJIQEIFRCNhuZCrFBIoSIhV5RzlaKQkVTH3IT5LRp8Cmj96bp02WqXoQPt04JUyMnTo7m3iF0aUfvLcokKIpCWWMfCds4gTcNmVb3m6AZMmFuHH3l7FNG703DshSNvchLCCJbBzqBfSZ0JcNNQN03hM4BA3JJ2MYpcuLkuNypZfxKlsoOW4ODzha9U/Pq9Xo9tm3bhu7ubojFYrz00ksIDg4e8Znnn38eZWVlEIttq7S9/vrrEAgEE57nToR8LrJjZFAMh0WYSseAHk09Q9h4XSLdUrySzGjbSpZljb2M3lqwbLjBQUbcOEdOvG0ly7NNGjB5v6naLj0C/emdJ+FUi/6jjz5CWloaPvzwQ6xatQqvv/76qM/U1NTgnXfeQXFxMYqLiyGVSid1nrvJSwxCtVrD6HHWZao+AEAuic87hUjIw+zoQJQ2MLtCL1P1QiTgYRaZJ+EUc+Lk4HCu5LOYirLTgJz4IFrnSThl9AqFAgUFBQCAJUuW4MSJEyOOW61WqFQqbN++HevXr8e+ffsmdZ4nyE8IhsliW0OGqZQ39kLI4yIzJpBuKV5LXkIwKpv7GL00dXljL7JjPbt3KJsI9BdgZoQUZxp66JZyTQb0Jqh6jbSH5yYM3ezduxfvv//+iPdCQkIgldpaIWKxGAMDI8eyDg4O4p577sH9998Pi8WCjRs3IjMzE1qtdtzz7CiVSqe+zGSQGmwt+a/PnIdUL5/0eXq93q26fsqPtWqkBAtQd/HChJ/1pK6pQqe2SL4OepMVXx2vxKyw0WuL0P3cDGbbPIm7Z8tH6aBb27Vgoq4ZMg4O1nVDlyNgnDYAKG8dAgUgFAO06pvQ6NesWYM1a9aMeG/Lli3Q6XQAAJ1Oh8DAkS1PkUiEjRs3QiSyjRi57rrrUFtbC4lEMu55dtLT06f+TaZAyvddUA3yp3QfpVLpdl0AYDRbcekfDdh4XcKk7ucpXc5Ap7aQGD12Hj6Ibo4M6enJo47T/dzONPTAQjVgRe4MpF+VR6Bb27Vgoq4VBjW+PF+BNj0Xd+QxSxsAHGi5CKAVdy7Khkzk3jH0CoXimsec6jPm5ubi8OHDAIAjR44gLy9vxPGGhgYUFRXBYrHAZDKhrKwMs2fPnvA8T5GfEAyFqpeRG1TUtGhgNFvJ+PlpEh7oj/jgAMbG6clEKdeQn2gbzFHTzsxN4UtVvYiXCdxu8hPh1KibwsJCPPnkkygsLIRAIMDu3bsBAHv27EF8fDxuvPFG3HnnnVi7di0EAgF+/vOfIzU1FbGxsWOe52nyEoNQUtqEy51apEYwKxFmH4lBErHTJz8xCEcudIKiKMYNUz3T0IukUDGZKDVNYuQiRMv8UdPBPKO3WCmUqXpRkED/XBinjF4kEuHVV18d9f7999/v+Hvz5s3YvHnzpM7zNPnDJlqq6mWc0StUPYiRixAR6Pk1q9lGfkIwPitTQ9U9iMRQz27GPB5WK4UzDT24dXYk3VJYQV5iMI5fbGdcha5s7ceAwYzMCPqN3ifT/UmhYoSIhYzL1lMUhdP1PZif5Lm5BWwmP9FWoTOtnC90DEAzZCLl7CLmJQahe9DCuKVNTtfbfneZEfQ32nzS6DkcDvISgqBQMSt+e7lThy6tEQuIAbiEGWESyEQCxpXzmWEDIEbvGuz5LKblY8402HrnYWL693v2SaMHbK09VfcgOgeYsx7KaWIALoXLtVXoTGvRn6rvQZTMH7FB9Hfp2cCsyECIBBxGlTNF2cJzTGm0+azR5yXYCqCUQT+OU/XdCJP6IYlB8WRvJz8xCJc7dejRGemWAmBkeI5J8WRvhsflICPMn1E9t7ouW+98HjF6esmKkUEk4OFkXTfdUgDYDOBUna0FQAzAdcwbHn5n7y3Rjap7EB0DBocugmvICPfH+fYBxuwJbQ/PMaWcfdbohXwu8hODcIIhRt/UM4S2fj1junpsYU6snFEV+unhHiQpZ9eSHekPirL1ipnA6foehIiFSAljRu/cZ40eABamhOBCu5YRcfqTwz/QBclMXofP+7BX6Mcvd9EtBYDNAILFQswIp28lQzYyM9Qf/gIujl9miNE39GBeInN65z5t9NenhAIAI1p7p+psBpBKDMDlXJ8SypgK/XR9D/LJPgMuR8DjYF5iME4wwOjVfUNo7h1iTHwe8HGjz4wOhMSPz4jwzan6bsxLJAbgDham2HpJdFfoTT2DaOwZxHWk1+YWrk8Jxfn2Ador9GOXbL3HRTOYU84+bfR8Hhfzk4JxkuZWgL0FsCCJOT8MNpEZHQgpAyp0e/ioIDWUVh1s5XqGVOjHLnUhVCLETAbNuvdpoweAhckhqOvSoU1D31oZRy92AgAWzSAG4A7sFTrd3fqjl7oRLvUj8Xk3MTs6EFJ/Pq1xeoqicOxSN65PCWVU75wY/XAr4EQdfcm6Hy92IVzqh7QIYgDuYmFKCOq7dGjV0DNN3mqlcPxSFxbNYJYBsAk+j4sFSSE4QWPi/UK7Fl1aAxYzrNHm80afHhUImUhAW2vPaqVw7FIXFqcSA3An9sQ7XeV8vn0A3Toj6bW5metTQtDQPUjbujeO+DzDwnM+b/Q8LgfXJQfj2KVuUJTn16evaelH76CJxG3dzKxIKYICBDh2iR6jZ2KCjo1cP/x8j1+ip1V/7FIXEkMCECNn1vIWPm/0ALA0LRzqviFc6tB6/N4/XiLxeU/A5XKwaEYojlzspGXDmaOXupASJkaUjFkGwDbSwqUIlQhxlAajN1msOFXfw8j/ZWL0AG6YGQYA+OF8p8fvffRiF2ZFShEupX8pU7Zzw8xwdA4YcK6136P3NZqtOFXHTANgG1wuB0vSwnD4QicsHq7QK5v6oDWYGVnOxOgBRMtFSIuQ4IcLHR6975DRgtKGXhK28RBL0+wVumfL+XR9D4ZMFhSkhnn0vr7Kspnh6Bs0obK5z6P3/b62A7zhniPTcGqhZL1ej23btqG7uxtisRgvvfQSgoOvzAJTKpXYuXOn43VFRQX+93//FwUFBViyZAkSExMBADk5Odi6dev0voGLWJoWhvePq6AzmCH288z60Sfru2G0WLGYGIBHCJP6IStGhh/Od+LGKM9t1fh9bQeEfC6Jz3uIJalh4HKAH2o7kBvvuXI+dL4T+QlBtO8POxZOteg/+ugjpKWl4cMPP8SqVavw+uuvjzienp6O4uJiFBcXo6ioCDfffDOWLFmCxsZGzJ4923GMKSYP2Lr1RovVo2NwDyrbESDkkQWuPMgNM8NQ1tiLAYPFY/c8dL4DC5NDECCkfwMKX0AWIEBufBAOeTAU26oZgrK1H8tmhXvsnlPBKaNXKBQoKCgAACxZsgQnTpwY83ODg4N47bXX8N///d8AgJqaGrS3t2PDhg3YvHkz6urqnJTtevITgxAg5HmsW09RFA4qO1CQGgp/Ac8j9yTYjN5KAWUtnhl+V9+lQ32XDssZagBsZdmscFSpNegY8MxESHt+j6nlPGETY+/evXj//fdHvBcSEgKp1Da9VywWY2BgYMxz9+3bh1tvvdUR1gkLC8ODDz6IlStXorS0FNu2bcOnn3466jylUjnlL+IKsiP88F11C345kzdqTLter3eprkvdBrRq9Fg/Wzqt67palythojY/KwWJkIuTjQNY6gFt/zynAQDE8/sn/SyY+NwA5uoCRmtLENrWu/n4h7NYMcP9SxF8caYN4WI+zN1NUPa41zucYUKjX7NmDdasWTPivS1btkCn0wEAdDodAgMDxzx3//79ePXVVx2vMzMzwePZWq/5+flobx975/b09PSpfQsX8QutGE9+WgVKHoOMaNmIY0ql0qW6/nPgIjgc4JfL5yBU4uf0dVyty5UwVdvydAN+qG1DatpM8HnuHY/wwrFTmBEuwbL52ZM+h6nPjam6gNHaZlEUdvzQBaWGh0fcrFlvsqDyIxXuyo1BRkbGhNrchUKhuOYxp37lubm5OHz4MADgyJEjyMvLG/WZgYEBGI1GREVFOd77y1/+4ugd1NbWIjo6mlGzQW9KjwCXA3xb3eb2ex2sbUdOnHxaJk9wjlszI9FvsDo2AXEXmkETTtZ148Z0Znbn2QyHw8HNsyPww/lODBndm485erELg0YLVmREuvU+08Epoy8sLMTFixdRWFiIkpISbNmyBQCwZ88eHDx4EABQX1+PmJiYEec9+OCDOHPmDO655x7s2rULu3btmqZ81xIi8cP8pGB8U+Neo2/VDOFsswY3pUe49T6EsblhZhj8eBy3V+j/OdcGs5XCbZlRE3+Y4HJWZkZhyGTB4QvuTcr+u7oVMpHAsXomE3FqGIBIJBoRkrFz//33O/7Ozs4eNRpHJpPhrbfecuaWHuPW2ZF4bv85XO7UIiXMPYuMfXW2FQBwWxYxADoIEPKRFyPCNzVtePbO2eBy3dOr/Lq6DTFyEbJjZRN/mOByFiQFIyhAgK+rW3Frpnta20azFd+da8ctsyMhcHMYcDowVxlN3Dzb9oP41o2t+i/PtmJ2dCCSQpmxn6QvsihejPZ+A8qb+txyfc2QCT9e7MTKzEhGhSd9CT6Pi5szInFQ2QGD2T3hm2OXuzCgN2OlmyoSV0GM/iqi5SLMiZM7Wt2upqlnEBVNfbgjO9ot1ydMjvlxARDwOPi6yj3lfFDZDpOFwm3ZpNdGJyuzIqE1mHH0onvWvvm6qhUSPz4WM3x2OzH6MViVE42aln5caB972Oh0+GrYWO4gBkArEiEPS9PC8EVlC8wWq8uv/9XZVkTJ/JETK3f5tQmT5/qUUMhEAnxe0eLyaxvMFnxb046b0sPhx2f2XBhi9GNw55xo8LgcfFamdvm1vzzbgjlxcsQFB7j82oSpcXduLDoGDDjm4tnQnQMG/HChEz/LiXZb/J8wOYR8Ln6eE41va9rQrze59NoHlR3QDJnwi9xYl17XHRCjH4NQiR9uSAvDv8rVLl0B71xLP6rV/ViVQ8I2TGB5ejhkIgE+VTS79LqfV9h+N6u9wAB8gbtzY2EwW10ejv1U0YyIQD/G7SY1FsTor8FdubFo69e7dEeiT0qbIORxsSonZuIPE9yOH5+HO+dEubS1R1EU9pY2Y06cHKkM2hzal8mOlSE1XIJ9LqzQ7b22VXNjwPOCXhsx+mtwY3o4Av35+PhMo0uupzdZ8M9yNW7JjESQWOiSaxKmj72192Wla1p7VWoNzrcPYHUeac0zBQ6Hg9V5sVCoenG50zWbC31W1uxVvTZi9NfAX8DDmvw4fFPdhvb+6S+M9E11GzRDJqzLj3OBOoKryImTY1akFH8/0eCSrSTfP65CgJCHn80h4Tkm8Yu5MeBzOfjgpGra17JYKRSfVGFBUrDX9NqI0Y/DxoUJsFAU/nFqeq16iqLwztE6JIeKGT17zhfhcDi4f1EiatsGcKJuemG6Lq0B+ytbcHduLCPXJPdlwgP9cUd2FPaWNmNgmmG6g8p2NPcO4b7rE10jzgMQox+HhBAxls0Mx4enGmG0ON/aO1nXg2p1PzYVJJFRGAzk5zkxCAoQYM+xhmld56NTjTBarLhvUaJLdBFcywOLk6A1mLG3dHqx+veONyBa5o8VGd6zhAkx+gnYtDgJXVoDvr3o/D6j7/xYh2CxEHd7STzP1/AX8PDLBQk4oGx3eu6EzmDGe8cbcMPMMLctnUGYHtmxcuQnBOHdo/Uwmp2bO6FQ9eL45W5svD7R7SufuhLvUUoT16eEYF5iEErO9kFvmvo06sqmPhys7cC9CxPJBiMMZtPiJIiFfLxy4IJT5793vAHdOiMevTHVxcoIrmTL8hlQ9w2hpLTJqfNfOXABwWIhNlyX4GJl7oUY/QRwOBw8dlMauocs+NCJWP2fvz2PoAABHlic6HpxBJcRJBbigUWJ+HdVG2paNFM6t19vwltH6rB8VjjmenCPUsLUWZoWhnmJQfjL9xen3HArbejBjxe78F9Lkj22r7SrIEY/CRamhCAnSoT/d/AiurWGSZ93UNmOo5e68NCyGZD6k+Qc09lUkAyZSIAd+89NaQTOK99dRL/ehMdXpLlRHcEVcDgc/J+bZ6K934C/Hp78VqYWK4Vnv6hBRKAfNiz0rtY8QIx+UnA4HPxmfgh0BjNe+PfktgTTGcx45l/VSIuQYOPCRPcKJLgEmUiAp1bOwqn6nklPrqlWa/De8XoUzY9HZgxZjtgbWJAcgjvnRON/D12a9Lj6v59oQE1LP565I8MrN3knRj9J4uVC/OaGFHxWpsYXleMvkERRFJ75VzVa+/XYdVcWhHzymL2FdflxmJcYhP9v/znUd+nG/azWYMYjH5cjROKHJ26Z5SGFBFew/Y4MiIQ8PPpx+YQhnHMt/dj1dS1umBmG2710D4lpOdB3332HrVu3jnnsk08+wV133YW1a9fi0KFDAGyb5D788MMoKirC5s2b0dPj3q3cXM0jN6YiLyEIT396FhXjrGP+zo/1+Kxcjd/dmIa8hGDPCSRMGy6Xg1fWzwWfx8F/FZeiV2cc83MmixWPl1SgoUuHV9fPhSyAhOa8iTCpH3avmYNqdT+e+vTsNde0au/X49cfKCAXCfB/18zx2r0FnDb6559/Hrt374bVOnqYUmdnJ4qLi/Hxxx/j3Xffxcsvvwyj0YiPPvoIaWlp+PDDD7Fq1apRO1AxHQGPi9d/mYsQiR82vnsKxy+NXOPaaqXw2sGLeOHfStw6OxIPL59Bk1LCdIiRi/B6US4augdR+PZJNPUMjjg+oDfht/8ow3/OtWP7HRlYSCbBeSU3ZUTg/9ychn9VtOCxkgoMGs0jjl/u1KLw7ZPo1hrw1w15Xr2/s9PBptzcXNx0000oKSkZdezs2bOYO3cuhEIhhEIh4uPjUVtbC4VCgV/96lcAgCVLlnid0QNARKA/Pty8APf+7TR++e4p3JYZhYLUUOiMFvyzvBnV6n78PCca/3fNHDI5you5fkYo3r03H7/9oAy3vHIEa/PjkBUjQ2PPID463YhunRHP3ZmB+xYl0S2VMA22LE8Fl8vBn789D4WqF+vmxSFaLkJFUy/2ljYjQMjDew/M9/rRVBMa/d69e/H++++PeG/nzp247bbbcOrUqTHP0Wq1kEqvrAEhFouh1WpHvC8WizEwMPbkFKVycglPT6LX60fo+vOKMJRUCfDV+XbHZiJxMgG2LQ7DsmQ/XLpwnhZdTMLbtYUCePX2KLxX1oMPTjbAbAU4ALIj/fFUQRTSg93z/Zj63JiqC5ietmURQMjNUfh7eS9e/s42j0LI42BJohgb5wZDMtQOpbKdFm2uYkKjX7NmDdasWTOli0okEuh0VxJZOp0OUql0xPs6nQ6BgYFjnp+enj6l+3kCpVI5SlduNrDTSqGlbwh+fC7CA/0ZoYspsEFbOoBl822rj7Zp9AgSC92+jg1TnxtTdQHT15aeDqxbBmgGTegbMiIi0N9lExw99dwUCsU1j7llOEh2djYUCgUMBgMGBgZw+fJlpKWlITc3F4cPHwYAHDlyBHl5ee64vUfhcTmICw6gxeQJnsNfwENiqJgsVsZyZAECJISIWTeL3aUDQvfs2YP4+HjceOON2LBhA4qKikBRFB577DH4+fmhsLAQTz75JAoLCyEQCLB7925X3p5AIBAIYzAto1+wYAEWLFjgeH3//fc7/l67di3Wrl074vMikQivvvrqdG5JIBAIhClCZvIQCAQCyyFGTyAQCCyHGD2BQCCwHGL0BAKBwHKI0RMIBALL4VBTWXjbA4w36J9AIBAI1+Zac5MYZ/QEAoFAcC0kdEMgEAgshxg9gUAgsBxi9JNgYGAAv/71r3HPPfdg3bp1KC8vp1sSrFYrtm/fjnXr1mHDhg1QqVR0S3JgMpmwbds2FBUVYfXq1Th48CDdkkbQ3d2NpUuX4vLly3RLGcFf//pXrFu3DnfddRf27t1LtxwHJpMJW7duxfr161FUVMSY51ZZWYkNGzYAAFQqFQoLC1FUVIRnn312zH0y6NClVCpRVFSEDRs2YNOmTejq6prgbPdAjH4S7NmzB9dddx0++OAD7Nq1Czt27KBbEg4cOACj0YiSkhJs3boVL774It2SHHzxxReQy+X48MMP8fbbb+OPf/wj3ZIcmEwmbN++Hf7+zFqE7tSpUygvL8dHH32E4uJitLW10S3JweHDh2E2m/Hxxx/joYcewiuvvEK3JLz99tv4wx/+AIPBAADYtWsXfve73+HDDz8ERVG0NS6u1vXCCy/gmWeeQXFxMVasWIG3336bFl3E6CfBfffdh/Xr1wMALBYL/Pzo32lGoVCgoKAAAJCTk4Pq6mqaFV3h1ltvxaOPPup4zeMxZyXAl156CevXr0d4eDjdUkZw9OhRpKWl4aGHHsKvf/1r3HDDDXRLcpCUlASLxQKr1QqtVgs+n/7NsePj4/Haa685XtfU1GD+/PkAbJsaHT9+nBG6Xn75ZccSxXR6B/0lxjCutdFKdnY2Ojs7sW3bNvz+97+nSd0VtFotJBKJ4zWPx4PZbGbEP6FYLAZg0/jII4/gd7/7Hb2Chvnss88QHByMgoICvPXWW3TLGUFvby9aWlrw5ptvorm5Gb/5zW/wzTffMGKP0oCAAKjVaqxcuRK9vb1488036ZaEW265Bc3NzY7XFEU5ntV4mxp5Wpe9QVFWVoYPPvgA//jHP2jRRb8rMIxrbbRy/vx5PP7443jiiSccLQc6uXpzF6vVygiTt9Pa2oqHHnoIRUVFuPPOO+mWAwD49NNPweFwcOLECSiVSjz55JN44403EBYWRrc0yOVyJCcnQygUIjk5GX5+fujp6UFICP370b733ntYvHgxtm7ditbWVtx7773Yv38/I3q2drjcK8GJ8TY1ooN///vfeOONN/DWW28hODiYFg0kdDMJLl26hEcffRS7d+/G0qVL6ZYDwLZn75EjRwAAFRUVSEtLo1nRFbq6uvDAAw9g27ZtWL16Nd1yHPzjH//ABx98gOLiYqSnp+Oll15ihMkDtokuP/74IyiKQnt7O4aGhiCXy+mWBQAIDAx0bAEqk8lgNpthsVhoVjWSjIwMx9amR44cQX5+Ps2KbHz++eeO31xcXBxtOpjTBGQwu3fvhtFoxAsvvADA1pp+4403aNW0YsUKHDt2DOvXrwdFUdi5cyeten7Km2++if7+frz++uuODeDffvttxiVAmcSyZctw5swZrF69GhRFYfv27YzJbdx33334/e9/j6KiIphMJjz22GMICAigW9YInnzySTzzzDN4+eWXkZycjFtuuYVuSbBYLHjhhRcQFRWFhx9+GAAwb948PPLIIx7XQmbGEggEAsshoRsCgUBgOcToCQQCgeUQoycQCASWQ4yeQCAQWA4xegKBQGA5xOgJBAKB5RCjJxAIBJZDjJ5AIBBYzv8PJqaqU3Wia2IAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"x = np.linspace(-np.pi, 4*np.pi, 1000)\n",
"y = np.sin(x)\n",
"\n",
"# 最初に0.25を超えるindexを求める\n",
"first_cross = np.argmax(y > 0.25)\n",
"\n",
"plt.plot(x,y)\n",
"plt.plot(x[first_cross], y[first_cross], 'x', color='red')\n",
"plt.show()"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"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.9.6"
},
"toc": {
"base_numbering": 1,
"nav_menu": {},
"number_sections": true,
"sideBar": true,
"skip_h1_title": false,
"title_cell": "Table of Contents",
"title_sidebar": "Contents",
"toc_cell": false,
"toc_position": {},
"toc_section_display": true,
"toc_window_display": false
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment