IPDOS Curves.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "import matplotlib.pyplot as plt\n",
    "import pandas as pd\n",
    "import numpy as np"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(-1.0, 0.07)"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAuEAAADgCAYAAACtg7geAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzt3Xm4JVV57/HvD1rQMCugiLagoIgDqO08IYMXjdqoTEYFDMjFCROHB25IvA5JbscpatR4uaI2ahRUFExwgFY0DiiNMgjKICIiCCqKoMj43j+qDmxOn737cPY+tc/w/TzPfk7VqrVrvdWnWf2yatWqVBWSJEmSurPOuAOQJEmSFhuTcEmSJKljJuGSJElSx0zCJUmSpI6ZhEuSJEkdMwmXJEmSOmYSLkmSJHXMJFyLQpLrez63JbmhZ//FSd6c5OZJ9X7f8/3lSc5K8ockv0myKsk2ST7UU/+mSef4Uvvdo5Nc0LZ70Nj+ECRpjkhy6aR++Pok953G9z7W9rW931u35/huSX6S5E9Jvp7kAX3O85Ukb52ifHmSXyVZ0uffhQf21N05yZltW2cm2Xmmfx5anEzCtShU1YYTH+Ay4Lk9ZZ9sqx3XW6+qNgVIsh1wLPB6YBNgW+CDwG1VdVjPef950jme1Z73bOCVwA+6u2JJmvOeO6nPvWKa33v7pO/dCpBkc+AE4B+AewKrgeP6nONjwEuTZFL5S4FPVtUt7f7kfxcuadtaDzgR+ASwGbASOLEtl6bFJFxau52Bn1XVqmpcV1Wfq6rLpvPlqvpAVa0C/jy7YUrS/JfkCUm+k+T3Sc5Osss0v/oC4Lyq+kxV/Rl4M7BTkh2mqPsFmkT9qT3tbgY8h2bQZW12AZYA76mqG6vqfUCAXacZq2QSLk3DD4Adkvxrkmck2XDcAUnSQpRka+C/gH+kSZLfAHwuyRY91V6Z5Jp2CsgLe8ofRnPnEYCq+iPw07b8TqrqBuB44ICe4n2Bn1TV2T1lz23bOi/JKya1dU5VVU/ZOVO1JfVjEi7dYd925GXi83WA9vbjLsDWNJ32b9p5iSbjkjRzX+jpb7/Qlr0EOLmqTq6q26rqFJppJc9uj78P2B7YkmbayceSPLk9tiFw7aQ2rgU26tP+SmCfJPdo9w9oyyYcDzwU2AJ4OfCmJC+aYVvSGkzCpTscX1Wb9nyeMXGgqk6vqn2ragua25dPA44aW6SSNP/t1dPf7tWWPYAmMb59QAR4CrAVQFX9oKp+W1W3VNXJwCdppqEAXA9sPKmNjYHrpmq8qr4F/BpY3j5w+VjgP3qOn19VV1TVrVX1HeC9wN4zaUuaikm4dBdV1Rk0D/88fNyxSNIC8wvg45MGRDaoqhV96hfNXGyA84CdJg4k2QB4UFvez7E0I+AvBb5aVVcNqDu5rUdOerDzkWtpS7oTk3BpLZI8JcnLk2zZ7u8APA84fZrfXy/J3Wk677sluXsS/9uTpDV9gmYe9v9Ism7bX+6S5H4ASfZOsmGSdZI8k2b6ykntdz8PPDzJC9s+900087Z/MqC9Y4Hdaaab9E5FmViucLM0HgccTrMiCsBpwK3A4UnWT/Lqtvxrw/4BaPEwEZDusN+k9WCvbxPv39Mk3ecmuR74Mk1n//ZpnverwA3Ak4Cj2+2njT58SZrfquoXwHLg72imivwCeCN35CuvBX5J0y+/A3h5VZ3WfvfXwAuBfwJ+Bzwe2H8t7V0KfAfYgDuS+Qn7AxfTTDE5FviXqlrZfu8mYC+aUfTfA39NM73mphlduBal3PnBXkmSJEmzzZFwSZIkqWMm4ZIkSVLHTMIlSZKkjpmES5IkSR0zCZckSZI6tmTcAYza5ptvXttss824w5CkGTnzzDN/076ZddGw35Y0Xw3TZy+4JHybbbZh9erV4w5DkmYkyc/HHUPX7LclzVfD9NlOR5EkSZI6ZhIuSZIkdWysSXiSPZNckOTiJEdOcXz9JMe1x7+XZJvuo5QkSZJGa2xJeJJ1gQ8AzwJ2BF6UZMdJ1Q4GfldV2wH/CvxLt1FKkiRJozfOBzMfB1xcVZcAJPk0sBw4v6fOcuDN7fZngfcnSVVV37PeeB1cctpsxCu1Mu4AJEnSPDfOJHxr4Bc9+5cDj+9Xp6puSXItcC/gN72VkhwKHArwmK3WgWOXz1bMkrSoJdkTeC+wLvDhqlox6fj6wLHAY4DfAvtV1aVdxylJc904k/CphhMnj3BPpw5VdTRwNMCyR+5QvOwjw0cnTWXATRhpeAVveeq4g+irZxrhHjQDJ2ckOamqeu9g3j6NMMn+NNMI9+s+Wkma28aZhF8O3L9n/37AFX3qXJ5kCbAJcM3As663ITzgSSMMU5LUmp1phH+6Bs761KwELE0xdifNCeNMws8Atk+yLfBLYH/grybVOQk4EPgusDfwtYEduSRpNs3eNMIvHDZbMUvSnDS2JLztnF8NfIVmbuFHquq8JG8FVlfVScAxwMeTXEwzAr7/uOKVJM3SNMKdH1EcftLw0Un9xAfqNUvesu2MvzrW19ZX1cnAyZPK3tSz/Wdgn67jkiRNaXamES5ZH+4583/IJGk+8o2ZkqTpun0aYZL1aO5OTh7CnphGCE4jlKS+xjoSLkmaP5xGKEmjYxIuSZo2pxFK0mg4HUWSJEnqmEm4JEmS1DGTcEmSJKljJuGSJElSx0zCJUmSpI6ZhEuSJEkdMwmXJEmSOmYSLkmSJHXMJFySJEnqmEm4JEmS1DGTcEmSJKljJuGSJElSx0zCJUmSpI6ZhEuSJEkdMwmXJEmSOmYSLkmSJHWsbxKe5AFJNunZf0aS9yZ5XZL1uglPkiRJWngGjYQfD2wAkGRn4DPAZcBOwAdnPzRJkiRpYVoy4Ng9quqKdvslwEeq6l1J1gHOmv3QJEmSpIVp0Eh4erZ3BVYBVNVtwzaa5J5JTklyUftzsynq7Jzku0nOS3JOkv2GbVeSJEmaCwYl4V9LcnyS9wKbAV8DSLIVcNOQ7R4JrKqq7WmS+yOnqPMn4ICqehiwJ/CeJJsO2a4kSZI0doOS8L8BTgAuBZ5SVTe35fcBjhqy3eXAynZ7JbDX5ApVdWFVXdRuXwFcDWwxZLuSpBnwDqYkjVbfJLwanwa+ADwqyV8meWBV/bCqvjJku/euqivbdq4EthxUOcnjgPWAnw7ZriRpZryDKUkj1PfBzCQbAx8GHgOcTTNHfKckZwIHV9UfBp04yak0o+aT3aVR9Hb6y8eBA/vNR09yKHAowNKlS+/K6SVJ07Mc2KXdXgmcBhzRW6GqLuzZviLJxB3M33cToiTNH4NWR3kfcD6w/0TymyTAPwDvBw4YdOKq2r3fsSRXJdmqqq5sk+yr+9TbGPgv4O+r6vQBbR0NHA2wbNmyGhSXJGlG7nQHM4l3MCVpCIOS8CdX1UG9BVVVwFuTXDRkuycBBwIr2p8nTq7QvhDo88CxVfWZIduTJK2FdzAlqTuDkvAMODasFcDxSQ6meQHQPgBJlgGHVdUhwL7A04B7JTmo/d5BVeUa5ZI0C7yDKUndGZSEfzvJm4C3tSPgACT5B6BvxzodVfVbYLcpylcDh7TbnwA+MUw7kqSR8Q6mJI3QoCUKXwM8Arg4yeeSfDbJJTSvrX91J9FJkuaKFcAe7XTEPdp9kixL8uG2zsQdzIOSnNV+dh5PuJI0t/UdCW9XP9knyYOAHWmmpxxRVT5kI0mLjHcwJWm0Bk1HIckSmgR8h7botiQ/r6pbZj0ySZIkaYHqOx0lyX2B84DXA/cFtgbeCJzXHpMkSZI0A4NGwv8Z+Peqek9vYZLDgf9D82COJEmSpLtoUBL+hMnrhANU1fuSXDB7IUmSJEkL26DVUW4YcOxPow5EkiRJWiwGjYRvkuQFU5QH2HiW4pEkSZIWvEFJ+DeA5/Y59s1ZiEWSNENJ/gK4uapubvcfAjwb+HlVnTDW4CRJaxi0TvjLugxEkjSULwMHAxcl2Q74LvBJ4DlJHltV/2us0UmS7mTQEoWPT3J2kuuTfDfJQ7sMTJJ0l2xWVRe12wcCn6qq1wDPAp4zvrAkSVMZ9GDmB4A3APcC3g28Z0BdSdJ4Vc/2rsApAFV1E3DbWCKSJPU1KAlfp6pOqaobq+ozwBZdBSVJusvOSfLOJH8LbAd8FSDJpuMNS5I0lUEPZm46aXWUO+37oI8kzSkvB14LbAM8s6omlpLdEXjnuIKSJE3trqyO0rtfgEm4JM0RVXUDsCLJ3YHtkjwM+GlVfQf4znijkyRN5uookrQAJFkC/DPwMuAymumG90vyUeCoiaULJUlzQ98kPMnrBn2xqt49+nAkSTP0DmAj4IFVdR1Ako1ppqK8k2aqiiRpjhg0HWWjzqKQJA3rOcCDq+r2VVKq6g9JXgH8BJNwSZpTBk1HeUuXgUiShlK9CXhP4a1J1iiXJI3XoCUKJUnzx/lJDphcmOQlNCPhkqQ5ZNB0FEnS/PEq4IQkfw2cSbOK1WOBewDPH2dgkqQ1mYRL0gJQVb8EHp9kV+BhQIAvVdWq8UYmSZrKoNVR1rit2auqjp1po0nuCRxH81KJS4F9q+p3fepuDPwY+HxVvXqmbUrSIvFr4Fft9lXjDESS1N+gkfDHTlEWmhf2bA3MOAkHjgRWVdWKJEe2+0f0qfs2mhcFSZL6SLIJcCJwf+Acmv76EUkuA5ZX1R/GGZ8k6c76PphZVa+Z+ACHA98Dng6cDjx6yHaXAyvb7ZXAXlNVSvIY4N7AV4dsT5IWurcBq4Htq+r5VbUXsD1wBvBPw548yT2TnJLkovbnZgPqbpzkl0neP2y7krRQDVwdJcmSJIcA5wO7A3tX1X5Vdc6Q7d67qq4EaH9uOUXb6wDvAt44ZFuStBjsDhxZVbdNFLTbf9ceG9bEHcztgVXtfj/ewZSktRg0J/xVNC93WAXsWVU/vysnTnIqcJ8pDh01zVO8Eji5qn6RZG1tHQocCrB06dK7EqYkLRQ3VdUtkwur6pYkN47g/MuBXdrtlcBpTDGNsOcO5peBZSNoV5IWpEFzwv8NuBp4CvDFnkQ4NC+FeOSgE1dV35GXJFcl2aqqrkyyVdvOZE8EnprklcCGwHpJrq+qNUZfqupo4GiAZcuW+VIKSYvR3ZM8iqaP7hVg/RGc/053MJMMuoP5UmC3QSdz8ETSYjcoCd92Fts9CTgQWNH+PHFyhap68cR2koOAZVMl4JIkAK4E3t3n2K/6lN9Jl3cwHTyRtNgNem39zwGSbEuz5mwBP66qS0bQ7grg+CQHA5cB+7RtLQMOq6pDRtCGJC0aVfWMEZyjszuYkrTYDZoTvjHwYZo5fWfR3NLcKcmZwMHDLHdVVb9liluVVbUaWCMBr6qPAR+baXuSpKF5B1OSRmjQ6ijvo1kVZbuqekFVPR94EHAu4LJTkrS4rAD2SHIRsEe7T5JlST481sgkaR4aNCf8yVV1UG9BVRXw1rYTliQtEt7BlKTRGpSED36qRpI0JyXZGngAPX18VX1zfBFJkiYblIR/O8mbgLe1I+AAJPkHmrdmSpLmmCT/AuxHM53w1ra4AJNwSZpDBiXhrwGOAS5OchZNJ/4o4IdMcetRkjQn7AU8pKpG8YIeSdIsGbRE4R+AfZI8CNiRZnrKEVX1066CkyTdZZcAdwNMwiVpDhu0ROGje3Z/2f7cZKK8qn4wm4FJkmbkT8BZSVbRk4hX1eHjC0mSNNmg6SjvGnCsgF1HHIskaXgntR9J0hw2KAn/u6r6bmeRSJKGVlUrk6wHPLgtuqCqbh5nTJKkNQ16Wc8HOotCkjQSSXYBLqLpwz8IXJjkaWMNSpK0BtcJl6SF5V3AM6vqAoAkDwY+BTxmrFFJku5kUBK+bZK+8wqr6nmzEI8kaTh3m0jAAarqwiR3G2dAkqQ1DUrCf83ghzMlSXPP6iTHAB9v918MnDnGeCRJUxiUhF9fVd/oLBJJ0ii8AngVcDjNtMJv0swNlyTNIYOS8J91FoUkaSTaN2W+u/1IkuaoQaujfD7JSycXJnl5kr+axZgkSTOU5MlJTklyYZJLJj7jjkuSdGeDRsJfB0y1rNWngdOA/5iNgCRJQzkG+FuaeeC3jjkWSVIfg5LwdavqusmFVXWdT9pL0px1bVV9adxBSJIGG5SE3y3JBlX1x97CJBsB681uWJKkGfp6kncAJwA3ThRW1Q/GF5IkabJBSfgxwGeTvKKqLgVIsg3NW9iOmfXIJEkz8fj257KesgJ2HUMskqQ++ibhVfXOJNcD30iyIU0n/kdgRVX9e1cBSpKmr6qeMe4YJElrN2gknKr6EPChNgnPxBzxJPeuqqu6CFCSNH1JNgH+N3c8WP8N4K1Vde34opIkTTZoicLbVdX1wDpJ/jrJqcBQcwuT3LNdQuui9udmfeotTfLVJD9Ocn47HUaS1N9HgOuAfdvPH4CPjjUiSdIaBibhSe6RZL8kJwI/onn5wz8C9x+y3SOBVVW1PbCq3Z/KscA7quqhwOOAq4dsV5IWugdV1f+uqkvaz1uABw57UgdPJGm0+ibhST4JXAg8E3g/sA3wu6o6rapuG7Ld5cDKdnslsNcU7e8ILKmqU6AZja+qPw3ZriQtdDckecrETpInAzeM4LwOnkjSCA2aE/5w4HfAj4GfVNWtSWpE7d67qq4EqKork2w5RZ0HA79PcgKwLXAqcGRV+fIJServMODYdm54gGuAg0Zw3uXALu32SpqXth3RW2GqwZMRtCtJC9Kg1VF2SrID8FfAqUmuBjZKcp+q+tXaTtzOHb/PFIeOuguxPRV4FHAZcBzNPyRrLI+Y5FDgUIClS5dO8/SStPBU1dnATkk2bvf/MKJTO3giSSO0ttVRfgK8CXhTkmXAi4DvJ7m8qp60lu/u3u9YkquSbNV25Fsx9e3Ky4EfVtUl7Xe+ADyBKZLwqjoaOBpg2bJloxqtl6R5J8n6wAtpphAuSQJAVb11Gt918ESSOjIwCe9VVauB1UneADx7yHZPAg4EVrQ/T5yizhnAZkm2qKpf07xoYvWQ7UrSQncicC1wJj1vzJwOB08kqTsDk/AkWwNbAedU1U3t7ce/oRnZuO8Q7a4Ajk9yMM1oyT5te8uAw6rqkHYO+huAVWmGcs4E/t8QbUrSYnC/qtpzFs7r4IkkjVDfJDzJ39DcgrwYWD/Je2mWKDwWeMwwjVbVb4HdpihfDRzSs38K8Mhh2pKkReY7SR5RVeeO+LwOnkjSCA0aCT8UeEhVXZNkKU0y/rSqOr2b0CRJ05XkXKBo+vWXJbmEZjpKgKqqoQY0HDyRpNEalIT/uaquAaiqy5JcaAIuSXPWc8YdgCRp+gYl4fdL8r6e/S1796vq8NkLS5J0F11Fs0b4dsC5wDFVdct4Q5Ik9TMoCX/jpP0zZzMQSdJQVgI3A/8NPAvYEXjtWCOSJPU16GU9K/sdSzLtpQ0lSZ3YsaoeAZDkGOD7Y45HkjTAOv0OJPlWz/bHJx22c5ekueXmiQ2noUjS3DdoRHuDnu2HTTqWWYhFkjRzOyWZeEV9gHu0+xOro2w8vtAkSZMNSsIHvcHMt5tJ0hxSVeuOOwZJ0vQNSsI3TfJ8mikrmyZ5QVseYJNZj0ySJElaoAYl4d8Antez/dyeY9+ctYgkSZKkBW5QEv7Fqjqhs0gkSZKkRaLv6ijA33cWhSRJkrSIDErCJUmSJM2CQdNRdkhyzhTlE8tdPXKWYpIkSZIWtEFJ+M+488OYkiRJkkZgUBJ+U1X9vLNIJEmSpEVi0Jzwb3cWhSRJkrSIDBoJPyPJAf0OVtWxsxCPJEmStOANSsKXTVEWmnniWwMm4ZIkSdIM9E3Cq+o1E9tJArwYOAI4Hfin2Q9NkiRJWpgGjYSTZAlwEPB64HvA3lV1QQdxSZIkSQtW3yQ8yauA1wKrgD1dKUWSJEkajUEj4f8GXA08BfhiMyMFGMHLepLcEzgO2Aa4FNi3qn43Rb23A39Js4rLKcBrq6pm2q4kSZI0FwxKwredxXaPBFZV1YokR7b7R/RWSPIk4MnARLL/LeDpwGmzGJckaQoOnkjSaPVdJ7yqfj7oM2S7y4GV7fZKYK+pQgDuDqwHrA/cDbhqyHYlSTMzMXiyPc00xSMnV5g0ePJw4LE0gyeSpEkGzQm/jiYRXuMQzXSUjYdo995VdSXNia5MsuXkClX13SRfB65s23x/Vf14iDYlSTO3HNil3V5Jc1fyiEl1egdPgoMnktTXoCUKNxrmxElOBe4zxaGjpvn97YCHAvdri05J8rSq+uYUdQ8FDgVYunTpzAKWJA0y0sET+21Ji93AJQqHUVW79zuW5KokW7Ud+VY0D4BO9nzg9Kq6vv3Ol4AnAGsk4VV1NHA0wLJly5x7KEkz0OXgif22pMWu75zwWXYScGC7fSBw4hR1LgOenmRJkrvRzCt0OookzZKq2r2qHj7F50TgqnbQhOkMnrQDKBODJ5KkScaVhK8A9khyEbBHu0+SZUk+3Nb5LPBT4FzgbODsqvriOIKVJDl4IkmjNGvTUQapqt8Cu01Rvho4pN2+FfifHYcmSZraCuD4JAfTJNv7QDN4AhxWVYfQDJ7sSjN4UsCXHTyRpKmNJQmXJM0vDp5I0miNazqKJEmStGiZhEuSJEkdMwmXJEmSOmYSLkmSJHXMJFySJEnqmEm4JEmS1DGTcEmSJKljJuGSJElSx0zCJUmSpI6ZhEuSJEkdMwmXJEmSOmYSLkmSJHXMJFySJEnqmEm4JEmS1DGTcEmSJKljJuGSJElSx0zCJUmSpI6ZhEuSJEkdMwmXJEmSOmYSLkmSJHXMJFySJEnq2FiS8CT7JDkvyW1Jlg2ot2eSC5JcnOTILmOUJN3BfluSRmtcI+E/Al4AfLNfhSTrAh8AngXsCLwoyY7dhCdJmsR+W5JGaMk4Gq2qHwMkGVTtccDFVXVJW/fTwHLg/FkPUJJ0J/bbkjRac3lO+NbAL3r2L2/LJElzk/22JE3TrI2EJzkVuM8Uh46qqhOnc4opyqpPW4cCh7a7Nyb50fSiXDA2B34z7iA65jUvDovxmh8yrobttzu12P5uL7brBa95sZhxnz1rSXhV7T7kKS4H7t+zfz/gij5tHQ0cDZBkdVX1fWhoIfKaFweveXFIsnpcbdtvd2exXfNiu17wmheLYfrsuTwd5Qxg+yTbJlkP2B84acwxSZL6s9+WpGka1xKFz09yOfBE4L+SfKUtv2+SkwGq6hbg1cBXgB8Dx1fVeeOIV5IWO/ttSRqtca2O8nng81OUXwE8u2f/ZODku3j6o4eLbl7ymhcHr3lxmJPXbL89covtmhfb9YLXvFjM+JpTNeUzM5IkSZJmyVyeEy5JkiQtSPM2CV/bq5GTrJ/kuPb495Js032UozWNa35dkvOTnJNkVZIHjCPOUZruK7CT7J2kBr1Oe76YzjUn2bf9XZ+X5D+6jnHUpvF3e2mSryf5Yfv3+9lTnWe+SPKRJFf3W5Yvjfe1fx7nJHl01zGOmn22ffakevbZ85h99hrHZ9ZnV9W8+wDrAj8FHgisB5wN7DipziuBD7Xb+wPHjTvuDq75GcBftNuvWAzX3NbbiOZV2qcDy8Yddwe/5+2BHwKbtftbjjvuDq75aOAV7faOwKXjjnvIa34a8GjgR32OPxv4Es26208AvjfumDv4HdtnL4JrbuvZZ8+B2Gf5mu2zp3He+ToSfvurkavqJmDi1ci9lgMr2+3PArslg9+3PMet9Zqr6utV9ad293SaNXrns+n8ngHeBrwd+HOXwc2S6Vzzy4EPVNXvAKrq6o5jHLXpXHMBG7fbm9Bn7en5oqq+CVwzoMpy4NhqnA5smmSrbqKbFfbZ9tm97LPnN/vsNc2oz56vSfh0Xo18e51qls26FrhXJ9HNjrv6OuiDaf6vbD5b6zUneRRw/6r6zy4Dm0XT+T0/GHhwkm8nOT3Jnp1FNzumc81vBl7SLpF3MvCabkIbm4X2+nf7bPtswD7bPnvBmlGfPZYlCkdgOq9Gnvbrk+eJu/I66JcAy4Cnz2pEs2/gNSdZB/hX4KCuAurAdH7PS2hub+5CM3L230keXlW/n+XYZst0rvlFwMeq6l1Jngh8vL3m22Y/vLFYjP3XYrzmpqJ99nxmn92wz17TWvuv+ToSPp1XI99eJ8kSmtshg24lzHXTeh10kt2Bo4DnVdWNHcU2W9Z2zRsBDwdOS3IpzTysk+b5gz7T/bt9YlXdXFU/Ay6g6eDnq+lc88HA8QBV9V3g7sDmnUQ3HtN+/fs8YZ9tnw322fbZC9eM+uz5moRP59XIJwEHttt7A1+rdvb8PLXWa25v8/1fms58vs85g7Vcc1VdW1WbV9U2VbUNzZzK51XV6vGEOxLT+bv9BZoHukiyOc2tzks6jXK0pnPNlwG7ASR5KE2H/utOo+zWScAB7RP3TwCuraorxx3UEOyz7bPts7HP7jTKbs2szx73E6cz/dA8iXohzRO6R7Vlb6X5DxqaX/hngIuB7wMPHHfMHVzzqcBVwFnt56Rxxzzb1zyp7mnM8yftp/l7DvBu4HzgXGD/ccfcwTXvCHyb5in8s4BnjjvmIa/3U8CVwM00IygHA4cBh/X8jj/Q/nmcu0j+Xttn22fPy499tn32TPts35gpSZIkdWy+TkeRJEmS5i2TcEmSJKljJuGSJElSx0zCJUmSpI6ZhEuSJEkdm69vzJRmJMmtNMsHTfh0Va0YVzySpP7ss7WQuUShFpUk11fVhiM+55KqumWU55Qk2WdrYXM6igQkuTTJW5L8IMm5SXZoyzdI8pEkZyT5YZLlbflBST6T5IvAV5Osk+SDSc5L8p9JTk6yd5Ldkny+p509kpwwpsuUpAXBPlsLgUm4Fpt7JDmr57Nfz7HfVNWjgX8H3tCWHUXz+uzH0rx2+B1JNmiPPRE4sKp2BV4AbAM8AjikPQbwNeChSbZo918GfHSWrk2SFhr7bC1YzgnXYnNDVe3c59jEaMeZNB00wDOB5yWZ6ODvDixtt0+pqmva7acAn6mq24BfJfk6QFVVko8DL0nyUZqO/oDRXY4kLWj22VqwTMKlO9zY/ryVO/7bCPDCqrqgt2KSxwN/7C0acN6PAl8E/kzT6TsXUZKGZ5+tec3pKNJgXwFekyThA5JBAAAA3UlEQVQASR7Vp963gBe28wzvDewycaCqrgCuAP4e+NisRitJi5t9tuYNR8K12NwjyVk9+1+uqiMH1H8b8B7gnLZTvxR4zhT1PgfsBvwIuBD4HnBtz/FPAltU1flDxC5Ji419thYslyiURiTJhlV1fZJ7Ad8HnlxVv2qPvR/4YVUdM9YgJUmAfbbGz5FwaXT+M8mmwHrA23o68zNp5iK+fpzBSZLuxD5bY+VIuCRJktQxH8yUJEmSOmYSLkmSJHXMJFySJEnqmEm4JEmS1DGTcEmSJKljJuGSJElSx/4/xhi9myI4lJYAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 864x216 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "#Data sets for FeV3 and Fe3V alloys\n",
    "\n",
    "d1 = pd.read_csv('TEST11.csv')\n",
    "d2 = pd.read_csv('Fe50V50_measured_DOS.csv')\n",
    "\n",
    "fig1, axes = plt.subplots(nrows=1, ncols=2, figsize=(12,3))\n",
    "\n",
    "#FIRST GRAPH\n",
    "axes[0].plot(d1)\n",
    "axes[0].set_title('TEST1')\n",
    "axes[0].set_xlabel('Energy')\n",
    "axes[0].set_xlim([0.0,1.0])\n",
    "axes[0].set_ylabel('INTERACTION IPDOS')\n",
    "axes[0].set_ylim([-1.0,0.07])\n",
    "\n",
    "#THIRD GRAPH \n",
    "axes[1].plot(d2)\n",
    "axes[1].set_title('Fe50 V50')\n",
    "axes[1].set_xlabel('Energy')\n",
    "axes[1].set_xlim([0.0,1.0]) #Tried rescaling energy from 0-1... fail\n",
    "axes[1].set_ylabel('Phonon DOS')\n",
    "axes[1].set_ylim([-1.0,0.07])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}