Skip to content

Instantly share code, notes, and snippets.

@alekseygenerozov
Created March 4, 2014 01:16
Show Gist options
  • Save alekseygenerozov/9338349 to your computer and use it in GitHub Desktop.
Save alekseygenerozov/9338349 to your computer and use it in GitHub Desktop.
{
"metadata": {
"name": ""
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "code",
"collapsed": true,
"input": [
"import nuker\n",
"from astropy.io import ascii\n",
"\n",
"galaxies=nuker.nuker_params()\n",
"rad=np.logspace(0,5,100)\n",
"grad_phi1=nuker.get_grad_phi(galaxies['NGC4551'])\n",
"\n",
"loglog()\n",
"plot(rad, -np.array(map(grad_phi1, rad)))\n",
"\n",
"# phi=np.empty_like(rad)\n",
"# for i in range(len(phi)):\n",
"# phi[i]=nuker.get_phi(rad[i], galaxies['NGC4551'])"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": "*"
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"rad=np.logspace(0,5,100)\n",
"loglog()\n",
"plot(rad, map(nuker.get_rho(galaxies['NGC4551']), rad))"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 25,
"text": [
"[<matplotlib.lines.Line2D at 0x6afcad0>]"
]
},
{
"metadata": {},
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAX8AAAEHCAYAAABGNUbLAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xt0lNW5x/FvQNAoVKCmaEgUDSVIEZBbVZgGBEGkgKhA\nvABFKAInoB4vIBRDbERAlKVB5FYBAUMQgUCDCRCNoq1B4XDQhsiliZdwrF1UikYJJJnzxy5JIAmZ\nyVzed2Z+n7WyJNPJO4+7uLPf/T77ecKcTqcTEREJKQ2sDkBERPxPk7+ISAjS5C8iEoI0+YuIhCBN\n/iIiIUiTv4hICNLkLyISgjT5i4iEoIt8eXGn08kf/vAHvv/+e7p168bo0aN9+XEiIuIin678t2zZ\nQlFREY0bNyYqKsqXHyUiIm5we/J/8MEHadmyJTfccMM5r2dmZtKuXTt++ctfMm/ePAAOHTpEz549\nWbBgAa+++qp3IhYREY+5PfmPHTuWzMzMc14rKysjISGBzMxM8vLySE1N5eDBg0RFRdGsWTPzQQ30\neEFExC7cnpEdDgfNmzc/57U9e/bQpk0bWrduTaNGjYiPjyc9PZ277rqLrKwspk6dSu/evb0Vs4iI\neMgrD3yLioqIjo6u+D4qKorc3FzCw8NZsWLFBX+2TZs2HD161BthiIiEjJiYGI4cOVLvn/fKXkxY\nWFi9f/bo0aM4nU6vfiUmJnr1vbW9x9XXL/R9bX/WWATXWLjymsbC/2Ph7vXsNBaeLpobzp49e7a7\nP3TixAlSU1OZPHkyACdPnmTr1q088MADAGRkZHD55ZfTq1evOq+VlJRU8efWrVu7G0qt3LmWK++t\n7T2uvn6h78/+OScnxyfbYxqL2j/b0/e6MxauvKaxqPl7X46Fu/OO1WNRWFjIqlWreO+996jH9F3J\nWQ8FBQXODh06VHx/5swZ53XXXecsKChwlpSUODt16uTMy8tz6Vr1DCEoJSYmWh2CbWgsKmksKmks\nKnk6d7q97XPvvfdyyy23cOjQIaKjo1m5ciUXXXQRixYtYsCAAbRv356RI0dy/fXXu3zN2bNnk5OT\n424oQUcPxStpLCppLCppLMzdj0cr/v8I+89vEMuEhYVhcQgiIgHH07lTyfciIiFIk7+ISAiyxeSv\nPX8REddoz19EJIRpz19ERNymyV9EJATZYvLXro+IiH/ZYvLv1Qs+/NDqKEREQoctJv+IiNkMG5bD\nXXfB559bHY2IiH0FRLZPTk4Os2bNokOHDsTHxxMXF1c9gP88sf7pJ0hJgeefh3vugcREuPJKX0Um\nIhLYbJ3t06BBA5o2bUpJSUmdPXzDw+HJJyE/Hy69FH71K5g9G77/3pcRioiEJp/28HU4HGzfvp25\nc+eSmJjo0vV//nN44QXYuxeOHIG2beHVV+HMGXcjFRGR2vi0h+/ZJi/NmjWjpKTErc9p3RrWroWM\nDHjrLejQATZtUmaQiIg3uN3G0eFwUFhYeM5rVXv4AhU9fPPz88nKyuLEiRNMmTKlXgF26QI7d8KO\nHWZb6IUXYP586NmzXpcTERF83MN3+vTpDBs2rM6fr/rkunfv3tVqdoeFwYAB0K8frFsH990HXbvC\nc89BbKw3/g1EROwtJyfHqzXQvDL5e9LDF3A5balhQxg9GkaMMJlBvXopM0hEQsP5C+OqLXDrwyvZ\nPq1ateKrr76q+P6rr76qM7unKnerel5yCTzxRPXMoB9+cCNoEZEAZGmef2FhIYMHD+bTTz8FoLS0\nlNjYWLKzs4mMjKRHjx6kpqa61MrRG1U9CwvhD3+A7Gx4+mkYPx4aNfLokiIitub3PH9f9PD1VE2Z\nQZs3KzNIRKQ2tqjnn5iYWOOD3vpwOiszg5o0UWaQiASXsw9+k5KSPFr522Ly90UIZWXmbmDWrMrM\noHbtvP4xIiKWsHV5Bys1bAhjxphCcbfcAg4HTJwI33xjdWQiItazxeTvyx6+4eEmM+jzz+Gyy5QZ\nJCKBLSCqeroUgJ97+J7NDHrnHZMZNG6cMoNEJPB4OneG3OR/1r595qHwV1/B3Llw553mJLGISCDQ\n5O8BZQaJSKDSA18PnK0ZtG8fTJgA994Lw4aZk8MiIsEspCf/s6pmBt18szKDRCT42WLy92W2jzuq\ndhM7mxmUmKhuYiJiH8r28QPVDBIRu7L9nn9xcTHdu3cnIyPD1x/ldeomJiLByucr/8TERJo2bcr1\n11/PoEGDqgdg45V/VVUzgy67DJ5/XplBImIdv6/83WngvnPnTtq3b09ERES9A7SLqplBEyeabmJ3\n3qnMIBEJTG6v/Hfv3k2TJk0YPXp0RT3/srIyYmNj2bVrF61ataJ79+6kpqaybt06iouLycvLIzw8\nnM2bN1fr+hUoK//znToFL79s7gDuvtuUjFA3MRHxF0/nTp82cE9OTgZg9erVRERE1Nrusa4evnZ0\nySVmC2j8eJgzx2QGTZkCjz0GTZtaHZ2IBBtv9/D1SievjRs3kpWVxfLlywFYu3Ytubm5pKSk1B1A\ngK78z1dYaMpH79qlzCAR8T1bZPt4o4G7HfL8PdG6NaxZA9u3m4wgZQaJiC94K8/fNg3cA2GrxxU3\n3gg7d0JKCjzzjMkI+vBDq6MSkWDRu3dv+0z+3bp14/DhwxQWFnL69GnS0tIYMmSINy4dsPr3N5lB\nkyaZzCDVDBIRO7FFA/dg2PapSYMGMGqUqRnUs6dqBomI51TeIQD9618mM2jlSkhIgMcfV2aQiNSP\nLR74imtatIAFC2DvXjh6FNq2hcWL4cwZqyMTkVBji8k/WLd9anO2ZtD27bBlizkj8NZbygwSkbpp\n2yeInK0ZdOmlpptYr15WRyQidqc2jkGivBzWrTMHxTp3hueeAzeemYtIiNGef5A4mxmUn2+yguLi\nTGvJ//s/qyMTkWCkyd9mLrnE1Af6/HO4/HJzUvjpp9VNTES8S5O/TTVvbiqG7ttn6gb98pewaBGc\nPm11ZCISDDT529w118Drr0NmJmzbZjKDNm5UZpCIeMYWk3+opXrWR+fOkJVlzgU8+yzcfDO8/77V\nUYmIvwVEqmd+fj4vvfQSx48fZ8CAAYwbN656AMr2cVt5Obzxhmku37EjzJ0L7dtbHZWI+FNApHqW\nl5cTHx/Phg0bqgegyb/eTp0ydwJz58LQoZCUBJGRVkclIv5g6x6+ANu2bWPQoEHEx8fXO0ip2SWX\nwH//t8kMatECbrjB3A2cPGl1ZCJidz7t4Vu1sufQoUNJT0+vHoBW/l7z5ZcmLTQzE2bMMBVEGze2\nOioR8QVb9/D99ttv2bRpE6dOnaJPnz71DlJcc/XVsGoVHDgA06bBSy+ZKqLDh5tDZCIiZ7k9+dek\nqKiI6Ojoiu+joqLIzc0lLi6OuLi4On8+EBu421nHjvD22/DOO6Zm0IIFpmaQfv+KBC5vN3D3yuTv\njR6+4n233gp79kBaGowbB+3awbx55tmAiASW8xfGSUlJHl3PNj18lefvGw0awL33wsGDMGAA9OsH\nY8dClf+7RCSAWJrnX1hYyODBgyse+JaWlhIbG0t2djaRkZH06NGj2gPfWgPQA1+/+ve/zRbQkiXw\n+9/D9OnQrJnVUYmIu/ye6qkevoHt8svNCeEDB+D4cdNN7MUXoaTE6shExBUBccLXpQC08rdUXp5Z\n/R84AMnJcN99ygwSCQQBccL3ggFo8reF3bvhiSfMHcC8edC/v9URiciFaPIXr3E6YdMmeOopc2Zg\n/nzo0sXqqESkJkHRyUt7/vYQFgZ33w1/+5v556BBZhvo73+3OjIROUt7/uJzP/xgHga//DI88ADM\nnAkREVZHJSIQJCt/sacmTUytoLw8KCszDeWTk6G42OrIRMRTmvylTr/4BaSkwEcfwWefmZaSS5fC\nmTNWRyYi9aVtH3HbJ5+YwnFff20Kx911l3leICL+o2wfsYTTCTt2mDMCF19sMoN+8xuroxIJHUGx\n569sn8ATFmZqBe3dC1OmwJgx8Nvfwn8qfoiIjwREtk96ejoZGRmcPHmScePGcdttt1UPQCv/oFBS\nYuoFzZkDt98OzzwD11xjdVQiwSsgtn1OnDjB448/zooVK6oHoMk/qJw8afoHvPKKuRuYMQOuuMLq\nqESCj+17+AIkJyeTkJBQ7yAlcPzsZ2bV/7e/mbuBdu2UHipiR25P/mPHjiUzM/Oc18rKykhISCAz\nM5O8vDxSU1M5ePAgTqeTadOmMXDgQDp37uy1oMX+rrzSrP6rpocuXqz0UBG7cHvydzgcNG/e/JzX\nqvbwbdSoUUUP30WLFpGdnc3GjRtZunSp14KWwNGmDaxfD9u2wZYt5qDY+vVQXm51ZCKhzac9fFNS\nUpgyZUqdP68evsGva1eTGpqdbdJD58+H554z1UN1RkCkburhKwGtb1/TV3jTJnj4YbjqKvNL4Kab\nrI5MxN7Uw1cC3tnqoZ99BvffD8OHw513mofEInJh6uErQeOnn8zD4HnzYOBASEqC1q2tjkrE3tTD\nVwJeeDg89hgcPmwOhnXtak4Nf/ON1ZGJ2E9AnPB1KQCt/OU8335rngO8/jo89JBpL3legplIyAuK\n2j4iVf3iF7BwIfzP/5hfBG3bml8GOigm4j2a/MW2rr4aVqyADz6A/fvNmYGXXzYnh0XEM7aY/LXn\nLxcSGwtpafD22+asQNu28Kc/QWmp1ZGJ+J/2/CVk/eUvpp9wUZGpIzRiBDSwxTJGxH8CoqrnBQPQ\n5C/14HSa08IzZ5pU0eRkGDxYp4UldGjyl5DmdJq6QbNmwSWXmF8C/frpl4AEP03+IphCcRs2wNNP\nQ2QkPPss9OxpdVQivqNUTxHMnn98POTlwejRcN99MGiQSRcVkepsMfkr20e85aKL4MEH4dAhUypi\n0CDzQDg/3+rIRLxD2T4iLiguhkWLTGvJwYMhMVG9hSU42Hrbp6CggPHjxzN8+HBffoxIrS67DKZN\nM3WDIiOhSxd49FH45z+tjkzEWj6d/K+99toam7aL+FuzZiYTKC/PHA5r1w5mz4bvv7c6MhFr+KWB\nu4hdtGwJKSnw8cdw5Ig5LfzKK+otLKHHpw3cRezquutg7VrYvh3S06F9e3jrLXNuQCQUuN3G0eFw\nUFhYeM5rVRu4AxUN3Fu2bMmMGTPYv38/8+bNY9q0aTVeUz18xSo33mjqBe3cCY8/Di++CC+8oLaS\nYj/e7uHrlU5eGzduJCsri+XLlwOwdu3aigbudQagbB+xibIy00Ng1izo1ct0FlNmkNiVLbJ9vNHA\nXXn+YrWGDWHsWPj8c7j+epMZNGsW/PCD1ZGJVPJWnr9tGrhrq0fs4rLLzHmA/fvh7383mUHr1ul5\ngNhD79697TP5d+vWjcOHD1NYWMjp06dJS0tjyJAh3ri0iGWio82kv2GDeRbgcKhchAQPNXAXqcMt\nt8CePTBmjCkZMXkyfPed1VFJqFJ5BxELfPed6SGwebPpKzx6tBrJiDVU0lnEAp98Yu4AGjeGJUug\nQwerI5JQY4tsH09p20cCTbdu8Ne/wv33Q58+MH06/Pij1VFJKNC2j4hNfPONKRb30UeweLF5LiDi\na9r2EbGJHTtg4kRzOnjhQlNHSMRXgmLbRyQY9O8Pn31mUkRvuAFWrtTZALEvrfxFfGD/fhg/3pSS\nXrbMFJIT8Sat/EVsqHNn8wzg9tuhRw/TSay01OqoRCrZYvJXto8Eo4suMpVCc3NN6ehbbjHbQiKe\nULaPSABxOmHFCpgxAxIS4KmnzBkBkfqydbZPcXExkydP5uKLL6Z3797cd9991QPQ5C8h5OuvYdIk\n+OIL80C4a1erI5JAZes9/02bNjFixAiWLVvG1q1bfflRIgEhKgq2bjVN5e+4w9wJnDpldVQSinza\nw7eoqIjo6GgAGjZs6IVwRQJfWJg5GXzgABw6ZLqJffSR1VFJqPFpD9+oqKiKOv/l5eXeiVgkSLRs\nCRs3wjPPwLBh8MQT8NNPVkclocKnPXynTp1KQkICGRkZF6zvrx6+EsqGD4fevWHKFJMiunKlyQwS\nqUo9fEWC2FtvmWyg+++HP/4RwsOtjkjsyhYPfNXDV8Q77r4bPv3UZAV17gx/+YvVEYndqIevSJC6\n4gpYvx7mzDG/DPQsQKpSD1+RIHf33SYj6MsvlREk3qceviI2FhEBaWmQnGwygqZN07mAUKfyDiIh\n5ttvTevIvDxYvRq6d7c6IrGSrcs7uBSAJn8RlzmdsGEDPPwwPPggJCbCxRdbHZVYwRbZPiLiH2Fh\nMHKk6ReQl2d6Ce/bZ3VUEohsMflrz1/EPVdeCZs3m8bxAweaO4DTp62OSvxBe/4iAsCxYzBhAhQV\nwapV0KmT1RGJP2jbRyTERUbCtm3mOcBtt5nMoDNnrI5K7E4rf5Eg8vXXMG4cHD9uMoJ+9SurIxJf\n0cpfRCpERUFmptkGiouDefPUO1hqppW/SJAqLDTpoD/+aO4CYmOtjki8KShW/sr2EfG+1q1h1y54\n4AHo2RMWLgS11Qh8AZHtU1BQwLPPPsu///1v3nzzzZoD0MpfxOeOHIGxY805gZUrISbG6ojEU7Ze\n+V977bWsWLHClx8hIi5o0wZyckx9oJtugsWLdRcQ6lya/N3p2ysi9tSwITz6KOzeDa+/Dv37wxdf\nWB2VWMWlyd+dvr1r1qzh0Ucf5dixYz4JWEQ8064dfPAB9OtnykOsWGFqBklocWnydzgcNG/e/JzX\nqvbtbdSoUUXf3lGjRrFw4UIiIyP517/+xcSJE9m/f7/uDERs5KKLTGmId9+FV1+FO+4wZwQkdLjd\nwP2soqIioqOjK76PiooiNzf3nPe0aNGCJUuW1HktNXAXsUaHDqZJzHPPQZcusGABjBplHgyLvXi7\ngXu9J39P+/ZW5Y20JRGpn0aN4OmnYfBgGDPGNJFfutQUjxP7OH9hnJSU5NH16p3t42nf3qqU5y9i\nvRtvhI8/NncDnTqZDmJiP37P8y8sLGTw4MF8+umnAJSWlhIbG0t2djaRkZH06NGD1NRUt9o3gvL8\nRexozx5zF3DDDfDKK6adpNiLX/L8fdG3tyqt/EXspUcP0yTm6quhY0fTO0DsISBO+LoUgFb+Irb2\nwQfwu9+Zw2EpKXBe4p9YxNYnfEUk8PXqBf/7v9CihdkG2r7d6ojEG7TyFxGXvfuuqRR6663w4otw\n+eVWRxS6gmLlrz1/kcDQpw8cOGAOiXXsaKqGin9pz19ELLVjB4wfD7/9LcyfD02aWB1RaAmKlb+I\nBJ7+/c1dwI8/mnMB779vdUTiDq38RcRjW7fCxIkwciTMmQPh4VZHFPy08hcRyw0ZAp9+Ct98A507\nm3pBYm9a+YuIV23cCAkJ5mxAUhJcfLHVEQWnoFj5K9tHJHjcc495FnDoEHTtCnv3Wh1RcAmYbJ/0\n9HQyMjI4efIk48aN47bbbjs3AK38RYKS0wmpqaZ72MSJMHMmNG5sdVTBw9O502/bPidOnODxxx+v\n1tNXk79IcDt2DCZMgKIiWL3anA8Qz/lt28fTPr7JyckkJCTUO1ARCUyRkbBtG0ydCn37wrPPQmmp\n1VGJyyv/3bt306RJE0aPHl1R1rmsrIzY2Fh27dpFq1at6N69O6mpqXzyySfs27ePJ554gquuuorp\n06fTv39/+vbtWz0ArfxFQsaXX5qDYSdOmLuAehYCFvy48q9vH9+UlBSys7PZuHEjS5curXegIhL4\nrr4asrJMfSCHw7SNLCuzOqrQVO82juBaH9+pU6cyderUC15HPXxFQkdYmHkA3L8/jB0LW7bAqlXQ\npo3VkdmbbXr4gvf6+KqHr0joue46UyU0JQVuvhkSE2HyZGhgiwR0+7FND1/wXh9f5fmLhKYGDeDh\nh+HDD2HdOujXDwoLrY7K3izJ8/dFH1898BURMHv/L7wAzz9v6gONH2+2iKRmfnvg6+s+viIS2ho2\nhCefhJwcWLoU7rgDvv7a6qiCly1q+yQmJupBr4hUOHMG5s41zwMWLIBRo3QXcNbZB79JSUmBccK3\n1gC07SMitdi/H8aMgdatzd3AlVdaHZF9BEVhNxGRmnTuDB9/bBrHd+oEaWlWRxQ8bDH5K9tHRGrT\nuDEkJ8Of/2xKRI8YAf/8p9VRWSdgqnrWGYC2fUTERadOwaxZsHYtLF4Mw4ZZHZF1AqaqZ60BaPIX\nETd9+KFpFvPrX5uHwudVngkJ2vMXkZDTs6d5GNy8uXkesH271REFHq38RSSgvfuuKRTXty+8+CL8\n7GdWR+QfWvmLSEjr08e0jWzY0DSKyc62OqLAYIvJX9k+IuKJpk3NOYClS02l0MmT4YcfrI7KN5Tt\nIyJSgxMn4JFHYPduWLkSfvMbqyPyDVtn++Tn5/PSSy9x/PhxBgwYwLhx46oHoMlfRHxg61aYNMmc\nC5gzB8LDrY7Iu2w9+Z9VXl5OfHw8GzZsqB6AJn8R8ZHjx2HKFNi717SNvOkmqyPyHr888PWkefu2\nbdsYNGgQ8fHx9Q5SRKQ+fv5zeOMN0zT+zjth+nQoKbE6KntwaeVf3+btkZGRFdcYOnQo6enp1QPQ\nyl9E/OAf/zDbQIcOmbuArl2tjsgzns6dLrVxdDgcFJ7XXqdq83agonn79OnTGTVqFADvvfcemzZt\n4tSpU/Tp06fW66uHr4j4WsuW8NZb5k5g4EDzi2DmTFM7KBB4u4evy3v+53fx2rhxI1lZWSxfvhyA\ntWvXkpubS0pKinsBaOUvIn527Bj8/vfmn6+/bk4JBxrLDnl5q3k7KM9fRPwrMtJUCZ0yBW691WQD\nlZZaHZVrvJXnX+/J31vN28FM/trqERF/CgszZSH27jUlInr2hIMHrY6qbr1797Z28u/WrRuHDx+m\nsLCQ06dPk5aWxpAhQzwOSETEn66+GnbsMB3DHA7TRL6szOqofM+lyd/Xzdu17SMiVgoLMyUhcnMh\nPR3i4uDIEaujqpnKO4iI+EB5Obz0kjkbMHu2+aXQwBZV0M4VECd8LxiAJn8RsaHPPzdbQZdeCq+9\nZprI20lQlHTWto+I2E1sLHzwAfTvD927w4oVYId1qrZ9RET85LPPYPRouPJKWL4cWrWyOqIgWfmL\niNhZhw7mYfCvfw033ghr1tjjLsATWvmLiLhh3z7zLCAmxjSPadnSmji08hcR8aMuXeCTT6B9e+jU\nCWqoVB8QtPIXEamnPXvMXUDHjvDKK3DFFf77bK38RUQs0qOH2QaKijK/AGqoWm9btpj8leopIoEq\nPNyUhEhLg8ceM1lB333nu88LmFTP4uLiikJEgwYNqh6Atn1EJEgUF8O0abBli0kJHTjQd59l+22f\n+fPnM3LkSF9/jIiI5S67DBYtMp3CJk0yPQNOnrQ6qpr5tIfvzp07ad++PREREd6JVkQkAPTtCwcO\nmD937AjvvGNtPDXxaQ/fxYsXU1xcTF5eHuHh4WzevLlaExht+4hIMHv7bZgwAYYOhXnzzN2BN/it\nsNv5bRz/+te/kpSURGZmJgBz584FYPr06dV+dvXq1URERHDHHXd4/V9ARMTuvvsOHnkEPvwQVq2C\nXr08v6ZfGrjXpKioiOjo6Irvo6KiyM3NrfG9Y8aMueC11MBdRIJZ8+bmOUB6OowYAffeC8nJJlPI\nVd5u4F7vyd/bPXxFRILd0KGmXWRCgqkRtHq1qRfkivMXxklJSR7FYpsevsrzF5FQcMUVsH49PPMM\nDBkCTz0FJSWu/7zf8/zP3/MvLS0lNjaW7OxsIiMj6dGjB6mpqW63ctSev4iEqn/8Ax56CI4eNXcB\nXbq4/rN+yfNXD18REe9r2RI2bzYHw26/3bSNPHPmwj8TMCd86wxAK38REYqKzKGwb74xdwHnHauq\nxvYnfEVEpG6tWkFGBvzXf8Gtt8Jzz0Fpqe8+Tyt/ERGb+eILGDcOfvjBnAto1676e4Ji5a89fxGR\nStdcAzt2mAqhDgcsXAjl5eZ/056/iEgIOHoUxo41f1650rSPhCBZ+YuISM1iYiAnB4YNg5tugsWL\nK+8CPKGVv4hIgMjPh9/9Dpo0gexsrfxFREJCu3bwwQfQr5/n19LKX0QkAAXFnr+yfUREXBMQ2T45\nOTnMmjWLDh06EB8fT1xcXPUAtPIXEXGbrVf+DRo0oGnTppSUlNS74mco0d1PJY1FJY1FJY2F9/i0\nh6/D4WD79u3MnTuXxMRE70QcxPQXu5LGopLGopLGwntcmvzHjh1b0a7xrLKyMhISEsjMzCQvL4/U\n1FQOHjzImjVrePTRRzl27FhFw5dmzZpR4k7Bag+58xfElffW9h5XX7/Q977+y6yxqP2zPX2vO2Ph\nymsai5q/9+VYuHvtYBoLlyZ/h8NB8+bNz3ltz549tGnThtatW9OoUSPi4+NJT09n1KhRLFy4kMjI\nSDZv3szEiRMZPXo0U6ZM8VrQddGEV/tne/pejUXd77Hbf+Q10VjU79pBNRZOFxUUFDg7dOhQ8f2b\nb77pHD9+fMX3a9ascSYkJLh6uQoxMTFOQF/60pe+9OXGV0xMjNvzbVWW9/A9cuSIV64jIiKus0UP\nXxER8a96T/7dunXj8OHDFBYWcvr0adLS0hgyZIg3YxMRER+xRQ9fERHxL8tr+4iIiP/ZoraPiIj4\nl+0m/+LiYsaMGcOECRN44403rA7HUgUFBYwfP57hw4dbHYrl0tPTmTBhAvHx8ezcudPqcCyVn5/P\npEmTGDFiBH/605+sDsdyxcXFdO/enYyMDKtDsVROTg4Oh4NJkybx3nvv1fl+203+mzZtYsSIESxb\ntoytW7daHY6lrr32WlasWGF1GLYwdOhQli1bxpIlS0hLS7M6HEu1a9eOV199lfXr15OVlWV1OJab\nP38+I0eOtDoMy7lbS80vk787tYGKioqIjo4GoGHDhv4Iz6/qWycpGNVnLJKTk0lISPBnmH7h7lhs\n27aNQYMGER8f7+9Qfc6dsdi5cyft27cnIiLCilB9zp2xcLuWmkdHxFz0/vvvO/ft23fOCeHS0lJn\nTEyMs6CgwHn69Glnp06dnHl5ec41a9Y4//znPzudTqczPj7eH+H5lTtjcdY999xjRag+585YlJeX\nO5988kmEaJJJAAAB20lEQVTnrl27LIzYd+rz98LpdDqHDBni71B9zp2xmDlzpvORRx5x9u/f3zl0\n6FBneXm5hZF7X33+XpSUlLg0Z9T7hK87HA4HhYWF57xWtTYQUFEbaOrUqSQkJJCRkRGU5wbcGYuW\nLVsyY8YM9u/fz7x585g2bZr/A/Yhd8Zi165dZGdnc/LkSY4cOcJDDz3k/4B9yJ2x+Pbbb9m0aROn\nTp2iT58+/g/Wx9wZi+TkZABWr15NRESE1yoP2IU7Y5Gfn09WVhYnTpxwqZaaXyb/mlTd3gGIiooi\nNzeXSy+9lNdee82qsCxR21i0aNGCJUuWWBiZ/9U2FikpKX4tDmgHtY1FXFxcjY2RglltY3HWmDFj\nrAjLErWNxfTp0xk2bJjL17HsgW+w/Yb2hMaiksaiksaiksaikrfGwrLJX7WBKmksKmksKmksKmks\nKnlrLCyb/FUbqJLGopLGopLGopLGopLXxsInj6jPEx8f77zqqqucjRs3dkZFRTlfe+01p9PpdG7f\nvt3Ztm1bZ0xMjHPOnDn+CMVyGotKGotKGotKGotKvhwL1fYREQlBtjvhKyIivqfJX0QkBGnyFxEJ\nQZr8RURCkCZ/EZEQpMlfRCQEafIXEQlBmvxFRELQ/wN5xaB3DSRS4wAAAABJRU5ErkJggg==\n",
"text": [
"<matplotlib.figure.Figure at 0x6afcd50>"
]
}
],
"prompt_number": 25
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"from scipy import integrate\n",
"\n",
"rho=nuker.get_rho(galaxies['NGC4551'])\n",
"r=np.logspace(0,3,100)\n",
"\n",
"f=lambda r1: 4.*np.pi*r1**2.*rho(r1)\n",
"menc=np.empty_like(r)\n",
"\n",
"for i in range(len(rad)):\n",
" menc[i]=integrate.quad(f, 0, r[i])[0]\n",
" \n",
"plot(r, menc)"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Warning: The algorithm does not converge. Roundoff error is detected\n",
" in the extrapolation table. It is assumed that the requested tolerance\n",
" cannot be achieved, and that the returned result (if full_output = 1) is \n",
" the best which can be obtained.\n"
]
},
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 15,
"text": [
"[<matplotlib.lines.Line2D at 0x5326510>]"
]
},
{
"metadata": {},
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAXIAAAEGCAYAAAB4lx7eAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl0lFWexvFvUI4jm0rLZgLChCVkIQkGUA5gAQIDArYN\n2IqCAoJIYwuDNI22ksYl0KLDYuMCgooKOHpQZEkr2iUimxA2BcHGxAYCDFuiIUCWuvPHlSAuSSWp\nqrcq9XzOybFIwlu/vEef/LzvXSKMMQYREQlZ1ZwuQEREKkdBLiIS4hTkIiIhTkEuIhLiFOQiIiFO\nQS4iEuIqHeTDhw+nQYMGJCQklPm93377Ld27dycxMZGuXbty6NChyr69iEjYq3SQDxs2jPT0dK++\n96GHHuKee+5hx44dPPbYY0yePLmyby8iEvYqHeSdO3fmqquuuuhz+/fvp3fv3qSkpNClSxf27t0L\nwJ49e+jWrRsALpeL9957r7JvLyIS9vwyRj5q1CjmzJnDli1bePrppxkzZgwAiYmJvPPOOwAsW7aM\n77//nlOnTvmjBBGRsHGpry+Yl5fHhg0bGDRoUMnnCgoKAJgxYwZjx47llVdeoUuXLkRGRnLJJZf4\nugQRkbDi8yD3eDxceeWVbNu27Wdfa9SoUUlHnpeXxzvvvEOdOnV8XYKISFgpc2hl1qxZJCQkEB8f\nz6xZs8q8YJ06dWjWrBlvv/02AMYYdu7cCcCJEyfweDwApKWlMWLEiMrULiIilBHkX3zxBfPnz+fz\nzz9nx44drFixgv3791/0PXfccQcdO3Zk7969NG7cmIULF/LGG2/w8ssvk5SURHx8PMuXLwfgn//8\nJzExMbRq1Ypjx47xyCOP+O8nExEJExGlbWP79ttvk56ezvz58wF44oknuOyyy5g4cWLAChQRkdKV\n2pHHx8fz6aefcvLkSfLz81m5ciUHDx4MVG0iIuKFUh92xsTEMGnSJHr27EnNmjVJTk6mWjWt6hcR\nCSalDq381MMPP0yTJk0YPXp0yeeaN2/+s3FzEREpXXR0NP/61798czFThqNHjxpjjPn2229NTEyM\nyc3NvejrXlwibEyZMsXpEoKG7sUFuhcX6F5c4MvsLHMe+cCBAzlx4gTVq1dn7ty5mvctIhJkygzy\ntWvXBqIOERGpID259CGXy+V0CUFD9+IC3YsLdC/8o1wPO3/xAhERVPISIiJhx5fZqY5cRCTEKchF\nRALo5Enw9TZTCnIRkQAwBhYvhrg4qFHDt9fWGLmIiJ9lZcGYMXDwIMybBx06aIxcRCQkFBXBs89C\nSgp07gxbt9oQ9zWfHywhIiKQkQEjR8KVV8KGDdCihf/eSx25iIgPnT4NDz0EvXvDH/8Ia9b4N8RB\nQS4i4jOrV0N8PBw9Cl98AXffDRER/n9fDa2IiFTS0aMwfjxs3Agvvgg9ewb2/dWRi4hUkDGwYAEk\nJEDjxrYLD3SIgzpyEZEK2bcP7rsP8vLggw8gKcm5WtSRi4iUQ0EBPPkkdOwIt9xih1OcDHFQRy4i\n4rUNG+yUwqZN7Zzwa691uiJLQS4iUobcXHj4YVi2DGbOhEGDAjMbxVtlDq2kpaURFxdHQkICgwcP\n5ty5c4GoS0QkKCxbZvdHKSyEL7+E224LrhCHMoI8KyuLefPmkZGRwa5duyguLmbJkiWBqk1ExDHZ\n2fC738HkyfDmm/DSS3DVVU5X9ctKDfI6depQvXp18vPzKSoqIj8/n8jIyEDVJiIScB6PDe3ERLu4\nZ/t26NLF6apKV+oYed26dZkwYQJNmjTh8ssvp1evXtx0002Bqk1EJKD27YNRo+DsWfj4Yzs/PBSU\nGuT79+9n5syZZGVlccUVVzBo0CDeeOMN7rzzzou+LzU1teS1y+XSuXwiElIKC2HGDHjmGXj0URg7\nFi65xLfv4Xa7cbvdvr3oD0rdj3zp0qV8+OGHzJ8/H4BFixaxceNG/v73v1+4gPYjF5EQtmUL3Hsv\nNGwIL7xgpxYGQsD2I4+JiWHjxo2cOXMGYwxr1qwhNjbWJ28sIuKk/Hy7S+HNN9t/rl4duBD3tVKD\nPDExkaFDh5KSkkKbNm0AGDVqVEAKExHxlzVr7Pj34cN2f5S77gq+KYXloaPeRCRsnDxpu++PPoLn\nn4c+fZyrRUe9iYiUgzHw1lt2OmGtWrYLdzLEfU1L9EWkSjt4EP7wB/j6a3jnHbjhBqcr8j115CJS\nJXk8dvgkORnatoVt26pmiIM6chGpgr76yu5SWFQEbrfdK6UqU0cuIlXG+b3CO3Wym1utW1f1QxzU\nkYtIFfH55zBiBERFQUYGNGnidEWBo45cRELa+YU9/frBn/8MK1eGV4iDglxEQpjbDW3a2C1nd+2C\nwYNDe2FPRWloRURCTm4uTJoEK1bYmSn9+jldkbPUkYtISFm50i7s8XjsiT3hHuKgjlxEQsTx4zBu\nnD0A+dVXoVs3pysKHurIRSSoGQNLl9pNrurXh507FeI/pY5cRIJWdjaMGWOX1y9bBtdf73RFwUkd\nuYgEHWPg5ZchKcnOSsnIUIiXRh25iASVzEy7vD4nx+4b/sNRCFIKdeQiEhSKi2HWLGjXDnr1go0b\nFeLeKjPI9+7dS3JycsnHFVdcwezZswNRm4iEiT17oHNnu83s+vUwcSJcqvECr5XrhCCPx0NkZCSb\nN2+mcePG9gI6IUhEKqiwEP72N5g5E6ZOhfvug2phMk7gy+ws1++8NWvWEB0dXRLiIiIVlZEBw4dD\no0awdWv47Y/iS+X63bdkyRIGDx7sr1pEJAycOwd/+Qv07g3//d+wapVCvLK87sgLCgp4//33mT59\n+s++lpqaWvLa5XLhcrl8UZuIVDGbN9suvEUL2L7dduPhwu1243a7/XJtr8fI33vvPZ5//nnS09Mv\nvoDGyEWkDGfOwJQp8Nprdjz8978Pz10Kf8yRMfLFixdzxx13+ORNRSR8rF9vu/A2bezy+vr1na6o\n6vGqIz99+jTXXnstmZmZ1K5d++ILqCMXkV+Qnw+PPGL3SZkzBwYMcLqi4OLL7PTqYWfNmjU5fvz4\nz0JcROSXfPKJ7cCPHbMHPijE/UtT7kXEZ/Ly7HFr774Lc+dC//5OVxQewmTqvYj420cf2a1mT5+2\nXbhCPHDUkYtIpXz3nV1Sv3o1vPiinR8ugaWOXEQqLD3dHrtmjO3CFeLOUEcuIuV26pRdlfnPf8KC\nBXDTTU5XFN7UkYtIubz/vh0Lr1HDduEKceepIxcRr5w6BQ8+CJ99Bq+/DtqJI3ioIxeRMq1YYbvw\nK66wqzMV4sFFHbmI/KpTp2DcOFi3Tl14MFNHLiK/aOVK24XXrg07dijEg5k6chG5SE6O7cLXroVF\ni6BrV6crkrKoIxeREqtW2S68Zk07Fq4QDw3qyEWEnBwYPx7cbnj1VejWzemKpDzUkYuEudWrbRd+\n+eV2XrhCPPSoIxcJUzk5F1ZnqgsPberIRcJQerrtwi+7zI6FK8RDW5lBnpOTw8CBA2ndujWxsbFs\n3LgxEHWJiB/k5sKIETB6NLzyCjz/vJ1eKKGtzCB/8MEH6dOnD3v27GHnzp20bt06EHWJiI998IHt\nwqtXt2Ph3bs7XZH4Sqlndubm5pKcnMw333zz6xfQmZ0iQe377+1+4atWwcsvQ48eTlckEMAzOzMz\nM6lXrx7Dhg2jbdu2jBw5kvz8fJ+8sYj4n9sNiYlQUGC7cIV41VTqrJWioiIyMjJ47rnnaNeuHePG\njWPatGlMnTr1ou9LTU0tee1yuXBpLa+Io/LzYfJkePtte2pP375OVyRutxu32+2Xa5c6tHLkyBFu\nuOEGMjMzAVi3bh3Tpk1jxYoVFy6goRWRoLJ+PdxzD7RrB3PmQN26TlckvyRgQysNGzakcePG7Nu3\nD4A1a9YQFxfnkzcWEd86exYmTYIBA2DaNHjjDYV4uChzQdCcOXO48847KSgoIDo6moULFwaiLhEp\nh61bYehQaN3a7lRYv77TFUkglTq04tUFNLQi4piCAnjiCTsOPnMm3H47REQ4XZV4w5fZqSX6IiFq\n5064+26IioLt26FRI6crEqdoib5IiCkqgieftAt6/vhHWL5cIR7u1JGLhJA9e2wXfuWVkJEBjRs7\nXZEEA3XkIiHA44Fnn4UuXWD4cPjHPxTicoE6cpEgl5lp54V7PLBxI0RHO12RBBt15CJByhiYNw/a\nt4d+/exye4W4/BJ15CJBKDsbRo6EI0dsgGsdnpRGHblIkFmyBJKTISXFDqUoxKUs6shFgsSJEzBm\njJ0fvmKF3StFxBvqyEWCwMqV0KYNREbaaYUKcSkPdeQiDvruO3sA8kcfwZtvwo03Ol2RhCJ15CIO\nOX/oQ0SEHU5RiEtFqSMXCbAzZ+Dhh+Gtt+Cll+Dmm52uSEKdOnKRANqyBdq2tdMLd+5UiItvqCMX\nCYDCQkhLg+eeg9mz7XazIr7iVZA3bdqUOnXqcMkll1C9enU2b97s77pEqoy9e2HIELjqKti2zc5M\nEfElr4I8IiICt9tNXZ0bJeI1jwfmzoXUVJg6Fe6/X4c+iH94PbSiU4BEvHfwoN2lMDfXHobcsqXT\nFUlV5tXDzoiICG666SZSUlKYN2+ev2sSCVnG2PngbdtC587w2WcKcfE/rzryzz77jEaNGnHs2DF6\n9OhBTEwMnTt39ndtIiHl/BL7Xbtg9Wq47jqnK5Jw4VWQN/rhHKl69epx6623snnz5ouCPDU1teS1\ny+XC5XL5tEiRYLd6Ndx7L9x2G7zyClx+udMVSbBxu9243W6/XDvClDH4nZ+fT3FxMbVr1+b06dP0\n7NmTKVOm0LNnT3sBH54ELRJq8vJg4kRYtQoWLoRu3ZyuSEKFL7OzzI786NGj3HrrrQAUFRVx5513\nloS4SDjbsAGGDoWOHWHHDnuOpogTyuzIy7yAOnIJM4WFdjrhSy/Z6YUDBjhdkYSigHbkInLBV1/B\nXXdB/fqwfTv88PhIxFHaa0XEC8bY5fWdOsGIEXb/cIW4BAt15CJlyM6GYcPg1Ckt7pHgpI5cpBT/\n+7/2/MwbbtDiHgle6shFfkFuLowdC5s2wfLl0KGD0xWJ/Dp15CI/4Xbb8zNr17a7FSrEJdipIxf5\nwblz8Mgjdq+U+fOhTx+nKxLxjoJcBHtaz113QfPm9vXVVztdkYj3NLQiYc3jgWefhe7d7Wn277yj\nEJfQo45cwtaBA3DPPXZIZdMm+M//dLoikYpRRy5hackSu81s9+7wyScKcQlt6sglrOTk2GmFW7bY\nHQtTUpyuSKTy1JFL2PjkE0hKgiuugIwMhbhUHerIpco7dw4eewwWLdK0QqmaFORSpX35pZ1WeO21\nds/wevWcrkjE9zS0IlWSxwOzZ8ONN8If/gDLlinEperyqiMvLi4mJSWFqKgo3n//fX/XJFIphw/b\n3QpzcmDjRrvIR6Qq86ojnzVrFrGxsURERPi7HpFKefddu1vh9dfDunUKcQkPZQb5wYMHWbVqFffe\ne6+OdJOglZcHI0fChAl2GCU1FS7VEyAJE2UG+fjx43n66aepVk3D6RKcNm2yXXhxsT1+7YYbnK5I\nJLBKTecVK1ZQv359kpOT1Y1L0Ckqgscfh/79IS0NFiywW8+KhJtS/+dz/fr1LF++nFWrVnH27Fm+\n++47hg4dymuvvXbR96Wmppa8drlcuFwuf9QqUuKbb2DIELj8cru4JzLS6YpESud2u3G73X65doTx\nstX+5JNPmDFjxs9mrURERKhbl4AxBl57DR56CCZPhnHjQKN+Eop8mZ3lehykWSvipJMnYfRo2L0b\n1qyBxESnKxIJDl535L96AXXkEgAffwx33w0DBsC0afAf/+F0RSKV41hHLhJo587Bo4/CG2/Yh5m9\nejldkUjwUZBL0NqzBwYPtvukbN+uJfYiv0aPiSToGANz50KXLnD//donRaQs6sglqBw9CiNGwJEj\ndol9q1ZOVyQS/NSRS9BYtcoe/NCmDaxfrxAX8ZY6cnHcmTMwcSK8/749S/PGG52uSCS0qCMXR23f\nbo9cO37cvlaIi5Sfglwc4fHAM89Ajx7w5z/D4sVw1VVOVyUSmjS0IgF36JBd3HPmDGzeDM2aOV2R\nSGhTRy4BtWwZtG1rpxZ+8olCXMQX1JFLQOTlwfjxdqn9e+/ZE3xExDfUkYvfff657cKLiuwDTYW4\niG+pIxe/KS6Gv/0N/ud/4Lnn4LbbnK5IpGpSkItf/Pvf9uCHatVg61Zo3NjpikSqLg2tiM8tXWrn\nhvfpY/cNV4iL+Jc6cvGZ776DBx6AjRth9Wq47jqnKxIJD2V25GfPnqVDhw4kJSURGxvL5MmTA1GX\nhJj16+0+KefP0FSIiwSOVycE5efnU6NGDYqKiujUqRMzZsygU6dO9gI6ISisFRXBE0/ACy/Aiy/C\nLbc4XZFIaAj4CUE1atQAoKCggOLiYurWreuTN5fQ9s03cNddUKuW7cKvucbpikTCk1cPOz0eD0lJ\nSTRo0ICuXbsSGxvr77okiJ0/yb5DBxg0CNLTFeIiTvKqI69WrRrbt28nNzeXXr164Xa7cblcJV9P\nTU0tee1yuS76mlQtp07ZU3t27dJJ9iLl4Xa7cbvdfrm2V2PkP/b4449z+eWX89BDD9kLaIw8bLjd\ndrOrW26B6dPtg00RqRhfZmeZQyvHjx8nJycHgDNnzvDhhx+SnJzskzeX0FBQAJMn24OQX3gBZs9W\niIsEkzKHVg4fPszdd9+Nx+PB4/EwZMgQunfvHojaJAjs3WsD/Jpr7D4p9es7XZGI/FS5h1Z+dgEN\nrVRJxsBLL8Ff/gJTp8Lo0RAR4XRVIlVHwKcfSng5dgzuvRcOHIC1a6F1a6crEpHSaK8VuUh6ul2h\n2aqVXWqvEBcJfurIBbDHrv3pT/bQh9dfh65dna5IRLyljlxKTrI/dgx27FCIi4QaBXkY83hgxgyd\nZC8S6jS0EqYOHLCLewoL7VFsTZs6XZGIVJQ68jD01lt2m9nu3e1qTYW4SGhTRx5GcnNh7FjYvBlW\nroR27ZyuSER8QR15mFi71m5wdX7LWYW4SNWhjryKKyiAxx6z287Omwc33+x0RSLiawryKmz3brjz\nTmjSRPukiFRlGlqpgjwemDMHbrwRxoyBd99ViItUZerIq5hDh2DYMHui/fr10KKF0xWJiL+pI69C\nli6Ftm2hc2dYt04hLhIu1JFXAadO2WmFW7fCihWakSISbtSRh7iPPrLTCuvW1bRCkXBVZpAfOHCA\nrl27EhcXR3x8PLNnzw5EXVKGM2dg/Hi7zH7ePPtws0YNp6sSESeUeULQkSNHOHLkCElJSeTl5XHd\nddfx7rvv0vqHjap1QlDgbdkCQ4bYTvzvf4ff/MbpikSkvAJ6+HLDhg1JSkoCoFatWrRu3Zrs7Gyf\nvLmUT2Eh/PWvdlHPlCmwZIlCXETK+bAzKyuLbdu20aFDB3/VI79izx4YOhSuvhq2bbOHIYuIQDmC\nPC8vj4EDBzJr1ixq1ap10ddSU1NLXrtcLlwul6/qC3vnF/c88YT9GDVKhyCLhCK3243b7fbLtcsc\nIwcoLCykb9++9O7dm3Hjxl18AY2R+01mJgwfbvdLefVVaN7c6YpExFcCOkZujGHEiBHExsb+LMTF\nP4yBl16C9u3tePjatQpxEfl1ZXbk69ato0uXLrRp04aIH/6fPi0tjf/6r/+yF1BH7lMHD8K998Lx\n47YLj4tzuiIR8QdfZqdXQyuBKiacGWNPr58wwa7SnDwZqld3uioR8RdfZqeW6AeBw4dh9Gj45htI\nT7f7pYiIeEtL9B1kDLzxBiQlQUKCXeijEBeR8lJH7pAjR2wX/q9/2fMzU1KcrkhEQpU68gAzBt58\n0y6vj4uzOxYqxEWkMtSRB1B2Ntx/P+zfr+1mRcR31JEHgDHwyit2LDwpyXbhCnER8RV15H7273/D\nfffB0aPwwQc2yEVEfEkduZ94PDB3Llx3nT16bdMmhbiI+Ic6cj/Yu9euzvR47PL6H7ZuFxHxC3Xk\nPlRYCE89BZ06we9/D59+qhAXEf9TR+4jn38OI0dCo0Z2Yc+11zpdkYiEC3XklfT99zBuHPTvD3/6\nE6xapRAXkcBSkFfCihUQHw+5ufDFFzB4sA59EJHA09BKBWRnw4MPwvbtsHAhdOvmdEUiEs7UkZdD\ncbE9di0xEWJiYOdOhbiIOE8duZcyMuzCnho1NKVQRIJLmR358OHDadCgAQkJCYGoJ+jk5MADD0Cf\nPjBmDLjdCnERCS5lBvmwYcNIT08PRC1B5fyJPa1b28OPv/wShg3Tw0wRCT5lDq107tyZrKysAJQS\nPL74wh63lpsLy5bB9dc7XZGIyK/Tw84fyc2F8eOha1cYNMgu7FGIi0iw88nDztTU1JLXLpcLl8vl\ni8sGjMcDixbZA4/79IHdu6FePaerEpGqxO1243a7/XLtCOPFMc5ZWVn069ePXbt2/fwCPjwJ2gmf\nf24fZhpjpxa2b+90RSISDnyZnWE7tHLkCAwfbpfWjx4NGzYoxEUkNJUZ5HfccQcdO3Zk3759NG7c\nmIULFwaiLr85examTbNL66++Gr76Cu65B6qF7a80EQl1Xg2tlHqBEBlaMQbefhsmTYI2bWDGDGje\n3OmqRCRc+TI7w2Jl5+bNdjbK6dMwf76W1YtI1VKlBxQyM+H22+HWW2HECHvosUJcRKqaKhnkJ07A\nhAn2pPq4ONi3zz7YvOQSpysTEfG9KhXk+fn2QWZMDJw5Y1doPvoo1KzpdGUiIv5TJcbICwthwQKY\nOhU6doTPPoOWLZ2uSkQkMEI6yD0eWLIEHnsMmjaFd9+1wykiIuEkJIPcGLuZ1ZQpUKsWvPSSHmKK\nSPgKqSA3xp6TmZpqX6elwc03a2tZEQlvIRHk5wN86lS7N3hqKvz2twpwEREI8iD3eOwQyhNP2DB/\n9FE7J1zL6UVELgjKIC8qsg8x09LsGZl//Sv07asAFxH5JUEV5Pn5sHCh3Qfl2mth5ky46SYNoYiI\nlCYogvz//g/mzoXnn7fzwN98E264wemqRERCg6ODFbt3w6hR0KoVZGfD2rV2TFwhLiLivYB35MXF\nsGqVPY1n5064/37Yuxfq1w90JSIiVUOZHXl6ejoxMTG0aNGC6dOnV/iNjh61+6A0bw6PPw5DhsC3\n39pFPQpxEZGKKzXIi4uLGTt2LOnp6ezevZvFixezZ88ery9eXAz/+AfcdpsdPtm3D956y+4PPmQI\nXHZZpesPKv46WDUU6V5coHtxge6Ff5Qa5Js3b6Z58+Y0bdqU6tWrc/vtt/Pee++VedFdu+DPf7Yz\nTx55BLp2td33ggVVey8U/Ut6ge7FBboXF+he+EepY+SHDh2icePGJX+Oiopi06ZNP/s+jwcyMmD5\ncnuc2unT9kCH9HR7NqaIiPhPqUEe4eUE7vr1oV49u2jn5ZehQwct3hERCRhTig0bNphevXqV/Pmp\np54y06ZNu+h7oqOjDaAPfehDH/oox0d0dHRp8VsuEaaUY5yLiopo1aoVH330Eddccw3t27dn8eLF\ntG7d+tf+ioiIBFipQyuXXnopzz33HL169aK4uJgRI0YoxEVEgkypHbmIiAS/Sj2S9NVioVBx4MAB\nunbtSlxcHPHx8cyePRuAkydP0qNHD1q2bEnPnj3Jyckp+TtpaWm0aNGCmJgYPvjgA6dK94vi4mKS\nk5Pp168fEL73ASAnJ4eBAwfSunVrYmNj2bRpU1jej7S0NOLi4khISGDw4MGcO3cubO7D8OHDadCg\nAQkJCSWfq8jPvnXrVhISEmjRogUPPvigd29e0cH1oqIiEx0dbTIzM01BQYFJTEw0u3fv9tXYfVA6\nfPiw2bZtmzHGmO+//960bNnS7N6920ycONFMnz7dGGPMtGnTzKRJk4wxxnz55ZcmMTHRFBQUmMzM\nTBMdHW2Ki4sdq9/XnnnmGTN48GDTr18/Y4wJ2/tgjDFDhw41L7/8sjHGmMLCQpOTkxN29yMzM9M0\na9bMnD171hhjzG233WZeeeWVsLkPa9euNRkZGSY+Pr7kc+X52T0ejzHGmHbt2plNmzYZY4zp3bu3\nWb16dZnvXeEgX79+/UUzWtLS0kxaWlpFLxeSbrnlFvPhhx+aVq1amSNHjhhjbNi3atXKGPPzWT69\nevUyGzZscKRWXztw4IDp3r27+fjjj03fvn2NMSYs74MxxuTk5JhmzZr97PPhdj9OnDhhWrZsaU6e\nPGkKCwtN3759zQcffBBW9yEzM/OiIC/vz56dnW1iYmJKPr948WJz3333lfm+FR5a+aXFQocOHaro\n5UJOVlYW27Zto0OHDhw9epQGDRoA0KBBA44ePQpAdnY2UVFRJX+nKt2j8ePH8/TTT1PtRwsGwvE+\nAGRmZlKvXj2GDRtG27ZtGTlyJKdPnw67+1G3bl0mTJhAkyZNuOaaa7jyyivp0aNH2N2HHyvvz/7T\nz0dGRnp1Tyoc5N4uFqqK8vLyGDBgALNmzaJ27doXfS0iIqLUe1MV7tuKFSuoX78+ycnJmF95Vh4O\n9+G8oqIiMjIyGDNmDBkZGdSsWZNp06Zd9D3hcD/279/PzJkzycrKIjs7m7y8PF5//fWLvicc7sOv\nKetnr4wKB3lkZCQHDhwo+fOBAwcu+k1SVRUWFjJgwACGDBnCb3/7W8D+pj1y5AgAhw8fpv4P2zn+\n9B4dPHiQyMjIwBftY+vXr2f58uU0a9aMO+64g48//pghQ4aE3X04LyoqiqioKNr9sJHQwIEDycjI\noGHDhmF1P7Zs2ULHjh35zW9+w6WXXsrvfvc7NmzYEHb34cfK899EVFQUkZGRHDx48KLPe3NPKhzk\nKSkpfP3112RlZVFQUMDSpUvp379/RS8XEowxjBgxgtjYWMaNG1fy+f79+/Pqq68C8Oqrr5YEfP/+\n/VmyZAkFBQVkZmby9ddf0759e0dq96WnnnqKAwcOkJmZyZIlS+jWrRuLFi0Ku/twXsOGDWncuDH7\n9u0DYM2aNcTFxdGvX7+wuh8xMTFs3LiRM2fOYIxhzZo1xMbGht19+LHy/jfRsGFD6tSpw6ZNmzDG\nsGjRopJJlD2mAAAA70lEQVS/U6rKDOyvWrXKtGzZ0kRHR5unnnqqMpcKCZ9++qmJiIgwiYmJJikp\nySQlJZnVq1ebEydOmO7du5sWLVqYHj16mFOnTpX8nSeffNJER0ebVq1amfT0dAer9w+3210yayWc\n78P27dtNSkqKadOmjbn11ltNTk5OWN6P6dOnm9jYWBMfH2+GDh1qCgoKwuY+3H777aZRo0amevXq\nJioqyixYsKBCP/uWLVtMfHy8iY6ONg888IBX760FQSIiIU57FIqIhDgFuYhIiFOQi4iEOAW5iEiI\nU5CLiIQ4BbmISIhTkIuIhDgFuYhIiPt/3kz9PiXPDNMAAAAASUVORK5CYII=\n",
"text": [
"<matplotlib.figure.Figure at 0x56aed90>"
]
}
],
"prompt_number": 15
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"loglog()\n",
"plot(r, menc)\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 16,
"text": [
"[<matplotlib.lines.Line2D at 0x4ef2a50>]"
]
},
{
"metadata": {},
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAYEAAAEHCAYAAABIsPrhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAG1VJREFUeJzt3X10VPWdx/FPBOyWdrfCVrOVyRFKkAQJtFsOrrApsbSL\nxcoqCAy6yiY8CDQotRSra02CVk3rHndrVEoRVJSATzU+jiu4g8W1REWeMqGAm6whRxflMYKAGe7+\n8dukEyVkHm7m3jv3/Ton53Qm8/C1l9zv/X3v9/f7ZVmWZQkA4EtnOB0AAMA5JAEA8DGSAAD4GEkA\nAHyMJAAAPkYSAAAfIwkAgI+RBADAx2xPAg0NDZo5c6YmT57c/tyRI0c0ffp0zZ49W6tWrbL7KwEA\nSbI9CQwYMEDLli3r8NwzzzyjKVOmaOnSpXruuefs/koAQJLiSgIlJSXKzs5WQUFBh+dDoZDy8vI0\naNAgVVZWdvr+5uZm5eTkSJJ69OiRQrgAADvFlQSKi4sVCoU6PBeNRlVaWqpQKKRIJKLq6mrV19ef\n8v2BQEBNTU2SpJMnT6YYMgDALnElgcLCQvXp06fDc7W1tcrNzVX//v3Vq1cvBYNB1dTUaP/+/Zoz\nZ442b97cPjqYOHGinn76ac2bN08TJkyw/78CAJCUnsm+MbbEI5mr/Y0bN6pv375asmRJh9f27t1b\ny5cvP+3n5ebm6r333ks2HADwpYEDB2r37t1Jvz/pG8NZWVlJf+mpvPfee7Isy/GfsrIyV3xevO+L\n53Vdvaaz3yfyvN3/v7nh+HX3sbPj+CX6u0SPtVePnVuOXzr+9lK9eE46CfTr16+9zi9JTU1NCgQC\nKQVTXl6ucDic0mekqqioyBWfF+/74nldV6/p7PeJPu8GdsbW3ccu3tee7jWJ/s4vxy6Vz7Pz+HXn\n3144HNbCheWn/fy4WHFqaGiwhg4d2v74s88+s775zW9aDQ0N1vHjx63hw4dbkUgk3o/7ggRCgQuV\nlZU5HQJSwPHzluZmy7rhBsvq0yf1c2dcI4Fp06Zp1KhR2rlzp3JycrRixQr17NlTVVVVGjdunIYM\nGaKpU6cqPz8/9awET3LzFSa6xvHzhuZmaf58aehQ6YwzpLq61D8zy7IsV2wvmZWVpbKyMhUVFfEP\nEgBi7Nkj3X23tGqVNGOGtHChVF8fVjgcVkVFhVI5jbsqCbgkFABwheZm6a67/nzy/9nPpHPO6fia\nVM+dLCAHAC7zwQfSDTdIBQXSX/yFtGOH9OtffzEB2MFVScAN3UEA4JSPPjKlngsuMDX/SES6555T\nn/zD4bDKy8tT/k7KQQDgsAMHzMl+yRIpGJRuuUXq1y++91IOAgCP+uQT6Ze/lM4/X/rwQ2nTJun+\n++NPAHYgCQBAmh0/Lv37v0uDBpk2zzfekB56SDrvvPTHkvTaQd2hvLycFlEAGSsalVaulMrLzU3f\nV16Rhg1L7rPC4bAt91C5JwAA3cyypOeeM7X+vn1Nz//o0fZ8dqrnTleNBAAg02zYIN10k3T4sPSr\nX0njx0s2r7+ZEpIAAHSD+nrp5z+XNm+Wbr9duvpqyY0bK3JjGABs9MEH0nXXSd/9rlRYKP3pT9K1\n17ozAUguSwJMFgPgVUeOSIsXm8Xd/vIvzcl/4UIz47c7MFkMAFwgGpUeeUT6xS+kMWNM3/+AAen7\nfm4MA4BDXntNuvFG6atflX7/e2nkSKcjShxJAAAStHOnKfXU1ZmOn4kT3dXxkwhX3RMAADc7eNBc\n+Y8ebW76RiLSpEneTQCSy5IAN4YBuFE0ahZ3GzzY3ACuqzNr+3/pS87FxI1hAEiDcNis7d+nj1nv\nZ/hwpyPqiBvDANAN/ud/TN3/rbfMMs9eL/t0xlXlIABw2qefmn7/v/1b0/NfXy9deWVmJgCJkQAA\nSPrzIm8LFkjf+Y5Z29+JpZ3TjSQAwPd27ZKuv15qbJR+9zvp+993OqL0oRwEwLeOHJH+5V+kiy6S\nxo6VtmzxVwKQXDYSYFMZAOlgWVJNjSn9jBplTv7p3NLRDmwqAwBJeO89U/r57/82+/l+73tOR5Qa\nNpoHgDgcP27W9b/wQrPM85Yt3k8AdnBVOQgAusPatdK8edIFF0jvvOOPrp94kQQAZKwPPjBr/fzx\nj9J990k/+pHTEbkP5SAAGScaNfX+YcPM2v51dSSAzjASAJBRNm0y2zt++cvS+vXSkCFOR+RujAQA\nZISWFtPy+cMfmvo/CSA+JAEAnvfss+aEf/iwKf0UF2fuWj92oxwEwLPef1+aP9/s9PXYY2aPXyTG\nVSMBNpUBEI/WVunee81KnyNGSJs3+y8BsKkMAF965x1p9mzprLOkBx+Uzj/f6YicxYxhAL7Q0iL9\n5CfSpZeanb7WriUB2IEkAMD1nnvOzPY9eFDavl269lpu/NqFG8MAXKu52Sz2tm2b9Mgj0sUXOx1R\n5mEkAMB12mb8Dh9uRgBbt5IAugsjAQCusnWrufHbq5f0+utM+OpujAQAuMLRo9LPf2529poxgxm/\n6UISAOC4//gPqaDA7PG7das0a5Z0BmentKAcBMAxe/eats//+i9zD2D8eKcj8h9yLYC0syzpoYek\noUPN3r7bt5MAnMJIAEBa7dhhlnr+9FNTBvrWt5yOyN/SMhKIRCKaOnWq5s2bp6effjodXwnAZY4d\nk8rLpcJCafJk6c03SQBukJYkEAqFNH/+fD3wwAN69NFH0/GVAFzkP//T9Pxv3Sq9+65UWir16OF0\nVJBSSAIlJSXKzs5WQUFBh+dDoZDy8vI0aNAgVVZWSpKuueYarV69WosWLdK+fftSixiAZ3z8sfTP\n/yxNny796lfSM89IgYDTUSFW0kmguLhYoVCow3PRaFSlpaUKhUKKRCKqrq5WfX29zj77bFVVVemu\nu+7S17/+9ZSDBuBulmWWebjgArPaZ12d9I//6HRUOJWkbwwXFhaqsbGxw3O1tbXKzc1V//79JUnB\nYFA1NTXq3bu37rzzTh05ckSLFi1KJV4ALvenP0lz5phdvl56SfrOd5yOCKdja3dQc3OzcnJy2h8H\nAgFt3LhR5513nn772992+f7YDRKKiopUVFRkZ3gAutGxY9Ldd0tVVdKtt5q6f0/6D20XDodt3XzL\n1kOUleLarnbskgMg/V57TZo715R/Nm+m7t+dPn+BXFFRkdLn2ZoE+vXrp6ampvbHTU1NCiTwr6G8\nvJwRAOAhH30k/fSnUjgs3Xcfdf90smtEkNL2ko2Njbrsssu0bds2SVJra6sGDx6sdevW6dxzz9XI\nkSNVXV2t/Pz8rgNhe0nAM06elJYvl265RbrmGqmiQvrqV52Oyp9SPXcmPRKYNm2a1q9fr3379ikn\nJ0eLFy9WcXGxqqqqNG7cOEWjUc2YMSOuBADAO+rqzI3fEyeY8ZsJXLXRfFlZGeUgwKWOHpVuv11a\ntkxavNis+c+EL+e0lYMqKipSGgm4Kgm4JBQAn/PSS6bb58ILpXvvlf7mb5yOCG0cKwcByHx79kgL\nFkhbtkhLlkj/8A9ORwS7uWop6fLyclv7XwEkp7VV+rd/k779bdP2uW0bCcBtwuGwLW31lIMAdLBx\no7nx27ev9MAD0uDBTkeE00n13OmqkQAA5xw4YCZ8XX65tHChtHYtCcAPSAKAz1mW9NhjZlP3rCwp\nEpGuvtr8b2Q+V90YZsYwkF47dkjz5kkHD0o1NdLIkU5HhHi5YsawnbgnAKTP0aPSL38pLV0q/eIX\nJhGw2Js30SIKICEvvijNn2+u+rdskc491+mI4CSSAOATTU3SDTeYdk96/tHGVTeGmScA2O+zz6R7\n7jE9/9/6Fj3/mYJ5AgC69Ic/mLbPQMBs9pKb63REsBv3BAB8wd690qJF0rp1ZubvxIm0fOLUXFUO\nApCaaNTU+4cOlf76r03P/6RJJAB0jpEAkCHeeceUfs4804wACgqcjghewEgA8LiDB80yz5deapLA\n66+TABA/VyUBuoOA+FmWtHKllJ9vVv2MRKTiYukMV/1Vo7vQHQT4WF2dmeX7ySdmpc8LL3Q6IjiF\nVUQBH/nkE9P1U1QkTZki1daSAJAakgDgAZYlPfWUWenzww+l7dulH/+YPX6ROrqDAJfbtcvc+G1u\nNks+f/e7TkeETMJIAHCpTz+VbrtNuugi6Qc/kN59lwQA+zESAFyobaXPESOkzZvNsg9Ad3BVEmBT\nGfhdY6O0YIFp92SlT5wOm8oAGeT4cbPS5733Sj/5idnj90tfcjoqeAELyAEe9+qr5sbv4MHSW29J\nAwY4HRH8hCQAOGTPHunGG6W335Z+8xvpRz9yOiL4Ed1BQJqdOCH9+tdmg5e8PDP7lwQApzASANIo\nHDbLPZx3nvTHP7LJC5xHEgDS4IMPzM3eDRvMJi+XX84a/3AHykFAN2ptNSf9YcPM1X8kIl1xBQkA\n7sFIAOgmGzaY9X3OPtvs9ZuX53REwBeRBACb/e//mpU+X3tN+td/lSZP5sof7uWqchCbysDLWlul\nqiqzv+8555jSz5QpJAB0DzaVAVzkzTdN189ZZ5lEcMEFTkcEv2DGMOCgjz6Sbr5Zevll0/s/bRpX\n/vAWV5WDAK+IRqUHHzRX/H/1V1J9vXTVVSQAeA8jASBBtbWm9NO7t7RunVRQ4HREQPIYCQBx+vhj\nafZsM9FrwQJp/XoSALyPJAB0IRqVli41pZ8vf9l0/fzTP1H6QWagHAScxttvm9JPr17SK6+YRd+A\nTMJIADiF/fuluXOlyy4zs37/8AcSADITSQCIcfKktGyZlJ8v9expSj/Tp0tn8JeCDEU5CPh/mzaZ\n0k9WlhQKSd/+ttMRAd2P6xv43v795uQ/frw0a5b0xhskAPhHWpLAnj17NHHiRM2YMUOVlZXp+Eqg\nSydPSitWSEOGmMeRiDRjBqUf+EtaykHbtm3TpEmTdPXVVysYDKbjK4HTevddc8M3GpVeeEEaMcLp\niABnJH3NU1JSouzsbBV8brZMKBRSXl6eBg0a1H7VP2rUKC1dulRjx47VJZdcklrEQAoOHJBKS6VL\nLjFX/W++SQKAvyWdBIqLixUKhTo8F41GVVpaqlAopEgkourqatXX12vFihW64447tG7dOr344osp\nBw0k6uRJ6eGHTeknGqX0A7RJuhxUWFioxsbGDs/V1tYqNzdX/fv3lyQFg0HV1NRo/PjxWrx4sVat\nWqUBAwZ0+pmxa2MXFRWpqKgo2fCAdps3m9LPZ59Jzz/PlT+8LRwO27rviq33BJqbm5WTk9P+OBAI\naOPGjRo2bJieeuqpLt9vxwYJQJuDB6XbbpPWrJFuv12aOZMrf3jf5y+QKyoqUvo8W/8kslhMBS5g\nWdKjj5oJX8ePm9LP7NkkAOBUbB0J9OvXT01NTe2Pm5qaFAgE4n5/eXk5ZSCkZMsWc+P32DGppkYa\nOdLpiIDuYVdZKKXtJRsbG3XZZZdp27ZtkqTW1lYNHjxY69at07nnnquRI0equrpa+fn5XQfC9pJI\nwaFDpvSzerUp/cyYIfXo4XRUQPdL9dyZ9AB52rRpGjVqlHbu3KmcnBytWLFCPXv2VFVVlcaNG6ch\nQ4Zo6tSpcSUAIFmxpZ9PP5Xq6kzphwQAxMdVG82XlZVRDkLctm41XT/Hjkn330/pB/7SVg6qqKhI\naSTgqiTgklDgcocOSWVlUnW1tHix6frhyh9+5Vg5CEg3y5JWrjSln6NHTennuutIAEAqXLWUNN1B\n6My2bab0c/So9PvfSxde6HREgLNc0R1kJ8pBOJVDh6Tycunxx03pZ9YsrvyBWJSDkJEsS3rsMbPW\nT0uLKf3MmUMCAOzmqnIQIJnST2mp9Mkn0tNPS3/3d05HBGQuV40EysvLbV0YCd5y+LB0443S2LFS\nMCjV1pIAgM6Ew2Fb1lvjngAcZ1nSqlXSokXSD38o3XWXdPbZTkcFeEOq507KQXDU9u2m66elhdIP\n4ARXlYPgH4cPSz/9qfS975nSz1tvkQAAJ7gqCXBPIPNZlmn3zM836/3X1Ulz59L1AySKewLwnNjS\nz/33Sxdd5HREgPcxTwCu19b1c/HF0pQppvRDAgDcgSSAbtPW9RNb+vnxjyn9AG5CdxC6xfbtZsLX\noUPSk09Ko0Y5HRGAU2EkAFvFdv1MmSK9/TYJAHAzVyUBuoO8K7b0c+CAGQnMm0fpB+gudAfBNdpq\n/YcOSQ88wE1fIJ3oDoJj2ko/RUXSlVea0g8JAPAWkgASZllma8e20k9dnbkJTOkH8B66g5CQttLP\nwYN0/QCZgJEA4tLSIi1c2LH0QwIAvI8kgNNqK/3k5Un79v259NOTMSSQEVz1p8xG8+4S2/VD6Qdw\nFzaaR7c5fFiqqJAefVQqK2OVT8DNaBGFbWK7fvbvp+sH8ANXlYPgnLYTPl0/gL8wEvA5un4AfyMJ\n+FRs18/HH7PMM+BXlIN8qK30c+CA9MQT0ujRTkcEwCmMBHwktvQzaZIp/ZAAAH8jCfiAZUmrV5uu\nn337/rzhCxO+ALjqNMBkMfvV1Unz55uWzyee4KYvkCmYLIbTamkxE74eecRM+Jozhyt/IBMxWQwd\nxJZ+Pv6Y0g+A0+PUkEEiEXPC379fWrOGm74AusZIIAO0df2MGSNNnEjXD4D4kQQ8LHatH7p+ACSD\n04VHxU74ovQDIFmMBDyGCV8A7EQS8Ai6fgB0B04hHhA74YvSDwA7MRJwsZYW6Wc/M6Ufun4AdAeS\ngAtR+gGQLmk5rWzYsEGPP/64WltbFYlE9MYbb6Tjaz2J0g+AdErr2kE1NTXau3evZs2a9cVAfL52\nEGv9AEiGY2sHlZSUKDs7WwUFBR2eD4VCysvL06BBg1RZWdnhd6tWrdJVV12V7FdmJEo/AJyUdBIo\nLi5WKBTq8Fw0GlVpaalCoZAikYiqq6tVX18vSXr//ff1ta99TV/5yldSiziD1NVJY8dKd99tSj8P\nPyxlZzsdFQA/SToJFBYWqk+fPh2eq62tVW5urvr3769evXopGAyqpqZGkrR8+XKVlJSkFm2GoOsH\ngFvYWnRobm5WTk5O++NAIKCNGzdKMhvGdCX2NZm4uYxlmSv+hQulH/zAlH648geQCLs2k2ljaxLI\nyspK6f3xJAqvYplnAHb4/AVyRUVFSp9n6zyBfv36qampqf1xU1OTAoGAnV/hOW2lnzFjpCuuoPQD\nwF1sTQIjRozQrl271NjYqBMnTmjNmjWaMGFC3O8vLy+3dZjjpNiun48+MqWf+fPp+gFgj3A4bEv1\nJOl5AtOmTdP69eu1b98+nXPOOVq8eLGKi4v18ssva8GCBYpGo5oxY4Zuvvnm+ALJoHkCsRO+7r+f\nK38A3SfVcycbzduICV8A0i2jNpr3ajmIHb4ApJvj5SC7eXUkELvDF6UfAOmWUSMBL2HCF4BM4Kok\n4IVy0KnW+qHrB0C6UQ5yAF0/ANyGclAaxG7uTukHQCYhCZwGyzwDyHSuOp2Vl5e7ZuE41voB4GZ2\nLSTHPYHPiZ3wddtt0ty5XPkDcC/uCdgktvTTNuGLrh8AmY5TnCj9APAvV40E0j1PIHaZZ7p+AHgJ\n8wRS8Pkdvu6+mx2+AHhTqudO35WDYid8UfoB4HeuKgd1p9i1ftjhCwCMjE8C7PAFAJ1z1anQ7sli\nlH4AZComi51GS4u0eLH08MPs8AUgszFZLMapSj+s9QMAncuY02MkYko/+/ZR+gGAeHl+JBA74evy\ny+n6AYBEeDYJsMMXAKTOk6fMtq6ffftMIvj7v3c6IgDwJleNBLpaOyh2h68rrpDeeYcEAMCffLV2\nUOxaP9//vlRZyVo/ACD5YO2gtmWeDxyg6wcA7OaqclCs2K6fK66Q3nqLBAAAdnNdEqDrBwDSx1Wn\nVtb6AYD0ctVIoKiIHb4AIJ1c1R304YcWXT8AkIBUu4NclQRcEgoAeAariAIAkuaqJNDVjGEAgOGr\nGcMAgFOjHAQASBpJAAB8jCQAAD5GEgAAHyMJAICPkQQAwMdIAgDgYyQBAPAxkgAA+BhJAAB8LC2b\nyliWpVtvvVUtLS0aMWKErr322nR8LQCgC2kZCTz77LNqbm7WmWeeqUAgkI6vRJqx8J+3cfz8K+kk\nUFJSouzsbBUUFHR4PhQKKS8vT4MGDVJlZaUkaefOnRo9erTuuecePfjgg6lFDFfiJOJtHD//SjoJ\nFBcXKxQKdXguGo2qtLRUoVBIkUhE1dXVqq+vVyAQ0FlnnWW+8Ax334aw+48h2c+L933xvK6r13T2\n+0SfdwM7Y+vuYxfva0/3mkR/55djl8rn2Xn8vPC3l/QZubCwUH369OnwXG1trXJzc9W/f3/16tVL\nwWBQNTU1mjhxol555RVdf/31KioqSjXmbuW1f4gkgY5IAqf/nV+OXSqf57ckICsFDQ0N1tChQ9sf\nP/nkk9bMmTPbH69cudIqLS2N67MGDhxoSeKHH3744SeBn4EDB6ZyGrds7Q7KyspK+r27d++2MRIA\nQDxsLdD369dPTU1N7Y+bmproBgIAF7M1CYwYMUK7du1SY2OjTpw4oTVr1mjChAl2fgUAwEZJJ4Fp\n06Zp1KhR2rlzp3JycrRixQr17NlTVVVVGjdunIYMGaKpU6cqPz/fzngBADZyzUbzAID0c3fTPgCg\nW7k2CRw5ckTTp0/X7NmztWrVKqfDQYIaGho0c+ZMTZ482elQkKCamhrNnj1bwWBQr776qtPhIEE7\nduzQ3LlzNWXKFD300ENdvt615aCVK1eqb9++uvTSSxUMBrV69WqnQ0ISJk+erCeffNLpMJCEgwcP\nauHChVq2bJnToSAJJ0+eVDAY1BNPPHHa16V1JJDIekPNzc3KycmRJPXo0SOdYaITiRw/uEsyx+6O\nO+5QaWlpOsNEJxI9fs8//3z7BXSXUppqlqDXX3/d2rRpU4dZxq2trdbAgQOthoYG68SJE9bw4cOt\nSCRirVy50nrhhRcsy7KsYDCYzjDRiUSOX5srr7zSiVDxOYkcu5MnT1qLFi2y1q5d62DEiJXM355l\nWdaECRO6/Oy07CfQprCwUI2NjR2ei11vSFL7ekPXX3+9SktL9eKLLzLXwCUSOX7Z2dm65ZZbtHnz\nZlVWVuqmm25Kf8Bol8ixW7t2rdatW6fDhw9r9+7duu6669IfMDpI5Pjt3btXzzzzjI4dO6aLL764\ny89OaxI4ldiyjyQFAgFt3LhRvXv31vLlyx2MDPHo7Pj17dtXS5YscTAydKWzY3ffffdp/vz5DkaG\neHR2/MaMGaMxY8bE/TmOdwelst4QnMfx8y6OnbfZdfwcTwKsN+RtHD/v4th5m13Hz/EkwHpD3sbx\n8y6OnbfZdvxsv419GsFg0PrGN75hnXnmmVYgELCWL19uWZZlvfTSS9b5559vDRw40LrzzjvTGRIS\nwPHzLo6dt3Xn8XPtZDEAQPdzvBwEAHAOSQAAfIwkAAA+RhIAAB8jCQCAj5EEAMDHSAIA4GMkAQDw\nsf8DqLUiU8tNdKIAAAAASUVORK5CYII=\n",
"text": [
"<matplotlib.figure.Figure at 0x4ef2750>"
]
}
],
"prompt_number": 16
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"f=lambda x:x**2.\n",
"\n",
"test=np.arange(0,10)\n",
"-np.array(map(f, test))"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 23,
"text": [
"array([ -0., -1., -4., -9., -16., -25., -36., -49., -64., -81.])"
]
}
],
"prompt_number": 23
},
{
"cell_type": "code",
"collapsed": false,
"input": [],
"language": "python",
"metadata": {},
"outputs": []
}
],
"metadata": {}
}
]
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment