School/phy-4660
1
0
Fork 0
mirror of https://asciireactor.com/otho/phy-4660.git synced 2024-11-22 15:25:06 +00:00
Code Issues Packages Projects Releases Wiki Activity
4ba19d4c00
phy-4660/lag/data/clag_analysis-origbins-noLFerrors.ipynb

888 lines
186 KiB
Plaintext
Raw Normal View History

added lag lab
2017-03-16 04:12:46 +00:00
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Populating the interactive namespace from numpy and matplotlib\n"
]
},
{
"name": "stderr",
"output_type": "stream",
"text": [
"/usr/local/lib/python2.7/dist-packages/IPython/core/magics/pylab.py:161: UserWarning: pylab import has clobbered these variables: ['norm']\n",
"`%matplotlib` prevents importing * from pylab and numpy\n",
" \"\\n`%matplotlib` prevents importing * from pylab and numpy\"\n"
]
},
{
"data": {
"text/plain": [
"<Container object of 3 artists>"
]
},
"execution_count": 1,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAg8AAAFkCAYAAACn/timAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzt3X18VPWd9/9XSAKSRBRwbFEoxFBYg1pNKhD6s+qCqFBy\n1UI35PGzbbjorle73V2vSzK22u7P31W07aRre117tbXXlpL21zWmWnYXoevderNU7mro9sb4kxID\nyO2MiDcQhNyc64/vHOcmM8lM5szMOTPv5+ORR2BmMufkmzPnfM73+/l+viAiIiIiIiIiIiIiIiIi\nIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIjJm+4GhBF//K4/7JCIiIi42Fbg46msxJnj4\neD53SkRERLzju8DefO+EiIiIeMN44A3gy/neEREREXFOWRbf+5PABUD7CK+ZFv4SERGR9BwNf+Vc\nSRbf+0ngPeA/JXl+2iWXXHLkyJEjWdwFERGRgnUYuJY8BBDZ6nmYiUmWvG2E10w7cuQIP/vZz7j8\n8suztBsS78477+S73/1uvnejqKjNc09tnntq89x65ZVXuP322y/F9N4XTPCwBjgObB3thZdffjl1\ndXVZ2g2Jd+GFF6q9c0xtnntq89xTmxeXcVl6zzXATzDTNEVERKSAZCN4WAJMB36chfcWERGRPMvG\nsMVTQGkW3ldERERcIBs9D+Jizc3N+d6FoqM2zz21ee6pzYtLNqdqjqYO6Orq6lKSjYiISBr27NlD\nfX09QD2wJ9fbV8+DiIiIpEXBg4iIiKRFwYOIiIikRcGDiIiIpEXBg4iIiKRFwYOIiIikRcGDiIiI\npEXBg4iIiKRFwYOIiIikRcGDiIiIpEXBg4iIiKRFwYOIiIikRcGDiIiIpEXBg4iIiKRFwYOIiIik\nRcGDiIiIpEXBg4iIiKRFwYOIiIikRcGDiIiIpKUs3zsgIlJIOjrMF8B778GBAzBzJpx3nnmsudl8\niXiZeh5ERBzU3AybN5vvZ8/C3r3wzjvw6qsmmOjogMbGSIAh4kXqeRARcVgoFOKppwK8/no3UMqh\nQ4McPlzL97/vZ/FiX753TyRj6nkQEXFQMBikoaGJ9vaV9PZuBOZy+LAF/Ae33HIDzc1/RSgUyvdu\nimREwYOIiIPuvruNnp4HgMuA1cBKYCvwNAMDv+eRR5ppaGhSACGepuBBRMRBu3d3AwuANuABYCFQ\nEn52HLCInp778fsDedpDkcwpeBARcdDAQCkmWLCDiEQWhIMMEW9S8CAi4qCyskHAAuwgIpFx4SBD\nxJsUPIiIOGj+/FpgF2AHEYkMhYMMEW9S8CAi4qBAwE9NzT3ARcDOJK/aFQ4yRLxJwYOIiIN8Ph87\ndnSyenUVZWX/GXgRGAo/OwTsoKbmXgIBf/52UiRDCh5ERBzm8/no6PgeTzzx78A/c+mlK4BGqqtX\n0NKyiR07OvH5VCxKvEsVJkVEHBS7toWPOXPamDwZDh+GGTNg6VJQ3CBep54HEREHRa9tcd55MHcu\nTJoEc+bAhAla20IKg3oeRESyQKtnSiFTz4OIiIikRcGDiIiIpEXDFiIiWRSbQAkHDsDMmSYfAjS8\nId6k4EFExCGjBQo33AD33mteU1eXt90UyZiCBxERh0T3IuzZA/X1JlCYMSOE3x/gRz/qBkpZtWqQ\n66+vJRDwq96DeJKCBxERB4VCJlB44QUTKNx2Wx8nT4Z4990fAgGghN7eIXp7d7NtW5MKRoknKWFS\nRMQhwWCQhoYm2ttX0tu7BdjMwYNX8+67DwELiayyOQ5YSE/P/fj9gbztr8hYKXgQEXHI3Xe30dPz\nALGBwivh/yeygN27u3OybyJOUvAgIuIQEwgsiHu0lEggEW8cAwOl2d0pkSxQ8CAi4hATCMQHCoOA\nleQnhigrG8zuTolkgYIHERGHmEAgPlCoBXYl+YldzJ9fm92dEskCBQ8iIg4xgUB8oOAH7gG2A0Ph\nx4aAHdTU3Esg4M/hHoo4Q8GDiIhDAgE/NTX3ADuIBApTgXVUVd3BzJm3Ao1UV6+gpWWTpmmKZ6nO\ng4iIQ3w+Hzt2dIbrPKynt7eU6mq7INSzvP66j/p6eOwxVZgUb1PwICLiIJ/Px8aNbe9XmPz852Hn\nTli71pSsnjMHvvxlrW0h3paN4OFS4FvALcBEYC+wFtiThW2JiLhG/NoWc+bA889HAoU1axQoSGFw\nOniYDLwI/BsmeAgCNcBbDm9HRMR11IsgxcLp4OFu4ACmp8F20OFtiIiISB45PduiEegCHgWOY4Yq\nPu/wNkRERCSPnA4eLgO+ALwKLAV+APxP4LMOb0dERETyxOlhi3HAbuCr4f//FrgC+C/ATxP9wJ13\n3smFF14Y81hzczPNGjgUERGho6ODDjsTN+ytt/KbSphstZax2g88BfxF1GNfAO4Fpse9tg7o6urq\nok4TnkVERFK2Z88e6uvrAerJw2xGp3seXgT+JO6xOZigQkSkaMVP4zxwAGbOVL0H8Sang4fvYAq4\nfwWTNDkf+PPwl4hI0WpuhiVLQuHqk9309pbS329Xn/SrTLV4itPBw0vAbcA3gL8FXgP+BugY6YdE\nRApdMBhk0aLV9PQ8AASAEnp7h+jt3c22bU1a50I8JRsLY20FrsJUl5wHbMjCNkREPOXuu9vCgcNC\nIulm44CF9PTcj98fyN/OiaRJq2qKiOTA7t3dwIIkzy4IPy/iDQoeRERyYGCglOQT3MaFnxfxBgUP\nIiI5UFY2CFhJnh0KPy/iDQoeRERyYP78WmBXkmd3hZ8X8QYFDyIiORAI+KmpuQfYAQyFHx0CdlBT\ncy+BgD9/OyeSJgUPIiI54PP52LGjk49/fBMVFSuARiZOXMGkSZuYPr2TtWt9NDZGCkmJuJnTdR5E\nRCQJn8/HCy+0sWcP1NfDP/wD3H47PPggqEq/eImCBxGRHIguT/3OOyEmTQrw53/eDZRy3XWDfPSj\ntTz2mCpNijdo2EJEJAeam2HzZvjRj4IcOtTEO++s5MyZLcBm+voe59//fSUNDU2EQqF876rIqBQ8\niIjkkCpNSiHQsIVIBrRSoqTLVJJMFiAsYPfu9bncHZExUfAgkoHo4MBOguvoUPKbJKdKk1IINGwh\nkqFQKMSaNa2sWrUcaGTVquWsWdOqsWtJSJUmpRCo50EkA1pmWdI1f34t3d27MDkP8VRpUrxBPQ8i\nGVDym6RLlSalEKjnQSQDSn4TSC9x1q406fcHeOGF9fT2llJdPcj119cSCKinSrxBwYNIBpT8JpBe\n4qwJNHxAG7NnQ3m5CTROnIC1azVDR7xBwxYiGVDym9hSTZy1i0U1N5ueiblzzeOvvmp6LTo60BoX\n4nrqeRDJgJLfBMaWONvcDEuWhMLDF9309pbS328PX6hMtbibeh5EMqDkN4GxJc4Gg0EaGppob19J\nb68pU93b+zjt7SpTLe6n4EEkA1pmWcBOnF2Q5NkF4edjaaaOeJmGLUTSEJ9V//vfw8CAj9LSNioq\nYNIkk/y2axdYlpLfCp19PBw4kH7irGbqiJcpeBBJQ6Ks+q4uk1UfCkXGr6GU118f5KmnalmyROPX\nhcrOW/jwh49gEmcTBRCJE2c1U0e8TMGDSJrig4RVqwaZP38mu3a9zP7930KVJouHnSj59tvzgJ1A\nQ4JXJU6cjczUST3gEHEL5TyIpKijA26+OcisWcOT3Do7T7N//zfQ+HVxieQtfBu4l+GJsy8mTZw1\nAcWuJO+smTribgoeRFLU3AxTpnydvr71DA8S3iDxXSckS5gT74skSvqATmATYBJn4RNccMFfJu11\nCgT8XHxx4pk6FRX3cuSIX4m24loKHkRSFAwG+cUvniVxkKDx62IUm7dgqkbCVmAz8EsGBj7E2rW+\nhEGAz+fjD3/oZPXqn1FVVQdcCXyUqqq/orFxHj/7mZJtxb2U8yCSorvvbqO//xISBwkavy5Go+Ut\nzJw5yObNyX/esix+/et
"text/plain": [
"<matplotlib.figure.Figure at 0x7f676e282ad0>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"import numpy as np\n",
"from scipy.stats import norm\n",
"from scipy.stats import lognorm\n",
"import sys\n",
"import getopt\n",
"sys.path.insert(1,\"/usr/local/science/clag/\")\n",
"import clag\n",
"%pylab inline\n",
"\n",
"ref_file=\"lc/1367A.lc\"\n",
"echo_file=\"lc/2246A.lc\"\n",
"\n",
"\n",
"dt = 0.01\n",
"t1, l1, l1e = np.loadtxt(ref_file).T\n",
"errorbar(t1, l1, yerr=l1e, fmt='o')"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"array([ 0.005 , 0.01861938, 0.04473305, 0.06933623, 0.10747115,\n",
" 0.16658029, 0.25819945, 0.40020915, 0.62032418])"
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"\n",
"\n",
"fqL = np.array([0.0049999999, 0.018619375, 0.044733049, 0.069336227, 0.10747115, 0.16658029, \n",
" 0.25819945, 0.40020915, 0.62032418])\n",
"# fqL = np.logspace(np.log10(0.0006),np.log10(1.2),11)\n",
"nfq = len(fqL) - 1\n",
"fqd = 10**(np.log10( (fqL[:-1]*fqL[1:]) )/2.)\n",
"\n",
"\n",
"fqL\n",
"\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
" 1 4.342e-01 5.077e+01 inf -- -5.530e+02 -- 1 1 1 1 1 1 1 1\n",
" 2 7.674e-01 5.065e+01 8.300e+01 -- -4.700e+02 -- 0.653018 0.587019 0.568277 0.567457 0.566281 0.566085 0.565773 0.566163\n",
" 3 3.298e+00 5.043e+01 8.075e+01 -- -3.893e+02 -- 0.414806 0.209135 0.141159 0.13784 0.133393 0.132728 0.131612 0.132761\n",
" 4 1.572e+00 5.010e+01 7.754e+01 -- -3.117e+02 -- 0.322539 -0.0834456 -0.273066 -0.284295 -0.297612 -0.299435 -0.302479 -0.300412\n",
" 5 5.908e-01 4.964e+01 7.386e+01 -- -2.379e+02 -- 0.302357 -0.214604 -0.654658 -0.688754 -0.723838 -0.729279 -0.736418 -0.733888\n",
" 6 3.713e-01 4.877e+01 6.953e+01 -- -1.683e+02 -- 0.284419 -0.200357 -0.96379 -1.05472 -1.13798 -1.15477 -1.17031 -1.16748\n",
" 7 2.709e-01 4.671e+01 6.269e+01 -- -1.056e+02 -- 0.277768 -0.185001 -1.13047 -1.34026 -1.52128 -1.56845 -1.6043 -1.60101\n",
" 8 2.135e-01 4.361e+01 5.281e+01 -- -5.282e+01 -- 0.277012 -0.185189 -1.16375 -1.49463 -1.83211 -1.9477 -2.03737 -2.03476\n",
" 9 1.764e-01 3.767e+01 4.019e+01 -- -1.264e+01 -- 0.282161 -0.184207 -1.17891 -1.53049 -2.01851 -2.24569 -2.46562 -2.46922\n",
" 10 1.508e-01 2.738e+01 2.645e+01 -- 1.382e+01 -- 0.289545 -0.182463 -1.18795 -1.52893 -2.08812 -2.41103 -2.87498 -2.90486\n",
" 11 1.349e-01 1.468e+01 1.390e+01 -- 2.772e+01 -- 0.293547 -0.180898 -1.1897 -1.52803 -2.10944 -2.46526 -3.22702 -3.34288\n",
" 12 1.358e-01 5.365e+00 5.378e+00 -- 3.309e+01 -- 0.295456 -0.179941 -1.19052 -1.52801 -2.11928 -2.48294 -3.46233 -3.79368\n",
" 13 1.868e-01 1.338e+00 1.633e+00 -- 3.473e+01 -- 0.297315 -0.17923 -1.19104 -1.52666 -2.12491 -2.48975 -3.55567 -4.30889\n",
" 14 6.248e-01 2.604e-01 4.517e-01 -- 3.518e+01 -- 0.299091 -0.178645 -1.191 -1.52463 -2.12805 -2.4915 -3.5672 -5.11363\n",
" 15 2.744e+02 2.611e-01 7.022e-02 -- 3.525e+01 -- 0.300248 -0.178286 -1.19075 -1.52307 -2.12961 -2.49161 -3.56337 -8\n",
" 16 2.745e+02 2.782e-01 4.368e-04 -- 3.525e+01 -- 0.300566 -0.178158 -1.19062 -1.52252 -2.13009 -2.49151 -3.5617 -8\n",
" 17 2.745e+02 2.805e-01 6.410e-05 -- 3.525e+01 -- 0.300599 -0.178134 -1.19059 -1.52242 -2.13017 -2.49148 -3.56148 -8\n",
"********************\n",
"0.300599 -0.178134 -1.19059 -1.52242 -2.13017 -2.49148 -3.56148 -8\n",
"0.238931 0.202434 0.232634 0.177249 0.153039 0.132988 0.297259 3285.23\n",
"-0.000915535 -0.00131061 -0.00195128 -0.00497899 -0.0226059 -0.0776175 -0.280541 -0.000206646\n",
"********************\n"
]
}
],
"source": [
"P1 = clag.clag('psd10r', [t1], [l1], [l1e], dt, fqL)\n",
"p1 = np.ones(nfq)\n",
"p1, p1e = clag.optimize(P1, p1)"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"\t### errors for param 0 ###\n",
"+++ 3.525e+01 3.480e+01 3.006e-01 5.395e-01 0.891 +++\n",
"+++ 3.525e+01 3.431e+01 3.006e-01 6.590e-01 1.87 +++\n",
"+++ 3.525e+01 3.458e+01 3.006e-01 5.993e-01 1.34 +++\n",
"+++ 3.525e+01 3.469e+01 3.006e-01 5.694e-01 1.11 +++\n",
"+++ 3.525e+01 3.475e+01 3.006e-01 5.545e-01 0.997 +++\n",
"\t### errors for param 1 ###\n",
"+++ 3.525e+01 3.476e+01 -1.781e-01 2.430e-02 0.973 +++\n",
"+++ 3.525e+01 3.421e+01 -1.781e-01 1.255e-01 2.07 +++\n",
"+++ 3.525e+01 3.451e+01 -1.781e-01 7.491e-02 1.48 +++\n",
"+++ 3.525e+01 3.464e+01 -1.781e-01 4.961e-02 1.21 +++\n",
"+++ 3.525e+01 3.470e+01 -1.781e-01 3.696e-02 1.09 +++\n",
"+++ 3.525e+01 3.473e+01 -1.781e-01 3.063e-02 1.03 +++\n",
"+++ 3.525e+01 3.475e+01 -1.781e-01 2.747e-02 1 +++\n",
"\t### errors for param 2 ###\n",
"+++ 3.525e+01 3.511e+01 -1.191e+00 -1.074e+00 0.276 +++\n",
"+++ 3.525e+01 3.495e+01 -1.191e+00 -1.016e+00 0.598 +++\n",
"+++ 3.525e+01 3.485e+01 -1.191e+00 -9.870e-01 0.8 +++\n",
"+++ 3.525e+01 3.479e+01 -1.191e+00 -9.725e-01 0.91 +++\n",
"+++ 3.525e+01 3.476e+01 -1.191e+00 -9.652e-01 0.968 +++\n",
"+++ 3.525e+01 3.475e+01 -1.191e+00 -9.616e-01 0.997 +++\n",
"\t### errors for param 3 ###\n",
"+++ 3.525e+01 3.482e+01 -1.522e+00 -1.345e+00 0.865 +++\n",
"+++ 3.525e+01 3.432e+01 -1.522e+00 -1.257e+00 1.86 +++\n",
"+++ 3.525e+01 3.459e+01 -1.522e+00 -1.301e+00 1.32 +++\n",
"+++ 3.525e+01 3.471e+01 -1.522e+00 -1.323e+00 1.08 +++\n",
"+++ 3.525e+01 3.476e+01 -1.522e+00 -1.334e+00 0.97 +++\n",
"+++ 3.525e+01 3.473e+01 -1.522e+00 -1.329e+00 1.03 +++\n",
"+++ 3.525e+01 3.475e+01 -1.522e+00 -1.331e+00 0.998 +++\n",
"\t### errors for param 4 ###\n",
"+++ 3.525e+01 3.481e+01 -2.130e+00 -1.977e+00 0.874 +++\n",
"+++ 3.525e+01 3.429e+01 -2.130e+00 -1.901e+00 1.91 +++\n",
"+++ 3.525e+01 3.457e+01 -2.130e+00 -1.939e+00 1.35 +++\n",
"+++ 3.525e+01 3.470e+01 -2.130e+00 -1.958e+00 1.1 +++\n",
"+++ 3.525e+01 3.476e+01 -2.130e+00 -1.968e+00 0.983 +++\n",
"+++ 3.525e+01 3.473e+01 -2.130e+00 -1.963e+00 1.04 +++\n",
"+++ 3.525e+01 3.474e+01 -2.130e+00 -1.965e+00 1.01 +++\n",
"+++ 3.525e+01 3.475e+01 -2.130e+00 -1.966e+00 0.997 +++\n",
"\t### errors for param 5 ###\n",
"+++ 3.525e+01 3.474e+01 -2.491e+00 -2.358e+00 1.01 +++\n",
"+++ 3.525e+01 3.512e+01 -2.491e+00 -2.425e+00 0.263 +++\n",
"+++ 3.525e+01 3.496e+01 -2.491e+00 -2.392e+00 0.579 +++\n",
"+++ 3.525e+01 3.486e+01 -2.491e+00 -2.375e+00 0.781 +++\n",
"+++ 3.525e+01 3.480e+01 -2.491e+00 -2.367e+00 0.893 +++\n",
"+++ 3.525e+01 3.477e+01 -2.491e+00 -2.363e+00 0.952 +++\n",
"+++ 3.525e+01 3.476e+01 -2.491e+00 -2.361e+00 0.982 +++\n",
"+++ 3.525e+01 3.475e+01 -2.491e+00 -2.360e+00 0.997 +++\n",
"\t### errors for param 6 ###\n",
"+++ 3.525e+01 3.507e+01 -3.561e+00 -3.413e+00 0.363 +++\n",
"+++ 3.525e+01 3.484e+01 -3.561e+00 -3.339e+00 0.814 +++\n",
"+++ 3.525e+01 3.468e+01 -3.561e+00 -3.301e+00 1.13 +++\n",
"+++ 3.525e+01 3.477e+01 -3.561e+00 -3.320e+00 0.962 +++\n",
"+++ 3.525e+01 3.473e+01 -3.561e+00 -3.311e+00 1.04 +++\n",
"+++ 3.525e+01 3.475e+01 -3.561e+00 -3.315e+00 1 +++\n",
"\t### errors for param 7 ###\n",
"+++ 3.525e+01 3.524e+01 -8.000e+00 -6.000e+00 0.0178 +++\n",
"+++ 3.525e+01 3.515e+01 -8.000e+00 -5.000e+00 0.187 +++\n",
"+++ 3.525e+01 3.494e+01 -8.000e+00 -4.500e+00 0.618 +++\n",
"+++ 3.525e+01 3.466e+01 -8.000e+00 -4.250e+00 1.18 +++\n",
"+++ 3.525e+01 3.482e+01 -8.000e+00 -4.375e+00 0.851 +++\n",
"+++ 3.525e+01 3.475e+01 -8.000e+00 -4.312e+00 1 +++\n",
"********************\n",
"0.300605 -0.178131 -1.19059 -1.52241 -2.13019 -2.49148 -3.56144 -8\n",
"0.253864 0.205597 0.228999 0.191096 0.1638 0.131949 0.246154 3.6875\n",
"********************\n"
]
}
],
"source": [
"p1, p1e = clag.errors(P1, p1, p1e)"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<Container object of 3 artists>"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgkAAAFrCAYAAABbtho0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAGp1JREFUeJzt3X9s4/d93/HnxZaj5dxOi12RZ+98bLhdaBhJDbLSHzrb\n5TIn2IY26dBOIbFsiLTCQdoOuHUrcMMgzZCAYS2G9uqua3FbD9kQjNQNaDYH2LXFUKW5UVqnimm7\nbmadUT+W5kxmTXvr6kSpVnt/UHeR1I9+8I7fL389HwAhivx8vp+3cJ+TXvx+P9/vFyRJkiRJkiRJ\nkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkvSA/iGwBvwh0AQ+A1zsakWSJKkn3AT+NvA08H7gs8AW\n8K4u1iRJknrQ48BbwHPdLkSSJJ3sHTGONbb39fdjHFOSJPW4M7QON/xqtwuRJEmn83BM4/xz4BmO\nP9Rwbu8hSZLa88beo6PiCAk/DXw38AJw+4g255544onbt28f9bYkSTrGl4EJOhwUogwJZ2gFhI8A\neWD7mLbnbt++zac//WmefvrpCEvqvMuXL3P16tW+HO9BttVu33ban6btSW2Oez/uf7NOca51vr1z\nLcy51vn2Uc611157jY997GNP0tob3zch4WeAIq2Q8CaQ3Hv9DrAT6vD000+TzWYjLKnzxsbGYq25\nk+M9yLba7dtO+9O0PanNce/H/W/WKc61zrd3roU51zrfPuq5FpWHItz2Z4F3AjPA39/3+CLwm4fa\nngM+8YlPfIJz5/pvWcL73ve+vh3vQbbVbt922p+m7Ultjnq/VCpRLBZPXUsvca51vr1zLcy51vn2\nUc21N954g2vXrgFco8N7Es50cmMPIAusr6+v92XqVn/58Ic/zKuvvtrtMjQEnGuKQ7VaJZfLAeSA\naie3Hed1EiRJUh8xJGjo9OvuX/Uf55r6nSFBQ8df3IqLc039zpAgSZKCDAmSJCnIkCBJkoIMCZIk\nKciQIEmSggwJkiQpyJAgSZKCDAmSJCnIkCBJkoIMCZIkKciQIEmSggwJkiQpyJAgSZKCDAmSJCnI\nkCBJkoIMCZIkKciQIEmSggwJkiQpyJAgSZKCDAmSJCnIkCBJkoIMCZIkKciQIEmSggwJkiQpyJAg\nSZKCDAmSJCnIkCBJkoIMCZIkKciQIEmSggwJkiQpyJAgSZKCDAmSJCnIkCBJkoIMCZIkKciQIEmS\nggwJkiQpyJAgSZKCDAmSJCkoypDwAvBZ4MvAW8BHIhxLkiR1WJQh4V3AF4Af2vv+7QjHkiRJHfZw\nhNv+xb2HJEnqQ65JkCRJQYYESZIUZEiQJElBUa5JaNvly5cZGxs78FqxWKRYLHapIkmSekepVKJU\nKh147c6dO5GNdyayLR/0FvC9wKtHvJ8F1tfX18lmszGVJElS/6tWq+RyOYAcUO3ktqPck3AW+Iv7\nvn8P8CzwVeBLEY4rSZI6IMqQMAH8yt7zt4Gf2Hv+KWA2wnElSVIHRBkSPocLIyVJ6lv+EZckSUGG\nBEmSFGRIkCRJQYYESZIUZEiQJElBhgRJkhRkSJAkSUGGBEmSFGRIkCRJQYYESZIUZEiQJElBhgRJ\nkhRkSJAkSUGGBEmSFGRIkCRJQYYESZIUZEiQJElBhgRJkhT0cLcLkKJQKpUolUoA7OzssL29zYUL\nFxgdHQWgWCxSLBa7WaIk9TxDggbS/hBQrVbJ5XKUSiWy2WyXK5Ok/uHhBkmSFGRIkCRJQYYESZIU\nZEiQJElBhgQNrK2tLWZnZ5mengZgenqa2dlZtra2uluYJPUJz27QwGk2mxQKBWq1Go1G497r9Xqd\ner3OzZs3yWQylMtlEolEFyuVpN5mSNBAaTabTE1NsbGxcWSbRqNBo9Hg0qVLVCoVg4IkHcHDDRoo\nhULh2ICwX71ep1AoRFyRJPUvQ4IGxubmJrVara0+tVrNNQqSdARDggbG4uLigTUIp9FoNFhYWIio\nIknqb4YEDYy1tbVY+0nSoDMkaGDs7u7G2k+SBp0hQQNjZGQk1n6SNOgMCRoYExMT99VvcnKyw5VI\n0mAwJGhgzM/Pk0wm2+qTTCaZm5uLqCJJ6m+GBA2MVCpFJpNpq08mkyGVSkVTkCT1OUOCBkq5XCad\nTp+qbTqdZmlpKeKKJKl/GRI0UBKJBJVKhXw+f+Shh2QyST6fZ2VlhfHx8ZgrlKT+YUjQwEkkEiwv\nL7O6usrMzMy9PQvpdJqZmRlWV1dZXl42IEjSCbzBkwZWKpXi+vXrVKtVcrkcN27cIJvNdrssSeob\nUe9J+EFgE/g68OvAcxGPJ0mSOiTKkPBR4CeBReBZ4BZwEzgf4ZiSJKlDogwJPwL8K+A68DvA3wO+\nBHwywjElSVKHRBUSHgGywC8fev2XgamIxpQkSR0U1cLFx4GHgOah178CtHdJPOk+lEolSqUSADs7\nO1y8eJErV64wOjoKQLFYpFgsdrNESep5nt2ggWQIOJ3DYWp7e5sLFy4YpiQBcCai7T4CvAl8P/Af\n9r3+U8D7gb90qH0WWH/++ecZGxs78Ia/pKR43D1VdH193VNFpR61P9jfdefOHW7dugWQA6qdHC+q\nPQl/DKwDH+JgSPgg8JmjOl29etVfTpIkHSH0wfluwI9ClGc3/ATwA8AM8DSt0yH/PPBzEY4pqU1b\nW1vMzs4yPT0NwPT0NLOzs2xtbXW3MEldF+WahBvAY8A8cA74b8Bfo3UapKQuazabFAoFarUajUbj\n3uv1ep16vc7NmzfJZDKUy2USiUQXK5XULVEvXPzZvYekHtJsNpmammJjY+PINo1Gg0ajwaVLl6hU\nKgYFaQh5gydpCBUKhWMDwn71ep1CoRBxRZJ6kSFBGjKbm5vUarW2+tRqNdcoSEPIkCANmcXFxQNr\nEE6j0WiwsLAQUUWSepUhQRoya2trsfaT1L8MCdKQ2d3djbWfpP5lSJCGzMjISKz9JPUvQ4I0ZCYm\nJu6r3+TkZIcrkdTrDAnSkJmfnyeZbO9mrMlkkrm5uYgqktSrDAnSkEmlUmQymbb6ZDIZUqlUNAVJ\n6lmGBGkIlctl0un0qdqm02mWlpYirkhSLzIkSEMokUhQqVTI5/NHHnpIJpPk83lWVlYYHx+PuUJJ\nvcCQIA2pRCLB8vIyq6urzMzM3NuzkE6nmZmZYXV1leXlZQOCNMSivsGTpB6XSqW4fv36vXvS37hx\ng2w22+2yJPUA9yRIkqQgQ4IkSQrycIM0xEqlEqVSCYCdnR0uXrzIlStXGB0dBaBYLFIsFrtZoqQu\nMiRIQ8wQIOk4Hm6QJElBhgRJkhRkSJAkSUGGBEmSFGRIkCRJQYYESZIUZEiQJElBhgRJkhRkSJAk\nSUGGBEmSFGRIkCRJQYYESZIUZEiQJElBhgRJkhRkSJAkSUGGBEmSFGRIkCRJQYYESZIUZEiQJElB\nhgRJkhRkSJAkSUGGBEmSFGRIkCRJQYYESZIUFFVI+EfACvA14A8iGkOSJEUoqpAwAiwB/yKi7UuS\npIg9HNF2X977+vGIti9JkiLmmgRJkhQU1Z4ESeq4UqlEqVQCYGdnh+3tbS5cuMDo6CgAxWKRYrHY\nzRKlgdJOSHgZmD+hzXcC1fuuRpKOsT8EVKtVcrkcpVKJbDbb5cqkwdROSPhp4N+e0Gb7AWrh8uXL\njI2NHXjNTwaSJLXs35t21507dyIbr52Q8NW9R2SuXr3qJwJJko4Q+uB8d69aFKJak/AU8O69rw8B\n3wGcAb4IvBnRmJIkqYOiOrthgdbahJeBs8AXgHUgmqgjaWhsbW0xOzvL9PQ0ANPT08zOzrK1tdXd\nwqQBFNWehI/jNRIkdVCz2aRQKFCr1Wg0Gvder9fr1Ot1bt68SSaToVwuk0gkulipNDg8BVJSz2s2\nm0xNTbGxsXFkm0ajQaPR4NKlS1QqFYOC1AFeTElSzysUCscGhP3q9TqFQiHiiqThYEiQ1NM2Nzep\n1Wpt9anVaq5RkDrAkCC
"text/plain": [
"<matplotlib.figure.Figure at 0x7f67900fa210>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"xscale('log'); ylim(-4,2)\n",
"errorbar(fqd, p1, yerr=p1e, fmt='o', ms=10, color=\"black\")\n"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<Container object of 3 artists>"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAg8AAAFkCAYAAACn/timAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzs3Xt4XPV97/v3WDdbkuUb8gXJ2I5sOcgxwTYVsg2BVL7g\nZJcAxQkDJJEPKUpKyaal284+3XQ7fXqe06O2NBuSbuy9va1QwiRgWuKGxDZKapCNwcGOwZWohYVs\nJGEj+W7drNucP5bWaEaa29Ksuerzeh49YM3MWmtGM7O+6/f7/r5fEBERERERERERERERERERERER\nERERERERERERERERERERERERERERERERERERGbNTwKCfnx/G8ZhEREQkgc0AZnr9lGMED1+I50GJ\niIhI8vgB0BDvgxAREZHkkAmcA74X7wMRERER+6RHcdv3AFOA6iD3mTP0IyIiItacGfqJOUcUt70X\n6AG+EuD2Oddff/0nn3zySRQPQUREJGW1Ar9HHAKIaI08zMNIlrw3yH3mfPLJJ7zwwgvceOONUToM\nGemJJ57gBz/4QbwPY1zRax57es1jT695bH3wwQc8/PDDBRij9ykTPGwCPgVeC3XHG2+8keXLl0fp\nMGSkqVOn6vWOMb3msafXPPb0mo8vE6K0zU3AjzGWaYqIiEgKiUbwsAYoBP5PFLYtIiIicRaNaYt9\nQFoUtisiIiIJIBojD5LAnE5nvA9h3NFrHnt6zWNPr/n4Es2lmqEsB44cOXJESTYiIiIWHD16lBUr\nVgCsAI7Gev8aeRARERFLFDyIiIiIJQoeRERExBIFDyIiImKJggcRERGxRMGDiIiIWKLgQURERCxR\n8CAiIiKWKHgQERERSxQ8iIiIiCUKHkRERMQSBQ8iIiJiiYIHERERsUTBg4iIiFii4EFEREQsSY/3\nAYiIpBLXcReuf3fReqWVj698THdfN1npWVzrv8akjEnckHcDBXkFOD/nxLnUGe/DFRkTBQ8iIjZy\nLnWyZvYaNm/dzLlD5zh35hzXuEZ/ej/XzbqOm269iaqtVeTn58f7UEXGTMGDiIiN2traWLVhFY2f\nbYSLwJehv7AfHHB68DTVrdXU3lXLoT2HFEBI0lLOg4iIjbZ8fwuNyxrhFFAOzAUcQzdOMP7duKyR\nzVs3x+sQRSKm4EFExEaHjx2GQqAd47/+FAzdTyRJKXgQEbFRP8YUhefHnwlD9xNJUgoeRERslE46\nuBn+8Wdw6H4iSUrBg4iIjUpvLoUWIB/jv/60Dt1PJEkpeBARsVHV1iqKflcE84FfA83A4NCNg8a/\ni35XRNXWqngdokjEFDyIiNgoPz+fQ3sOUTG1ghum3QCvQfr2dPjfMO+1eVRkVWiZpiQ9TbqJiNjI\nU2Hypla65neR05fjqTDZmdHJ+3nv80jNI6owKUlNwYOIiI2cSxUUSOrTtIWIiIhYouBBRERELFHw\nICISRa7jLtZvW8/cDXPJXZJLZkkmuUtymbthLuu3rcd13BXvQxSxTDkPIiI2MZMlAXr6ezh9+TRz\nHHP47TO/pev2LrgVcEDfYB+drZ1kbc9izX1r4nvQImOgkQcREZs4lzrZsWYHMw7M4OQPT9LwbANH\nqo4YgYMaZEkK0ciDiIhNPO24lzXCBsABHS90BG+QVaMGWZJ8FDyIiNjE0457rtcvJ6AGWZJyNG0h\nImITTztub2qQJSlIwYOIiE087bi9qUGWpCAFDyIiNvG04/a2GqNB1seoQZakDAUPIiI28bTj9pYD\nbASOQO6Pc+FFWLBngRpkSVJT8CAiYhNPO+6RbbgvQHZ3Niv+ywqKHy9m4Z8s5Pxt53mk5hEViZKk\nFI1MnQLg/wPuAiYBDcAjwNEo7EtEJGHUnK2h6NEirr16jYuHLtLr7iXTkcm0G6ZR8mclVKyqUNMs\nSQl2Bw/TgIMYM3x3AW1AEXDJ5v2IiCQcT0fNyngfiUh02R08bAFOY4w0mD62eR8iIiISR3bnPNwN\nHAFeBj7FmKr4ls37EBERkTiyO3j4DPAd4ASwDvifwDPAN2zej4iIiMSJ3dMWE4DDwH8b+vd7wOeA\nbwPP+3vAE088wdSpU31+53Q6cTqVVCQiIuJyuXC5fFflXLoU31TCQBXXx+oUsA941Ot33wH+gtFF\nW5cDR44cOcLy5cttPgwRkcTir133vCnzmJg+EQDn55xaiSFhO3r0KCtWrABYQRxWM9o98nAQ+OyI\n3xVjBBUiIuOWc6mTNbPXsHnrZt549w2aLjXRN7WPO265g6qtVSoWJUnF7uDhH4C3gP+KkTRZCvzR\n0I+IyLjlr11302ATTa1N1N5Vq2qTklTsTph8F7gXcALHMaYr/jOgEmoiMq75tOs2J4wnAHOhcVkj\nm7dujuPRiVgTjQqTrw39iIjIkMPHDsPaADcWwOGawzE9HpFIqLeFiEgM+G3XbZowdLtIklDwICIS\nA37bdZsGh24XSRIKHkREYsBvu25T69DtIklCwYOISAwEbNfdDJPenMQnN33C3a671aJbkoKCBxGR\nGMjPz+fQnkNUZFVww7/eAM9B+vZ02Asz82Zy/fvXs2PNDhWKkqSgSTYRkRgwK0xeW3qN9tfb4cvQ\nX2gkUZ4ePE11a7XqPUjS0MiDiEgMOJc62e3czfXHr6f7C92q9yBJTSMPIhFQvwKxSvUeJBUoeBCJ\ngPoViFWq9yCpQMGDSATUr0Cs8tR78BdAqN6DJAnlPIhEQP0KxCrVe5BUoOBBJAKHjx2GwgA3Fgzd\nLuIlWL2Hot8VUbW1Ko5HJxIeBQ8iEdD8tYCROLt+23rmbphL7pJcMksyyV2Sy9wNc1m/bb1P4aea\nszUUPVpEYWshObtyyPhZBjm7cihsLaTo0SJqztbE8ZmIhEeTayIR0Py1AJTPKuep7U/RsqwFbgUc\n0DfYR2drJ1nbs1hz3xrPfZ1LjRU4rlUuqt+qpv7Vei5+fJFPT3/KxWcuUv9qPdX3VFOxqkIrdSRh\naeRBJAKavxYYW+5L+axyGrc30lLQQufGTvq+1kfn/Z20FLTQuL2RNbPXjHqMSKJQ8CASAc1fC4wt\n90XJtpLMFDyIRMCcv55+ajoTXpgAz2H87IWmK0181vnZUXPekjpcx13c7bqb01dOW859UbKtJDNN\nyIpY4DruO0/d6+4l05FJ3qw8siZk0f3lbuOE4IDBwUEutF4whqDv0xB0KjKLhC36i0WWc1+UbCvJ\nTCMPIhYEmqc+c+EM3XeoX8F409bWxsq7VnI577Ll3BdPsq0/SraVBKfgQSRMruMubn7wZv/z1F1o\nCHoc8uQtrAN+zejcl48D576ESra9kndF012SsBTaioSpfFY55+rOwW1+bnSgIehxyNPkygFsBA4C\nbw79exCm9E3h0G/9lyiv2lpF7V21NHY3QhNwbuhxvZDbn8u+ffu48cYbY/dkRCzQyINImLZ8fwt9\n2X3+gwQ3GoIeh3zyFnIwRiAeAh4EHob+yf08UvOI3xGE/Px8dv/TbibXToaSocc8CHwDOtZ38AcP\n/wHt7e0xeiYi1ih4EAnT4WOHIQ3/QUI+qvcwDoXKW5iXN4/dzt1+iz25jrtY/+R6rq67qlwZSToK\nHkTC1E9/4CBhNbAX+BjVe/BipWxzMoqkSJhzqZO8K3nKlZGkpOBBJEzppMMq/CfGnYf0a+k8wAMs\n2LMAXoQFexZQkVUxrttyp3oVxUiLhGm5piQrTcSKhKn05lLqL9aPToxzA9mQVpRGQ2kDC+9cSMbl\nDOZNmcf59PM8UvMIzs85x2WfAp8qiiZzWB5jWH7nj3bG6/DGzKz38e7L73Lp0iV4DSNoyIQJEycw\ndd5UT5MrZ37gv7t6o0iy0jtTJEye7PhljbAG4yQ4CLQaV5mHdo/fEYZAPKsR/CmAwzXJOSxvNsK6\nsOyCsfpmaHUFrbDgdws45ArvvVB6cyn1LfW+wZVJuTKSwDRtIRKm/Px8Du05REVWRcCpiVSf47cq\nVYfl7epLod4okqw08iA
"text/plain": [
"<matplotlib.figure.Figure at 0x7f676c054a10>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"t2, l2, l2e = np.loadtxt(echo_file).T\n",
"errorbar(t1, l1, yerr=l1e, fmt='o', color=\"green\")\n",
"errorbar(t2, l2, yerr=l2e, fmt='o', color=\"black\")"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
" 1 4.379e-01 6.748e+01 inf -- -3.302e+02 -- 1 1 1 1 1 1 1 1\n",
" 2 7.796e-01 6.679e+01 7.758e+01 -- -2.526e+02 -- 0.618991 0.579498 0.564066 0.564121 0.56476 0.563682 0.563867 0.56205\n",
" 3 3.540e+00 6.595e+01 7.576e+01 -- -1.769e+02 -- 0.314779 0.183356 0.127531 0.127909 0.129377 0.127333 0.127748 0.123875\n",
" 4 1.452e+00 6.481e+01 7.321e+01 -- -1.036e+02 -- 0.145822 -0.154058 -0.309225 -0.308235 -0.306082 -0.309022 -0.308575 -0.314603\n",
" 5 5.873e-01 6.313e+01 7.037e+01 -- -3.327e+01 -- 0.101594 -0.377738 -0.742409 -0.742747 -0.740971 -0.745083 -0.745333 -0.753741\n",
" 6 3.716e-01 6.044e+01 6.710e+01 -- 3.383e+01 -- 0.0846307 -0.469976 -1.15278 -1.17038 -1.1732 -1.18022 -1.18308 -1.19385\n",
" 7 3.359e-01 5.601e+01 6.209e+01 -- 9.591e+01 -- 0.0659456 -0.50867 -1.48798 -1.574 -1.5971 -1.61323 -1.62275 -1.63479\n",
" 8 5.149e-01 4.857e+01 5.399e+01 -- 1.499e+02 -- 0.0437973 -0.538092 -1.68322 -1.90776 -2.00001 -2.04069 -2.06646 -2.07709\n",
" 9 7.563e-01 3.652e+01 4.228e+01 -- 1.922e+02 -- 0.0212475 -0.55813 -1.75872 -2.09423 -2.35209 -2.45331 -2.51958 -2.52359\n",
" 10 1.974e+00 2.058e+01 2.756e+01 -- 2.197e+02 -- 0.00517809 -0.569653 -1.80128 -2.12815 -2.59485 -2.83321 -2.9937 -2.97691\n",
" 11 1.071e+00 7.198e+00 1.293e+01 -- 2.327e+02 -- -0.00504406 -0.575918 -1.82651 -2.12905 -2.68784 -3.15505 -3.51079 -3.42807\n",
" 12 2.091e-01 1.531e+00 3.715e+00 -- 2.364e+02 -- -0.0104477 -0.580461 -1.83547 -2.13193 -2.69807 -3.39738 -4.11475 -3.83333\n",
" 13 5.653e-01 3.932e-01 6.297e-01 -- 2.370e+02 -- -0.0124146 -0.583395 -1.83764 -2.13329 -2.69789 -3.54657 -4.97511 -4.11701\n",
" 14 1.852e+02 1.372e-01 6.354e-02 -- 2.371e+02 -- -0.0129899 -0.584709 -1.83799 -2.13471 -2.69697 -3.6137 -7.78741 -4.23122\n",
" 15 2.921e+02 7.159e-02 2.904e-03 -- 2.371e+02 -- -0.0130755 -0.585195 -1.83773 -2.13556 -2.696 -3.63655 -8 -4.24469\n",
" 16 2.917e+02 5.374e-02 4.141e-04 -- 2.371e+02 -- -0.0130579 -0.585335 -1.83758 -2.13596 -2.69543 -3.64347 -8 -4.24419\n",
" 17 2.916e+02 4.858e-02 9.106e-05 -- 2.371e+02 -- -0.0130506 -0.585376 -1.83754 -2.13607 -2.69524 -3.64553 -8 -4.24379\n",
"********************\n",
"-0.0130506 -0.585376 -1.83754 -2.13607 -2.69524 -3.64553 -8 -4.24379\n",
"0.269171 0.232445 0.302922 0.233996 0.234912 0.539942 6377.59 0.989373\n",
"-0.000636424 -0.00149059 -0.00259185 -0.00860622 -0.0213922 -0.0485818 -5.7966e-05 -0.0291693\n",
"********************\n"
]
}
],
"source": [
"P2 = clag.clag('psd10r', [t2], [l2], [l2e], dt, fqL)\n",
"p2 = np.ones(nfq)\n",
"p2, p2e = clag.optimize(P2, p2)"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {
"collapsed": false,
"scrolled": true
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"\t### errors for param 0 ###\n",
"+++ 2.371e+02 2.366e+02 -1.305e-02 2.561e-01 0.96 +++\n",
"+++ 2.371e+02 2.361e+02 -1.305e-02 3.907e-01 2 +++\n",
"+++ 2.371e+02 2.364e+02 -1.305e-02 3.234e-01 1.44 +++\n",
"+++ 2.371e+02 2.365e+02 -1.305e-02 2.898e-01 1.19 +++\n",
"+++ 2.371e+02 2.365e+02 -1.305e-02 2.729e-01 1.07 +++\n",
"+++ 2.371e+02 2.366e+02 -1.305e-02 2.645e-01 1.02 +++\n",
"+++ 2.371e+02 2.366e+02 -1.305e-02 2.603e-01 0.987 +++\n",
"+++ 2.371e+02 2.366e+02 -1.305e-02 2.624e-01 1 +++\n",
"\t### errors for param 1 ###\n",
"+++ 2.371e+02 2.366e+02 -5.854e-01 -3.529e-01 0.949 +++\n",
"+++ 2.371e+02 2.361e+02 -5.854e-01 -2.367e-01 1.99 +++\n",
"+++ 2.371e+02 2.364e+02 -5.854e-01 -2.948e-01 1.43 +++\n",
"+++ 2.371e+02 2.365e+02 -5.854e-01 -3.239e-01 1.18 +++\n",
"+++ 2.371e+02 2.365e+02 -5.854e-01 -3.384e-01 1.06 +++\n",
"+++ 2.371e+02 2.366e+02 -5.854e-01 -3.457e-01 1 +++\n",
"\t### errors for param 2 ###\n",
"+++ 2.371e+02 2.369e+02 -1.838e+00 -1.686e+00 0.273 +++\n",
"+++ 2.371e+02 2.368e+02 -1.838e+00 -1.610e+00 0.6 +++\n",
"+++ 2.371e+02 2.367e+02 -1.838e+00 -1.572e+00 0.806 +++\n",
"+++ 2.371e+02 2.366e+02 -1.838e+00 -1.554e+00 0.92 +++\n",
"+++ 2.371e+02 2.366e+02 -1.838e+00 -1.544e+00 0.979 +++\n",
"+++ 2.371e+02 2.366e+02 -1.838e+00 -1.539e+00 1.01 +++\n",
"\t### errors for param 3 ###\n",
"+++ 2.371e+02 2.369e+02 -2.136e+00 -2.019e+00 0.267 +++\n",
"+++ 2.371e+02 2.368e+02 -2.136e+00 -1.961e+00 0.589 +++\n",
"+++ 2.371e+02 2.367e+02 -2.136e+00 -1.931e+00 0.793 +++\n",
"+++ 2.371e+02 2.366e+02 -2.136e+00 -1.917e+00 0.905 +++\n",
"+++ 2.371e+02 2.366e+02 -2.136e+00 -1.909e+00 0.964 +++\n",
"+++ 2.371e+02 2.366e+02 -2.136e+00 -1.906e+00 0.994 +++\n",
"\t### errors for param 4 ###\n",
"+++ 2.371e+02 2.370e+02 -2.695e+00 -2.578e+00 0.256 +++\n",
"+++ 2.371e+02 2.368e+02 -2.695e+00 -2.519e+00 0.577 +++\n",
"+++ 2.371e+02 2.367e+02 -2.695e+00 -2.490e+00 0.787 +++\n",
"+++ 2.371e+02 2.366e+02 -2.695e+00 -2.475e+00 0.904 +++\n",
"+++ 2.371e+02 2.366e+02 -2.695e+00 -2.468e+00 0.965 +++\n",
"+++ 2.371e+02 2.366e+02 -2.695e+00 -2.464e+00 0.997 +++\n",
"\t### errors for param 5 ###\n",
"+++ 2.371e+02 2.369e+02 -3.646e+00 -3.376e+00 0.301 +++\n",
"+++ 2.371e+02 2.367e+02 -3.646e+00 -3.241e+00 0.783 +++\n",
"+++ 2.371e+02 2.365e+02 -3.646e+00 -3.173e+00 1.15 +++\n",
"+++ 2.371e+02 2.366e+02 -3.646e+00 -3.207e+00 0.953 +++\n",
"+++ 2.371e+02 2.366e+02 -3.646e+00 -3.190e+00 1.05 +++\n",
"+++ 2.371e+02 2.366e+02 -3.646e+00 -3.199e+00 1 +++\n",
"\t### errors for param 6 ###\n",
"+++ 2.371e+02 2.371e+02 -8.000e+00 -5.000e+00 0.0545 +++\n",
"+++ 2.371e+02 2.353e+02 -8.000e+00 -3.500e+00 3.58 +++\n",
"+++ 2.371e+02 2.369e+02 -8.000e+00 -4.250e+00 0.39 +++\n",
"+++ 2.371e+02 2.365e+02 -8.000e+00 -3.875e+00 1.18 +++\n",
"+++ 2.371e+02 2.367e+02 -8.000e+00 -4.062e+00 0.674 +++\n",
"+++ 2.371e+02 2.366e+02 -8.000e+00 -3.969e+00 0.891 +++\n",
"+++ 2.371e+02 2.366e+02 -8.000e+00 -3.922e+00 1.03 +++\n",
"+++ 2.371e+02 2.366e+02 -8.000e+00 -3.945e+00 0.956 +++\n",
"+++ 2.371e+02 2.366e+02 -8.000e+00 -3.934e+00 0.99 +++\n",
"\t### errors for param 7 ###\n",
"+++ 2.371e+02 2.367e+02 -4.244e+00 -3.749e+00 0.707 +++\n",
"+++ 2.371e+02 2.359e+02 -4.244e+00 -3.502e+00 2.45 +++\n",
"+++ 2.371e+02 2.364e+02 -4.244e+00 -3.625e+00 1.37 +++\n",
"+++ 2.371e+02 2.366e+02 -4.244e+00 -3.687e+00 0.998 +++\n",
"********************\n",
"-0.0130485 -0.585388 -1.83753 -2.1361 -2.69519 -3.64613 -8 -4.24367\n",
"0.275479 0.239709 0.298186 0.230343 0.23123 0.447579 2 0.556379\n",
"********************\n"
]
}
],
"source": [
"p2, p2e = clag.errors(P2, p2, p2e)"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<Container object of 3 artists>"
]
},
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgkAAAFrCAYAAABbtho0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzt3X9w3Hd95/Gn44i4/DQm590Eg7dsz6wpTnIScmOLuEqP\ndkraAgdXdbdlbpDIhdL0Or47Mudrz7qMfNOjJVNSoJRxE9H2gJV808KRmfiAUrlJZYUTEpC49lJu\nJZmYeFcNxoFSlBOx74/vKpbNV5ZW2u/+0Pf5mNmRvPv5fD8fT76RX/p+P9/PGyRJkiRJkiRJkiRJ\nkiRJkiRJkiRJkiRJkiRJkiRJkiRJkrRG/xkYB74LlIFPAzsaOiNJktQUjgL/BtgJ3AQ8BMwAL2zg\nnCRJUhO6HrgAvLHRE5EkScu7po5jba58PVfHMSVJUpPbQHC74W8aPRFJkrQy19ZpnI8AP8nVbzXc\nUHlJkqTqnK28aqoeIeHDwC8C+4Cnlmhzw4033vjUU08t9bEkSbqKbwGd1DgoRBkSNhAEhLcC3cDp\nq7S94amnnuITn/gEO3fujHBKtbd//37uv//+lhxvLceqtm817VfSdrk2V/u83v/NasVzrfbtPdfC\nea7Vvn2U59qpU6d45zvf+UqCq/EtExL+CMgRhITvA8nK++eBubAOO3fupL29PcIp1d7mzZvrOuda\njreWY1Xbt5r2K2m7XJurfV7v/2a14rlW+/aea+E812rfPupzLSobIzz2Q8B1QC/wHxe9vgF87Yq2\nNwDvec973sMNN7TesoRdu3a17HhrOVa1fatpv5K2y7VZ6vN8Pk8ul1vxXJqJ51rt23uuhfNcq337\nqM61s2fPcvjwYYDD1PhKwoZaHmwN2oGJiYmJlkzdai1vectb+OxnP9voaSgGPNdUD5OTk3R0dAB0\nAJO1PHY990mQJEktxJCg2GnVy79qPZ5ranWGBMWOP7hVL55ranWGBEmSFMqQIEmSQhkSJElSKEOC\nJEkKZUiQJEmhDAmSJCmUIUGSJIUyJEiSpFCGBEmSFMqQIEmSQhkSJElSKEOCJEkKZUiQJEmhDAmS\nJCmUIUGSJIUyJEiSpFCGBEmSFMqQIEmSQhkSJElSKEOCJEkKZUiQJEmhDAmSJCmUIUGSJIUyJEiS\npFCGBEmSFCrKkLAPeAj4FnABeGuEY0mSpBqLMiS8EPgKcHflzxcjHEuSJNXYtREe+39XXpIkqQW5\nJkGSJIUyJEiSpFCGBEmSFCrKNQlV279/P5s3b77svVwuRy6Xa9CMJElqHvl8nnw+f9l758+fj2y8\nDZEd+XIXgLcBn13i83ZgYmJigvb29jpNSZKk1jc5OUlHRwdABzBZy2NHeSXhRcA/X/Tn1wC3AN8G\nnoxwXEmSVANRhoRO4K8r318E/qDy/Z8CfRGOK0mSaiDKkHAMF0ZKktSy/EdckiSFMiRIkqRQhgRJ\nkhTKkCBJkkIZEiRJUihDgiRJCmVIkCRJoQwJkiQplCFBkiSFMiRIkqRQhgRJkhTKkCBJkkJFWeBJ\napj8E3nyJ/IAzP1wjtPPnGb7y7az6dpNAORenyO3K9fIKUpS0zMkaF3K7boUAibPTtJxuIP8O/K0\n39De4JlJUuvwdoMkSQplSNC6NTMzQ9/dffS8vQc+BT1v76Hv7j5mZmYaPTVJagnebtC6Uy6Xyd6Z\npXCuQOl1Jfj54P0iRYpnihz9taNktmQYemCIRCLR2MlKUhMzJGhdKZfL7L1jL1O3TsEbQhpsg9K2\nEqXZEl13dDH68KhBQZKW4O0GrSvZO7NBQNi6TMOtULy1SPbObF3mJUmtyJCgdWN6eprCucLyAWHB\nViicK7hGQZKWYEjQunHovkPBGoQqlHaWGLhvIKIZSVJrMyRo3Rh/fBy2VdlpG4x/bTyS+UhSqzMk\naN2Yf26++k4bYP7CKvpJUgz4dIPWjbaNbZf+8ETltZRdlddFaLum7SoNJSm+DAlaNzpv6uTEmRPB\nLYeFELCcM7D75t0Rz0ySWpO3G7Ru9N/TT/Jksqo+yVNJDr7vYEQzkqTWZkjQupFKpchsycDsCjvM\nQmZLhlQqFeW0JKllRR0SfgOYBn4AfBl4Y8TjKeaGHhgi/Vg6CArfAT4DfBT4cOXrZyrvz0L6sTTD\nDw43brKS1OSiXJPwK8AHgfcCo8CvA0eB1wFPRjiuYiyRSPCZP/sMb+x+I9/77ve4MH/h8gazcM3f\nXcNLXvoSPnPsM2zdutKdlyQpfqK8kvAfgAeAQeDrwL8nCAfvjXBMxVy5XOatb30rz3z7mR8NCBUX\n5i/wzLef4W1vexvlcrnOM5Sk1hFVSHgB0A58/or3Pw/sjWhMiWw2y9TU1IraFotFsllrN0jSUqIK\nCdcDG4Erf02bBapbfi6t0PT0NIVCoao+hYK1GyRpKe6ToHXj0KFDlEpV1m4olRgYGGBwcDCiWTW3\n/BN58ifyAMz9cI7Tz5xm+8u2s+naTQDkXp8jtyvXyClKaqCoQsLTwHNA4or3E8DZpTrt37+fzZs3\nX/ZeLpcjl/OHlJY3Pr66Ggyr7bce5Hbl2POSPQx8YIBHJh+heK7Ic1ueY1/7Pvrv6ffxUKnJ5PN5\n8vn8Ze+dP38+svE2RHZkeAyYAO5e9N5J4NPA71zRth2YmJiYoL29PcIpaT3LZDJ8/etfr7rfa1/7\n2qpvU6wH5XKZ7J1ZCucKQfXMxcWxzkDyZJLMlgxDDwyRSFyZ9yU1i8nJSTo6OgA6gMlaHjvK2w1/\nAPwPgv0RHgPuIvgx9LEIx1SMtbWtrgbDavu1snK5zN479jJ16xS8IaTBNihtK1GaLdF1RxejD48a\nFKQYivIRyCPAfqAf+ArBRkp34B4JikhnZ+eq+u3eHb/aDdk7s0FAWG6biK1QvLVI9k6fApHiKOod\nF/8Y+HFgE9AJ/G3E4ynG+vv7SSarrN2QTHLwYLxqN0xPT1M4V1g+ICzYCoVzPgUixZG1G7RupFIp\nMplMVX0ymfjVbjh036FgDUIVSjtLDNw3ENGMJDUrQ4LWlaGhIdLp9IraptNphofjV7th/PHxyxcp\nrsQ2GP9afJ8CkeLKkKB1JZFIMDo6Snd395K3HpLJJN3d3Rw/fjyWtRvmn5uvvtMGmL+win6SWpoh\nQetOIpFgZGSEsbExent7n7+ykE6n6e3tZWxsjJGRkVgGBIC2jat4muMitF0Tv6dApLhzx0WtW6lU\nisHBweefIT5y5Ij7cACdN3Vy4syJ6m45nIHdN8fvKRAp7rySIMVM/z39JE9W+RTIqSQH3xevp0Ak\nGRKk2EmlUmS2ZIJyaysxC5kt8XsKRJK3G7ROLd7ffG5ujh07dnDgwAE2baoULop5TZChB4bouqOL\n4q3Fq++XMAvpx9IMH43fUyCSDAlap+IeApaTSCQYfXg0qN3w1QKlnZXaDRuAiwS1G04FtRuGjw7H\ndpGnFHeGBCmmEokEd/3qXTz44IO0fbGNp7/zND949gf82HU/xvUvv54d6R28+1ffbUCQYsyQIMXZ\n6+GF734ht3ALs+dm+dLol7ip6ya2btn6/OeS4suFi1KM5Xbl+NCeD3H9F67n6Q8+DXl4+oNPc/0X\nrudDez5Ebpe3bKQ480qCFFPlcplsNkuhUKBUulTLoVgsUiwWOXr0KJlMhqGhIctESzFlSJBiqFwu\ns3fvXqamppZsUyqVKJVKdHV1MTo6alCQYsjbDVIMZbPZqwaExYrFItlsNuIZSWpGhgQpZqanpykU\nClX1KRQKzMzMRDMhSU3LkCDFzKFDhy5bg7ASpVKJgYGBiGYkqVkZEqSYGR8fr2s/Sa3LkCDFzPz8\nfF37SWpdhgQpZtra2uraT1LrMiRIMdPZ2bmqfrt3767xTCQ1O0OCFDP9/f0kk8mq+iSTSQ4ePBjR\njCQ1K0OCFDOpVIpMJlNVn0wmQyqVimZCkpqWIUGKoaGhIdLp9IraptNphoeHI56RpGZkSJBiKJFI\nMDo6Snd395K3HpLJJN3
"text/plain": [
"<matplotlib.figure.Figure at 0x7f676b9aaa90>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"xscale('log'); ylim(-6,2)\n",
"errorbar(fqd, p1, yerr=p1e, fmt='o', ms=10, color=\"green\")\n",
"errorbar(fqd, p2, yerr=p2e, fmt='o', ms=10, color=\"black\")"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
" 1 5.670e+03 9.675e+00 inf -- 2.919e+02 -- -0.156222 -0.68176 -1.81406 -2.12926 -2.71269 -3.36881 -6.08072 -6.42184 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1\n",
" 3 1.047e+03 1.151e+01 2.879e+00 -- 2.948e+02 -- -0.118762 -0.639838 -1.76688 -2.08895 -2.66889 -3.26353 -5.78072 -6.12184 0.0424515 0.105123 0.138683 0.0948799 0.140121 0.0393153 2.81332 0.248074\n",
" 5 4.807e+02 1.352e+01 2.643e+00 -- 2.974e+02 -- -0.0866187 -0.605762 -1.72887 -2.05583 -2.63238 -3.19092 -6.08072 -5.82184 0.000648562 0.108523 0.164491 0.0915061 0.170413 0.00913963 2.85855 1.07978\n",
" 7 1.906e+02 1.570e+01 2.338e+00 -- 2.998e+02 -- -0.0590998 -0.577513 -1.69772 -2.02818 -2.60158 -3.1372 -6.38072 -5.52184 -0.0305296 0.110846 0.182461 0.0892265 0.194212 -0.00858619 2.71513 -3.11054\n",
" 9 5.891e+01 1.810e+01 2.317e+00 -- 3.021e+02 -- -0.0354896 -0.55376 -1.67186 -2.0048 -2.57534 -3.09575 -6.30939 -5.82184 -0.0543126 0.112411 0.195339 0.0877287 0.21341 -0.0199574 1.23313 2.11694\n",
" 11 5.912e+01 2.072e+01 2.113e+00 -- 3.042e+02 -- -0.0151439 -0.533554 -1.65012 -1.98485 -2.55283 -3.06281 -6.00939 -5.52184 -0.0727701 0.113507 0.204816 0.0866854 0.229407 -0.0274598 -2.84324 2.21189\n",
" 13 5.318e+02 2.358e+01 2.081e+00 -- 3.063e+02 -- 0.00246155 -0.516217 -1.63169 -1.96771 -2.5334 -3.03613 -6.30939 -5.22184 -0.0873173 0.114241 0.211893 0.0859833 0.242956 -0.0325361 -0.803414 1.81103\n",
" 15 6.000e+01 2.669e+01 1.754e+00 -- 3.080e+02 -- 0.0177641 -0.501232 -1.61594 -1.95288 -2.51654 -3.01421 -6.60939 -4.92184 -0.0989192 0.114745 0.21724 0.085486 0.254699 -0.0358684 -2.05745 -0.439047\n",
" 17 2.760e+02 3.007e+01 2.029e+00 -- 3.101e+02 -- 0.0311238 -0.488205 -1.60242 -1.93999 -2.50187 -2.99581 -6.90939 -4.99492 -0.108184 0.115134 0.221508 0.0852356 0.26512 -0.0377338 1.31525 1.49865\n",
" 19 3.455e+02 3.375e+01 1.932e+00 -- 3.120e+02 -- 0.0428159 -0.476826 -1.59071 -1.92874 -2.48903 -2.98053 -6.60939 -4.69492 -0.115759 0.115347 0.224648 0.0849875 0.274357 -0.0388865 2.71661 0.824234\n",
" 21 5.627e+02 3.773e+01 2.106e+00 -- 3.141e+02 -- 0.0530961 -0.466839 -1.58055 -1.91888 -2.47782 -2.96759 -6.90939 -4.39492 -0.121894 0.115561 0.227171 0.0847372 0.282893 -0.0391746 3.09606 1.24474\n",
" 23 1.412e+03 4.204e+01 2.098e+00 -- 3.162e+02 -- 0.0621536 -0.458045 -1.57168 -1.9102 -2.468 -2.95676 -7.20939 -4.18138 -0.126953 0.115732 0.229039 0.084359 0.29081 -0.0390341 -1.48027 1.23885\n",
" 25 4.162e+04 4.667e+01 1.855e+00 -- 3.181e+02 -- 0.0701599 -0.450276 -1.56393 -1.90253 -2.45941 -2.94767 -7.50939 -4.06761 -0.131108 0.115932 0.230463 0.0838354 0.298357 -0.0385057 0.151626 1.24726\n",
" 27 1.368e+03 5.163e+01 1.694e+00 -- 3.197e+02 -- 0.0772509 -0.443397 -1.55713 -1.89576 -2.45187 -2.94001 -7.20939 -3.99039 -0.134522 0.11614 0.231552 0.0832518 0.30556 -0.0376963 2.85355 1.25729\n",
" 29 3.278e+03 5.692e+01 1.564e+00 -- 3.213e+02 -- 0.0835427 -0.437294 -1.55116 -1.88977 -2.44524 -2.93357 -7.50939 -3.9332 -0.137326 0.116351 0.232387 0.0826471 0.312468 -0.0366776 2.17754 1.26747\n",
" 31 7.767e+02 6.255e+01 1.446e+00 -- 3.228e+02 -- 0.0891349 -0.43187 -1.54591 -1.88447 -2.4394 -2.92816 -7.20939 -3.88882 -0.139623 0.116565 0.233032 0.0820438 0.319124 -0.0354977 -1.72243 1.27737\n",
" 33 6.146e+04 6.851e+01 1.343e+00 -- 3.241e+02 -- 0.0941131 -0.427041 -1.5413 -1.87977 -2.43426 -2.92364 -7.50939 -3.85335 -0.141498 0.116784 0.233532 0.0814561 0.325567 -0.0341934 0.110556 1.28681\n",
" 35 5.141e+02 7.481e+01 1.246e+00 -- 3.253e+02 -- 0.0985508 -0.422736 -1.53723 -1.87561 -2.42973 -2.91988 -7.20939 -3.82441 -0.143017 0.117008 0.233923 0.0808916 0.331826 -0.0327907 0.968192 1.29571\n",
" 37 2.561e+03 8.143e+01 1.157e+00 -- 3.265e+02 -- 0.102512 -0.418892 -1.53365 -1.87193 -2.42573 -2.9168 -6.90939 -3.80046 -0.144236 0.11724 0.234234 0.0803546 0.337926 -0.0313104 0.474092 1.30403\n",
" 39 4.175e+02 8.835e+01 1.071e+00 -- 3.276e+02 -- 0.106051 -0.415457 -1.53049 -1.86867 -2.4222 -2.91429 -6.60939 -3.78043 -0.145202 0.117482 0.234484 0.0798467 0.343885 -0.0297682 2.50013 1.31174\n",
" 41 2.484e+02 9.557e+01 9.861e-01 -- 3.286e+02 -- 0.109218 -0.412381 -1.52772 -1.86578 -2.41908 -2.91229 -6.30939 -3.76351 -0.14595 0.117736 0.234692 0.079367 0.349713 -0.0281813 -1.35794 1.31885\n",
" 43 1.664e+02 1.030e+02 9.286e-01 -- 3.295e+02 -- 0.112053 -0.409626 -1.52527 -1.86322 -2.41632 -2.91073 -6.60939 -3.74914 -0.146514 0.118004 0.234875 0.0789157 0.355418 -0.0265569 0.955065 1.32532\n",
" 45 3.715e+01 1.108e+02 8.360e-01 -- 3.303e+02 -- 0.114595 -0.407155 -1.52313 -1.86096 -2.41388 -2.90956 -6.30939 -3.73686 -0.146919 0.118288 0.235035 0.0784863 0.361005 -0.024906 -2.00648 1.33124\n",
" 47 1.757e+03 1.187e+02 7.899e-01 -- 3.311e+02 -- 0.116875 -0.404936 -1.52125 -1.85897 -2.41171 -2.90871 -6.60939 -3.7263 -0.147187 0.118587 0.235191 0.0780786 0.366467 -0.0232463 -0.606714 1.3366\n",
" 49 2.419e+03 1.268e+02 7.198e-01 -- 3.318e+02 -- 0.118922 -0.402941 -1.5196 -1.85722 -2.40977 -2.90815 -6.90939 -3.71719 -0.147336 0.118903 0.235343 0.0776863 0.371808 -0.0215787 -0.828333 1.34144\n",
" 50 3.758e+03 1.322e+04 4.674e+00 -- 3.272e+02 -- 0.137314 -0.384994 -1.50519 -1.84182 -2.39252 -2.9049 -8 -3.63812 -0.147803 0.122227 0.236936 0.0738852 0.423919 -0.00499483 -1.20255 1.3852\n",
" 52 2.909e+04 6.385e+02 1.073e+01 -- 3.379e+02 -- 0.136943 -0.384318 -1.50651 -1.84446 -2.39241 -2.90988 -8 -3.64036 -0.1468 0.121814 0.239523 0.0607648 0.447789 0.00582078 -0.683535 1.38625\n",
" 53 8.944e+04 2.504e+02 4.851e-01 -- 3.384e+02 -- 0.137108 -0.382857 -1.511 -1.85395 -2.38877 -2.93977 -8 -3.6534 -0.132666 0.131148 0.26174 0.0372861 0.519783 0.0387377 0.222008 1.38794\n",
" 55 1.907e+04 2.097e+02 2.143e-01 -- 3.386e+02 -- 0.137241 -0.382905 -1.51059 -1.85324 -2.38834 -2.93581 -7.7 -3.65133 -0.133058 0.131502 0.26506 0.0380769 0.517294 0.0298916 0.427458 1.38872\n",
" 57 2.701e+03 1.743e+02 1.707e-01 -- 3.388e+02 -- 0.137362 -0.382947 -1.51023 -1.85264 -2.38795 -2.93261 -7.4 -3.64964 -0.133373 0.131817 0.267849 0.0388293 0.515228 0.0226558 -1.3799 1.38928\n",
" 58 1.901e+04 6.807e+03 9.857e+00 -- 3.289e+02 -- 0.138463 -0.383316 -1.50699 -1.84755 -2.38449 -2.90718 -8 -3.63586 -0.135895 0.13461 0.291172 0.0458143 0.498444 -0.0363648 -0.545357 1.39317\n",
" 60 5.237e+04 4.383e+02 8.985e+00 -- 3.379e+02 -- 0.138201 -0.382988 -1.50809 -1.84815 -2.38423 -2.91079 -8 -3.63639 -0.138737 0.135181 0.288716 0.034039 0.510877 -0.0254366 -0.396913 1.39809\n",
" 61 8.766e+03 4.097e+02 8.172e-01 -- 3.387e+02 -- 0.137968 -0.382536 -1.51047 -1.85072 -2.38194 -2.93875 -8 -3.64468 -0.141543 0.135484 0.281935 0.0183923 0.556885 0.045315 1.79478 1.42693\n",
" 63 1.814e+04 3.493e+02 1.618e-01 -- 3.389e+02 -- 0.13803 -0.382587 -1.51013 -1.85062 -2.38179 -2.93566 -7.7 -3.64395 -0.140901 0.135271 0.284223 0.0205977 0.55377 0.0372067 -0.67593 1.42469\n",
" 64 4.667e+03 9.194e+03 7.138e+00 -- 3.317e+02 -- 0.138594 -0.383049 -1.50705 -1.84979 -2.38045 -2.91127 -8 -3.63796 -0.135166 0.133564 0.302461 0.0405415 0.526549 -0.0285933 1.72877 1.40601\n",
" 66 2.229e+03 8.468e+02 6.334e+00 -- 3.381e+02 -- 0.138375 -0.382821 -1.50788 -1.8493 -2.38027 -2.91526 -7.7 -3.63713 -0.138411 0.134331 0.300768 0.0296213 0.534063 -0.0172368 -2.02693 1.41805\n",
" 67 6.383e+03 4.433e+02 6.946e-01 -- 3.388e+02 -- 0.138019 -0.382446 -1.51012 -1.84832 -2.37854 -2.94579 -8 -3.63552 -0.145232 0.136059 0.290585 0.00900103 0.566205 0.0595408 -1.58605 1.49345\n",
" 69 2.791e+04 3.818e+02 1.832e-01 -- 3.389e+02 -- 0.13807 -0.382495 -1.50981 -1.84878 -2.37843 -2.94234 -8 -3.63605 -0.14425 0.135725 0.292234 0.0124456 0.563364 0.0504802 -0.811369 1.48719\n",
" 71 7.592e+03 3.289e+02 1.419e-01 -- 3.391e+02 -- 0.138118 -0.382539 -1.50952 -1.84915 -2.37833 -2.93967 -8 -3.63644 -0.143394 0.135448 0.293502 0.0154512 0.560893 0.0432337 1.73025 1.48211\n",
" 72 9.792e+00 1.238e+04 4.549e+00 -- 3.345e+02 -- 0.138556 -0.382946 -1.50702 -1.85213 -2.37736 -2.91947 -5 -3.63926 -0.135949 0.133151 0.302987 0.0415893 0.539598 -0.0143255 -2.01979 1.44199\n",
" 73 1.790e+03 2.036e+04 5.362e+00 -- 3.292e+02 -- 0.13686 -0.381484 -1.51285 -1.83573 -2.37645 -2.96346 -8 -3.61694 -0.169822 0.139603 0.306242 -0.0616507 0.578831 0.105014 -2.94898 1.60824\n",
" 75 2.000e+04 3.767e+03 7.143e+00 -- 3.363e+02 -- 0.136937 -0.381554 -1.51207 -1.84035 -2.37648 -2.95217 -8 -3.62182 -0.164447 0.137434 0.31164 -0.0445826 0.577957 0.0688955 -1.04621 1.58574\n",
" 77 1.670e+04 2.864e+03 1.308e+00 -- 3.376e+02 -- 0.137027 -0.381619 -1.51152 -1.84295 -2.37643 -2.94676 -8 -3.62459 -0.16153 0.135849 0.316596 -0.035261 0.576777 0.0522889 -1.10409 1.57318\n",
" 79 3.940e+03 2.269e+03 6.679e-01 -- 3.383e+02 -- 0.137114 -0.381682 -1.51105 -1.8448 -2.37637 -2.94335 -8 -3.62651 -0.159205 0.134652 0.320529 -0.028154 0.575644 0.041706 1.23968 1.56483\n",
" 81 3.136e+03 1.843e+03 4.148e-01 -- 3.387e+02 -- 0.137197 -0.381742 -1.51062 -1.84618 -2.3763 -2.9411 -7.7 -3.62789 -0.157267 0.133753 0.323527 -0.0224013 0.574634 0.0345912 2.94827 1.55921\n",
" 83 7.332e+03 1.524e+03 2.814e-01 -- 3.390e+02 -- 0.137277 -0.381798 -1.51025 -1.84723 -2.37623 -2.9396 -8 -3.6289 -0.155627 0.133083 0.325719 -0.0176294 0.573764 0.0297776 1.97421 1.55551\n",
" 84 4.352e+01 2.419e+05 8.521e+00 -- 3.305e+02 -- 0.13802 -0.382337 -1.50692 -1.85513 -2.37556 -2.92991 -5 -3.63632 -0.141645 0.128163 0.340891 0.0224931 0.566412 -0.00221606 2.81464 1.53266\n",
" 86 2.646e+01 1.360e+05 5.846e+00 -- 3.363e+02 -- 0.137542 -0.382135 -1.50734 -1.85175 -2.3755 -2.93273 -4.7 -3.63458 -0.154265 0.130343 0.343497 0.00406723 0.567925 0.0074277 1.76444 1.5434\n",
" 88 3.599e+00 1.247e+03 4.005e+00 -- 3.403e+02 -- 0.137281 -0.382038 -1.50761 -1.84969 -2.37546 -2.93454 -4.44346 -3.63355 -0.162752 0.1317 0.344853 -0.00669382 0.568721 0.0140179 2.19787 1.55056\n",
" 90 1.187e+00 3.038e+03 3.904e-01 -- 3.407e+02 -- 0.13727 -0.382043 -1.50767 -1.84906 -2.37542 -2.93544 -4.35974 -3.63317 -0.162904 0.131742 0.343848 -0.00910299 0.569009 0.0171214 2.09178 1.55417\n",
" 92 7.462e-01 5.052e+03 2.290e-01 -- 3.409e+02 -- 0.13728 -0.382056 -1.50771 -1.84873 -2.3754 -2.93603 -4.30688 -3.63296 -0.16231 0.131703 0.342609 -0.0100924 0.569167 0.0191537 2.00906 1.55662\n",
" 94 5.215e-01 6.434e+03 1.735e-01 -- 3.411e+02 -- 0.137299 -0.382072 -1.50773 -1.84856 -2.37538 -2.93645 -4.26782 -3.63284 -0.161418 0.131643 0.341287 -0.0103928 0.569258 0.0205828 1.9431 1.5584\n",
" 96 3.952e-01 7.618e+03 1.468e-01 -- 3.413e+02 -- 0.137323 -0.38209 -1.50774 -1.84848 -2.37536 -2.93676 -4.23675 -3.63278 -0.16042 0.131584 0.339946 -0.0103403 0.569314 0.0216563 1.89017 1.55979\n",
" 98 3.181e-01 8.762e+03 1.305e-01 -- 3.414e+02 -- 0.137349 -0.382109 -1.50775 -1.84843 -2.37535 -2.93701 -4.21105 -3.63274 -0.159411 0.131539 0.338623 -0.0101102 0.56935 0.0225121 1.84758 1.56093\n",
" 100 3.494e-01 9.944e+03 1.188e-01 -- 3.415e+02 -- 0.137376 -0.382128 -1.50775 -1.84841 -2.37533 -2.93721 -4.18933 -3.63272 -0.158443 0.131511 0.337339 -0.00979793 0.569375 0.0232282 1.81316 1.56189\n",
" 102 3.644e-01 1.121e+04 1.096e-01 -- 3.416e+02 -- 0.137403 -0.382146 -1.50776 -1.8484 -2.37532 -2.93738 -4.17073 -3.63271 -0.157537 0.131501 0.336106 -0.0094556 0.569394 0.023849 1.78523 1.56274\n",
" 104 3.648e-01 1.258e+04 1.021e-01 -- 3.417e+02 -- 0.13743 -0.382164 -1.50776 -1.84839 -2.37531 -2.93754 -4.15467 -3.6327 -0.156705 0.131508 0.33493 -0.00911106 0.569409 0.0244002 1.7624 1.56349\n",
" 106 3.569e-01 1.408e+04 9.589e-02 -- 3.418e+02 -- 0.137456 -0.382182 -1.50777 -1.84838 -2.3753 -2.93768 -4.14072 -3.6327 -0.155947 0.13153 0.333816 -0.00877869 0.569422 0.0248968 1.74363 1.56416\n",
" 108 3.443e-01 1.574e+04 9.062e-02 -- 3.419e+02 -- 0.137481 -0.382199 -1.50777 -1.84838 -2.3753 -2.9378 -4.12855 -3.63269 -0.155261 0.131565 0.332764 -0.00846538 0.569432 0.0253482 1.72809 1.56477\n",
" 110 3.290e-01 1.757e+04 8.614e-02 -- 3.420e+02 -- 0.137504 -0.382216 -1.50778 -1.84837 -2.37529 -2.93791 -4.11789 -3.63269 -0.154643 0.131611 0.331773 -0.00817392 0.569442 0.0257602 1.71513 1.56532\n",
" 112 3.124e-01 1.961e+04 8.229e-02 -- 3.421e+02 -- 0.137527 -0.382232 -1.50779 -1.84837 -2.37528 -2.93801 -4.10852 -3.63269 -0.154086 0.131666 0.330844 -0.00790497 0.569451 0.0261371 1.70424 1.56582\n",
" 114 2.950e-01 2.186e+04 7.896e-02 -- 3.422e+02 -- 0.137548 -0.382247 -1.50779 -1.84836 -2.37527 -2.9381 -4.10026 -3.63269 -0.153586 0.131728 0.329975 -0.00765804 0.569459 0.0264824 1.69505 1.56628\n",
" 116 2.774e-01 2.436e+04 7.607e-02 -- 3.422e+02 -- 0.137569 -0.382262 -1.5078 -1.84835 -2.37527 -2.93819 -4.09297 -3.63269 -0.153136 0.131795 0.329162 -0.00743212 0.569466 0.0267986 1.68724 1.5667\n",
" 118 2.597e-01 2.714e+04 7.354e-02 -- 3.423e+02 -- 0.137588 -0.382276 -1.50781 -1.84835 -2.37526 -2.93826 -4.08651 -3.63268 -0.152733 0.131866 0.328406 -0.00722598 0.569472 0.0270881 1.68057 1.56707\n",
" 120 2.423e-01 3.022e+04 7.133e-02 -- 3.424e+02 -- 0.137606 -0.38229 -1.50781 -1.84834 -2.37526 -2.93833 -4.08078 -3.63268 -0.15237 0.13194 0.327701 -0.0070383 0.569478 0.0273531 1.67485 1.56742\n",
" 122 2.253e-01 3.365e+04 6.937e-02 -- 3.424e+02 -- 0.137623 -0.382303 -1.50782 -1.84833 -2.37525 -2.93839 -4.07569 -3.63268 -0.152044 0.132015 0.327048 -0.00686775 0.569484 0.0275954 1.66993 1.56773\n",
" 124 2.087e-01 3.745e+04 6.764e-02 -- 3.425e+02 -- 0.137638 -0.382316 -1.50783 -1.84833 -2.37525 -2.93845 -4.07116 -3.63268 -0.151752 0.132092 0.326442 -0.00671305 0.569489 0.027817 1.66568 1.56802\n",
" 126 1.926e-01 4.168e+04 6.611e-02 -- 3.426e+02 -- 0.137653 -0.382328 -1.50783 -1.84832 -2.37525 -2.9385 -4.06713 -3.63268 -0.15149 0.132168 0.325881 -0.00657297 0.569493 0.0280195 1.66199 1.56827\n",
" 128 1.773e-01 4.637e+04 6.474e-02 -- 3.426e+02 -- 0.137667 -0.38234 -1.50784 -1.84831 -2.37524 -2.93855 -4.06353 -3.63268 -0.151255 0.132244 0.325363 -0.00644634 0.569497 0.0282044 1.65879 1.56851\n",
" 130 1.626e-01 5.158e+04 6.353e-02 -- 3.427e+02 -- 0.13768 -0.382351 -1.50784 -1.84831 -2.37524 -2.93859 -4.06031 -3.63268 -0.151045 0.132319 0.324886 -0.00633207 0.569501 0.0283732 1.656 1.56872\n",
" 132 1.487e-01 5.737e+04 6.244e-02 -- 3.428e+02 -- 0.137692 -0.382362 -1.50785 -1.8483 -2.37524 -2.93863 -4.05743 -3.63268 -0.150856 0.132392 0.324446 -0.00622911 0.569505 0.0285272 1.65356 1.56892\n",
" 134 1.356e-01 6.381e+04 6.147e-02 -- 3.428e+02 -- 0.137703 -0.382372 -1.50785 -1.84829 -2.37524 -2.93866 -4.05486 -3.63268 -0.150688 0.132464 0.324041 -0.00613648 0.569508 0.0286677 1.65143 1.56909\n",
" 136 1.233e-01 7.096e+04 6.059e-02 -- 3.429e+02 -- 0.137714 -0.382381 -1.50786 -1.84829 -2.37523 -2.9387 -4.05256 -3.63268 -0.150538 0.132535 0.32367 -0.00605327 0.569511 0.0287957 1.64957 1.56925\n",
" 138 1.118e-01 7.890e+04 5.981e-02 -- 3.429e+02 -- 0.137723 -0.382391 -1.50786 -1.84828 -2.37523 -2.93873 -4.05049 -3.63268 -0.150404 0.132603 0.32333 -0.00597864 0.569514 0.0289123 1.64794 1.5694\n",
" 140 1.012e-01 8.773e+04 5.910e-02 -- 3.430e+02 -- 0.137732 -0.382399 -1.50787 -1.84828 -2.37523 -2.93875 -4.04864 -3.63267 -0.150285 0.132669 0.323018 -0.00591178 0.569516 0.0290186 1.6465 1.56953\n",
" 142 9.132e-02 9.753e+04 5.847e-02 -- 3.431e+02 -- 0.137741 -0.382407 -1.50787 -1.84827 -2.37523 -2.93878 -4.04698 -3.63267 -0.150179 0.132732 0.322732 -0.00585198 0.569518 0.0291153 1.64524 1.56965\n",
" 144 8.225e-02 1.084e+05 5.790e-02 -- 3.431e+02 -- 0.137748 -0.382415 -1.50788 -1.84827 -2.37523 -2.9388 -4.04549 -3.63267 -0.150085 0.132794 0.322472 -0.00579854 0.56952 0.0292034 1.64413 1.56976\n",
" 146 7.394e-02 1.205e+05 5.739e-02 -- 3.432e+02 -- 0.137755 -0.382422 -1.50788 -1.84826 -2.37523 -2.93882 -4.04416 -3.63267 -0.150003 0.132852 0.322234 -0.00575085 0.569522 0.0292835 1.64315 1.56985\n",
" 148 6.635e-02 1.340e+05 5.692e-02 -- 3.432e+02 -- 0.137762 -0.382429 -1.50788 -1.84826 -2.37523 -2.93883 -4.04295 -3.63267 -0.14993 0.132909 0.322017 -0.00570833 0.569524 0.0293563 1.64229 1.56994\n",
" 150 5.944e-02 1.489e+05 5.650e-02 -- 3.433e+02 -- 0.137767 -0.382436 -1.50788 -1.84825 -2.37522 -2.93885 -4.04187 -3.63267 -0.149865 0.132963 0.321819 -0.00567045 0.569526 0.0294226 1.64153 1.57002\n",
" 152 5.317e-02 1.655e+05 5.613e-02 -- 3.434e+02 -- 0.137773 -0.382442 -1.50789 -1.84825 -2.37522 -2.93887 -4.0409 -3.63267 -0.149809 0.133014 0.32164 -0.00563675 0.569527 0.0294828 1.64086 1.5701\n",
" 154 4.750e-02 1.840e+05 5.579e-02 -- 3.434e+02 -- 0.137778 -0.382448 -1.50789 -1.84825 -2.37522 -2.93888 -4.04003 -3.63267 -0.14976 0.133063 0.321476 -0.00560678 0.569529 0.0295374 1.64027 1.57016\n",
" 156 4.239e-02 2.045e+05 5.548e-02 -- 3.435e+02 -- 0.137782 -0.382453 -1.50789 -1.84824 -2.37522 -2.93889 -4.03925 -3.63267 -0.149717 0.13311 0.321328 -0.00558015 0.56953 0.0295871 1.63974 1.57022\n",
" 158 3.778e-02 2.273e+05 5.521e-02 -- 3.435e+02 -- 0.137786 -0.382458 -1.50789 -1.84824 -2.37522 -2.9389 -4.03854 -3.63267 -0.149679 0.133154 0.321193 -0.00555649 0.569531 0.0296321 1.63928 1.57028\n",
" 160 3.366e-02 2.526e+05 5.496e-02 -- 3.436e+02 -- 0.13779 -0.382462 -1.5079 -1.84824 -2.37522 -2.93891 -4.03791 -3.63267 -0.149647 0.133196 0.321071 -0.0055355 0.569532 0.029673 1.63887 1.57033\n",
" 162 2.995e-02 2.807e+05 5.473e-02 -- 3.436e+02 -- 0.137794 -0.382467 -1.5079 -1.84823 -2.37522 -2.93892 -4.03734 -3.63267 -0.149619 0.133236 0.32096 -0.00551686 0.569533 0.0297101 1.63851 1.57037\n",
" 164 2.663e-02 3.119e+05 5.453e-02 -- 3.437e+02 -- 0.137797 -0.382471 -1.5079 -1.84823 -2.37522 -2.93893 -4.03683 -3.63267 -0.149594 0.133273 0.32086 -0.00550034 0.569534 0.0297438 1.63819 1.57041\n",
" 166 2.370e-02 3.467e+05 5.435e-02 -- 3.437e+02 -- 0.137799 -0.382475 -1.5079 -1.84823 -2.37522 -2.93894 -4.03636 -3.63266 -0.149574 0.133308 0.320769 -0.0054857 0.569535 0.0297742 1.6379 1.57045\n",
" 168 2.105e-02 3.852e+05 5.418e-02 -- 3.438e+02 -- 0.137802 -0.382478 -1.5079 -1.84823 -2.37522 -2.93894 -4.03595 -3.63266 -0.149556 0.133341 0.320686 -0.00547269 0.569536 0.0298019 1.63765 1.57048\n",
" 170 1.872e-02 4.281e+05 5.403e-02 -- 3.438e+02 -- 0.137804 -0.382481 -1.5079 -1.84823 -2.37522 -2.93895 -4.03557 -3.63266 -0.14954 0.133373 0.320612 -0.00546117 0.569536 0.0298269 1.63743 1.57051\n",
" 172 1.661e-02 4.757e+05 5.390e-02 -- 3.439e+02 -- 0.137806 -0.382484 -1.5079 -1.84822 -2.37522 -2.93895 -4.03524 -3.63266 -0.149527 0.133402 0.320545 -0.00545095 0.569537 0.0298495 1.63723 1.57054\n",
" 174 1.477e-02 5.286e+05 5.378e-02 -- 3.440e+02 -- 0.137808 -0.382487 -1.5079 -1.84822 -2.37522 -2.93896 -4.03493 -3.63266 -0.149516 0.133429 0.320484 -0.00544189 0.569537 0.02987 1.63705 1.57056\n",
" 176 1.314e-02 5.874e+05 5.367e-02 -- 3.440e+02 -- 0.13781 -0.382489 -1.50791 -1.84822 -2.37522 -2.93896 -4.03466 -3.63266 -0.149506 0.133454 0.320429 -0.00543385 0.569538 0.0298886 1.63689 1.57058\n",
" 178 1.169e-02 6.528e+05 5.357e-02 -- 3.441e+02 -- 0.137812 -0.382492 -1.50791 -1.84822 -2.37522 -2.93897 -4.03442 -3.63266 -0.149498 0.133478 0.320379 -0.00542671 0.569538 0.0299054 1.63675 1.5706\n",
" 180 1.036e-02 7.254e+05 5.348e-02 -- 3.441e+02 -- 0.137813 -0.382494 -1.50791 -1.84822 -2.37522 -2.93897 -4.0342 -3.63266 -0.149491 0.1335 0.320335 -0.00542037 0.569539 0.0299206 1.63663 1.57062\n",
" 182 9.239e-03 8.060e+05 5.340e-02 -- 3.442e+02 -- 0.137814 -0.382496 -1.50791 -1.84822 -2.37522 -2.93897 -4.034 -3.63266 -0.149485 0.133521 0.320294 -0.00541475 0.569539 0.0299343 1.63652 1.57064\n",
" 184 8.217e-03 8.957e+05 5.333e-02 -- 3.442e+02 -- 0.137815 -0.382498 -1.50791 -1.84822 -2.37522 -2.93898 -4.03382 -3.63266 -0.149481 0.13354 0.320258 -0.00540975 0.56954 0.0299467 1.63642 1.57065\n",
" 186 7.314e-03 9.952e+05 5.327e-02 -- 3.443e+02 -- 0.137816 -0.382499 -1.50791 -1.84821 -2.37522 -2.93898 -4.03366 -3.63266 -0.149476 0.133557 0.320225 -0.00540531 0.56954 0.0299579 1.63633 1.57066\n",
" 188 6.496e-03 1.106e+06 5.321e-02 -- 3.443e+02 -- 0.137817 -0.382501 -1.50791 -1.84821 -2.37522 -2.93898 -4.03351 -3.63266 -0.149472 0.133573 0.320195 -0.00540135 0.56954 0.029968 1.63625 1.57068\n",
" 190 5.793e-03 1.229e+06 5.315e-02 -- 3.444e+02 -- 0.137818 -0.382502 -1.50791 -1.84821 -2.37522 -2.93899 -4.03338 -3.63266 -0.14947 0.133589 0.320169 -0.00539784 0.56954 0.0299772 1.63618 1.57069\n",
" 192 5.217e-03 1.365e+06 5.311e-02 -- 3.444e+02 -- 0.137819 -0.382503 -1.50791 -1.84821 -2.37522 -2.93899 -4.03326 -3.63266 -0.149468 0.133603 0.320144 -0.00539472 0.569541 0.0299854 1.63612 1.5707\n",
" 194 4.529e-03 1.517e+06 5.306e-02 -- 3.445e+02 -- 0.13782 -0.382504 -1.50791 -1.84821 -2.37522 -2.93899 -4.03316 -3.63266 -0.149465 0.133615 0.320123 -0.0053919 0.569541 0.0299928 1.63606 1.57071\n",
" 196 4.139e-03 1.686e+06 5.303e-02 -- 3.445e+02 -- 0.13782 -0.382505 -1.50791 -1.84821 -2.37522 -2.93899 -4.03306 -3.63266 -0.149464 0.133627 0.320103 -0.00538946 0.569541 0.0299996 1.63601 1.57071\n",
" 198 3.573e-03 1.873e+06 5.299e-02 -- 3.446e+02 -- 0.137821 -0.382506 -1.50791 -1.84821 -2.37522 -2.93899 -4.03298 -3.63266 -0.149462 0.133636 0.320085 -0.00538723 0.569541 0.0300056 1.63597 1.57072\n",
" 200 3.302e-03 2.081e+06 5.296e-02 -- 3.446e+02 -- 0.137821 -0.382507 -1.50791 -1.84821 -2.37522 -2.93899 -4.0329 -3.63266 -0.14946 0.133646 0.32007 -0.00538531 0.569541 0.0300111 1.63593 1.57073\n",
" 202 2.976e-03 2.313e+06 5.293e-02 -- 3.447e+02 -- 0.137822 -0.382508 -1.50791 -1.84821 -2.37521 -2.93899 -4.03283 -3.63266 -0.149459 0.133656 0.320055 -0.00538353 0.569542 0.0300161 1.63589 1.57073\n",
" 204 2.400e-03 2.570e+06 5.291e-02 -- 3.448e+02 -- 0.137822 -0.382509 -1.50791 -1.84821 -2.37521 -2.939 -4.03277 -3.63266 -0.149457 0.133663 0.320042 -0.00538193 0.569542 0.0300205 1.63586 1.57074\n",
" 206 2.292e-03 2.855e+06 5.288e-02 -- 3.448e+02 -- 0.137823 -0.382509 -1.50791 -1.84821 -2.37521 -2.939 -4.03271 -3.63266 -0.149457 0.13367 0.320031 -0.00538063 0.569542 0.0300245 1.63583 1.57074\n",
" 208 2.150e-03 3.173e+06 5.287e-02 -- 3.449e+02 -- 0.137823 -0.38251 -1.50791 -1.84821 -2.37521 -2.939 -4.03266 -3.63266 -0.149456 0.13368 0.32002 -0.0053794 0.569542 0.0300281 1.6358 1.57075\n",
" 210 1.886e-03 3.525e+06 5.285e-02 -- 3.449e+02 -- 0.137823 -0.38251 -1.50791 -1.84821 -2.37521 -2.939 -4.03262 -3.63266 -0.149456 0.133687 0.320011 -0.00537824 0.569542 0.0300314 1.63578 1.57075\n",
" 212 1.670e-03 3.917e+06 5.283e-02 -- 3.450e+02 -- 0.137823 -0.382511 -1.50791 -1.84821 -2.37521 -2.939 -4.03258 -3.63266 -0.149454 0.133693 0.320002 -0.00537723 0.569542 0.0300343 1.63576 1.57075\n",
" 214 1.507e-03 4.352e+06 5.281e-02 -- 3.450e+02 -- 0.137824 -0.382511 -1.50791 -1.84821 -2.37521 -2.939 -4.03254 -3.63266 -0.149455 0.133696 0.319995 -0.00537633 0.569542 0.030037 1.63574 1.57076\n",
" 216 1.104e-03 4.836e+06 5.280e-02 -- 3.451e+02 -- 0.137824 -0.382512 -1.50791 -1.84821 -2.37521 -2.939 -4.03251 -3.63266 -0.149453 0.133698 0.319988 -0.00537552 0.569542 0.0300394 1.63572 1.57076\n",
" 218 1.248e-03 5.374e+06 5.279e-02 -- 3.451e+02 -- 0.137824 -0.382512 -1.50791 -1.84821 -2.37521 -2.939 -4.03247 -3.63266 -0.14945 0.133705 0.319982 -0.00537493 0.569542 0.0300416 1.63571 1.57076\n",
" 220 9.861e-04 5.971e+06 5.278e-02 -- 3.452e+02 -- 0.137824 -0.382512 -1.50791 -1.84821 -2.37521 -2.939 -4.03245 -3.63266 -0.14945 0.13371 0.319976 -0.00537426 0.569542 0.0300435 1.6357 1.57077\n",
" 222 1.025e-03 6.634e+06 5.277e-02 -- 3.452e+02 -- 0.137824 -0.382513 -1.50791 -1.84821 -2.37521 -2.939 -4.03242 -3.63266 -0.149451 0.133714 0.319971 -0.00537373 0.569542 0.0300453 1.63568 1.57077\n",
" 224 7.583e-04 7.371e+06 5.276e-02 -- 3.453e+02 -- 0.137825 -0.382513 -1.50791 -1.84821 -2.37521 -2.939 -4.0324 -3.63266 -0.14945 0.133716 0.319966 -0.00537318 0.569542 0.0300468 1.63567 1.57077\n",
" 226 6.529e-04 8.190e+06 5.275e-02 -- 3.453e+02 -- 0.137825 -0.382513 -1.50791 -1.84821 -2.37521 -2.939 -4.03238 -3.63266 -0.149445 0.133722 0.319962 -0.00537277 0.569542 0.0300482 1.63566 1.57077\n",
" 228 8.712e-04 9.101e+06 5.274e-02 -- 3.454e+02 -- 0.137825 -0.382513 -1.50791 -1.84821 -2.37521 -2.939 -4.03236 -3.63266 -0.14945 0.133724 0.319959 -0.00537242 0.569543 0.0300495 1.63566 1.57077\n",
" 230 7.297e-04 1.011e+07 5.274e-02 -- 3.454e+02 -- 0.137825 -0.382514 -1.50791 -1.84821 -2.37521 -2.939 -4.03235 -3.63266 -0.149455 0.133729 0.319955 -0.00537195 0.569543 0.0300506 1.63565 1.57077\n",
" 232 7.713e-04 1.124e+07 5.273e-02 -- 3.455e+02 -- 0.137825 -0.382514 -1.50791 -1.84821 -2.37521 -2.939 -4.03233 -3.63266 -0.149461 0.133726 0.319952 -0.00537156 0.569543 0.0300516 1.63564 1.57077\n",
" 234 1.232e-03 1.248e+07 5.273e-02 -- 3.455e+02 -- 0.137825 -0.382514 -1.50791 -1.84821 -2.37521 -2.939 -4.03232 -3.63266 -0.14946 0.133736 0.319949 -0.00537133 0.569543 0.0300525 1.63563 1.57078\n",
" 236 8.332e-04 1.387e+07 5.272e-02 -- 3.456e+02 -- 0.137825 -0.382514 -1.50791 -1.84821 -2.37521 -2.939 -4.03231 -3.63266 -0.149466 0.133721 0.319947 -0.00537067 0.569543 0.0300532 1.63562 1.57078\n",
" 238 1.904e-03 1.541e+07 5.271e-02 -- 3.456e+02 -- 0.137825 -0.382514 -1.50791 -1.84821 -2.37521 -2.939 -4.0323 -3.63266 -0.149478 0.133724 0.319945 -0.00537055 0.569543 0.0300541 1.63562 1.57078\n",
" 240 1.685e-03 1.712e+07 5.271e-02 -- 3.457e+02 -- 0.137825 -0.382514 -1.50791 -1.84821 -2.37521 -2.939 -4.03229 -3.63266 -0.149472 0.13372 0.319943 -0.00536953 0.569543 0.0300545 1.63561 1.57078\n",
" 242 2.904e-03 1.903e+07 5.268e-02 -- 3.458e+02 -- 0.137825 -0.382514 -1.50791 -1.84821 -2.37521 -2.939 -4.03228 -3.63266 -0.149485 0.133743 0.319941 -0.00537003 0.569543 0.0300553 1.63561 1.57078\n",
" 244 2.822e-03 2.113e+07 5.268e-02 -- 3.458e+02 -- 0.137826 -0.382515 -1.50791 -1.84821 -2.37521 -2.939 -4.03227 -3.63266 -0.149462 0.133738 0.319939 -0.00536847 0.569543 0.0300552 1.63559 1.57078\n",
" 246 3.296e-03 2.349e+07 5.271e-02 -- 3.459e+02 -- 0.137826 -0.382514 -1.50791 -1.84821 -2.37521 -2.939 -4.03227 -3.63266 -0.149478 0.133747 0.319938 -0.00536999 0.569543 0.0300566 1.63561 1.57078\n",
" 248 8.922e-04 2.610e+07 5.276e-02 -- 3.459e+02 -- 0.137826 -0.382515 -1.50791 -1.84821 -2.37521 -2.939 -4.03226 -3.63266 -0.149469 0.133756 0.319936 -0.00536822 0.569543 0.0300563 1.63559 1.57078\n",
" 250 1.022e-03 2.900e+07 5.270e-02 -- 3.460e+02 -- 0.137826 -0.382515 -1.50791 -1.84821 -2.37521 -2.939 -4.03226 -3.63266 -0.149478 0.133751 0.319935 -0.0053687 0.569543 0.0300571 1.63559 1.57078\n",
" 252 1.275e-03 3.223e+07 5.272e-02 -- 3.460e+02 -- 0.137826 -0.382515 -1.50791 -1.84821 -2.37521 -2.939 -4.03225 -3.63266 -0.149493 0.133745 0.319934 -0.00536846 0.569543 0.0300572 1.63559 1.57078\n",
" 254 2.477e-03 3.581e+07 5.270e-02 -- 3.461e+02 -- 0.137826 -0.382515 -1.50791 -1.84821 -2.37521 -2.939 -4.03225 -3.63266 -0.149483 0.133762 0.319933 -0.00536868 0.569543 0.0300576 1.63558 1.57078\n",
" 256 2.786e-03 3.976e+07 5.253e-02 -- 3.461e+02 -- 0.137826 -0.382515 -1.50791 -1.84821 -2.37521 -2.939 -4.03224 -3.63266 -0.149446 0.133745 0.319932 -0.00536896 0.569543 0.0300578 1.63558 1.57078\n",
" 258 2.532e-03 4.421e+07 5.284e-02 -- 3.462e+02 -- 0.137827 -0.382514 -1.50791 -1.84821 -2.37521 -2.939 -4.03224 -3.63266 -0.149451 0.133732 0.319933 -0.00537046 0.569543 0.0300589 1.6356 1.57078\n",
" 260 5.032e-03 4.907e+07 5.240e-02 -- 3.462e+02 -- 0.137827 -0.382515 -1.50791 -1.84821 -2.37521 -2.939 -4.03223 -3.63266 -0.149413 0.133704 0.319932 -0.00536945 0.569543 0.0300589 1.63559 1.57078\n",
" 262 3.782e-03 5.456e+07 5.281e-02 -- 3.463e+02 -- 0.137827 -0.382514 -1.50791 -1.84821 -2.37521 -2.939 -4.03223 -3.63266 -0.149414 0.133736 0.319934 -0.00537215 0.569543 0.0300602 1.63561 1.57078\n",
" 264 4.099e-03 6.056e+07 5.240e-02 -- 3.463e+02 -- 0.137827 -0.382515 -1.50791 -1.84821 -2.37521 -2.939 -4.03223 -3.63266 -0.149357 0.133751 0.319931 -0.00537099 0.569543 0.0300595 1.63559 1.57078\n",
" 266 3.335e-03 6.736e+07 5.298e-02 -- 3.464e+02 -- 0.137827 -0.382514 -1.50791 -1.8482 -2.37521 -2.939 -4.03222 -3.63266 -0.14935 0.133787 0.319931 -0.00537319 0.569543 0.030061 1.63561 1.57079\n",
" 268 2.128e-03 7.485e+07 5.274e-02 -- 3.464e+02 -- 0.137827 -0.382515 -1.50791 -1.84821 -2.37521 -2.939 -4.03222 -3.63266 -0.149373 0.133832 0.319929 -0.00537226 0.569543 0.0300605 1.6356 1.57079\n",
" 270 4.362e-03 8.314e+07 5.273e-02 -- 3.465e+02 -- 0.137827 -0.382516 -1.50791 -1.84821 -2.37521 -2.939 -4.03222 -3.63266 -0.149357 0.133842 0.319928 -0.00537111 0.569543 0.0300599 1.63558 1.57078\n",
" 272 7.401e-03 9.241e+07 5.281e-02 -- 3.465e+02 -- 0.137827 -0.382515 -1.50791 -1.84821 -2.37521 -2.939 -4.03222 -3.63266 -0.149297 0.133867 0.319929 -0.00537346 0.569543 0.0300607 1.6356 1.57079\n",
" 273 2.449e-01 9.778e+11 3.957e+01 -- 3.070e+02 -- 0.137829 -0.382514 -1.50792 -1.8482 -2.37521 -2.93901 -4.03219 -3.63266 -0.149811 0.134858 0.319919 -0.00539408 0.569543 0.0300672 1.63566 1.57079\n",
" 275 2.085e-01 1.216e+12 5.713e+01 -- 2.498e+02 -- 0.137838 -0.382547 -1.50791 -1.84822 -2.37521 -2.93901 -4.03229 -3.63267 -0.146142 0.136356 0.319715 -0.00536682 0.569539 0.0300298 1.63519 1.57077\n",
" 278 1.768e-01 1.553e+11 6.527e+01 -- 3.151e+02 -- 0.13783 -0.382542 -1.50791 -1.84822 -2.37521 -2.93901 -4.03228 -3.63266 -0.146083 0.136071 0.319727 -0.00537179 0.56954 0.0300328 1.63531 1.57077\n",
" 280 1.401e-01 2.293e+10 2.877e+01 -- 3.439e+02 -- 0.13781 -0.382487 -1.50792 -1.8482 -2.37521 -2.93902 -4.03218 -3.63266 -0.147503 0.133665 0.319814 -0.00542418 0.569545 0.0300764 1.63632 1.57081\n",
" 282 3.005e-01 1.584e+12 6.015e+01 -- 2.837e+02 -- 0.137818 -0.38251 -1.50792 -1.84821 -2.37521 -2.93902 -4.03219 -3.63266 -0.149569 0.134696 0.319799 -0.00539672 0.56954 0.0300658 1.63602 1.5708\n",
" 285 6.887e-01 3.338e+12 5.965e+01 -- 2.241e+02 -- 0.137808 -0.382511 -1.50792 -1.84821 -2.37521 -2.93901 -4.03218 -3.63266 -0.149948 0.13503 0.319788 -0.0053805 0.569539 0.0300609 1.63597 1.5708\n",
" 286 2.694e+01 1.474e+08 6.253e+01 -- 2.866e+02 -- 0.13812 -0.382806 -1.50771 -1.84841 -2.37523 -2.93893 -4.03632 -3.63262 -0.0891567 0.0420396 0.319568 -0.00184545 0.56953 0.029434 1.62857 1.57017\n",
" 289 1.835e+01 2.710e+07 3.352e+01 -- 3.201e+02 -- 0.137998 -0.382649 -1.50775 -1.84833 -2.37522 -2.93901 -4.03602 -3.63259 -0.0918629 0.0445307 0.319829 -0.00234263 0.56956 0.0297005 1.63399 1.57042\n",
" 291 1.677e+00 4.770e+04 1.799e+01 -- 3.381e+02 -- 0.13705 -0.381771 -1.50801 -1.84766 -2.3752 -2.93967 -4.03323 -3.63229 -0.113651 0.0634674 0.321859 -0.00664201 0.569811 0.0320009 1.68021 1.57254\n",
" 293 8.980e-01 2.894e+05 2.261e+00 -- 3.404e+02 -- 0.137302 -0.381985 -1.50798 -1.84785 -2.37521 -2.93948 -4.03313 -3.63237 -0.110751 0.0693232 0.32106 -0.00552825 0.569729 0.0313773 1.66678 1.57201\n",
" 295 7.368e-01 3.177e+05 9.101e-01 -- 3.413e+02 -- 0.137338 -0.382099 -1.50797 -1.84786 -2.37521 -2.93948 -4.03253 -3.63238 -0.112228 0.0755485 0.320805 -0.00556173 0.569719 0.0313963 1.66607 1.57205\n",
" 297 5.900e-01 3.573e+05 6.175e-01 -- 3.419e+02 -- 0.137345 -0.382178 -1.50797 -1.84784 -2.37521 -2.93949 -4.03194 -3.63237 -0.114034 0.0811151 0.320645 -0.00570323 0.56972 0.0314796 1.66672 1.57213\n",
" 299 4.790e-01 3.999e+05 4.486e-01 -- 3.424e+02 -- 0.137357 -0.382242 -1.50797 -1.84783 -2.37521 -2.9395 -4.03146 -3.63237 -0.11548 0.0859006 0.320493 -0.00578771 0.569719 0.0315348 1.66691 1.57219\n",
" 301 3.950e-01 4.458e+05 3.396e-01 -- 3.427e+02 -- 0.137372 -0.382298 -1.50797 -1.84783 -2.37521 -2.93951 -4.03105 -3.63237 -0.116614 0.0900153 0.320346 -0.00582755 0.569717 0.0315672 1.66675 1.57223\n",
" 303 3.299e-01 4.961e+05 2.661e-01 -- 3.430e+02 -- 0.137388 -0.382345 -1.50797 -1.84783 -2.37521 -2.93951 -4.03071 -3.63237 -0.117524 0.0935712 0.320206 -0.00584073 0.569714 0.0315858 1.6664 1.57225\n",
" 305 2.783e-01 5.518e+05 2.149e-01 -- 3.432e+02 -- 0.137405 -0.382386 -1.50797 -1.84783 -2.37521 -2.93951 -4.03042 -3.63237 -0.118268 0.0966579 0.320075 -0.00583822 0.569712 0.0315961 1.66596 1.57227\n",
" 307 2.366e-01 6.135e+05 1.782e-01 -- 3.434e+02 -- 0.137421 -0.382421 -1.50797 -1.84783 -2.37521 -2.93951 -4.03017 -3.63237 -0.118886 0.0993476 0.319954 -0.00582657 0.569709 0.0316013 1.66548 1.57227\n",
" 309 2.026e-01 6.820e+05 1.512e-01 -- 3.435e+02 -- 0.137436 -0.382452 -1.50797 -1.84784 -2.37521 -2.93951 -4.02995 -3.63238 -0.119407 0.101698 0.319843 -0.00580983 0.569707 0.0316033 1.665 1.57228\n",
" 311 1.745e-01 7.581e+05 1.310e-01 -- 3.436e+02 -- 0.137449 -0.382478 -1.50797 -1.84784 -2.37521 -2.9395 -4.02976 -3.63238 -0.11985 0.103759 0.319742 -0.00579053 0.569705 0.0316034 1.66454 1.57228\n",
" 313 1.511e-01 8.426e+05 1.157e-01 -- 3.438e+02 -- 0.137462 -0.382502 -1.50797 -1.84784 -2.37521 -2.9395 -4.02959 -3.63238 -0.120231 0.10557 0.319649 -0.00577023 0.569703 0.0316024 1.6641 1.57228\n",
" 315 1.314e-01 9.364e+05 1.039e-01 -- 3.439e+02 -- 0.137474 -0.382522 -1.50797 -1.84785 -2.37521 -2.9395 -4.02944 -3.63238 -0.120561 0.107165 0.319566 -0.00574994 0.569701 0.0316007 1.66369 1.57228\n",
" 317 1.146e-01 1.041e+06 9.470e-02 -- 3.440e+02 -- 0.137484 -0.38254 -1.50797 -1.84785 -2.37521 -2.9395 -4.0293 -3.63238 -0.120849 0.108572 0.31949 -0.00573024 0.569699 0.0315987 1.66331 1.57228\n",
" 319 1.002e-01 1.157e+06 8.741e-02 -- 3.440e+02 -- 0.137494 -0.382556 -1.50797 -1.84785 -2.37521 -2.9395 -4.02918 -3.63238 -0.121102 0.109816 0.319421 -0.00571149 0.569698 0.0315965 1.66296 1.57228\n",
" 321 8.813e-02 1.285e+06 8.158e-02 -- 3.441e+02 -- 0.137503 -0.382571 -1.50797 -1.84786 -2.37521 -2.9395 -4.02907 -3.63239 -0.121323 0.110917 0.319359 -0.0056939 0.569696 0.0315944 1.66265 1.57228\n",
" 323 7.747e-02 1.428e+06 7.694e-02 -- 3.442e+02 -- 0.13751 -0.382583 -1.50797 -1.84786 -2.37521 -2.9395 -4.02898 -3.63239 -0.12152 0.111894 0.319303 -0.00567758 0.569695 0.0315923 1.66236 1.57228\n",
" 325 6.841e-02 1.587e+06 7.313e-02 -- 3.443e+02 -- 0.137518 -0.382594 -1.50797 -1.84786 -2.37521 -2.93949 -4.02889 -3.63239 -0.121694 0.112761 0.319252 -0.00566248 0.569694 0.0315904 1.6621 1.57228\n",
" 327 6.037e-02 1.764e+06 7.004e-02 -- 3.443e+02 -- 0.137524 -0.382604 -1.50797 -1.84786 -2.37521 -2.93949 -4.02881 -3.63239 -0.121849 0.113533 0.319207 -0.00564862 0.569693 0.0315885 1.66187 1.57228\n",
" 329 5.354e-02 1.960e+06 6.747e-02 -- 3.444e+02 -- 0.13753 -0.382613 -1.50797 -1.84786 -2.37521 -2.93949 -4.02875 -3.63239 -0.121987 0.114218 0.319166 -0.00563591 0.569692 0.0315868 1.66165 1.57228\n",
" 331 4.752e-02 2.178e+06 6.536e-02 -- 3.445e+02 -- 0.137535 -0.382621 -1.50797 -1.84787 -2.37521 -2.93949 -4.02868 -3.63239 -0.12211 0.114829 0.319129 -0.00562439 0.569691 0.0315853 1.66146 1.57228\n",
" 333 4.222e-02 2.420e+06 6.360e-02 -- 3.445e+02 -- 0.13754 -0.382628 -1.50797 -1.84787 -2.37521 -2.93949 -4.02863 -3.63239 -0.122219 0.115375 0.319095 -0.00561382 0.56969 0.0315839 1.66129 1.57228\n",
" 335 3.769e-02 2.689e+06 6.211e-02 -- 3.446e+02 -- 0.137544 -0.382634 -1.50797 -1.84787 -2.37521 -2.93949 -4.02858 -3.63239 -0.122317 0.115862 0.319065 -0.00560428 0.56969 0.0315826 1.66113 1.57228\n",
" 337 3.311e-02 2.989e+06 6.088e-02 -- 3.447e+02 -- 0.137547 -0.38264 -1.50797 -1.84787 -2.37521 -2.93949 -4.02853 -3.63239 -0.122404 0.116299 0.319038 -0.00559562 0.569689 0.0315815 1.66099 1.57228\n",
" 339 2.982e-02 3.321e+06 5.975e-02 -- 3.447e+02 -- 0.137551 -0.382645 -1.50797 -1.84787 -2.37521 -2.93949 -4.02849 -3.63239 -0.122484 0.116684 0.319014 -0.00558779 0.569689 0.0315804 1.66087 1.57228\n",
" 341 2.639e-02 3.690e+06 5.889e-02 -- 3.448e+02 -- 0.137554 -0.382649 -1.50797 -1.84787 -2.37521 -2.93949 -4.02845 -3.63239 -0.122555 0.117032 0.318992 -0.00558071 0.569688 0.0315795 1.66076 1.57228\n",
" 343 2.409e-02 4.100e+06 5.809e-02 -- 3.448e+02 -- 0.137557 -0.382653 -1.50797 -1.84787 -2.37521 -2.93949 -4.02842 -3.63239 -0.122619 0.117341 0.318972 -0.00557437 0.569688 0.0315786 1.66065 1.57228\n",
" 345 2.129e-02 4.556e+06 5.751e-02 -- 3.449e+02 -- 0.137559 -0.382657 -1.50797 -1.84787 -2.37521 -2.93949 -4.02839 -3.63239 -0.122676 0.117624 0.318955 -0.00556861 0.569687 0.0315779 1.66056 1.57227\n",
" 347 1.889e-02 5.062e+06 5.690e-02 -- 3.450e+02 -- 0.137561 -0.38266 -1.50797 -1.84787 -2.37521 -2.93949 -4.02837 -3.6324 -0.12273 0.117874 0.318939 -0.0055633 0.569687 0.0315772 1.66048 1.57227\n",
" 349 1.666e-02 5.625e+06 5.637e-02 -- 3.450e+02 -- 0.137563 -0.382663 -1.50797 -1.84788 -2.37521 -2.93949 -4.02834 -3.6324 -0.122775 0.118097 0.318924 -0.00555856 0.569687 0.0315765 1.6604 1.57227\n",
" 351 1.548e-02 6.250e+06 5.589e-02 -- 3.451e+02 -- 0.137565 -0.382666 -1.50797 -1.84788 -2.37521 -2.93949 -4.02832 -3.6324 -0.122815 0.118293 0.318911 -0.00555438 0.569686 0.031576 1.66034 1.57227\n",
" 353 1.367e-02 6.945e+06 5.559e-02 -- 3.451e+02 -- 0.137567 -0.382668 -1.50797 -1.84788 -2.37521 -2.93949 -4.0283 -3.6324 -0.122851 0.118477 0.3189 -0.00555074 0.569686 0.0315756 1.66028 1.57227\n",
" 355 1.218e-02 7.717e+06 5.527e-02 -- 3.452e+02 -- 0.137568 -0.38267 -1.50797 -1.84788 -2.37521 -2.93949 -4.02828 -3.6324 -0.122888 0.118639 0.318889 -0.00554714 0.569686 0.0315751 1.66022 1.57227\n",
" 357 1.125e-02 8.574e+06 5.496e-02 -- 3.452e+02 -- 0.13757 -0.382672 -1.50797 -1.84788 -2.37521 -2.93949 -4.02827 -3.6324 -0.122919 0.118783 0.31888 -0.00554391 0.569686 0.0315747 1.66017 1.57227\n",
" 359 9.239e-03 9.527e+06 5.475e-02 -- 3.453e+02 -- 0.137571 -0.382674 -1.50797 -1.84788 -2.37521 -2.93949 -4.02825 -3.6324 -0.122951 0.118917 0.318871 -0.00554137 0.569686 0.0315743 1.66013 1.57227\n",
" 361 8.671e-03 1.059e+07 5.442e-02 -- 3.453e+02 -- 0.137572 -0.382675 -1.50797 -1.84788 -2.37521 -2.93949 -4.02824 -3.6324 -0.122974 0.119027 0.318864 -0.00553843 0.569685 0.0315738 1.66009 1.57227\n",
" 363 7.118e-03 1.176e+07 5.421e-02 -- 3.454e+02 -- 0.137573 -0.382677 -1.50797 -1.84788 -2.37521 -2.93949 -4.02823 -3.6324 -0.122989 0.11913 0.318857 -0.00553662 0.569685 0.0315737 1.66005 1.57227\n",
" 365 6.502e-03 1.307e+07 5.400e-02 -- 3.454e+02 -- 0.137574 -0.382678 -1.50797 -1.84788 -2.37521 -2.93948 -4.02822 -3.6324 -0.123017 0.119215 0.318851 -0.0055347 0.569685 0.0315735 1.66002 1.57227\n",
" 367 5.241e-03 1.452e+07 5.386e-02 -- 3.455e+02 -- 0.137574 -0.382679 -1.50797 -1.84788 -2.37521 -2.93948 -4.02821 -3.6324 -0.123025 0.119292 0.318845 -0.00553265 0.569685 0.0315732 1.65999 1.57227\n",
" 369 4.797e-03 1.614e+07 5.362e-02 -- 3.456e+02 -- 0.137575 -0.382679 -1.50797 -1.84788 -2.37521 -2.93948 -4.0282 -3.6324 -0.123045 0.119355 0.318841 -0.00553154 0.569685 0.0315732 1.65997 1.57227\n",
" 371 4.242e-03 1.793e+07 5.357e-02 -- 3.456e+02 -- 0.137575 -0.38268 -1.50797 -1.84788 -2.37521 -2.93948 -4.0282 -3.6324 -0.123059 0.119412 0.318836 -0.00553001 0.569685 0.031573 1.65995 1.57227\n",
" 373 5.841e-03 1.992e+07 5.347e-02 -- 3.457e+02 -- 0.137576 -0.382681 -1.50797 -1.84788 -2.37521 -2.93948 -4.02819 -3.6324 -0.123078 0.119463 0.318832 -0.00552879 0.569685 0.0315729 1.65993 1.57227\n",
" 375 3.219e-03 2.213e+07 5.366e-02 -- 3.457e+02 -- 0.137577 -0.382682 -1.50797 -1.84788 -2.37521 -2.93948 -4.02818 -3.6324 -0.123089 0.119532 0.318828 -0.00552759 0.569685 0.0315728 1.65991 1.57227\n",
" 377 2.633e-03 2.459e+07 5.334e-02 -- 3.458e+02 -- 0.137577 -0.382683 -1.50797 -1.84788 -2.37521 -2.93948 -4.02818 -3.6324 -0.123107 0.119571 0.318825 -0.00552604 0.569685 0.0315725 1.65989 1.57227\n",
" 379 4.116e-03 2.733e+07 5.319e-02 -- 3.458e+02 -- 0.137578 -0.382683 -1.50797 -1.84788 -2.37521 -2.93948 -4.02817 -3.6324 -0.12312 0.119602 0.318822 -0.00552505 0.569685 0.0315725 1.65988 1.57227\n",
" 380 5.533e+01 8.512e+13 3.580e+01 -- 3.100e+02 -- 0.13758 -0.382687 -1.50797 -1.84788 -2.37521 -2.93948 -4.02813 -3.6324 -0.123171 0.120095 0.318799 -0.0055189 0.569684 0.0315729 1.65977 1.57227\n",
" 382 7.194e+00 3.296e+06 1.119e+02 -- 1.981e+02 -- 0.136564 -0.384498 -1.50899 -1.84727 -2.37514 -2.93869 -4.00752 -3.63339 -0.76763 0.784575 0.315418 -0.0211879 0.568207 0.0342865 1.65043 1.57386\n",
" 385 6.393e+00 7.364e+05 4.355e+01 -- 2.417e+02 -- 0.134399 -0.385404 -1.5086 -1.84755 -2.37516 -2.93833 -4.01279 -3.63358 -0.756411 0.778329 0.313579 -0.0196636 0.5681 0.0330514 1.62964 1.57277\n",
" 387 6.586e+00 2.759e+05 6.486e+01 -- 3.066e+02 -- 0.116998 -0.394901 -1.50618 -1.84965 -2.37527 -2.93545 -4.04221 -3.6351 -0.661407 0.711541 0.301895 -0.00709193 0.567185 0.023737 1.46306 1.56458\n",
" 389 8.793e-01 6.706e+03 1.268e+01 -- 3.192e+02 -- 0.115989 -0.400915 -1.50697 -1.8488 -2.37519 -2.93697 -4.03138 -3.63433 -0.618922 0.64982 0.305913 -0.0117627 0.567838 0.0280683 1.55938 1.56838\n",
" 391 8.721e-01 1.576e+04 4.935e+00 -- 3.242e+02 -- 0.112264 -0.405631 -1.5073 -1.84862 -2.37517 -2.93742 -4.02343 -3.63415 -0.565235 0.59268 0.305403 -0.0120062 0.568 0.0294253 1.57977 1.56959\n",
" 393 8.609e-01 2.268e+04 2.879e+00 -- 3.271e+02 -- 0.110008 -0.408683 -1.50757 -1.84848 -2.37515 -2.93776 -4.01681 -3.63402 -0.51694 0.540991 0.304308 -0.0120189 0.568108 0.0304757 1.59266 1.57052\n",
" 395 8.448e-01 2.811e+04 1.940e+00 -- 3.290e+02 -- 0.108886 -0.410427 -1.50779 -1.8484 -2.37513 -2.938 -4.01122 -3.63393 -0.473875 0.494416 0.302852 -0.011813 0.568178 0.0312724 1.60029 1.57123\n",
" 397 8.241e-01 3.293e+04 1.442e+00 -- 3.304e+02 -- 0.108628 -0.411163 -1.50798 -1.84834 -2.37512 -2.93818 -4.00644 -3.63387 -0.435812 0.452647 0.30122 -0.0114726 0.568224 0.0318892 1.60467 1.57177\n",
" 399 7.993e-01 3.764e+04 1.149e+00 -- 3.316e+02 -- 0.108995 -0.411144 -1.50813 -1.84831 -2.37512 -2.93832 -4.00231 -3.63383 -0.402379 0.415346 0.299533 -0.0110612 0.568256 0.0323797 1.60708 1.5722\n",
" 401 7.711e-01 4.253e+04 9.635e-01 -- 3.325e+02 -- 0.109791 -0.410578 -1.50827 -1.84829 -2.37511 -2.93843 -3.99871 -3.6338 -0.37314 0.382148 0.297866 -0.01062 0.568277 0.0327796 1.60828 1.57255\n",
" 403 7.403e-01 4.778e+04 8.370e-01 -- 3.334e+02 -- 0.110862 -0.409634 -1.50838 -1.84828 -2.3751 -2.93852 -3.99556 -3.63378 -0.347639 0.35268 0.296263 -0.0101744 0.568292 0.0331127 1.60875 1.57284\n",
" 405 7.072e-01 5.351e+04 7.456e-01 -- 3.341e+02 -- 0.11209 -0.408441 -1.50848 -1.84827 -2.3751 -2.93859 -3.9928 -3.63377 -0.325436 0.326572 0.29475 -0.00973915 0.568303 0.0333948 1.60877 1.57308\n",
" 407 6.726e-01 5.984e+04 6.759e-01 -- 3.348e+02 -- 0.11339 -0.407097 -1.50857 -1.84826 -2.3751 -2.93866 -3.99036 -3.63375 -0.306122 0.303476 0.29334 -0.00932269 0.568311 0.0336367 1.60852 1.57328\n",
" 409 6.368e-01 6.686e+04 6.201e-01 -- 3.354e+02 -- 0.114702 -0.405678 -1.50865 -1.84826 -2.37509 -2.93871 -3.98821 -3.63374 -0.289327 0.283064 0.292038 -0.00892928 0.568317 0.033846 1.60811 1.57346\n",
" 411 6.003e-01 7.469e+04 5.735e-01 -- 3.360e+02 -- 0.115982 -0.404238 -1.50873 -1.84826 -2.37509 -2.93876 -3.98632 -3.63374 -0.274719 0.265039 0.290845 -0.00856069 0.568322 0.0340282 1.60761 1.57361\n",
" 413 5.635e-01 8.343e+04 5.332e-01 -- 3.365e+02 -- 0.117205 -0.402814 -1.50879 -1.84826 -2.37509 -2.9388 -3.98465 -3.63373 -0.262009 0.249128 0.289757 -0.0082171 0.568325 0.0341871 1.60706 1.57375\n",
" 415 5.268e-01 9.320e+04 4.976e-01 -- 3.370e+02 -- 0.118353 -0.401435 -1.50885 -1.84826 -2.37509 -2.93884 -3.98317 -3.63373 -0.25094 0.235089 0.288769 -0.00789773 0.568328 0.034326 1.60648 1.57386\n",
" 417 4.906e-01 1.041e+05 4.652e-01 -- 3.375e+02 -- 0.119416 -0.400117 -1.5089 -1.84826 -2.37508 -2.93887 -3.98187 -3.63372 -0.241292 0.222703 0.287874 -0.00760125 0.56833 0.0344471 1.60589 1.57396\n",
" 419 4.551e-01 1.164e+05 4.356e-01 -- 3.379e+02 -- 0.120392 -0.398872 -1.50895 -1.84826 -2.37508 -2.93889 -3.98073 -3.63372 -0.23287 0.211778 0.287065 -0.00732599 0.568331 0.0345524 1.6053 1.57405\n",
" 421 4.207e-01 1.301e+05 4.080e-01 -- 3.383e+02 -- 0.12128 -0.397707 -1.50899 -1.84826 -2.37508 -2.93892 -3.97972 -3.63372 -0.225507 0.202139 0.286336 -0.00707017 0.568332 0.0346434 1.6047 1.57413\n",
" 423 3.877e-01 1.454e+05 3.822e-01 -- 3.387e+02 -- 0.122083 -0.396622 -1.50902 -1.84826 -2.37508 -2.93894 -3.97883 -3.63372 -0.219056 0.193634 0.285679 -0.00683198 0.568332 0.0347216 1.6041 1.57419\n",
" 425 3.562e-01 1.626e+05 3.580e-01 -- 3.391e+02 -- 0.122806 -0.39562 -1.50906 -1.84827 -2.37508 -2.93895 -3.97805 -3.63372 -0.213392 0.186128 0.285086 -0.0066097 0.568332 0.034788 1.60349 1.57425\n",
" 427 3.264e-01 1.818e+05 3.352e-01 -- 3.394e+02 -- 0.123454 -0.394697 -1.50909 -1.84827 -2.37508 -2.93897 -3.97737 -3.63372 -0.208405 0.179498 0.284553 -0.00640172 0.568331 0.0348437 1.60289 1.57429\n",
" 429 2.984e-01 2.032e+05 3.137e-01 -- 3.397e+02 -- 0.124033 -0.39385 -1.50911 -1.84827 -2.37508 -2.93898 -3.97677 -3.63372 -0.204002 0.173639 0.284072 -0.0062066 0.56833 0.0348898 1.60227 1.57433\n",
" 431 2.872e-01 2.272e+05 2.934e-01 -- 3.400e+02 -- 0.124549 -0.393077 -1.50913 -1.84828 -2.37508 -2.93899 -3.97624 -3.63372 -0.200103 0.168457 0.283638 -0.00602309 0.568329 0.034927 1.60166 1.57436\n",
" 433 2.792e-01 2.540e+05 2.744e-01 -- 3.403e+02 -- 0.125009 -0.392371 -1.50915 -1.84828 -2.37508 -2.939 -3.97578 -3.63372 -0.196637 0.163868 0.283246 -0.00585012 0.568328 0.0349564 1.60104 1.57439\n",
" 435 2.715e-01 2.839e+05 2.565e-01 -- 3.406e+02 -- 0.125417 -0.39173 -1.50917 -1.84828 -2.37508 -2.939 -3.97538 -3.63372 -0.193545 0.159801 0.282892 -0.00568681 0.568326 0.0349787 1.60041 1.57441\n",
" 437 2.639e-01 3.173e+05 2.396e-01 -- 3.408e+02 -- 0.12578 -0.391148 -1.50919 -1.84829 -2.37508 -2.939 -3.97503 -3.63373 -0.190778 0.156191 0.282571 -0.00553244 0.568324 0.0349948 1.59979 1.57442\n",
" 439 2.564e-01 3.545e+05 2.239e-01 -- 3.410e+02 -- 0.126102 -0.390621 -1.5092 -1.84829 -2.37508 -2.93901 -3.97472 -3.63373 -0.188292 0.152983 0.282281 -0.00538642 0.568322 0.0350054 1.59916 1.57443\n",
" 441 2.488e-01 3.960e+05 2.092e-01 -- 3.412e+02 -- 0.126388 -0.390144 -1.50921 -1.8483 -2.37508 -2.93901 -3.97446 -3.63373 -0.186054 0.150128 0.282017 -0.00524832 0.56832 0.0350114 1.59854 1.57443\n",
" 443 2.409e-01 4.422e+05 1.954e-01 -- 3.414e+02 -- 0.126643 -0.389713 -1.50922 -1.84831 -2.37508 -2.93901 -3.97422 -3.63373 -0.18403 0.147582 0.281778 -0.00511776 0.568317 0.0350133 1.59793 1.57443\n",
" 445 2.327e-01 4.937e+05 1.827e-01 -- 3.416e+02 -- 0.126869 -0.389324 -1.50923 -1.84831 -2.37508 -2.93901 -3.97402 -3.63374 -0.182195 0.14531 0.281562 -0.00499447 0.568315 0.0350119 1.59733 1.57443\n",
" 447 2.243e-01 5.511e+05 1.709e-01 -- 3.418e+02 -- 0.12707 -0.388974 -1.50923 -1.84832 -2.37508 -2.93901 -3.97384 -3.63374 -0.18053 0.143277 0.281365 -0.00487824 0.568313 0.0350078 1.59674 1.57443\n",
" 449 2.155e-01 6.150e+05 1.600e-01 -- 3.419e+02 -- 0.127249 -0.388658 -1.50924 -1.84833 -2.37508 -2.939 -3.97369 -3.63374 -0.179014 0.141458 0.281186 -0.00476884 0.56831 0.0350015 1.59617 1.57442\n",
" 451 2.064e-01 6.861e+05 1.499e-01 -- 3.421e+02 -- 0.127409 -0.388373 -1.50924 -1.84833 -2.37508 -2.939 -3.97356 -3.63375 -0.177636 0.139827 0.281023 -0.00466608 0.568308 0.0349935 1.59562 1.57441\n",
" 453 1.969e-01 7.653e+05 1.407e-01 -- 3.422e+02 -- 0.127551 -0.388118 -1.50925 -1.84834 -2.37508 -2.939 -3.97344 -3.63375 -0.176375 0.138362 0.280875 -0.00456977 0.568305 0.0349843 1.59509 1.5744\n",
" 455 1.874e-01 8.534e+05 1.322e-01 -- 3.424e+02 -- 0.127678 -0.387888 -1.50925 -1.84835 -2.37508 -2.939 -3.97334 -3.63375 -0.175226 0.137048 0.280741 -0.00447979 0.568303 0.0349741 1.59459 1.57439\n",
" 457 1.775e-01 9.514e+05 1.245e-01 -- 3.425e+02 -- 0.127792 -0.387681 -1.50925 -1.84835 -2.37508 -2.93899 -3.97325 -3.63376 -0.174176 0.135866 0.28062 -0.00439585 0.568301 0.0349634 1.59411 1.57438\n",
" 459 1.678e-01 1.060e+06 1.174e-01 -- 3.426e+02 -- 0.127894 -0.387495 -1.50926 -1.84836 -2.37508 -2.93899 -3.97317 -3.63376 -0.17322 0.134805 0.280509 -0.0043178 0.568299 0.0349524 1.59365 1.57437\n",
" 461 1.581e-01 1.182e+06 1.110e-01 -- 3.427e+02 -- 0.127985 -0.387328 -1.50926 -1.84836 -2.37508 -2.93899 -3.9731 -3.63376 -0.172347 0.13385 0.280409 -0.00424533 0.568297 0.0349413 1.59322 1.57436\n",
" 463 1.484e-01 1.317e+06 1.052e-01 -- 3.428e+02 -- 0.128066 -0.387178 -1.50926 -1.84837 -2.37508 -2.93898 -3.97304 -3.63376 -0.171549 0.132987 0.280318 -0.0041782 0.568295 0.0349303 1.59282 1.57435\n",
" 465 1.387e-01 1.467e+06 9.997e-02 -- 3.429e+02 -- 0.128139 -0.387043 -1.50926 -1.84837 -2.37508 -2.93898 -3.97299 -3.63377 -0.170819 0.132211 0.280236 -0.00411621 0.568293 0.0349196 1.59244 1.57434\n",
" 467 1.295e-01 1.634e+06 9.520e-02 -- 3.430e+02 -- 0.128204 -0.386922 -1.50926 -1.84838 -2.37508 -2.93898 -3.97295 -3.63377 -0.170155 0.131513 0.280162 -0.0040591 0.568292 0.0349092 1.59209 1.57433\n",
" 469 1.206e-01 1.820e+06 9.088e-02 -- 3.431e+02 -- 0.128263 -0.386813 -1.50926 -1.84838 -2.37508 -2.93898 -3.97291 -3.63377 -0.169555 0.130886 0.280094 -0.00400655 0.56829 0.0348992 1.59177 1.57432\n",
" 471 1.118e-01 2.026e+06 8.702e-02 -- 3.432e+02 -- 0.128316 -0.386716 -1.50926 -1.84839 -2.37508 -2.93897 -3.97287 -3.63377 -0.169004 0.130321 0.280033 -0.00395823 0.568289 0.0348897 1.59147 1.57431\n",
" 473 1.033e-01 2.256e+06 8.358e-02 -- 3.433e+02 -- 0.128363 -0.386628 -1.50926 -1.84839 -2.37508 -2.93897 -3.97284 -3.63377 -0.168497 0.129809 0.279978 -0.00391395 0.568287 0.0348807 1.59119 1.57431\n",
" 475 9.579e-02 2.511e+06 8.038e-02 -- 3.433e+02 -- 0.128405 -0.386549 -1.50926 -1.84839 -2.37508 -2.93897 -3.97282 -3.63378 -0.168047 0.129353 0.279928 -0.00387351 0.568286 0.0348723 1.59093 1.5743\n",
" 477 8.818e-02 2.794e+06 7.758e-02 -- 3.434e+02 -- 0.128443 -0.386478 -1.50926 -1.8484 -2.37508 -2.93897 -3.97279 -3.63378 -0.167639 0.128938 0.279883 -0.00383641 0.568285 0.0348643 1.59069 1.57429\n",
" 479 8.091e-02 3.109e+06 7.509e-02 -- 3.435e+02 -- 0.128478 -0.386415 -1.50926 -1.8484 -2.37508 -2.93896 -3.97277 -3.63378 -0.167267 0.128559 0.279842 -0.00380258 0.568284 0.0348569 1.59048 1.57428\n",
" 481 7.424e-02 3.460e+06 7.280e-02 -- 3.436e+02 -- 0.128509 -0.386357 -1.50926 -1.8484 -2.37508 -2.93896 -3.97275 -3.63378 -0.166936 0.128218 0.279805 -0.00377181 0.568283 0.0348501 1.59028 1.57428\n",
" 483 6.806e-02 3.849e+06 7.073e-02 -- 3.436e+02 -- 0.128537 -0.386306 -1.50926 -1.84841 -2.37508 -2.93896 -3.97274 -3.63378 -0.166628 0.127921 0.279772 -0.00374381 0.568282 0.0348438 1.5901 1.57427\n",
" 485 6.227e-02 4.282e+06 6.891e-02 -- 3.437e+02 -- 0.128562 -0.38626 -1.50926 -1.84841 -2.37508 -2.93896 -3.97272 -3.63378 -0.166346 0.127652 0.279742 -0.00371833 0.568281 0.034838 1.58993 1.57426\n",
" 487 5.690e-02 4.762e+06 6.723e-02 -- 3.438e+02 -- 0.128585 -0.386219 -1.50926 -1.84841 -2.37508 -2.93896 -3.97271 -3.63378 -0.166098 0.127408 0.279715 -0.00369517 0.56828 0.0348326 1.58978 1.57426\n",
" 489 5.246e-02 5.297e+06 6.572e-02 -- 3.438e+02 -- 0.128605 -0.386181 -1.50926 -1.84841 -2.37508 -2.93896 -3.9727 -3.63378 -0.165881 0.127186 0.27969 -0.00367415 0.56828 0.0348277 1.58964 1.57426\n",
" 491 4.668e-02 5.891e+06 6.451e-02 -- 3.439e+02 -- 0.128623 -0.386148 -1.50926 -1.84841 -2.37508 -2.93895 -3.97269 -3.63378 -0.165662 0.126983 0.279668 -0.00365488 0.568279 0.0348231 1.58952 1.57425\n",
" 493 4.336e-02 6.551e+06 6.330e-02 -- 3.440e+02 -- 0.12864 -0.386118 -1.50927 -1.84842 -2.37508 -2.93895 -3.97268 -3.63379 -0.165486 0.126802 0.279649 -0.00363782 0.568279 0.0348191 1.58941 1.57425\n",
" 495 3.896e-02 7.284e+06 6.215e-02 -- 3.440e+02 -- 0.128655 -0.386091 -1.50927 -1.84842 -2.37508 -2.93895 -3.97267 -3.63379 -0.165321 0.12665 0.279631 -0.00362204 0.568278 0.0348153 1.5893 1.57424\n",
" 497 3.577e-02 8.098e+06 6.116e-02 -- 3.441e+02 -- 0.128668 -0.386067 -1.50927 -1.84842 -2.37508 -2.93895 -3.97266 -3.63379 -0.165157 0.126529 0.279614 -0.00360793 0.568278 0.0348118 1.58921 1.57424\n",
" 499 3.106e-02 9.005e+06 6.042e-02 -- 3.442e+02 -- 0.12868 -0.386045 -1.50927 -1.84842 -2.37508 -2.93895 -3.97266 -3.63379 -0.164994 0.126396 0.2796 -0.00359503 0.568277 0.0348086 1.58912 1.57424\n",
"********************\n",
"0.128691 -0.386025 -1.50927 -1.84842 -2.37508 -2.93895 -3.97265 -3.63379 -0.164908 0.12629 0.279587 -0.00358386 0.568277 0.034806 1.58905 1.57423\n",
"0.00126105 0.000998709 0.00018594 0.000356512 2.30678e-05 0.000253636 0.00869718 0.000210956 0.0200512 0.0186844 0.00479074 0.00354031 0.000510855 0.00183101 0.0188838 0.0015161\n",
"-10903 199.23 -2.70075 -1.73798 -1.01372 0.00748608 0.0139392 -0.0830399 -358.669 233.027 -1.33222 -0.312077 -0.267079 -0.122273 0.0137976 -0.0193582\n",
"********************\n"
]
}
],
"source": [
"Cx = clag.clag('cxd10r', [[t1,t2]], [[l1,l2]], [[l1e,l2e]], dt, fqL, p1, p2)\n",
"p = np.concatenate( ((p1+p2)*0.5-0.3,p1*0+0.1) ) # a good starting point generally\n",
"p, pe = clag.optimize(Cx, p)"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"phi, phie = p[nfq:], pe[nfq:]\n",
"lag, lage = phi/(2*np.pi*fqd), phie/(2*np.pi*fqd) \n",
"cx, cxe = p[:nfq], pe[:nfq]"
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<Container object of 3 artists>"
]
},
"execution_count": 12,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAhIAAAFrCAYAAACE+GArAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAFI1JREFUeJzt3WtsZOddB+Bf2rgttMD2QjxJKXJqqnghC8WuA8lWxRE0\nQghSJBDYUhGsuURctRICpKKYsEFFQkCXD1wUpKRIFd4FBGoRhMuHpBK7LRg7XAxxgawNIbvjtKVb\neoHWasyH46Ver73reT0zxzPzPNJox+e858zfPu/O+c25vJMAAAAAAAAAAAAAAAAAAAAAAAAAAAAA\nQN94S5I/TvJckheSvG2XNg9tzf90kieSfGW3igMADu5FHVz3FyZ5KsmPbv28uWP+zyQ5uTV/Mkkz\nyV8meUUHawIAetALSe7f9vNNSS4l+alt016S5GNJfqiLdQEAB9DJIxLXc3uS4SR/sW3aZ5O8P8k9\ntVQEALSsriDR2Pp3fcf057fNAwAOuZvrLmAXO6+luOLWrQcA0JpLW4+2qytINLf+Hd72fLefr7j1\ntttuu3jx4sWOFwYAfei5VDc2tD1M1BUkVlMFhvuS/P3WtJck+YZcfQHmFbdevHgx73nPe3L06NEu\nldg+J0+ezOnTp3vytQ6yvlaX3W/7/bS7UZvrze/m9mo3fa297fW1velr7W3fyb729NNP5+1vf/tr\nUx3V76kg8fIkb9j28+uTvDHJR5M8m+R0knck+dck/7b1/JNJfnevFR49ejTj4+Odqrdjjhw50rW6\n2/1aB1lfq8vut/1+2t2ozfXmd3N7tZu+1t72+tre9LX2tu90X+ukF3dw3ceTnE/yQKrrHr556/kr\nk7w3ybkkL0vyc0l+IsnHk8wk2e38xa1JHnjggQdy6629eZnEsWPHeva1DrK+Vpfdb/v9tLtRm73m\nz8/PZ2ZmZl91HEb6Wnvb62t709fa275Tfe3SpUt55JFHkuSRdOCIxE3tXmGHjCdZXFxc7Nn0Tu+4\n//778773va/uMhgA+hrdsLS0lImJiSSZSLLU7vXXdfsnANAHBAnYoZcPNdNb9DX6gSABO3hzp1v0\nNfqBIAEAFBMkAIBiggQAUEyQAACKCRIAQDFBAgAoJkgAAMUECQCgmCABABQTJACAYoIEAFBMkAAA\nigkSAEAxQQIAKCZIAADFBAkAoJggAQAUEyQAgGKCBABQTJAAAIoJEgBAMUECACgmSAAAxQQJAKCY\nIAEAFBMkAIBiggQAUEyQAACKCRIAQDFBAgAoJkgAAMUECQCgmCABABQTJACAYoIEAFBMkAAAigkS\nAEAxQQIAKCZIAADFBAkAoJggAQAUEyQAgGKCBABQTJAAAIoJEgBAMUECACgmSAAAxQQJAKCYIAEA\nFBMkAIBiggQAUEyQAACKCRIAQDFBAgAoJkgAAMUECQCgmCABABQTJACAYoIEAFBMkAAAigkSAEAx\nQQIAKCZIAADF6gwSDyV5YcfjYo31AAAturnm119O8k3bfv5cXYUAAK2rO0h8LsnzNdcAABSq+xqJ\nNyR5LsmFJPNJbq+3HACgFXUGiQ8m+Z4k9yX5wSSNJOeTvKrGmgCAFtR5auPPtj3/pyQfSPJMku9N\n8q5aKgIAWlL3NRLbfTrJPyb5ir0anDx5MkeOHLlq2szMTGZmZjpcGgAcfvPz85mfn79q2uXLlzv6\nmjd1dO2teWmqIxK/leQXdswbT7K4uLiY8fHxrhcGAL1qaWkpExMTSTKRZKnd66/zGolfTvKWVBdY\nfl2SP0jyiiS/U2NNAEAL6jy18dpUd2q8JsmHU10j8fVJnq2xJgCgBXUGCRc2AECPq3scCQCghwkS\nAEAxQQIAKCZIAADFBAkAoJggAQAUEyQAgGKCBABQTJAAAIoJEgBAMUECACgmSAAAxQQJAKCYIAEA\nFBMkAIBiggQAUEyQAACKCRIAQDFBAgAoJkgAAMUECQCgmCABABQTJACAYoIEAFBMkAAAigkSAEAx\nQQIAKCZIAADFBAkAoJggAQAUEyQAgGKCBABQTJAAAIoJEgBAMUECACgmSAAAxQQJAKCYIAEAFBMk\nAIBiggQAUEyQAACKCRIAQDFBAgAoJkgAAMUECQCgmCABABQTJACAYoIEAFBMkAAAigkSAEAxQQIA\nKCZIAADFBAkAoJggAfSdtbW1zM7O5tixYxkbG8uxY8cyOzubtbW1ukuDvnNz3QUAtMv6+nqmp6ez\nsrKSZrN51bzl5eU8/vjjGRsby5kzZzI8PFxTlZRaW1vLqVOnsrCwkI2NjQwNDWVycjJzc3MZGRmp\nu7yBJUgAN9QLb+Dr6+u55557cuHChT3bNJvNNJvNHD9+POfOneuLMNEL2+agBETaYTzJ5uLi4ibQ\nPc1mc3Nqamqz0WhsJrnm0Wg0NqempjabzWbdpW5OTU3tWuNej6mpqbpLPpBe2jYH0Ww2N1//+tfv\na5uOjo72/O/bCYuLi1f+RuOd2EG7RgLY1ZVP+E8++eQ1nwKvaDabefLJJ3P8+PGsr693ucLPW11d\nzcrKSkvLrKys9Ow1E720bQ5qenr6ukeZtnvmmWcyPT3d4YrYSZAAdtVLb+APP/zwnjvUvTSbzZw6\ndapDFXVWL22bgzgsAdHFu9fnGgngGgd5A6/jvPzCwkJXl6tTr22bgzhIQHz00UcP/PquzdgfRySA\na/TaJ/yNjY2uLlenXts2B1FnQByk00cHJUgA1+i1T/hDQ0NdXa5OvbZtDqLOgDgop4/aQZAArtFr\nn/AnJyeLlrvrrrvaXEnn9dq2OYi6AuJhuTajVwgSwDV67RP+3NxcGo1GS8s0Go08+OCDHaqoc3pt\n2xxEXQFxkE4ftYMgAVyj1z7hj4yMZGxsrKVlxsbGeu7iw6T3ts1B1BUQB+n0UTsIEsA1evET/pkz\nZzI6OrqvtqOjozl79myHK+qMw7JtunFLZF0BcZBOH7WDIAFcoxc/4Q8PD+fcuXOZmprac0fbaDQy\nNTWV8+fP55Zbbml7Df28c71ifX099957b+6+++489thjWV5ezoc+9KEsLy/nsccey91335177723\nbXcx1BEQB+n0Ub/4kSSrSf4nyd8mefMubQyRDV3WbDY3R0dH9z008fr6et0l/7/V1dXNEydObN55\n552bd9xxx+add965eeLEic3V1dWOvF63h6uua9vUNVz1fv++7fo9T5w40dJw61ces7OzbXn9duv0\nENl1++4kn0kym+SOJO9K8okkr9vRTpCAGnT7DbwXDcrOdXOz/u8z6VZAXF1d3fPvutej0Wh0LKge\nVL8Hib9O8us7pv1zknfumCZIQI26/Qm/lwzKzvXChQt9tXO9kbq3azv1c5B4SZKNJG/bMf10kid3\nTBMkgENnkHaupYf7T5w4UXfpRXr51N5O/fztn69J8uIkO6/IeT5Ja5ckA9RgkMYbGLRbIg/Dxbu9\nwpd2ARQapJ3rIN4SOTw8nCeeeCJra2s5depUFhYWsrGxkaGhoUxOTmZubq4nxyJptzqDxEeSfC7J\nzq9MG05yabcFTp48mSNHjlw1bWZmJjMzMx0pEOB6BmnnOsi3RI6MjLTl20S7YX5+PvPz81dNu3z5\nckdfs84g8dkki0nuS/LebdPfmuSPdlvg9OnTGR/vu2tFgB41SDvXycnJLC8vt7xcL46o2ct2+3C9\ntLSUiYmJjr1m3QNS/WqSH0hyIsnRVLd/flmS36qzKID9MFz19dU92indUXeQ+L0kJ5PMJXkq1WBU\n35Lk2TqLAtiPQdq51j2iJodX3UEiSX4zye1JXpZkMslf1VsOwP4M2s51UL7PhNYchiAB0LMGaefq\nlkh24/ZPgAO4snOdnp7OysrKruNKNBqNjI2N5ezZsz2/c3VLJDsJEgAHNIg71166JZLOEiQA2sTO\nlUHkGgkAoJggAQAUEyQAgGKCBABQTJAAAIoJEgBAMUECACgmSAAAxQQJAKCYIAEAFBMkAIBiggQA\nUEyQAACKCRIAQDFBAgA
"text/plain": [
"<matplotlib.figure.Figure at 0x7f676b953610>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"\n",
"\n",
"xscale('log'); ylim(-10,10)\n",
"errorbar(fqd, lag, yerr=lage, fmt='o', ms=10,color=\"black\")\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 49,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"(0.095117985155199827, 1.1056585227168947)"
]
},
"execution_count": 49,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgkAAAFkCAYAAACq4KjhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAHwFJREFUeJzt3X98XHWd7/FX2sYWQagiTVpFQ+OWVItgYvjRIqassPcq\nCyiKmb0otrIg668o8lhQiNx2F68sAuLeqwvSqrBMio9d3eralcdiqkubqzXRK5EGsW24ru0E5aao\nQGulvX98E5u0J00mM2dmzpzX8/GYR6cz58en7bft+5zz/QGSJEmSJEmSJEmSJEmSJEmSJEmSJEmS\nJEmSJEkJVlPuAsaYP/KSJEn52TXyKqpKCQnzm5qadg4MDJS7DkmSkmgr8KcUOShUSkhoBnrvu+8+\nFi9eXO5a8tLR0cEdd9yRyPMVcqx8981n+6lsO9k2R/q+o6OD2267g9ZW+MQn4G1vm1JZZWdbG+8n\nP4EVK+Daa6G9fXrnKkVbK+WfWbHY1oq/fZxtbevWrVx22WUALUDflAqaolnFPFihFi9eTHNzc7nL\nyMvcuXNLWnMxz1fIsfLdN5/tp7LtZNsc6fu5c+fS2Bi+O/VUSEqTs62N19wMvb3w938PV1wBixbl\nf65StLWk/ZsGtrU4to+7rcVlZknPNrH5wFVXXXUV8+cnr1vCKaecktjzFXKsfPfNZ/upbDvZNhN9\nn81mWb48w2c/C1deCQsXTrmssrOtjffGN8L998ODD4a7CjNm5H+uuNtaJpOZci2VxLZW/O3jamu7\ndu3irrvuAriLan7c0Nvbm8jUrWS58MILuemm9bS0wJYt8PrXl7siFWLzZnjDG+Dmm+Gv/7rc1Yx3\n4YUXsn79+nKXoSrX19dHS0sLxPC4YUYxDyYlxfBw+PHFLy5vHSrc0qVwzTXQ2QmPPFLuaqTqYkhQ\n6mQyGXbvDu8NCdVh1Sp41avg8sth375yV3NQUh81SKMMCUqdTCbzxzsJxx1X3lpUHHPmwFe+EkY8\n/O3flruagwwJSjpDglJpeBiOPRZmVkrXXRWspQVuuAH+5m/CqAdJhTMkKJWGh33UUI0+8Ql47Wvh\n3e+GPXvKXY2UfIYEpdLu3YaEalRbGx47/PznoSOjpMIYEpRKw8Mwd265q1AcliyB1avh1lth06Zy\nVyMlmyFBqeTjhup2zTVw5plhtMMzz5S7Gim5DAlKJR83VLeZM+FLX4KdO+G668pdjZRchgSlkncS\nqt+iRfDpT4e1HR56qNzVSMlkSFAq2SchHd7/fli+HFauhKefLnc1UvIYEpQ6Bw54JyEtZsyANWvC\nn/dHP1ruaqTkMSQodZ59Fv7wB0NCWjQ0wO23h7DwzW+WuxopWQwJSp3RKZl93JAeK1fCm98Mf/mX\n8NRT5a5GSg5DglLHFSDTp6YG7r4b9u6FD36w3NVIyWFIUOq4AmQ6LVgQRjpks/DVr5a7GikZDAlK\nHe8kpFcmA5dcAldfDUND5a5GqnxxhoRzgG8AvwT2AxfFeC5pyuyTkF41NfD5z4dRD1deGUa6SJpY\nnCHhhcCPgPeP/Ny/jqoIw8Nw1FEwe3a5K1E5nHAC3HUXrF8P995b7mqkyjYrxmP/28hLqihOyayL\nL4Z3vQs+9CE491x4+cvLXZFUmeyToNRxtkUBfPazcMwx8N73+thBmoghQanjbIuC0AbuuQcefDA8\nfpB0OEOCUsfHDRr1Z38WOjBecw1s317uaqTKE2efhLx1dHQw95D7wJlMhkwmU6aKVI2Gh+Gkk8pd\nhSrFrbeGuwnveQ9s3BhGPkiVKpvNks1mx322e3TylxjUxHbk8fYDFwPrJ/i+Gejt7e2lubm5RCUp\nrZYsCZ3V7ryz3JWoUnz3u9DWBrfdBh/5SLmrkfLT19dHS0sLQAvQV8xjx5mZjwZOG3kBLBx5f2KM\n55QmZZ8EHeqNb4SODrj+ehgYKHc1UuWIMyS0EhJNH2GOhNtG3v/3GM8pTco+CYpy881hxcjLLw+r\nhEqKNyRsHDn+DGDmmPcrYzyndES//31YKtohkMk1ODjIypUrOeWUU2hqauKUU05h5cqVDA4OFnTc\no46CL38ZfvhDuOWW4tQqJV1FdVyU4ua6Dck1NDREe3s7AwMD5HK5cd/19/ezYcMGmpqa6Orqoq6u\nblrnOOMMuO46uOkmeMtb4NRTi1C4pmRwcJBVq1axZcsW9u3bR21tLa2trXR2dtLQ0FDu8lLLkKBU\ncQXIZBoaGmLp0qVsP8I4xVwuRy6XY9myZWzatGnaQaGzE77xDXj3u2HLFnjBC6ZbtaaiFOFP0+dg\nH6WKdxKSqb29/YgBYaxt27bR3t4+7XPNng1f+Qo8+iisWjXtw2gKRsPfxo0bDwsIo3K5HBs3bmTZ\nsmUMuXRnyRkSlCquAJk8O3bsYCDPIQcDAwMF9VE47TT45CfhU5+CH/xg2ofRJEoZ/jQ9hgSlio8b\nkmf16tUTXmVOJJfLsarA2wDXXQctLWG0w3PPFXQoRShH+FP+DAlKleFhmDULjj663JVUjrhGCxTL\nli1bSrrfqFmzwmiHHTvghhsKOpQilCv8KT92XFSqjK4AWVOquUYr2NDQEMuXL2f79u3s3bt33Hf9\n/f3cf//9LFy4kO7u7rJ2GNu3b19J9xtr8eIwf8LHPgYXXQTnnFPwIfNSzT3+yxX+lB9DglLF2RaD\nqYwW2Lt3L1u3bi14tEChamtrS7rfoT78Yfj618PaDj/5SVheOm5p6PFfzvCnqfNxg1LF2RaDJHUY\na21tndZ+p59+elHOP3MmrF0LTz4J115blEMeUVp6/Jc7/GlqDAlKFe8kJK/DWGdnJ/X19XntU19f\nz4033li0Ghob4e/+Dr7whbBiZJySFOAKUe7wp6kxJChVRvskpFnSOow1NDTQ1NSU1z5NTU1Ff2b/\nvvfBeefBe997cJRMsSUtwBWiEsKfJmdIUKr4uCGZHca6urpobGyc0raNjY2sW7eu6DXU1MA998Bv\nfxumb37ooaKfInEBrhCVEv50ZIYEpYqPG5LZYayuro5NmzbR1tY24dVnfX09bW1tbN68mXnz5sVS\nx4knwubNUF8Pb3oT/MVfwK5dxTt+EgNcISoh/OnIDAlKFR83JLfDWF1dHd3d3fT09LBixQqWLFnC\nySefzJIlS1ixYgU9PT10d3fHFhBGvfrVsHFjmEPh3/8dmprgc5+D558v/NhJDHCFqJTwp4k5BFKp\n8fzz8PTT3klobW2lv78/7/0qpcNYQ0MDa9asKWsNNTVhAagLLoCPfzwMk1y7NnRsLOS3KakBrhCj\n4a+a54RIMu8kKDV+85vwY9pDgh3GiuclLwnBoKcHDhyAM8+Eq68+uEZIvtLc4380/D3yyCMMDAzw\nyCOPsGbNGgNCmRkSlBquABnYYaz4zjgjLCt9xx3wj/8YHkHce28IDvkwwJVepU9LXm6GBKWGK0Ae\nZIex4ps1Cz70IXjsMTj33PA4YvnysOT0VBngSmd0WvKzzjqLtWvX0t/fz2OPPUZ/fz9r167lrLPO\nYvny5YmdrKpYDAlKDVeAPGi0w9jixYuZPXt25DazZ89m8eLFdhjL0/z5kM2GSZd27oRTT4Xrr4dn\nnpna/ga4+KVlVstiMCQoNXzcMF5dXR2PPvooAwMDkaMFBgYGePTRRw0I03TeeWGth85OuP12eM1r\nYP36yfezx3/80jKrZTE4ukGpMRoSjj22vHVUmkoYLVCt5syBG28M8yl84ANhJckLL4Q774RXvnLi\n/ezxH59CZrVM4++5IUGpMTwMxx0XFuyRSqmxEb71Lfinf4KOjrAEdWcnfPSj8IIXTLyfAa74CpnV\nMo1/Fj5uUGo4JbPKqaYG3v522Lo1DJO84QZ43evgu98td2XpkrZZLQtlSFBqOCWzKsGLXgSf+Qz0\n9oY7W21tcPnlYSlqxS9ts1oWypCg1HBKZlWSU0+Fhx+Gu++Gb34TTj45TMxUjOmdNbE0zmpZCEOC\nUsPHDao0M2bAFVeEuRX
"text/plain": [
"<matplotlib.figure.Figure at 0x7f676881b110>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"s, loc, scale = lognorm.fit(lag,loc=.01)\n",
"xscale('log'); ylim(-4,1.5)\n",
"errorbar(fqd, lag, yerr=lage, fmt='o', ms=10,color=\"black\")\n",
"#plot(fqd,norm.pdf(fqd,mu,sigma))\n",
"plot(fqd,lognorm.pdf(fqd,s,loc,scale))\n",
"mu,sigma\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 50,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"/usr/lib/python2.7/dist-packages/numpy/core/numeric.py:460: ComplexWarning: Casting complex values to real discards the imaginary part\n",
" return array(a, dtype, copy=False, order=order)\n"
]
},
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x7f676847e7d0>]"
]
},
"execution_count": 50,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAg0AAAFkCAYAAACjCwibAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzt3XucleP+//HXjEolNREdsEOKkU3NsCuHCpEoClumZtoi\nx1B03vv7JWyaksohp2x0HJEdQkmiQsKMECX7m9oSFRodqan1++Oz+s00zUxrZta9rnXf6/18PO7H\nzNyzDh9LM/Ne93VdnwtERERERERERERERERERERERERERERERERERERERERERERERERERERERERE\nRCQB3QJ8B+wAPgXOLuO2DYBpwDfAbmCs59WJiIhIxJI9fOzu2B/++4AWwCJgNnBMKbc/GNgQvv3n\nQMjD2kRERCSOLAHGFzv3NfBABPd9FxgT9YpERESkwry60lANSAPmFjs/FzjTo+cUERERD1Xx6HHr\nAQcB64ud34DNXYiWhuFDREREyufH8BExr0JDLDRs1KjRunXr1rmuQ0RExI9+AM6gHMHBq9DwM7YC\non6x8/UpZ6opQ8N169YxZcoUUlNTo/SQ/tS/f3/GjRvnuoy4oNfC6HUopNfC6HUopNcCli9fTmZm\n5lHY1XrnoWEnkAtcCLxa5PwFwMxoPlFqaippaWnRfEjfSUlJSfjXYC+9FkavQyG9FkavQyG9FhXn\n5fDEGGAy1p/hI+AG4GjgyfD3RwCNgL8VuU+L8MdDgSPDX+/EVl2UaNu2qNYsIiIipfAyNLwIHA7c\nhV3++BK4GPg+/P0G7N+zIS/8MYStvugBrAaOL+1J3n0XzjknajWLiIhIKbxs7gTwBHAcUB2bbPF+\nke/1Bs4roZ5kbOXF3s9LDQwAb74ZrVJFRESkLF6HBs99/DEk+gKKjIwM1yXEDb0WRq9DIb0WRq9D\nIb0WFZfkuoBKSANyq1bNZcSINAYMcF2OiIiIP+Tl5ZGeng6QTuHUgAPy/ZWGtm1h8mTXVYiIiASf\n70PDJZfA55/Dl1+6rkRERCTYfB8azjwTDj8cpk51XYmIiEiw+T40VK0K3btbaNizx3U1IiIiweX7\n0ACQmQlr18KCBa4rERERCa5AhIbWraFJE02IFBER8VIgQkNSkl1tmDEDduxwXY2IiEgwBSI0APTs\nCVu2wGuvua5EREQkmAITGpo2tWGKKVNcVyIiIhJMgQkNYEMUc+bAxo2uKxEREQmeQIWG7t3t4/Tp\nbusQEREJokCFhnr1oFMnDVGIiIh4IVChASArC5YsgZUrXVciIiISLIELDZ07Q+3aaistIiISbYEL\nDTVqwJVX2hBFKOS6GhERkeAIXGgAG6JYtQoWL3ZdiYiISHAEMjS0bQvHHKO20iIiItEUyNCQnGwd\nIqdPh507XVcjIiISDIEMDWCNnjZtgjffdF2JiIhIMAQ2NDRvDi1bqmeDiIhItAQ2NIBdbZg1C/Lz\nXVciIiLif4EODRkZUFAAL73kuhIRERH/C3RoaNgQOnTQEIWIiEg0BDo0gA1RLFwIq1e7rkRERMTf\nAh8aunWDmjVh2jTXlYiIiPhb4ENDrVoWHCZPVltpERGRygh8aABrK71iBeTlua5ERETEvxIiNJx/\nPtSvrwmRIiIilZEQoaFKFejRA3JybAmmiIiIlF9ChAawVRTr18O8ea4rERER8aeECQ0tW0Jqqna+\nFBERqSivQ8MtwHfADuBT4OwD3L4dkBu+/f8BN0arkKQkmxA5cyZs2RKtRxUREUkcXoaG7sBY4D6g\nBbAImA0cU8rtjwPeBBaEb/8A8AhwebQK6tEDduyw4CAiIiLl42VouBN4BngW+Aa4A/geuLmU298E\nrA7f7xvgX+H7DoxWQY0bQ7t2WkUhIiJSEV6FhmpAGjC32Pm5wJml3KdNKbc/HTgoWoVlZsI778C6\nddF6RBERkcTgVWioh/2hX1/s/AagQSn3qV/C7dcDVcKPFxVXXglVq9rySxEREYlcwqye2CslBbp0\n0RCFiIhIeVXx6HF/BnZjVw+Kqg/8WMp9fmL/qxD1gYLw45Wof//+pKSk7HMuIyODjIyMUovLyoLL\nLoNly+CUU0q9mYiIiO/l5OSQU+zyen5+foUeKykaBZXiI2z5ZN8i574GZgL/KOH22UAXoHmRc08A\npwJnlXD7NCA3NzeXtLS0chW2cyc0bAjXXw/Z2eW6q4iIiO/l5eWRnp4OkA5EvDOTl8MTY4A+QG8g\nFVt+eTTwZPj7I4CJRW7/JNAYeCh8+2vDx+hoF1atGnTvDlOnwp490X50ERGRYPIyNLwI9AfuAj7D\nGjtdjC27BBuKKNqzYXX4++3Dt/8HcBt2ZSLqsrJg7VpYsMCLRxcREQker+Y07PVE+ChJ7xLOLcQu\nlXiudWto0sTaSp97biyeUURExN8SbvXEXklJ1rNhxgzrEikiIiJlS9jQANCzp+1D8dprrisRERGJ\nfwkdGpo2tWEK9WwQERE5sIQODWBDFHPmwMaNrisRERGJbwkfGrp3t4/Tp7utQ0REJN4lfGioVw86\nddIQhYiIyIEkfGgAG6JYsgRWrnRdiYiISPxSaMA2sKpd2zpEioiISMkUGoAaNWzL7ClTIBRyXY2I\niEh8UmgIy8qCVatg8WLXlYiIiMQnhYawtm3hmGOsrbSIiIjsT6EhLDnZOkROn25bZ4uIiMi+FBqK\nyMyETZtg9mzXlYiIiMQfhYYimjeHFi00RCEiIlIShYZisrJg1izIz3ddiYiISHxRaCgmIwMKCuCl\nl1xXIiIiEl8UGopp2BA6dFBbaRERkeIUGkqQmQkLF8Lq1a4rERERiR8KDSXo1g1q1oRp01xXIiIi\nEj8UGkpQq5YFh8mT1VZaRERkL4WGUmRlwYoVkJfnuhIREZH4oNBQivPPh/r1NSFSRERkL4WGUlSp\nYssvc3JsCaaIiEiiU2goQ1YWrF8P8+a5rkRERMQ9hYYytGwJqalqKy0iIgIKDWVKSrKrDTNnwpYt\nrqsRERFxS6HhAHr0gB07LDiIiIgkMoWGA2jcGNq10yoKERERhYYIZGbCO+/AunWuKxEREXFHoSEC\nV14JVava8ksREZFEpdAQgZQU6NJFQxQiIpLYFBoilJkJS5fCsmWuKxEREXFDoSFCnTrBYYfpaoOI\niCQuhYYIVasG3bvD1KmwZ4/rakRERGLPq9BQF5gM5IePSUCdA9zncmAu8AuwBzjVo9oqLCsL1q6F\nBQtcVyIiIhJ7XoWGadgf/Y7ARUALLESUpSawEBjsUU2V1ro1NGmittIiIpKYvAgNqVhY6AMsAT4C\nrgc6A83KuN8U4J/AOx7UFBVJSTYhcsYM6xIpIiKSSLwIDW2A34BPipxbEj7XxoPni6mePW0filmz\nXFciIiISW16EhgbAhhLObwh/z9eaNrVhCg1RiIhIoqlSjtsOB+46wG3OqHgpFdO/f39SUlL2OZeR\nkUFGRoZnz5mZCf37w8aNcMQRnj2NiIhIpeXk5JBTrKVxfn5+hR4rqRy3PTx8lGUN0BN4CFtBUdQm\noD8w8QCPcSywCps8+UUZt0sDcnNzc0lLSzvAQ0bXzz9Dw4YwdizcemtMn1pERKTS8vLySE9PB0gH\n8iK9X3muNPwSPg5kMba88gwK5zW0Cp/7sBzPF7fq1bNmT1OmKDSIiEji8GJOw3JgDjABCwutw5/P\nAr4tcrsVQNciX9fFri6cHP76pPDX9T2osdIyM2HJEli50nUlIiIiseFVn4YewJdYs6a3gKVAVrHb\nNANqF/n6MuwSyetACHgh/PWNHtVYKV26QO3a1iFSREQkEZRneKI88tk/JBRXPLA8Hz58oUYN2zJ7\nyhQYPtx6OIiIiASZ9p6ohKwsWLUKFi92XYmIiIj3FBoqoW1bOPpo7XwpIiKJQaGhEpKTrUPk9Omw\nc6frakRERLyl0FBJWVnw668we7brSkRERLyl0FBJzZtDixZqKy0iIsGn0BAFWVm2gVUFu3KKiIj4\ngkJDFGRkQEEBvPSS60p
"text/plain": [
"<matplotlib.figure.Figure at 0x7f6768614950>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plot(ifft(lag))"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 2",
"language": "python",
"name": "python2"
},
"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": 2
}
Reference in New Issue Copy Permalink