{
"cells": [
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"44"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"# A\n",
"sum([6, 6, 1, 5, 3, 6, 3, 9, 5])"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"31"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"# indep\n",
"sum([5, 2, 2, 3, 4, 3, 3, 5, 4])"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"sum([2,4,2])"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"10-element Vector{Int64}:\n",
" 4\n",
" 3\n",
" 3\n",
" 2\n",
" 5\n",
" 1\n",
" 4\n",
" 2\n",
" 5\n",
" 4"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"rand(1:5, 10)"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"154//75"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"(88*7) // (100*3)"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"0.5366666666666666"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"float((7 *23)//(300))"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"ename": "ErrorException",
"evalue": "syntax: invalid identifier name \"?\"",
"output_type": "error",
"traceback": [
"syntax: invalid identifier name \"?\"\n",
"\n",
"Stacktrace:\n",
" [1] top-level scope\n",
" @ ~/sci/projetos/zugzwang/code/drafts/drafts.ipynb:1"
]
}
],
"source": [
"?Rational"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"-61//100"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"1-(21*23)//300"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"1000"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"35+644+68+211+7+35"
]
},
{
"cell_type": "code",
"execution_count": 16,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"391//750"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"θ = (68 * 23)//(3000)\n"
]
},
{
"cell_type": "code",
"execution_count": 17,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"0.5213333333333333"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"float(θ)"
]
},
{
"cell_type": "code",
"execution_count": 19,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"4853//3000"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"θ_1 = (211 * 23) // 3000"
]
},
{
"cell_type": "code",
"execution_count": 20,
"metadata": {},
"outputs": [],
"source": [
"using StatsBase"
]
},
{
"cell_type": "code",
"execution_count": 24,
"metadata": {},
"outputs": [],
"source": [
"ds = 0.001 .* [0, 35, 644, 68, 211, 0, 0, 7, 35];\n",
"de = (1/23) .* [0, 0, 7, 3*θ, 3*(1-θ), 0, 0, 3, 10];"
]
},
{
"cell_type": "code",
"execution_count": 25,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"0.18670312228405667"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"rmsd(ds, de)"
]
},
{
"cell_type": "code",
"execution_count": 27,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"0.5213333333333333"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"float(θ)"
]
},
{
"cell_type": "code",
"execution_count": 28,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"3.1029411764705883"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"211/68"
]
},
{
"cell_type": "code",
"execution_count": 29,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"0.3222748815165877"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"68/211"
]
},
{
"cell_type": "code",
"execution_count": 30,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"7.666666666666667"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"23/3"
]
},
{
"cell_type": "code",
"execution_count": 36,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"ab: 253//500, 0.506\n",
"ac: -2543//3000, -0.8476666666666667\n"
]
}
],
"source": [
"n_ab = 66\n",
"n_ac = 241\n",
"\n",
"t_ab = (23*n_ab) // 3000\n",
"println(\"ab: $t_ab, $(float(t_ab))\")\n",
"\n",
"t_ac = 1 - (23 * n_ac // 3000) \n",
"println(\"ac: $t_ac, $(float(t_ac))\")"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"a: 0.03402, b: -0.03298, c: 0.30903\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAAGQCAIAAAD9V4nPAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nOzdeUAU5f8H8M/MLDcoICCXIaKIoKgcXngnlppWXl/vDkstv13mkXaZaVZ0WmZ5VJalpmnlkaJi3nkfKCgIKnIJCnKz7M7z/P7Yfvs1RV10l9nj/fprd3l29rPPLvPeZ+aZGYFzTgAAALZKVLoAAAAAJSEIAQDApiEIAQDApiEIAQDApiEIAQDApiEIAQDApiEIAQDApiEIAQDApiEIAQDApiEIAQDAppl7EF68eHHBggWGt9doNKYrBu4APa8UrVardAm2iHOOnlcE51yWZeMu09yDMD09fePGjYa3r66uNl0xcAfoeUVwztHzimCM1dTUKF2FLZJl2eg/u809CAEAAEwKQQgAADYNQQgAADYNQQgAADYNQQgAADYNQQgAADYNQQgAADbN0CBkjM2dO7d9+/Y9e/bctGnTrQ327t37+OOPt23bNi4ubs6cOWq1Wvf49u3bX3vtteHDh//22283tt+zZ098fHzbtm1nzJiBw3EAAEAphgbhl19+uWrVquXLl7/yyiujRo06d+7cTQ1UKtVTTz21cuXK+fPnr1u37u2339Y9npiYSETp6empqan6xgUFBQMHDhwzZsyqVasOHDgwe/ZsI7wVAACwdudK6OBVwbjLVBnY7quvvpo7d25kZGRkZOSQIUOWLFny0Ucf3digU6dO+ttPP/3077//rrv94YcfElFaWtqNjZcvX965c+cnnniCiObPn//YY4/NmTNHpTK0GAAAsE1zT1JnL+oRaMxlGjQirKmpSUtLi42N1d2NiYlJTk6+tZlarc7IyEhKSvr+++9Hjx59hwWePn06JiZGdzs2Nvbq1au5ubl1rBwAAGzL9RramsOHNzXyYg0ahF29epVz3rBhQ91dDw+PgoKCW5tlZGSMGzfu8uXLYWFh/fv3v8MCCwoKIiMjdbft7e1dXFwKCgoeeOCBW1tmZ2fv3bs3ODhY/8gHH3xwh4VXVFQIgpFHzWAI9LwiOOeVlZVKV2GLZFlWq9WMMaULsS3L0qUHfUhVU1FebujpRp2dnUXxLkM+g4JQF4EVFRXu7u5EVF5e7uHhcWuz8PDwI0eOMMamTZs2cuTInTt33mGB5eXlutuyLFdVVdW6QCIKDAyMjo5evny5/pEmTZrY2dndbsmcc1dXV0PeFBgXel4RnHNBENDz9U+WZTs7O2dnZ6ULsS0/X9LOjyZXV1dHR0cjLtagIHRxcfH29k5LSwsICCCitLS0pk2b3q6xKIpDhgxZtmzZHRYYHBycnp6uu33+/Hk7Ozt/f//bNXZycmrWrJkhdQIAgLU6VcSL1NTTl5iRr8Jk8KzRMWPGLFiwgHNeUFCwcuXKMWPGEFFpaenMmTNLSkqI6NChQ7rLwVRUVCxevFg/d6aysrK4uLimpqa6urq4uFh3WMWYMWM2btyYmZlJRJ9//vngwYOdnJyM/M4AAMCKLD3HngoVRRPsgTE0CN98883S0lJfX9+WLVuOGTOmd+/eRFRRUbFw4cKysjIiWrNmjZeXl7+/v4+PT2Fh4ZIlS/RPDAkJ2b9//xdffBESErJ27VoiioiIeOONN6KiogICAo4cOZKQkGD8dwYAANZCLdOqDPZEC5NMRBA454a3LioqcnJyut3oTavVFhcXe3h4GHggRHV1dXl5uZeX1x3abNu2LSEhQXcwoiHKysrc3NwMbAxGhJ5XBOe8oqIC+wjrn26yDPYR1ptVGez7dLblYZVWq9VqtQrsI9Tz9PS807JUKm9vb8OX5ujoaNw3AwAAVmlZGpsQZqpzguJcowAAYNYulvGT1/igBxCEAABgk75NY2Oaiw6SqZaPIAQAAPPFOC1P50+GmjCtEIQAAGC+tmZzP2eK9DThiasQhAAAYL6WnmPjTTkcJAQhAACYrStVtDOP/ScEQQgAADbp+zQ2pKnY4LanlzYOBCEAAJip79PZ+JYmzykEIQAAmKNdeVwlUCcfk1/fDUEIAADmaOk59ozph4OEIAQAADNUrKZNl9mYFghCAACwSSvOs36BYiOH+ngtBCEAAJidb9PYsyY7y/ZNEIQAAGBeDhbwMg318DP5NBkdBCEAAJiXJefYhDBTXIu+dnW7HiEAAIBJlWto3UWWMtTER9HfACNCAAAwIz9lsAf9RV+n+ntFBCEAAJiRxWfrb5qMDoIQAADMxZGrvFhNffzrbf8gEYIQAADMh244WH/zZIgIk2UAAMBMlGvo1wsseUh9BxNGhAAAYBZ+ymC9/EV/5/odDyIIAQDATCw+y56tl7Ns3wRBCAAAytNNk4kPqO/hICEIAQDAHHydyibU+zQZHUyWAQAAhZVpaN1FdmZI/Z1N5kYYEQIAgMJ+SGfxAaKfszKvjiAEAACFLT3HJtbv2WRuhCAEAAAl7b/CyzTUq37PJnMjqwrCGkY78iWlqwAAgDr45ix7rpUis2T+YVVBqJbpqQOqgiql6wAAAMNcU9OGLPZkqJJhZFVB6GZHAwPYd2lM6UIAAMAg36exgQ+IjRyUrMGqgpCIxjeXF59ljCtdBwAA3A0nWnyWTWqlcBJZWxBGeTIPB0rMQRICAJi7HTncUaLOPgruHySyviAkokmtxK9TsXUUAMDcfX2WPaf0cJCsMghHhoh78tnlCgwKAQDMV14lJeWy0c2VjyHlKzA6FxWNbi4uPYtBIQCA+Vp6jo1oJropc1a1f7HCICSi51qJS84xDaIQAMAsaRktNo/tomStQdjKXQhrKPx2CUkIAGCO/shiwW7UxlPhaTI61hmERPRcuLgoBUEIAGCOFqWay3CQrDgIHw8S00rpdDGmzAAAmJfzpfxUER/c1FwCyFzqMDqVSONDhcWYMgMAYGYWpbLxoaKD2ZwZ2mqDkIieDRN/Ps/KNUrXAQAA/69SSz+ks4lms12UrDsIA12EHn7iTxkYFAIAmIuVGSyusRjkahbTZHSsOQiJaHK4+MUZnHkUAMBcLEplz4ebV/SYVzVG19tfEIj25iMKAQCUd6CAX6+hPspdg7dWVh6ERDSxlbgQx1EAAJiBr1LY5HBRyYvw1sb6g3BcC3FbDsvH1XoBABRVWE2bLrMnW5hd7phdQUbXwI6GNxNxHAUAgLKWnmODm4oeil6Dt1bWH4RE9EKE+E0qTj0KAKAYmZvRyUVvYo41GV24uxDakNZdRBICACjj90ss0IWivcxs9yAR2UgQEtHkcEyZAQBQzMIUNtksh4NkO0H4WJB4oYyOXcVxFAAA9S31Ok+9zgcHm2nimGlZRqcSaVIr8WtMmQEAqHcLU9izLUV7cw0cc63LBJ5tKa69wIrVStcBAGBLSjW0MsO8Ti56E/OtzOh8nOiRJuKyNAwKAQDqz3fn2EOBor+zOU6T0bGhICSil1uLX55hMnYUAgDUC060KJW9EGHWWWPWxRldlJfg60ybsjAoBACoD1suc2cVdfYx3+Eg2VoQEtEL4eIXOI4CAKBefJEiv9Ta3IPG3OszumHNxDPFPPU6No8CAJhWegk/epX/p5m5B42512d09iJNCBO/OINBIQCAaX2ZwsaHio6S0nXcjc0FIRFNDJNWZrAiHEcBAGAyZRr66bzZXYO3VhZQotH5OdOAB8TvcBwFAIDJfHuO9QkQA13MepqMji0GIRG90lr8AsdRAACYhkUcNaGnMrzpb7/9tmbNGmdn58mTJ7dr1+6mv549e3b58uUXL150dHSMj48fMWKEKIpExDlfsmRJUlKSj4/Pq6++GhQUpGt/4MCB77//vqysrE2bNpMnT27QoIGx3pIhor0EX2famMUeDbKMzwkAwIJsvsydVRTX2AKGg2T4iHDdunUTJ0585JFHwsLCevbsmZWVdVODrKwsNze34cOHd+3adebMmR9++KHu8Q8++OCzzz4bMWKEnZ1d9+7dq6qqiOjYsWPx8fHt2rV74okndu/ePXLkSCO+JQO9GCEuwJQZAAAT+PKMBRw18T/cMHFxcYsWLdLdHjly5Ouvv36HxgsXLuzWrRvnXKPR+Pn5JSUl6R6Piopavnw55/zjjz9++OGHdQ+ePHnS3t6eMVbrohITE+Pj4w0sknNeWlpqYEu1zAN+1iQX1f66UFeG9zwYEWOsrKxM6SpskVarraioULoKM5VazPx+qqnWmmThGo2mqqrKuMs0KLE554cPH+7evbvubvfu3Q8ePHi7xnl5eZs2berRowcRZWVlXblypVu3bjc9sWPHjqdPn87Ozuacb9iwIS4uThDqewRtL9JzrTAoBAAwss/PsIlhkoPZHzWhZ9A+wqKiopqamkaNGunuenl55eXl3drs9OnTXbt2LSkpiYuLmzFjBhFduXKlQYMGKpVK/8SjR48SUVxc3HPPPRcUFOTo6Ojj45OUlHS7l87Pzz969Gjv3r11dyVJmjJlij5Zb1VRUWF4po57QGi70W5mWLW3I6bN3K869TwYC+e8srJS6SpskSzLarWaMfySvlmJRlidYXe4v6a83CTrVe3/M7C9s7OzbsLKHRgUhM7OzkSkVv9z5F11dbWLi8utzVq3bn39+vXi4uL//ve/o0eP/v33352cnKqrq/UN9E9cs2bNt99+e+rUqWbNmn3++ecDBgw4ceKEvb39rcv08vIKDg6eNWuW/pGoqChXV9fblco5v8Nfb+JKNCRYXpmtmtXOcrZlm6s69TwYC+dcEAT0fP2TZdnOzk63boQbfXWKDQriIV61rM+NQpeCjo6ORlymQQHg5OTUqFGjS5cu6e5eunSpSZMmt2vs4eHxyiuvbNu2jYgCAgLUanVBQYH+iYGBgUS0fv364cOHR0REODk5TZ8+/fz586mpqbUuTaVSeXp69rmBp6dnnd7hnb3cWvwqldXgVx0AwH2TOX2dyv5rIUdN6Bla7tChQ7///nsiqq6uXrVq1dChQ4lIrVYvX768rKyMiHJzc/WNN23aFBYWRkTe3t49evTQPbGwsHDz5s26JwYGBh47doxzTkQnT55kjPn7+xv1fRkqwkMIa0i/XkASAgDcr/UXWaALxXhZ2F4SQ48jfP3113v16hUXF1dYWBgaGjpkyBAiKi0tffLJJ+Pi4tzc3GbOnLl3796goKDs7GzG2KpVq3RPTEhIGDhw4JYtW9LS0kaNGhUdHU1EU6dO7du3b5s2bYKDg/fv3//+++97e3ub6B3e1UutxXePs5EhFvYTBgDA3Hx2mk1pY3nrUkODsEmTJqmpqUePHnV1dW3durXuwUaNGuXm5vr4+BDR8uXLMzIycnNzvby8QkNDJemfCUMxMTHnz58/fvx4QEBAcHCw7kEfH59jx46dO3euqKgoLCxMPw1HEY88IE4/xPZd4ZZy7CcAgBk6cpXnVJIlnqWkDmeWsbOz69Sp042PiKLo5+envxsSEhISEnLrE11cXLp27XrTg6IotmrVqi6lmopA9N9w8bPTLK6x5cz2BQAwMx8ns5ciRMkCBxSWF92m8FSouCuPZZbhIAoAgHuRU8G3ZbOnQi0yUyyyaKNzVtGToeJCXLkeAOCefJnCxrYQG5rqoAnTQhD+44UIcXkaK9MoXQcAgKWp0NKycxZzrYlbWWrdRtfERYgPFJeew6AQAKBuvktjPfzEZm4WuHuQiBCEN3q1jbjgDNMiCgEADMY4fX6avWJB15q4hQWXbnQxXkITF1p/CUkIAGCoP7KYhwN1seTDzxCE//JKazHhFIIQAMBQnyazqRZ4EP2NLLt6oxsUJBap6UABjqMAALi7I1f5pXIa3NSyo8Syqzc6SaCXIsSPkzEoBAC4u0+S2YsRosrCk8TCyzeB8S3F3XksoxSDQgCAO8mu4InZbHxLi88Ri38DRuesomdaip/jyvUAAHf0STJ7KtRSD6K/EYKwFv+NEFecZ9fUStcBAGCuSjX0Q7oFH0R/I2t4D0bn7ywMekBcfBaDQgCA2n2Tyh5uIj7gasFHTeghCGs3pY345RlcuR4AoBYaRl+msCmWfBD9jazkbRhdpKfQxpNWZiAJAQBu9ksma9GAoiztSvS3gyC8raltpI9OMUweBQC4ycfJbFqk9VzAFUF4W30CBHuJtlxGFAIA/M+2HK7l1DfQSoaDhCC8s1fbiB8ly0pXAQBgRhJOyVPbiNYTgwjCOxseLGaU0tGrGBQCABARnSriqddpRIhVZYdVvRmjU4n0Ek7DDQDw/xJOsRcjRHvrig7rejcmMCFM3IkzrgEAEGVX8D8vswlh1hYc1vZ+jM5FRc+0FD87jUEhANi6j5PZM2HWcE61myAI7+6lCGllBiusVroOAADlFKvpx3T2olWcU+0mVviWjM7HiYYGiwtTMH0UAGzXwhT2WFPR39mapov+A0FokFfbiItSWaVW6ToAAJRQpaWFKfIUC78S/e1Y57syuhYNhe6+4tJz2FMIALbo+3TWwUcMd7fC4SAhCA03PVL86BTTIAoBwMbInD5JZjMirTYvrPaNGV2stxDakFZlIgkBwLasyWT+ztSlsXUOBwlBWCcz2krzTzCchxsAbEpCMpvR1npOsX0rBGEdxAcITirajNNwA4DNSMzhGkb9mljtcJAQhHU1I1J8/ySOowAAWzH/hDw90qpOsX0rBGHdDAkWC6ppTz4GhQBg/Q4V8kvlNKKZlSeFlb89o5MEmtZGnI9BIQDYgHkn2NQ2osrag8La358JjGshniqiY7g2EwBYtdTr/FABeyrU+mPC+t+h0TlI9Epr8UNcmwkArNr8E+zF1pKTSuk6TA9BeC8mhIlJuSy9BINCALBOF8r45svs+VY2kRE28SaNzs2Ong8XE5IxKAQA65Rwik1qZYVXXKoVgvAevRgh/XqBZZVjUAgA1ia/in7JZC9GWPNB9DdCEN4jTwd6KlT8FBfsBQCr8/EpeUwL0cdJ6TrqC4Lw3k2NlH5MxwV7AcCqFKnp2zT2SmsbSgcbeqtG5+tEw5qJn5/GMYUAYD2+OMMebyoGuVr3yWT+BUF4X6ZHit+cZSU1StcBAGAMZRpamCK/1ta2osG23q3RBbsJA5qIX6ZgTyEAWIOvUljfQLF5AxsaDhKC8P690V5ccEYu1yhdBwDA/amWacEZZmvDQUIQ3r/mDYRefuI3ZzEoBADLtvgs6+QjtPawreEgIQiNYlY78ZNkVo1JMwBgsdQyfXSKzWpni6Fgi+/Z6CI9hVhvYek5DAoBwFJ9n87aeFK0l80NBwlBaCxvR4kfnmRqDAoBwAJpGH1wks1qZyunkrkJgtA42jcS2njS8nQMCgHA8vx4njVvQHGNbXE4SAhCI3qrvTTvBKtBFAKARZE5fXCSvdHeRoeDhCA0oo4+QsuGtOI8khAALMnKDObnTN19bXQ4SAhC43qzvfTeCaZFFAKAhZA5zT3O3rTh4SAhCI2rm6/QxIV+ykASAoBlWJXBGjnSg/62OxwkBKHRvR0lvXscg0IAsACM0/sn2TtRNj0cJASh0fX0EwJdaGUmkhAAzN3qTOZmR30CbHo4SAhCU3g7Spp7nMm4dj0AmDHGad4J9rbNDwcJQWgKvfwEXydaiT2FAGDG1lxgDe3poUBbHw4SgtBEZkdLc7CnEADMlW44OBvDQSJCEJpILz8hwJlWYU8hAJilXzKZq4ribX7voA6C0FTejZbeOYZBIQCYHcbpvZNsTjSGg/9AEJpKV1+hiQv9jD2FAGBmVmeyBpgsegMEoQnNjsYxhQBgXmROc46zdzAcvAGC0IS6+wpBrvQDzj4KAGbjp/OskYOtn0rmJghC05oTLc05hktSAIBZkDm9d4LNjcFw8F9Uhjc9fPjwxo0bPTw8xo0b5+npedNfs7Ozt27deunSJVdX1/79+7du3Vr/p8TExN27dwcFBY0dO9bR0VH/+I4dO3bv3u3o6Pjwww+3b9/+Pt+JeerSWAhzp+/T2IQw/OYAAIUtT2eBLtTTD8PBfzF07bx58+aHHnrIwcHhyJEjHTt2rKiouKnB1q1b9+/f7+jomJub27lz599++033+IIFC5599tmGDRuuXbu2X79++vYvvfTSxIkTtVpteXn577//bpQ3Y57eiZLmncDF6wFAYRpGc7F3sFbcMF26dFm0aJH+9uLFi+/Q+K233nrsscc45zU1NX5+fklJSZzz6upqPz+/Xbt2cc537drl6el59erVu75uYmJifHy8gUVyzktLSw1vXG8GbNF8cUZWugrTMs+et3qMsbKyMqWrsEVarbaiokLpKurm61T5oT81SldxvzQaTVVVlXGXadCIsLq6ev/+/frx3EMPPZSUlHS7xhqN5syZMy1btiSi1NTUkpKS7t27E5GDg0Pv3r11T9ywYcPQoUMzMjIWLVq0b98+I+S5eZsTLc0/wSq1StcBALaqWqZ5x3GhidoZtI8wLy+PiHx8fHR3fX19d+zYcWuzgwcPjho1qqioqGPHjnPmzNE90cvLS5L+6frGjRvn5uYSUWZmZmpqamZmZocOHT788MNRo0bNmzev1pcuKipKS0ubOnWq/pFhw4a1bdv2dqVWV1fb2dkZ8qbqU7grdfASFyTXvNzKaqfNmGfPWz3OeXV1tUpVh539YBSyLKvValG0mH3/X6YK7TzFtg3U1dVKl3J/tFqtVluHUYW9vf1dPyaD/n90ScbYPytxWZZrXeVFRkZu27YtJydnypQp77zzzrx580RR1D/rpidyzrdu3SqK4ujRoyMjI1999dVbJ+AQkYODg729vf5PgiB4eHjok7XWUu/wVwXNbk99E2lCS8nNSsPCbHveunHO0fNKsaCer9LSZ6n024NkKQXfgW5jpuFvRBDuPjPIoCD09fUVRTE3N7dFixZElJub6+fnd2szJyenZs2aNWvWbO7cuePHj583b15AQMDVq1dramrs7e11T4yMjCSigIAAQRB0KR0eHi6KYlZWVq1B6OLi0rRp01mzZhlSJxHZ2dmZ57ikrTfFB8hfnRPeaG8xPyHrxGx73rpxztHzitD9yreUnv8slXVpzGMaW3wKEpEgCIIgGLfnDVop29vbx8fHr1+/noi0Wu0ff/wxYMAAIlKr1QcPHtSNUmtqavTtk5OT/f39iSgsLMzf33/z5s1EVFpaumPHjv79+xPRwIEDT548qXviqVOnBEEIDg424rsyT7OjxAVn5GK10nUAgC0p1dDHyfKcaOv8CW4Uhu5aePvttwcMGHD+/Pn09HRnZ+fBgwcTUVZWVqdOnQoLC728vB5//HFRFAMCAi5dunTo0KFff/2ViERRnDt37oQJE7Zt27Zv376+fftGRUURUZ8+fUJDQ3v06NGxY8c1a9bMnTu3YcOGpnuTZqJ5A+HxpuJHyfI8HM0KAPXlk2S5fxMxzB3HDt6WwLmhV1K/fPnyjh073N3d+/fvr9vUWVVV9ffff3ft2tXOzq64uHjfvn25ubm+vr7du3d3d3fXPzE1NXX//v1NmjSJj4/Xb65ljCUmJubn58fGxkZERNzuRbdt25aQkJCYmGhgkWVlZW5ubgY2rn+5lbztOu3pIXaNnZQuxdjMvOetFee8oqLC1dVV6UJsjm6yjLOzs9KF3MXVamq1VnPwUVUzNysJQt1kmRvPzXL/6hCEirCyICSiKX/LjOizTtY2KDT/nrdKCEKlWEoQTj0o1zBa0Nl6VjimCEJsNa5vM9tJP51nl8rN+vcHAFiB7Ar+fRqb1c56UtBEEIT1zduRnmslvnPMag8oBAAz8e5x9myY6Gt1O2KMDkGogGmR0pZslnIdg0IAMJX0Er7uIpsWieHg3SEIFeBmR6+2kd48gkEhAJjKG0fZq20kTwel67AECEJlTA4XDxfyQ4UYFAKA8R2/xvfm8xcisIY3CLpJGY4SvRUlzjyMizMBgPHNPCy/0V50wTloDYMgVMxToWJ+JSXmYFAIAMa0K4+nl9D4lli9Gwo9pRhJoHdjxOkHZYYoBAAj4USvHZbnx4r2WLsbDF2lpMFNRWcV/ZKJWTMAYBy/XmBVWhoajHV7HaCzFPZerPTGUaZBFALAfdMyeuMIS+goiVZyPrV6giBUWE8/IbQhLT6LJASA+7X0HGviSvEBiMG6QRAq74NYae5xuVSjdB0AYMkqtPTucfZ+LI6grzMEofLaeAoPBYofncKhFABw7z5OZj39hGgvDAfrDEFoFubFiF+lsOwKzB8FgHtRWE1fnMHVd+8Res0sBLgI41uKc49jTyEA3Is5x+TRzcWQBhgO3gsEobmY2U5afwln4gaAOksr4asz2ZvtsXfwHiEIzYW7Pb3WVppxCHsKAaBuXjvMpkVKjXB+7XuFIDQjk8PFs9dpRy4GhQBgqL8L+JFC/t9wrMzvHfrOjNiLNC9GnHFIRhICgCE40ZS/5fdiRSecX/s+IAjNy7Bmokqkn89j1gwA3N2aTFbDaFQI1uT3Bd1nXgSizzpJMw+zSq3SpQCAeathNOsI+6ADTqh2vxCEZqeTj9DRR/jsNAaFAHAnn59mbTyEB/0Rg/cLQWiOPuwgfpIs51cpXQcAmKsiNX2cLH/QAetwI0AnmqNgN+HJUHH2URxKAQC1m31MHt5MDG2I4aARIAjN1OvtpN8usdPFmEAKADc7e52vzmRvR+EIeuNAEJopDwd6s7300gEMCgHgZq8elGe1xRH0RoMgNF+TWokFVbTpMgaFAPA/Sbn8XAk9hyPojQddab4kgT7sIE09KOP69QCgI3N6+W/5446iPVbexoO+NGv9mghNXWlRKpIQAIiIFp9l3o70aBBW3caE3jR3n3SS5h6Xr1YrXQcAKK1YTe8ckz/thDkyRoYgNHet3IURIeLbxzBrBsDWzT4mD24qRnrikAkjQxBagHeipF8vsOQizJoBsF1nr/OVGWxONIaDxocgtAAeDvRGe2nKQQwKAWzXK3/Ls9pJXo5K12GNEISWYVKYmFdJ6y9i1gyALdqYxS+V02QcMmEa6FbLoBJpQWfp5b9xVQoAm1PD6NWD8iedJDussE0D/WoxevsLMT03KJcAACAASURBVF7CJ8kYFALYlk+TWSt34eFAzJExFQShJfmoo/jZaflyBWbNANiKvEr6KFn+qCPW1SaEzrUkwW7C5HBx+iEMCgFsxYxD8oQwsXkDDAdNCEFoYV5rKx0s4DvzMCgEsH77r/C/8vjMtjhkwrQQhBbGSUUfdRRf2I8TkAJYOcbppQPyRx1FVzulS7F2CELLM7ipGOBMX6UgCQGs2depzM2OhjfDWtrk0MUWaUEXad4JOb9K6ToAwDQKq+md4/KCLtgoWh8QhBapZUPh6ZbidJxrBsBKvXZIHttcbO2BOTL1AUFoqd5qL+25wv/CrBkAq3O4kG/J5m+2x3CwniAILZWzihI6iC8ekLXYVwhgRWROE/fKCR3FhvZKl2IzEIQWbGiw6OdEX2DWDIAVWZTK3O1pZAhWzvUHfW3ZFsZJ80/IOTjXDIBVuFJF7x6XF8ZJ2DdYnxCElq15A+G5VuKUgxgUAliDKX/Lz7QUW7kjB+sVgtDizWonnbjGN1/GoBDAsu3O5/uu8FntMEemviEILZ6DRJ93ll46IFfjYAoAi6WWaeJe+YsuootK6VJsD4LQGjwcKLRtJLx3AkkIYKk+OMVauQsDH8A6WQHodCvxRWfp61SWeh0bSAEsz/lSvuC0/AmutaQQ9LuV8HOmN9pLk/bKSEIAi/P8Pvn19lJTN8yRUQaC0HpMDhcrtfRDOmaQAliSnzPY1Wp6IRxrY8Wg662HJNDibtJrh+Sr1UqXAgCGKVLT1IPyN10lFVbGykHfW5X2jYTRzcVXcTJuAAsx7aA8vJkY642NokpCEFqbOdHS3ny+LQf7CgHM3e58vj2XvxuNAwcVhiC0Ns4q+rKL9Pw+uVKrdCkAcHvVMk3YIy/oLLrhAvRKQxBaoX5NhFhv4Z1j2EAKYL7mHpfbeAqPBmElrDx8BtZpQWfpx/Ps6FVsIAUwR8lF/Juz7LNOWAObBXwM1snLkebHSuN3yxocTAFgZhiniXvl92OlABfMkTELCEKr9UQLsbETfXoaSQhgXj4/wxwlerolVr/mAp+ENVvcTUo4JeO8awDm42IZn39CXtwNVxw0IwhCaxbkKrzVXhq/W2aIQgAzwIkm7JVntpOaN0AOmpE6BGFWVtaKFSt27NjBWC1b2yoqKnbt2rVixYqNGzeWlZXd+KczZ86sWLHi77//vvVZ2dnZR48erWvRYLjJ4aIo0KJUbCAFUN7Ss6xUQy9GYARiXgz9PLZv396+ffsdO3ZMnTr18ccf5/zmIca77777zjvvJCUlffzxxyEhIWfOnNE9vmTJkt69e+/atWv06NHTpk278SklJSVdu3aNjY29/7cBtyMKtKy79M4x+WIZRoUASrpcwV8/Ii/tJmGrqNnhhunQocPXX3/NOS8vLw8MDNy5c+cdGo8ePfr555/nnFdXV/v4+OgaX7582cnJKSsrS99s/Pjx06ZNEwThDotKTEyMj483sEjOeWlpqeGNbcf7J+TemzTMlC+BnlcEY6ysrEzpKmyRVqutqKio01P6bdHMPS6bqB7bodFoqqqqjLtMg0aEV65cOXTo0LBhw4jIxcWlX79+GzZsuEN7tVrt6elJRAcPHiSiHj16EFFgYGBsbOzmzZt1bZKSkjIzM5944on7DHIwxNRIsVxDS85iAymAMr5NY3mVND0SG0XNkcqQRtnZ2c7OzrpsI6LAwMCUlJRbm/39998JCQnp6emhoaGvvfYaEeXk5AQEBAiCoH9iTk4OEVVUVLz00kvr16+vqam580tXVlZmZ2d//fXX+kfi4+ObNm16u/ayLMsyzqhSi6VdhQf/lPv6syamOXQJPa8Izjl6XhHy/zOkcU4Fn3mYb3lIFDk+q/tVp54nIlEU9Rl0OwYFoSzLkvS/08KqVCqNRnNrs2bNmk2cODE9PX3+/Pl//vnn0KFDZVkWRfHWJ86YMWPcuHHNmzevNVBvVF1dXVpaeuTIEd1dSZLatWsXEBBwu/YajabW2qCFC/03TJywl//em5kiCdHzitBtKULP1z9Zlg3v+Un7xImh1MoVH5QRaLVarVZ7YyTdmb29vXGC0M/Pr6ysrKqqysnJiYiuXLni7+9/azMfH5++ffv27dvXwcHhgw8+GDp0qK+vb2Fhob7BlStXYmJiSktLv/nmm3Hjxk2cOLG4uJhzPnHixOnTp4eEhNy6TE9Pz/Dw8KVLlxpSJxFpNBpHR0cDG9uaWVEUt0H740XVhDDjb59BzytCNyJEz9c/WZYFQTCk579NY/nV7M1olR02ixqDLgiN+5036JMJCAgICQnZtm0bETHGtm/frtvtJ8tySUnJre2LiopcXV2JKCYm5tq1a2fPniWiioqK/fv39+jRw8HBYeHChR07doyOjo6IiBAEITo62sXFxYjvCmqlEml5D+n1I3JGKWaQAtSHnAo+67C8vIeEFDRnBo0IRVGcNWvW888/f/HixQMHDoii+OijjxLRwYMH4+LiOOdENGzYsBYtWjRu3DgtLe3HH39cvXo1Ebm7u0+ePHn48OETJ05cv359r169IiMjiWjChAm6JaekpLzzzjv6u2BqYe7C9Ejp2T3yjgEqTOEGMClO9PRu+cUIqY0n/tvMmqG/Up5++unvvvvu8uXLMTExe/bssbOzI6KQkJBFixbpGkybNq1BgwaXL19u1qzZiRMn+vXrp3v8/fffnzVr1sWLF4cOHbpmzZqbFuvn53fjRBioB1PaiGpGX5zBDFIA0/o6lV2vwUxRCyDwWw6NNyvbtm1LSEhITEw0sH1ZWZmbm5tJS7ICF8p4h9+1fw1QRXgY7Ycqel4RnPOKigrdngioT7Isq9VqZ2fn2zXILOOd/9DufkTVsiGGg8ak2D5CsDLBbsKcaGncLlykCcAkGKcnd8lvtJOQghYBQWijJrUSvR3p/ZNIQgDj+/AUsxfpvzinqIXA52SjBKJvu0sLU+TDhWa9bRzA4hy/xj89LX/XA6cUtRgIQtvl7ywsipNG7pTLcZAvgJGoZXpil/xZJ8lEp3ACU0AQ2rTHm4pdfIQZh3HSJwDjeO2wHOEhjAzBqtWS4NOydV90kTZf5psvYwMpwP3alsPXXeRfdTH07F9gJhCEtq6hPa3oKT2zR5tfpXQpAJbsajU9tVv+rrvk4aB0KVBHCEKguMbC06HiU7u0GBUC3LPn98ljmgu9/bFr0PIgCIGI6O0oqUhNi1JwNAXAvVh8lmWU8TnR2ChqkRCEQERkJ9KKntLsY/LpYgwLAeom5Tp/44j8U0/JHitUy4TPDf7RoqHwYQfpPzvkSq3SpQBYDrVMo3fK82OlMHdsFLVUCEL4nydDxXaNhKkHcTQFgKGmH5JbNBDGt8S61ILhw4N/WdRVSszh6y9iZyHA3W3IYn9k8cXdsGvQsiEI4V8a2NGq3tKkffLFMuwsBLiT3EqasEde0VNyt1e6FLg/CEK4WYyXMD1S+k8Srk0BcFuM0zMHpJdaS3GNsWvQ4iEIoRZT2ohejvTWUewsBKjdnBNcEnDRXSuBTxFqIRAt76H66Tz/E6deA7jF9hz+fTot68JEjAatAoIQauflSL88KD21W3upHFkI8D9XqujJ3fKyroKPI/41rASCEG6rk4/wSmtpzF+yFjsLAYiISOY0Ikk7MUx80F/pUsB4EIRwJ9Pbig3saCau0wRARERvHpHtRHq9HdacVgUfJ9yJQLSip+rXi/zXCxgVgq3bdJn/lMF/6qXCrkErgyCEu/BwoFW9pef3yxml2CMCtutiGR+/W7uyl+TtqHQpYGwIQri7Dt7CW+2loTvkKpyGFGxStUzDdsivtZW64KhBa4QgBINMDhfbegrP7sXOQrBF/90vhzQQXm6NFaZ1wucKhloUJ6UU80Wp2FkItmXpObY3HycUtWYIQjCUk4p+eVCafUz+uwA7C8FWHLnKXz8i/xYvNbBTuhQwGQQh1EHzBsIPPVRDd8i5lchCsH7X1DR8h7ywC641aOUQhFA3DwUKz7YUh26Xa7CJFKyaltGw7doxzYWhwVhPWjl8wFBnb0WJPk7CK39j4gxYs2mHZEcVzY7CrkHrhyCEOhOIfugp7czlS89hVAjWaXk623yZ/4xj520DghDuRQM72vSQ9OYReVcedhaCtTl2lU8/JK/rgyvu2goEIdyjYDfh+x6qUTvlyxXIQrAeORX80W3ysm6qCA8MBm0FghDu3UOBwpQ24qOJcgXOOANWoUpLj22TJ4eLjzyAFLQhCEK4L6+2EWO9hXF/yQzDQrBwnGj8HjnYTZjRFitG24LPG+7XF12kgmr+3hmV0oUA3Je5x1lmGf+hp4TBoK1BEML9shdpXR/V6ovSivOYRAqWanUmW3aO/RavcsThErYHQQhG4O1Ia7vXvHpQ/guTSMECHbnKXzog/95X8nVSuhRQAoIQjKNlA/5DD9WondrMMmQhWJKLZfzRRHlZd1VbT2wTtVEIQjCahwKFt9pL/bfIRWqlSwEwzPUaGrBVntlOHNAEKWi7EIRgTJNaiQMfEB7bplXj/Gtg9jSMhu3QPhQo/Dcca0Kbho8fjOyDDlJjJ+Hp3TK2kII540TP7JFdVcJHHTE9xtYhCMHIRIF+6CFdKuevHcKoEMzXm0fktBL+Uy8JZxMFBCEYn5OK/uir+iOLf34aB1SAOVpylq3K5L/Hq5xx+CsgCMFEPB3oz4ekhGT26wVkIZiXPy6x2cfY1oclHxwsAUSEIATTaeombOwrPb8fBxeCGdmbz5/dK//RVwppgE2i8A8EIZhQu0bC6t6q/yRpTxYhC0F5Z4r50B3aFT1V0V5IQfgfBCGYVk8/YWEXacBWOaMUWQhKuljG+22RP+0kxQcgBeFfsKcYTG5osFihpQc3y3sGSk1csA4CBRRW08Nb5CltxJEh+PUPN8N3AurDEy3E51qJOOkMKOJ6DcVv1o5uLr7cGms8qAW+FlBPZrQV+zcR+m3RlmmULgVsSbmG+m/R9vYX3myP1R3UDt8MqD/vd5CivYSHt2hxRXuoHzWMhu3QtnQXPu6E08fAbSEIof4IRF92kULchKHbcTJSMLkaRkO2az0chGXdcK1duBMEIdQrUaBvu0tudsLwJFmDQ+3BZLSMRibJdqLwQw+cRA3uAkEI9U0l0k+9JJFo5E5ZiywEE5A5jd0lq2W+qrekwkoO7gbfEVCAnUirH5SqtXz0X8hCMDKZ0xO75KJqvraPyh5rODAAviagDHuR1sWrqmWMC8GYZE5P7ZILqvhv8SpHzI8BwyAIQTH2Iv3SW6rQ8jEYF4IxyJye3CXnV/Hf41VOOFkIGAxBCEpykGh9H1W5ho/YibkzcF+0jEbvlAurkYJQZwhCUJiDROviVTKjoTtkHFMB96aG0YidcrkGKQj3AkEIyrMX6ZcHJQeRHtumrcSx9lBH1TIN3qZlnNbFqxywXxDqDkEIZsFOpJW9JX9n4eEt2lKcgw0MVqGlQYlaN3thdW8Jc0Th3uCLA+ZCEmhJN6m1hxC/WXsN5+YGAxSrKX6ztqmb8FNPyQ4rM7hX+O6AGREFWhgn9fYXum/Q5lTg+oVwJ3mV1GOTNq6x8E1XnDsG7guCEMyLQDQ/Vnqqpdh1o5xWgiyE2p0v5V03aEeFiAkdcR5RuF91m1+VlZV17dq18PBwBweHW/+q0WjS09NramqaNm3q7u5+45/S0tLUanVERIQo/hO9nPMLFy6Ul5c3a9bM1dX1nt8AWKWpbcRGDtRzk/a3eFUHb6zo4F+OXuWDEuXZUeKzYfgpD0ZQhyB8+eWXV61aFRwcnJOT8+eff0ZERNz416qqqsaNGwcFBbm4uKSkpLz55pvTpk0jIo1GM3To0FOnTjVs2JAxtn37dh8fn4yMjN69ewuC4OXllZGR8emnnz755JPGfWNg6Z4KFb0dhYGJ2u+7q/o1QRbCP7bn8NF/aRd3lR4NQgqCkXDDHDp0yMvLKz8/n3P+1ltv9e/f/6YGsiyXlpbqbh8+fFgQhMLCQs75ihUrWrduXVlZyTkfMWLEK6+8wjnPy8s7dOiQrvGmTZvs7e3Lyspqfd3ExMT4+HgDi+Sc62uAemainj9whfmuqFl2TjbFwq0AY+x2/ztWaXma7LuiZk8eU7oQrtVqKyoqlK7CFmk0mqqqKuMu09CfVL/88ssjjzzSuHFjIho/fvyWLVtKSkpubCCKopubm+52kyZNRFHUarW6J44aNcrJyUn3xNWrVxORr69vbGysrnGXLl1qamqKioqMlOxgVTr5CPsGqT48yV46IGOHoY37/DR75xjbOUDV1RdbCMCYDN00mpWVpd8W2qRJE0mScnJyGjZseFOzBQsWZGZm7tu37/PPP/f19SWiS5cujRo1SvfXZs2a5eXlaTQaOzs7/VMWL14cHR3dpEmTWl9Xq9UWFRVt375d/0hMTMxNOyDBujVzE/YMVA1K1I77S17aTcIR0zaohtHEvfKZYr5/kKqxk9LVgNUxNAgrKyv1E2QEQXBwcCgvL7+1WWBgYHV1tb29/e7duydNmiRJUlVVlf6Jjo6OnPPKykp9gm7duvWjjz5KSkoShNp/4hUWFl64cOG9997TPzJlypSePXvers6KiorbLQpMyqQ970S0oQdNPmTX/Q/tqm5ab0cMDv+h+4dSugrTul4jjN6ramhPG3tonGWqbcWjAFmW1Wo1YzhDbn3T/j8D2zs7O+snad6OoUHo6+ur33qpVqvLy8v9/PxubTZ48GAieumll4KDgxMTE/v169e4cWP9E69du+bs7KxPwZ07d44dO/a3335r3br17V7Xz88vOjo6MTHRwDo555iDqghT97wr0ap4euOIHJ8k/dFXCnfHzx0iIs65IAhW/J0/e50P2iE/HiTMj5VEoZbJ6kqRZdnOzs7Z2VnpQmyOLgUdHR2NuExD9xFGR0fv27dPd3vfvn3+/v7+/v63a2xnZ2dvby/LMhHFxMTc+MSYmBjd7f37948YMWL16tVdunS59/LBlghE82Kkt6PEXpu0my5jUGj9tmTzHpu0M9uKH3TAIfNgQoaOCEeNGjV79uzZs2d36tRp2rRpL7zwgiRJRDR27NjWrVvPmDFj3759W7ZsiY2NZYz9/PPPdnZ2PXr0IKJJkybFxsbGxsb6+vrOnj170aJFRJSRkfHwww8PGDAgPT09PT2diAYNGqTbpwhwZ2Obi6ENhCE75BfCxeltsXq0Wh8ls0+T2bo+qrjG+JDBtAwNwgYNGuzatevjjz9euHDhc889N2nSJN3j3bt3DwgIIKKmTZsS0XfffadSqSIjI7/66ivdJNLQ0NAtW7YsXLiwqqrqyy+/fPzxx3VPfP7554koMzNTd1etxsklwVAdfYSDg6TB2+WjV/m33SVXu7s/BSxIpZae2SOnl/C/H5WauCAFweQEzs16E9O2bdsSEhIM30dYVlamP4oD6lP997xapsn75UOF/NcHpRYNbXR1yTmvqKiwpn2EGaV8yHa5XSPh666SoxnPENZNlsE+wvqn5D5CAHPjINHSbtJ/w8WuG7XrL2LynjXYkMXiNmgnhInf9zDrFAQrg2s5g2WbECZGewnDdsj7rvD5sbgWj6XSMnrjqLwqg//RF2eXhfqGIASLF+0lHHlM9cQubbeN2tW9pSBXrEYtTHYFH7lTbmBHRx5TeRlzixcQEX3xxRfLly9XugojCA0N/fnnn02xZAQhWANPB/qjr+qTZNbxd+1XcdLgphgYWozfL7FJe+VX2kjTIjEH2CROnjzZv3//Rx99VOlC7ktBQcGECRNMtHAEIVgJgejVNmJ3X2HkTnlrNv+0k+SMb7d5q9LS1EPyn5f5+nhVJx+EoAkFBQVFR0crXcV9yc7ONt3C8cMZrEqst3D8cVWVlmJ+0x69atYzom3c8Ws85jdtsZqOP44UBIUhCMHauNnRDz2lt6PEAVu1804wXLTC3Mic3j/JHt6ifaO9+HMvqaG90gWBzUMQgnX6TzPx6GOq3XksboM29TrC0FycK+FdN2i357Ajj6lGhmD9A2YBX0SwWgEuwpZ+qqdDxR4btR+eYlocaqgomdPHyazbBu3Y5uK2/iqcMgbMB4IQrJlANCFMPPyYKimXdfpDe7IIQ0NlJBfxLn9o/7zM/n5U9Xw4ZoeCeUEQgvULchW2PKx6IUJ86E/tjENypaEXMgMjqNTSzMPyg5u1E8LEbf1VzdwQgmB2EIRgK55oIZ4abJdbSRG/ajdkYTtpfdh8mbf5VXupnE4NsRvfEgNBuJODBw/Omzdv1KhRP/74Yz2/NI60Ahvi40Q/9pSScvnk/fLis+yzTlJIA6ycTeJCGX/5b3b2Ol/UVeobgE6Gu9uwYcP169cvXLhw6tSpen5pjAjB5vT2F04NVnX3FTv9oZ11WC7TKF2QdSnX0OtH5A6/azv5CMlDVEhBuMmmTZv++OMP/d0NGzZs3LiRiObOnfvll1+2bt26/ktCEIItshNpWqR4crAqr4rC1miXnsPhhkbAOH2XxsLWai+X04nHVTPbivZYwcAtGjZs+PLLLzPGiIhzPmXKFMWvI4ZNo2C7/J2F77pLR6/yaQflz06z92OlRx7A8OUe/XmZzzgse9jTuj4SLh9h5t49zr5Pq6fd5G08hd/i/3VJra5du7q6um7fvr1v377btm2TJKlHjx71U8ztIAjB1kV7CUkDVJsv89cOy/NP0rwYqacf1uN1sCefzzoiF1XT/FhxUBDGgBbgxQhxdPN6+pK7qGp5oUmTJi1evLhv375Llix59tlnBUHh/zgEIQARUf8mwsOBqp8z2LN75CBXejtK6uaLOLyLfVf4O8fk86U0O0oc3VyU0GEWoqE9NbRX8tMaN27cm2++efz48S1btnz11VcKVqKDIAT4hyjQmObif5qJP55nT+2Wm7jQG+2lB/2xdq9FUi5/74ScWUaz2olPtBBxPWSoE1dX1+HDhw8ZMmTQoEHe3t5Kl4PJMgD/ZifS06Hi2aGqp0LFF/fLHX7Xrr2AqTT/kDn9eoF1/F07eb88prl4bpjqmZZIQbgXEydOvHDhwo2XGHzjjTcEQVi6dOlHH30kCEJCQkK9FYMRIUAtVCKNayGOaS5uyGIJp9j0Q+yFCPHpUNFmL5VQqqHv09jnp5mvM81oKz4WJOLweLgfWVlZ4eHh3bt31z8yd+7cuXPnKlIMghDgtkSBHg0SHw0SDxfyT0+zucc1w5uJz7USIz1tKAROF/NFqWxVBosPEH/qJeHagXCfqqurExISvv322/fff1/xaTI6CEKAu4v1Fn7uJeVXSUvPske2yv4u9ExL8T/NRDc7pSszmXINrbnAlpxllytofEsheYjK39ks1llgBRwdHb/++uuHHnpI6UL+gSAEMJSvE73RXpzZTtyazZeeY9MOagY8II5tLj7oL6isZT+ZzGlHLl+RzjZksR5+4qx2Ur8mAqaDghE5OjpOmzZN6Sr+BUEIUDeSQP2bCP2bSNfU0qoMNvuYPPYvPiRYHBos9vC11ETUMtqVz9ddZGsvsGA3YVSI+HEnO29HpcsCqBcIQoB71MiBJoeLk8PFS+X8l0z++hH5fAnv30R85AEhPkD0cFC6PgOU1FBiDtuYxTdfZsFuwuCm4t+DVMG4UhLYGAQhwP0KchWmRQrTIsWcCv5HFv/xPHt2j9zaU+jjL/QJEDt4Cw7S3RdSb9QyHS7kO3L59lx28hrv5is88oA4NwaXjAfbhSAEMJoAF+G5VsJzrcRqmfbm8x25bOpBOeU6j2okdGksdPQROngLikw5yaukw4XsYCHfd4Ufu8pbuQs9/YS32ktdGwtOWAeAzcM/AYDxOUrUJ0DoEyARUbmG/i7g+wv4krPs2T1cEqhdI6Gtp9DKXQj3EEIaCI2MvRG1SE3nS3nqdZ5SzE8V8RPXuIZRBx+ho7c4s60Y11hwtd7JrnArBweHmTNnvvfee0oXcl+0Wq2Dg6n2NyAIAUzL1U4XioLuRE7ZFfzENTpdzJNy+cIUdr6UCwIFuwmBLkKQK3k5Co2dyMuBGtoLHg7kakd2IjW0F/RTcLScyjS8RqYKLV2voeJqllcmXWessJpfKqfsCn6xjMucmjcQwtyF1h7C8+FCu0YCNnvasoSEhFdffVXpKoygUaNGJloyghCgXgW6CIEudOP1nq5W06Vynl3BL5VTYTU/UkhFarpew0pqqExDWkbXa7j+FG8qgdzsBDuRXO2ooT252wsNJKFJA96igdDbn5q4iA+4CpjtCTdydnZu1qyZ0lWYNQQhgMK8HMnLUYj2updBG+e8okLr6oroA7h3lnnQEwAAgJEgCAEAwKYhCAEAwKZZVRCWlZUtWrRI6SpsUWFh4bfffqt0FbYoKyvr559/VroKW3T27Nn169crXYUt0l3X3rjLtKogzMnJWbZsmdJV2KLz589jdayI06dPr1u3TukqbNGxY8c2bNigdBW26MCBA4mJicZdplUFIQAAQF0hCAEAwKYhCAEAwKaZ+wH1V65c2bdvX0hIiCGNtVptfn6+gY3BiNRq9bVr19Dz9a+qqqqkpAQ9X/8qKioqKyvR8/WvrKyspqZm69atBrYfNWrUu+++e+c2Av/fyZvM1Llz5+zsDD1JsFqtNt2JWeEO0PNKQc8rgjEmy7LhqyYwFlmWOecqlaGjOD8/Pycnpzu3sYAgBAAAMB3sIwQAAJuGIAQAAJuGIAQAAJuGIAQAAJtm7odP3OTYsWOrV692dHR88skng4ODb21QWFi4bNmygoKCAQMGPPjgg/rHt2/fvnnz5saNG48fP97Ly6seS7YSycnJK1euFEVx7NixLVu2vOmvJSUlmzZtOnnypJOT06BBg6KionSPb9++PTMzU3fb0dFx3Lhx9Vq0VUhKStq0aZO3t/f48eO9vb1v+uvZs2d3796tvzt48GD913vNmjUHDhwICgp69tlnnZ2d669ia7Fv377169e7u7s//fTT/v7+210evwAACYVJREFUN/01KSnp/Pnz+rtOTk5jx44lovXr1xcWFuoe9Pb2fvzxx+utYOtw4cKFo0ePFhUVDR8+3N3dvdY2GzZs2Llzp7+//4QJExo0aKB78OrVq8uWLbty5Ur//v379OlTpxe1pBHh/v37e/Xq1ahRo6qqqtjY2Ozs7JsaVFVVde7c+ezZs02bNh09evSqVat0j69YsWLs2LFNmzZNSUnp0qVLdXV1vddu2U6dOhUXF+fs7CyKYqdOndLT029qMGfOnJ9++snT01Oj0fTs2XP16tW6x5csWbJu3brMzMzMzMyLFy/Wd92Wb/Xq1SNHjgwKCkpLS+vcuXNlZeVNDfbs2fPJJ59k/j+1Wq17/N13333rrbeaN2++bdu2/v3713vhFm/z5s0DBw4MCAjIy8vr0KFDcXHxTQ0KCgr03f7NN9989913usffe++9HTt26B7Pycmp98It27Vr16Kior7++uuJEyfm5eXV2uaLL7548cUXQ0JCDh482LNnT1mWiai6urpLly4pKSnBwcHjxo1bsWJF3V6YW47HHntszpw5utsjRox4/fXXb2rw3XffxcTEMMY45ytXrmzTpo3u8fDw8NWrV3POGWNRUVE//PBDPVZtDZ544olXX31Vd3vixIkvvPDCTQ2qqqr0t+fPn9+rVy/d7eHDhy9atKh+irRKkZGRK1as0N3u0KHDsmXLbmqwePHiIUOG3PRgZWWlh4fHoUOHOOdqtdrb23vv3r31UK016dat28KFC3W34+PjP/3009u1ZIy1aNFC/zHFxMRs3bq1Pkq0Rowx3QpcEISUlJRbG2g0msDAQF0Pa7XakJCQDRs2cM5/+OGH9u3b6567Zs2a8PDwOr2uJY0Id+/e3bdvX93t+Pj4Xbt23dRg165d8fHxgiDoGiQnJxcVFRUWFqakpOhGyoIg9OnT59Ynwp3dtecdHR31t9Vqtaurq/7u3r17ExIS1q1bp9Vq66FUa3L9+vVTp07Fx8fr7tba80R06dKlhISE77777tq1a7pHTp06RUQxMTFEZG9v36NHD3zn60Sr1e7fv/+uPa+ze/fugoKCGzeBbty48ZNPPtmyZQvHUdp1JAiCbgV+O5mZmfn5+b169SIiSZJ69+6t+2h2797dp08f/co/JSWloKDA8Ne1mCBUq9XFxcX6fSQ+Pj63Dpzz8vL0DTw9PVUqVV5eXl5enkql8vDw0D3euHHj3NzceivbOtzYsbX2vF5mZubnn38+ffp03d3g4GAPD4+ioqLZs2fHxcXpN9yBIfLy8kRR1O/zq/Wr6+7uHh4eXlJSsnbt2pYtW6amphJRfn6+t7e3foWC73xdFRQUyLLs4+Oju9u4ceM7fOe//fbbUaNG6ffCtmnTRqVS5efnT5o0STdYr4+KbUZ+fr6Hh4f+hD767/aN66iGDRs6ODjc4SO7lcVMlpEkSRAE/ahCq9Xa29vf1MbOzk7fQDfEtre312q1+rE2EWk0GpyPqq5UKtWde16noKDgkUcemTVrVteuXXWPvP/++7obs2fPbtOmzYoVK8aPH18PBVsHOzs73Xm8RFGk23x1hw0bNmzYMN3t8ePHz5kzZ+XKlTd+Xron3jhGh7vSrWf1fajRaG73nS8tLV27du2N85X0V6iePn16ixYt/vrrL93wBYzi1u+27p/ixsc557Is3+4jq5XFjAhVKpWPj49+53NOTs6t87gCAgL0v3xzc3M5535+fv7+/owx/a+DnJwcPz+/eivbOgQEBNy554mooKCgV69eo0ePnjp16q1/dXBwiImJuXDhgmkLtS5+fn6CIOi/0nf96nbp0kU3R9ff3//KlSv69QK+83Xl5eXl4OBw43f+dh24cuXKkJCQ6OjoWhfSsmVLfOeNy9/fv6SkpLy8XHdX/9HcuPLPz8/XarV1+s5bTBAS0cCBA9esWUNEnPO1a9cOHDiQiBhjSUlJFRUVugYbNmzQTQpdu3Ztr169XF1dGzZs2L1797Vr1xJRVVXVxo0bBw0apOj7sDz6nieiNWvW6HqeiP7666+SkhIiKi4u7tev37Bhw15//XX9szjn+gm6xcXFe/bsiYiIqN/CLZuLi0vv3r11Pa9Wqzds2KD76lZUVCQlJekny+kaM8Y2bdrUunVrIoqMjPTy8vrzzz+JqLCwcNeuXfqPDAwhCMKAAQN0Pa/VatevX6/rebVanZSUVFNTo2/57bffPvPMM/q7Go1G97kQ0YULF06fPo3vvFGkp6efOXOGiJo2bRoREbFu3ToiKi0tTUxM1H00AwcO3LhxY1VVFRGtXbu2e/futzv0onb3M8Onnp0/f97Pz2/o0KG9e/eOjIwsKSnhnOt+Gpw8eZJzrtVq+/btGx0dPXbs2EaNGu3bt0/3xN27dzdq1Gjs2LFRUVH9+vXTnbwcDJednf3AAw8MGjTo4YcfDg0NLSws1D0uSdLu3bs5588995xKpYr+fwMHDuScV1dXe3p6Dho0aOTIkY0bNx4yZIhWq1XybVig/fv3N2rUaMyYMTExMX379tV1oG4uTFlZGee8T58+ffr0GTNmTERERFhYWHZ2tu6JK1eu9Pb2fvLJJ0NDQ5977jkl34NlOnHihJeX18iRI7t06RIXF1ddXc05v3TpEhHpOzk5Odne3l7/78A5T0lJCQgIGDJkyLBhw9zd3V955RVlqrdkunU4EUVERERHR+vW8y+++OKIESN0DTZt2tSoUaMnnngiIiJC/6Asy/369YuKitKt/HXrJcNZ2NUnrl+/vn37dicnpwcffFA3U5ExduTIkTZt2ugutCHL8l9//XX16tXu3bvfODTOzc3ds2ePt7d3z549dXtcoE7Kysq2b98uSVKfPn308wIOHz7cqlUrV1fXS5cuXb16Vd/YwcFBNzTJyMhITk6uqakJCwuLjIxUpnQLl5+fv2vXLi8vr549e0qSRERVVVXJyckxMTGiKF67du3QoUPXr18PDAzs3LnzjdemOX/+/OHDh4OCgrp06aJc+Rbs6tWrSUlJ7u7uvXr10u01rKmpOXHiRPv27XV3CwsLCwsLw8PD9U/hnKekpJw9e5aIIiMjW7RooVTxluvUqVMajUZ/t23btiqVKisrS6PR6K/+mJWVtW/fvoCAgG7duuknhTHG/vrrr8L/a88OagAAgRgIKjv/YjCBAcKLD9kZE5uma83M8b65+CyEAPCWbQRAmhACkCaEAKQJIQBpQghAmhACkCaEAKQJIQBpQghAmhACkCaEAKRtl2+enzOqmsoAAAAASUVORK5CYII=",
"image/svg+xml": [
"\n",
"\n"
],
"text/html": [
"\n",
"\n"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"using Plots\n",
"\n",
"t = 0:0.005:1\n",
"a = 0.03402\n",
"b = -0.03298\n",
"c = 0.30903\n",
"e = (xi -> a*(xi^2) + b*xi + c).(t)\n",
"\n",
"println(\"a: $a, b: $b, c: $c\")\n",
"plot(t,e)"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"0.48471487360376253"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"0.03298/(0.03402 * 2)"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"1.2167205406993828"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"0.08281/(0.03403 * 2)"
]
},
{
"cell_type": "code",
"execution_count": 15,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"b: -0.08280907372400756; 2a: 0.06805293005671077; argmin: 1.2168333333333334 => 0.30699680718336486\n"
]
}
],
"source": [
"a = 18 /529\n",
"b = -21903/264500\n",
"c = 188207311/529000000\n",
"\n",
"err(t) = c + b*t + a*(t^2)\n",
"m = -b / (2 * a)\n",
"\n",
"println(\"b: $b; 2a: $(2*a); argmin: $m => $(err(1))\")"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Julia 1.9.3",
"language": "julia",
"name": "julia-1.9"
},
"language_info": {
"file_extension": ".jl",
"mimetype": "application/julia",
"name": "julia",
"version": "1.9.3"
},
"orig_nbformat": 4
},
"nbformat": 4,
"nbformat_minor": 2
}