{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "5eb52d01-f75e-47dc-bdf1-da61f5841c56",
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "7d33e159-d051-4fb8-a14f-04ca8c7db66f",
   "metadata": {},
   "source": [
    "Diesel cycle efficiency term. Greater than zero implies Otto cycle is more efficiency (for a given compression ratio)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "388b373e-e002-435f-a2b0-3f278bacbb61",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAD4CAYAAAD8Zh1EAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAfwElEQVR4nO3deVxVdeL/8dcHBUVESXFBAXHfcEPcmjLbsyxLW7TFyoqWmZZpmek709RUM+3bqJWZmVmpNVamtqilZatLBooLizuoASIgIrLcz+8P8PvzayIol3vuvbyfj4ePLp4T5/346H17OPdzPsdYaxEREd8X4HQAERFxDxW6iIifUKGLiPgJFbqIiJ9QoYuI+ImGTh04PDzcxsTEOHV4ERGf9Msvv+RYa1sdb5tjhR4TE8OaNWucOryIiE8yxuyoapsuuYiI+AkVuoiIn1Chi4j4CRW6iIifUKGLiPgJFbqIiJ9QoYuI+AkVuoiIhxQeLuONb7ewaltunXx/x24sEhGpL/YfLOHtH7fzzo/byT9Uyp0jOjO4Ywu3H0eFLiJSR7IKipn+/Tbe+3kHRSXlnN+rDX88uwv9o8Lq5HgqdBERN9uVW8QbK7bw4ZoMyspdXNqvHXeN6EL3tqF1elwVuoiIm2zNLuS1b7Yw/9dMjIErB0Zy+/DOxISHeOT4KnQRkVratKeAV5en89n6PTRqGMANwzqQMLwTEc2DPZqj2kI3xswARgFZ1trY42wfAXwKbKv8rY+ttU+4MaOIiFdK2pXHlOXpLN34GyFBDbh9eGduPbMj4U0bOZKnJmfoM4EpwKwT7POdtXaUWxKJiHi5NdtzmbQsnRWp2TRr3JD7zuvKTafHENYkyNFc1Ra6tXaFMSbGA1lERLyWtZaftu5j0tdp/Lw1l5YhQfz1oh5cPzSa0MaBTscD3HcNfZgxJgnYDTxord1wvJ2MMQlAAkB0dLSbDi0iUnestXybms2UZems2bGf1qGNeOSSnlw7JJomQd71MaQ70qwFOlhrC40xFwPzga7H29FaOw2YBhAfH2/dcGwRkTphreXrTVlMXpZGUkY+7Zo35snRvbkqPorGgQ2cjndctS50a23BUa8/N8a8ZowJt9bm1PZ7i4h4mstlWbJxL5O+TmfjngKiWgTz9Jg+jI2LJKihd6+WUutCN8a0BX6z1lpjzGAq1ofZV+tkIiIe5HJZvkjey+RlaWzee4CO4SE8f2VfLh/QnsAG3l3kR9Rk2uIcYAQQbozJAB4DAgGstVOBK4E7jTFlwCFgnLVWl1NExCeUuyyfrd/D5K/TSMsqpHOrEF65pj+j+kbQ0EeK/IiazHIZX832KVRMaxQR8RnlLsuidbuZ9HUaW7IP0rV1UyaNH8AlfSJoEGCcjndKvOsjWhGROlZW7mLhut1MXpbO1uyDdG8TyqvXxjEyti0BPlrkR6jQRaRe+N8i/zqdrTkH6dE2lNevi+PC3r5f5Eeo0EXErx2vyKdeP5ALerXxmyI/QoUuIn6p3GVZkJR5TJHHcUEv/zkjP5YKXUT8SrnLsjCp4sNOfz8jP5YKXUT8wrGzVurDGfmxVOgi4tOO3BD0yleppGUV0q1NU167Lo6L/OjDzppSoYuITzpyi/4rX1Xc2dmldVOmXDuAi2Mj6l2RH6FCFxGfcmTRrJeWprJxTwGdwkP4z7j+jOrbzmdvCHIXFbqI+IQjy9i+vDSVpIx8ols04cWr+jG6fzufu0W/rqjQRcTr/Ziew4tLU/llx37ahwXz7Ng+jImL9JlFszxFhS4iXmvN9lxeXJLKT1v30bZZY568PJZr4qO8fhlbp6jQRcTrrMvI48UlqXybmk1400Y8OqoX1w6J9toHS3gLFbqIeI3Newt4aUkqSzb+RliTQB4e2YMJwzp43aPevJVGSUQcty3nIC8vTWXhut00DWrIn8/rxsQzYrzm4cu+QoUuIo7JzDvEpK/SmLc2g6AGAdxxVmduH96JsCZBTkfzSSp0EfG4rAPFvLZ8C7NX7gRgwrAO3DWiC61CGzmczLep0EXEY/KLSnljxRbe/mE7JeUurhoYyd3ndqV9WLDT0fyCCl1E6tzBw2XM/HE7U7/dwoHiMi7r144/n9+NjuEhTkfzKyp0Eakzh8vKmbNyJ1OWp5NTWMJ5PVvzwAXd6RnRzOlofkmFLiJuV+6yfLw2g1e+SiMz7xBDOrbgjRu6M7BDC6ej+TUVuoi4jbWWJRt/44XFKaRlFdKnfXOeHtOHM7uGY0z9XjjLE1ToIuIWP23Zx7NfbiZxVx6dwkN49do4Lu7TVkXuQSp0EamV5Mx8nlucworUbCKaN+bZsX0YGxepFRAdoEIXkVOyPecgLy5NZWHSbsKaBPL3i3tyw7AOWm/FQSp0ETkpWQeKmfx1OnNW7SSwQQB/PLszt5/VmWa6Td9xKnQRqZEDxaW8uWIrb363jdJyF+MGR3HPOV1p3ayx09GkkgpdRE7ocFk5s1fuZPKydHIPljCqbwQPXNBdNwV5IRW6iByXy2VZuG43LyxJYVfuIU7v3JKHR/agb2SY09GkCtUWujFmBjAKyLLWxp5gv0HAT8A4a+0890UUEU/7Li2bZ77YzIbdBfSKaMasiZpL7gtqcoY+E5gCzKpqB2NMA+BZYIl7YomIEzbszueZLzbzXVoO7cOCefmafozu156AABW5L6i20K21K4wxMdXsdjfwETDIHaFExLMy9hfx0pJUPknMpHlwII9cUjEFsVFDTUH0JbW+hm6MaQ9cAZxNNYVujEkAEgCio6Nre2gRqaX8olJe+yadt3/cDkDC8E7cNaILzYM1BdEXueND0VeAv1prXdVdX7PWTgOmAcTHx1s3HFtETsHhsnLe/WkHk5elU1BcypgBkdx/QTetS+7j3FHo8cDcyjIPBy42xpRZa+e74XuLiBtZa1m0bg/PLd7MrtxDnNk1nP8Z2ZNe7bScrT+odaFbazseeW2MmQksUpmLeJ9V23L59+ebSNqVR4+2ocyaOJjh3Vo5HUvcqCbTFucAI4BwY0wG8BgQCGCtnVqn6USk1rZmF/Lsl5tZvOE32jZrzPNX9mVMXCQNNHPF79Rklsv4mn4za+1NtUojIm6Te7CESV+n8d7PO2jUMIAHzu/GrWd2IjhIM1f8le4UFfEzxaXlvPPjdqYsT6eopJxxg6K477xutApt5HQ0qWMqdBE/Ya3ls/V7eOaLzWTsP8Q5PVrzPyN70LVNqNPRxENU6CJ+YO3O/fxr0UbW7qz4wPO9W4ZwRtdwp2OJh6nQRXxYxv4invsyhQVJu2kV2ojnxvZl7EB94FlfqdBFfNCB4lJe/2YL07/fRoCBe87pwu1ndSakkd7S9Zn+9EV8SLnL8uGaXby4JIWcwhKuGNCehy7sTjvd4Smo0EV8xo/pOTyxaCOb9x5gYIfTmH7jIPpHhTkdS7yICl3Ey23POchTn29iycbfaB8WzJRrB3BJnwitTS6/o0IX8VIFxaVMWZbO2z9sI6hBAA9d2J1bzuhI40DdGCTHp0IX8TJHrpO/sDiF3KISroyL5KELu+thzFItFbqIF1m5dR+PL9zIxj0FDIo5jZmjBtMnsrnTscRHqNBFvEDG/iKe/nwzn63fQ/uwYCaPH8CovrpOLidHhS7ioKKSMqZ+s4U3VmzFGPjzed1IGK4FtOTUqNBFHGCtZeG6PTz9+Sb25BdzWb92PDyyh+aTS62o0EU8bMPufB5fsJFV23Pp3a4Zk8YPYFBMC6djiR9QoYt4SO7BEl5cksKcVTsJaxLE02P6cHV8lNZdEbdRoYvUsbJyF7NX7eTFJakUHi5jwrAY/nxeN5o3CXQ6mvgZFbpIHVq5dR+PLdjA5r0HOL1zS/55WW+6aX1yqSMqdJE6sDe/mKc+38SCpN20Dwvm9eviuCi2raYhSp1SoYu4UUmZixk/bGPS12mUuSz3nNuVO8/qrGmI4hEqdBE3WZGazT8XbGBrzkHO69mGR0f1IrplE6djST2iQheppcy8Qzy5cCNfbthLTMsmvH3zIM7u3trpWFIPqdBFTtHhsnLeXLGVKcvTAXjowu7cemZHGjXU5RVxhgpd5BR8m5rNY58ms31fESNj2/LIqF60112e4jAVushJOPrySsfwEN6ZOJizurVyOpYIoEIXqZGSMhdvfV8xe8VidXlFvJIKXaQaP23Zxz8+TSY9q5Dze1XMXolqodkr4n1U6CJVyDpQzFOfbWJ+4m6iWgTz1o3xnNuzjdOxRKqkQhc5RrnL8v7KHTy/OIXDpS7uOacLd53dRc/yFK9XbaEbY2YAo4Asa23scbaPBp4EXEAZcJ+19nt3BxXxhHUZefz9k2TWZ+ZzRpdwHh/dm86tmjodS6RGanKGPhOYAsyqYvvXwAJrrTXG9AU+BHq4J56IZxQUl/LC4hTe/XkH4U0bMWn8AC7VI+DEx1Rb6NbaFcaYmBNsLzzqyxDAuiGXiEcceXLQk4s2sq/wMBOGduCBC7vTrLGWthXf45Zr6MaYK4CngdbAJSfYLwFIAIiOjnbHoUVO2Y59B3lkfjLfpeXQp31z3roxnr6RYU7HEjllbil0a+0nwCfGmOFUXE8/r4r9pgHTAOLj43UmL44oKXMxbcUWJi9LJ7BBAP+8tBc3DIvRk4PE57l1lkvl5ZlOxphwa22OO7+3iDus3p7L3z5eT1pWIRf3acujo3rTtnljp2OJuEWtC90Y0wXYUvmhaBzQCNhX62QibpRXVMIzX2xm7updtA8LZsZN8ZzTQ3PKxb/UZNriHGAEEG6MyQAeAwIBrLVTgbHABGNMKXAIuMZaq8sp4hWstSxI2s2Tizayv6iU24d34t7zutIkSLdgiP+pySyX8dVsfxZ41m2JRNxkV24Rf5+fzIrUbPpFNuediYPp3a6507FE6oxOU8TvlJVXPAbupaWpNDCGxy7txQR96Cn1gApd/Mr6jHwe/ngdG3YXcF7P1jwxOpZ2Wqdc6gkVuviFopIyXl6aylvfb6Nl00a8fl0cF8W21Z2eUq+o0MXnfZeWzd8+Wc+u3EOMHxzNwyN70DxYd3pK/aNCF5+VV1TCk4s28dHaDDqFh/BBwlCGdGrpdCwRx6jQxedYa/ls/R7+uWADeUWl/OnsLvzpHC1vK6JCF5+yN7+YR+Yn89Wm3+gb2ZxZE4fQq10zp2OJeAUVuvgEay1zV+/iqc82Uepy8beLezDxDx1p2CDA6WgiXkOFLl5vx76DPPzRen7auo+hnVrwzJi+xISHOB1LxOuo0MVrlbssb/+wjReWpBAYEMBTV/Rh/OAoTUUUqYIKXbxSelYhf5mXxNqdeZzTozX/viKWiOa6QUjkRFTo4lXKyl1M+24rr3yVRpOgBrx8TT8u799eZ+UiNaBCF6+RsvcAD81LYl1GPiNj2/LE6FhahTZyOpaIz1Chi+NKy11M/WYLk5al0axxIK9eG8clfSOcjiXic1To4qhNewp4aF4SyZkFjOobweOX9aZlU52Vi5wKFbo4orTcxevfbGHysjSaBwfy+nVxjOyjs3KR2lChi8el7D3AA/9NJDmzgEv7tePxy3rTIiTI6VgiPk+FLh5TVu7ijRVb+c9XaYQ2bsjU6+O4KFZn5SLuokIXj0jPKuSB/yaRtCuPS/pE8MRoXSsXcTcVutQpl8sy44dtPL84heCgBkweP4BL+7VzOpaIX1KhS53Zua+IB+clsWpbLuf1bMNTY2JpHdrY6VgifkuFLm5nrWXOql3867ONNDCG56/sy5UDI3W3p0gdU6GLW/1WUMxfP1rHNynZ/KFLS567sh/t9ZBmEY9QoYvbLEzazSPzkzlcVs7jl/XmhqEdCAjQWbmIp6jQpdbyi0r5x6fJLEjaTb+oMF6+uh+dWjV1OpZIvaNCl1r5Pi2HB/+bRE7hYe4/vxt3jeispwiJOESFLqekuLScZ77YzMwft9O5VQhvTvgDfSKbOx1LpF5ToctJS87M574PEknPKuSm02N4eGQPGgc2cDqWSL2nQpcaK3dZ3lixhZeXptIiJIh3bxnMmV1bOR1LRCpVW+jGmBnAKCDLWht7nO3XAX8FDHAAuNNam+TuoOKsjP1F3P9BEqu253JJnwj+fUUsYU20oJaIN6nJGfpMYAowq4rt24CzrLX7jTEjgWnAEPfEE28w/9dM/jE/GQu8eFU/xsTpkXAi3qjaQrfWrjDGxJxg+49HffkzEOmGXOIFCopL+cf8ZD5N3E18h9N4+Zr+RLVo4nQsEamCu6+h3wJ8UdVGY0wCkAAQHR3t5kOLO63enst9cxPZW1DMA+d3405NRxTxem4rdGPM2VQU+hlV7WOtnUbFJRni4+Otu44t7lNW7mLS12lMWZ5OVIsmzLtjGAOiT3M6lojUgFsK3RjTF5gOjLTW7nPH9xTP25VbxL1zf2XtzjzGxkXy+OjeNG2kiVAivqLW71ZjTDTwMXCDtTa19pHECZ8mZvL3T5IxBq1ZLuKjajJtcQ4wAgg3xmQAjwGBANbaqcCjQEvgtcqZD2XW2vi6CizuVXi4jEc/TebjtZnEdziNV8b1J/I0ffAp4otqMstlfDXbbwVudVsi8Zj1GfncPWctO3OLuPfcrtx9Thd98Cniw3SBtB5yuSxvfb+N5xZvJrxpI+YmDGNwxxZOxxKRWlKh1zM5hYd58L9JfJOSzYW92/Ds2L6641PET6jQ65Ef03O474NE8g6V8uTlsVw/JFp3fIr4ERV6PVBW7uI/lXPLO4WH8M7EwfSMaOZ0LBFxMxW6n9uTf4h75ySyansuVw2smFveJEh/7CL+SO9sP7Z8cxb3f5jI4TIXL1/TjysGaJkdEX+mQvdDpeUuXlicwhsrttIzohmvXjtAz/gUqQdU6H4mM+8Qd89ey9qdeVw3JJp/jOqlpwmJ1BMqdD+ybPNv3P9hEmXlVrfvi9RDKnQ/UFbu4oUlqUz9dgu9Iprx6nVxdAwPcTqWiHiYCt3H7c0v5p45v7Jqey7XDonmUV1iEam3VOg+7If0HO6Z8yuHSst55Zr+XD6gvdORRMRBKnQf5HJZXl2ezktfpdKlVVNevz6OLq1DnY4lIg5TofuY/QdL+POHiXyTks3l/dvx1Jg+ulFIRAAVuk9J2pXHXe+vJfvAYf51eSzXaS0WETmKCt0HWGuZvWonjy/YSKvQRvz3jmH0iwpzOpaIeBkVupcrLi3n758k89HaDIZ3a8V/runPaSFa7lZEfk+F7sV27ivijvd+YdPeAu49tyv3nNuVBgG6xCIix6dC91LLU7K4d86vAMy4cRBn92jtcCIR8XYqdC/jclmmLE/n5a9S6dm2GVOvH0h0Sz20WUSqp0L3IgXFpdz/QSJfbcpizID2/PuKPgQH6a5PEakZFbqXSPvtAAnv/sKu3CIev6w3E4Z10JREETkpKnQv8GXyHh74MIngoAbMvm0ogzu2cDqSiPggFbqDXC7LS0tTmbI8nf5RYUy9fiBtmzd2OpaI+CgVukMKiku5b24iyzZncXV8JE9eHkujhrpeLiKnToXugPSsQhJmrWFnbhFPju7N9UN1vVxEak+F7mHLNv/GvXMSCWoYwPu3DmFIp5ZORxIRP6FC9xBrLa9/u4XnF6fQK6IZ0ybE0z4s2OlYIuJHVOgeUFxazl/mrWNB0m4u7deO58b21fxyEXG7gOp2MMbMMMZkGWOSq9jewxjzkzHmsDHmQfdH9G1784u5aupPLFy3m4cu7M6kcf1V5iJSJ2pyhj4TmALMqmJ7LnAPcLl7IvmPxF153DZrDUWHy5h2Qzzn92rjdCQR8WPVnqFba1dQUdpVbc+y1q4GSt0ZzNd9mpjJ1W/8ROPAAD6+6w8qcxGpcx69hm6MSQASAKKjoz15aI85+mahwR1bMPX6gbTQ+uUi4gEeLXRr7TRgGkB8fLz15LE9oaikjPs/SOLLDXsZNyiKJ0bHEtSw2h+CRETcQrNc3GRvfjG3zlrNht0FPHJJT245o6NuFhIRj1Khu0FyZj63vLOawuIypk+I59yeul4uIp5XbaEbY+YAI4BwY0wG8BgQCGCtnWqMaQusAZoBLmPMfUAva21BXYX2Jos37OW+uYm0CAli3p2n0zOimdORRKSeqrbQrbXjq9m+F4h0WyIfYa1l+nfbeOqLTfSNDOPNCQNpHaqVEkXEObrkcgpKy108+ukG5qzaySV9Injx6n40DtTNQiLiLBX6STpQXMpd76/lu7Qc7hrRmQcv6E5AgD78FBHnqdBPwu68Q0ycuZr0rEKeHduHawb551x6EfFNKvQaSs7MZ+LM1RwqKWfmzYM5o2u405FERP4PFXoNLE/J4k/vr6V5cCDz7jyd7m1DnY4kIvI7KvRqzFm1k0fmJ9OjbSgzbhpEm2aaySIi3kmFXgVrK9ZkmbwsnRHdW/HqtXGENNJwiYj3UkMdR0mZi4c/XsfHazMZNyiKf10eS8MGWpNFRLybCv0YhYfLuPO9X/guLYf7z+/G3ed00ZosIuITVOhHyTpQzMSZq9m05wDPXdmXq+OjnI4kIlJjKvRK23IOMmHGSnIOlDD9xnjO7t7a6UgiIidFhQ4k7crj5pmrMcDchKH0iwpzOpKIyEmr94X+bWo2d773Cy2bBjFr4hA6hoc4HUlE5JTU60L/NDGTBz5MomubUN65eRCtNcdcRHxYvS30t3/YxuMLNzKkYwvevDGeZo0DnY4kIlIr9a7Qj75h6MLebfjPuAFa+lZE/EK9KvRyl+XRT5N5f+VOromP4t9X6IYhEfEf9abQS8pc3P9hIovW7eGOszrz14u664YhEfEr9aLQD5WUc+f7v/BNSjYPj+zBHWd1djqSiIjb+X2hFxSXcsvM1azZsZ+nx/Rh/GA9lEJE/JNfF/q+wsPc+PYqUvYeYPL4AYzq287pSCIidcZvC31vfjHXv7WSXblFTJugW/lFxP/5ZaHvyi3i2uk/k1tYwjsTBzO0U0unI4mI1Dm/K/St2YVcN30lRSXlvH/bUPprXRYRqSf8qtBT9h7guukrsdYyN2EoPSOaOR1JRMRj/KbQkzPzueGtlQQ1DOD9W4fRpXVTpyOJiHiUXxR64q48Jry1ktDGgcy+bQgdWmrFRBGpf3y+0Ndsz+Wmt1fTIiSI2bcNIfK0Jk5HEhFxhE8X+sqt+7h55mraNmvM7NuG0ra5lr8Vkfqr2pWpjDEzjDFZxpjkKrYbY8wkY0y6MWadMSbO/TF/76ct+7jp7dW0CwtmboLKXESkJksNzgQuOsH2kUDXyl8JwOu1j3ViP6TncPPMVUS1CGbObUP1YAoREWpQ6NbaFUDuCXYZDcyyFX4GwowxEe4KeKwf03OYOHM1MS1DmH3bUFqFNqqrQ4mI+BR3LAbeHth11NcZlb/3O8aYBGPMGmPMmuzs7FM6WOtmjRncsQXv3zqE8KYqcxGRIzz6dAdr7TRrbby1Nr5Vq1an9D26tG7Ku7cMoaXKXETk/3BHoWcCUUd9HVn5eyIi4kHuKPQFwITK2S5DgXxr7R43fF8RETkJ1c5DN8bMAUYA4caYDOAxIBDAWjsV+By4GEgHioCb6yqsiIhUrdpCt9aOr2a7Bf7otkQiInJK9Mh7ERE/oUIXEfETKnQRET+hQhcR8ROm4jNNBw5sTDaw4xT/93Agx41x3MVbc4H3ZlOuk6NcJ8cfc3Ww1h73zkzHCr02jDFrrLXxTuc4lrfmAu/NplwnR7lOTn3LpUsuIiJ+QoUuIuInfLXQpzkdoAremgu8N5tynRzlOjn1KpdPXkMXEZHf89UzdBEROYYKXUTET3h1oXvrA6prkGuEMSbfGJNY+etRD2SKMsYsN8ZsNMZsMMbce5x9PD5eNczlxHg1NsasMsYkVeZ6/Dj7NDLGfFA5XiuNMTFekusmY0z2UeN1a13nOurYDYwxvxpjFh1nm8fHq4a5nByv7caY9ZXHXXOc7e59T1prvfYXMByIA5Kr2H4x8AVggKHASi/JNQJY5OGxigDiKl+HAqlAL6fHq4a5nBgvAzStfB0IrASGHrPPXcDUytfjgA+8JNdNwBRPjtdRx74fmH28Py8nxquGuZwcr+1A+Am2u/U96dVn6NbLHlB9Erk8zlq7x1q7tvL1AWATv3+2q8fHq4a5PK5yDAorvwys/HXsDIHRwDuVr+cB5xpjjBfkcoQxJhK4BJhexS4eH68a5vJmbn1PenWh10CNH1DtgGGVPzZ/YYzp7ckDV/6oO4CKs7ujOTpeJ8gFDoxX5Y/piUAWsNRaW+V4WWvLgHygpRfkAhhb+SP6PGNM1HG214VXgL8Ariq2OzJeNcgFzowXVPxjvMQY84sxJuE42936nvT1QvdWa6lYb6EfMBmY76kDG2OaAh8B91lrCzx13OpUk8uR8bLWlltr+1PxHNzBxphYTxy3OjXItRCIsdb2BZby/8+K64wxZhSQZa39pa6PdTJqmMvj43WUM6y1ccBI4I/GmOF1eTBfL3SvfEC1tbbgyI/N1trPgUBjTHhdH9cYE0hFab5vrf34OLs4Ml7V5XJqvI46fh6wHLjomE3/O17GmIZAc2Cf07mstfustYcrv5wODPRAnD8AlxljtgNzgXOMMe8ds48T41VtLofG68ixMyv/mwV8Agw+Zhe3vid9vdC98gHVxpi2R64dGmMGUzHOdfoXu/J4bwGbrLUvVbGbx8erJrkcGq9WxpiwytfBwPnA5mN2WwDcWPn6SmCZrfwky8lcx1xjvYyKzyXqlLX2f6y1kdbaGCo+8Fxmrb3+mN08Pl41yeXEeFUeN8QYE3rkNXABcOzMOLe+J6t9pqiTjJc+oLoGua4E7jTGlAGHgHF1/RebijOVG4D1lddfAf4GRB+Vy4nxqkkuJ8YrAnjHGNOAin9APrTWLjLGPAGssdYuoOIfoneNMelUfAg+ro4z1TTXPcaYy4Cyylw3eSDXcXnBeNUkl1Pj1Qb4pPJcpSEw21r7pTHmDqib96Ru/RcR8RO+fslFREQqqdBFRPyECl1ExE+o0EVE/IQKXUTET6jQRUT8hApdRMRP/D8O+txYknoW3QAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "γ = 1.4\n",
    "α = np.linspace(1.0001,5,100)\n",
    "\n",
    "y = (α**γ - 1)/γ/(α-1)\n",
    "\n",
    "plt.plot(α,y);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "80fadbf7-c059-42c9-9148-bfc0ae6ea032",
   "metadata": {},
   "source": [
    "Carnot efficiency"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "752bbe8b-f2dd-4725-9c94-00765e299220",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZsAAAELCAYAAAAP/iu7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAo/UlEQVR4nO3deXxcdb3/8dcnW7M0S9O9abqX0o1SSGsVkLAooNdWRIWrInBlEa+KV72KP7xeLou7cFXuFYo7oqgoglyVCiXsBcrWNqUrSZsuaZtm37fv749z0o5p0kzSmTmZzPv5eMzjZE7OOfP5Tqbz7lm+32POOURERKIpKegCRERk5FPYiIhI1ClsREQk6hQ2IiISdQobERGJOoWNiIhEXSBhY2ZTzeyHZvaCmTWbmTOzGWGum2RmXzGzcjNrNbM3zOySKJcsIiInIKg9mznAh4Ea4JlBrnsrcDNwF3ARsA74vZm9J5IFiohI5FgQnTrNLMk51+3/fDVwLzDTOVc+wHoTgArgm865/wyZ/wQw3jl3SvSqFhGRoQpkz6YnaIbgAiAN+FWv+b8CFpvZzBMqTEREoiLeLhBYCLQBO3rNL/WnC2JbjoiIhCMl6AIGKR+odcce+6sO+f0xzOxa4FqAjIyM0wsLC6NWYHd3N0lJ8ZbhQ6f2jmxq78i2bdu2Kufc+Fi8VryFzZA451YDqwGKiorc+vXro/ZaJSUlFBcXR237w43aO7KpvSObme2K1WvFW4TXAHlmZr3m9+zRVCMiIsNOvIVNKTAKmN1rfs+5ms2xLUdERMIRb2HzN6AD+Giv+R8DNjnnymJfkoiIDCSwczZm9kH/x9P96UVmdgg45Jx7yl+mE/iFc+4TAM65g2Z2B/AVM2sAXgUuBc4FVsa0ASIiErYgLxD4fa/n/+tPnwKK/Z+T/Ueom4BG4AZgErAV+LBz7tHolCkiIicqsLBxzvU+yR/WMs65LuA2/yEiInEg3s7ZiIhIHFLYiIhI1ClsREQk6hQ2IiISdQobERGJOoWNiIhEncJGRESiTmEjIiJRp7AREZGoU9iIiEjUKWxERCTqEuJOnSIiclRXt6OxrTOmr6mwERGJI845Gto6qW/poKG10390HJnW+/PqQ+b947Qz5kEDChsRkZjr7j4aGHUDPHovU9/SQbc7/vbTkpPIyUghOz2V7PQUstNTGD96tP/z0XnXfCs27QWFjYjIkHV3OxpaO6lubqemuZ2apnZqmjuoaWqnurndC4jmYwOkofX4gZGSZORmpJKbkUpORipjMtOYMTbryDxvfgo56an/EB49P6en9r4NWN+uidD7EA6FjYgIXnA0tjt2Hmqktrmd6qaOYwKkxg+V6qZ2apu93/cXGilJRl6mFxa5GamMHZ3GrPG9AyP1H573PDLTkjEb8JZfcUVhIyIjVnN7J1UN7VQ1tXG4sZ2qxjYON7ZRdeTndg43ec9re4Jj7VPHbCctOYm8zFTys9LIy0xl3qRs8jLTyM9MY0xWGmMyU/2pNy8vK5XsUSkjLjBOhMJGROKGc95VVAfq2zjY0MqhhpDAaDgaHD1B0tLR1ed2stNTGDd6FGOz0pg5LotlM/LJz0rj8L7dLF+y4EiwjPHDJGsE7mnEmsJGRALnnKO+pZODDa0cbGjjQL03PVjfxoGGVg754XKgvq3PAElJMvKz0rwAGe0FyLjRaYz1A2Xc6FFHfpefldbvOY2Skv0ULy2IdnMTksJGRKKqs6ubAw1t7K9tYX9dK/vrWthX23okPA42tHKwvo22zu5j1s1KS2ZiTjrjs0exeGoe52ePYkLOKCZkp/tTL0Ry0lNJStKex3CmsBGRIevqdhxsaPVCpPZokFTWe9P9dS0camg75iR6VloyE3PTmZidzmnTxjAxJ50J2aOY0DP1fx49Sl9RI4X+kiLSr65uR0V1MxU1zeypbmFPTTN7alqoqGlmb00LBxra6OqVJBmpyUzOS2dKbgbvnDueybnpTM7L8Ka5GUzOSycnPTWgFklQFDYiCayr23GgvpWK6qMhsqem5cjz/XUtdK958sjySQaTczOYOiaDFbPHMsUPjym5GUzK9aY5GboKS46lsBEZ4bq6HftqWyg/3ER5VRPlh5v9aRMV1S20dx09V2IGE7PTKczPYPnMfDrrD3LWqfOZOiaDwvxMJuWmk5qs8Xtl8BQ2IiNAd7djb20LZX6IlFc1e9PDTVRUN9PRdfRQV0ZqMtPHZjJ3QjbvWjCJafmZFOZnMHVMJlPy0hmVcvRKrZKSEoqXFQbRJBlhFDYicaSts4vyqmZ2HGz0Hoe86VuHGv/haq6eQJk3MZt3L5jEzHGZTB+bxcxxWUzIHqXDXBJzChuRYai5vZNtBxrZfqCBHYca2emHy+7q5iNXdplBQV4GcyaM5ozZY5k9YTSzxmUxQ4Eiw5DCRiRA3d2OXdXNbK2s5839DWyprGdrZQO7qptxfqikJhszx2WxYEoOK5dMYfaE0cwe7z0y0sIbcFEkaIGEjZkVAncC7wIMeBz4nHNudxjrTgNuBc4BxgMVwO+AbzjnmqJWtMgJqm/tYPO+et7c7wXKm5UNbKtsONIj3gxmjvVC5eKlU5k3KZuTJo5mWn4mKTopL3Eu5mFjZpnAWqANuAJwwG3Ak2Z2yvECw8yy8IIpFfgPYDewDPgvYC5waXSrFwlPXUsHpXvr2Og/Nu2to/xw85Hfj8lM5eRJOVy2vJD5k3L8YMnWnoqMWEHs2VwDzALmOed2AJjZBmA7cB1wx3HWPQMvVC5wzq3x5z1pZvnAF80s0znX3P/qIpFX39rBpj11bAgJll0hwTIlN51FBblcctpUFk3NZeHkHMbrnIokmCDCZiWwridoAJxzZWb2HLCK44dNmj+t7zW/FkjCOyQnEjXdzrH9QAOv7a7l1d01vLq7hu0HG4+cXynIy2BxQS4fLipkUUEui6bkMHb0qGCLFhkGggibhcDDfcwvBT40wLqP4+0BfcvMrsc7jLYcuAG4W+dsJNLqWjp4vaKWV3fV8FpFLevfaqb5sacByM1IZem0PN67eAqnTstjcUEu+VlpA2xRJDEFETb5QE0f86uBMcdb0TnXamZnAn/AC6cePwY+3d96ZnYtcC3AxIkTKSkpGWTJ4WtsbIzq9oebkdbehnbH1uouttZ0sbW6m4qGbhzeLvPU7CSWjnOcPG4Uc/KSmZhlJFkz0Izbt5cN+wIuPgpG2t93IInW3liKq0ufzSwd+C0wAbico3s2XwM6gev7Ws85txpYDVBUVOSKi4ujVmNJSQnR3P5wE+/tPdjQyktl1bz4VjUvlh1m2wFv5zg9NYnTpo3hA2/LZ9mMfJYU5jF6VErct3ew1F6JlCDCpoa+92D62+MJ9QmgGJjjnNvpz3vazOqA1WZ2t3PujYhVKiNOU1sn6946zDPbq3hm+yF2HvLCJSstmdNn5LPq1ALeNjOfU6bmkZaiy41FIiWIsCnFO2/T2wJg8wDrLgZqQoKmx0v+dD6gsJEjurodG/fW8cy2Qzyzo4pXd9XQ2e1IT01i+cyxfLiokLfNGsuiKTnqyyISRUGEzSPAd81slnPuLQAzm4F3WfONA6xbCYwxszmhV7MBb/OneyNdrMSf6qZ2ntxykLVbDvLsjirqWjoAWDglh6vPmsVZc8dx+vQx/d4aWEQiL4iwuRfvZP7DZvZVvE6dt+KNBHBPz0JmNh3YCdzinLvFn/1z4PPAX8zsdrxzNkV4HTxfAZ6LURtkGHHOsf1gI4+/eYC1bx7k1d01dDsYnz2Kdy2YyFlzx3HmnHG6BFkkQDEPG+dck5mdizdczX14F/o8gTdcTWPIogYk4/Wf6Vm33MxWADfjjTowDi+kVgO3O+eOvYm5jEidXd28VFbNms0HeGLLASqqWwBYVJDDp8+dy/nzJ7BoSq7uSy8yTARyNZo/BtolAyxTTh+dNJ1zm4EPR6cyGc46u7p54a3D/GVjJWtKKznc1M6olCTOnDOOT549m/NOnsik3PSgyxSRPsTVpc+SeDq6unl+52H+smE/azZXUtPcQWZaMufNn8h7Fk2ieN4EjScmEgcUNjLsOOd4vaKWh17by5/f2EdNcwejR6Vw3vwJvGfxZM4+abxO7ovEGYWNDBsV1c089NpeHnptL2VVTYxKSeJdCyay6tQCzpo7TgEjEscUNhKoprZOHt2wjwdf2cPL5V6f3hWz8rn+7NlcuHgSOempAVcoIpGgsJFAlO6r4zcv7eZPr+2jsa2T2eOz+PcL5vH+pQUU5GUEXZ6IRJjCRmKmZy/m1y/u5o09dYxKSeK9p0zmI8uncfr0Mbq/i8gIprCRqKuobuaXL5TzwMsVNLR2ctLE0dz8vgVcvHQquZk6TCaSCBQ2EhXOOdbvquGnz5bxWGklZsZ7Fk/mirdP116MSAJS2EhEdXZ1838b93PvM2+xaW89uRmpXHf2bC5fMZ0pOhcjkrAUNhIRbZ1d/PHVvdz91E52HW5m9vgsbr94ERcvLSAzTR8zkUSnbwE5IW2djp88W8a9T79FZX0riwtyuftjp/PuBRM1LpmIHKGwkSFp7ejivhd28YOnmmno2Mzymfl8+4OncNbccTofIyLHUNjIoHR0dfPblyv44drtHKhvY+HYJG7+0NtYNiM/6NJEZBhT2EhYurodj7yxlzv/vp3d1c2cPn0M/33pUtoqNipoRGRAChsZ0Lq3DnPro5sp3VfPgsk5/OzKZRTPG4+ZUVIRdHUiEg8UNtKv3Yeb+cZf3+SvmyqZkpvO9y87lfedMkUn/kVk0BQ2coymtk5+uHYHP322jOQk4wvvOolr3jlLoy6LyJApbOQfrCmt5OZHStlX18oHTivgSxecrLtfisgJU9gIAPtqW7j5kVLWbD7AvInZ/OEjSzl9uk78i0hkKGwSXFe34+fPl3PHmq10OceNF53MJ86cSWpyUtClicgIorBJYOVVTXzx92+wflcN58wbzy2rFlGYnxl0WSIyAilsElB3t+P+F3fx9b9sISXZuPPSJbz/1AL1/BeRqFHYJJh9tS186cENPLujineeNJ5vXbKYybkajVlEokthk0DWlFby7w9uoKOrm9svXsRHlk/T3oyIxITCJgG0dXbxzb9u4WfPlbOoIIe7/vk0ZozLCrosEUkgCpsRrryqiU//5lU27a3nqjNmcONFJzMqRZ0zRSS2FDYj2JNbDvLZ37yGGdxz+elcsHBS0CWJSIIKpDOFmRWa2YNmVmdm9Wb2RzObNoj155vZ782sysxazGyrmd0QzZrjiXOO/3lyB//yi5eZNjaT//vsWQoaEQlUzPdszCwTWAu0AVcADrgNeNLMTnHONQ2wfpG/fglwNVAHzAVGR7HsuNHc3smXHtzAoxv2874lU/j2JaeQkabDZiISrCAOo10DzALmOed2AJjZBmA7cB1wR38rmlkS8EvgCefcxSG/ejJ65caP/XUtfOLn63mzsp4bLzqZ6945S1ebiciwEETYrATW9QQNgHOuzMyeA1ZxnLABioH5eKEkIbZU1nPlT1+msa2Tn165jHPmTQi6JBGRI4I4Z7MQ2NTH/FJgwQDrnulP081snZl1mNlBM/uBmSVsz8Tnd1TxoR+9gMPxu+verqARkWEniD2bfKCmj/nVwJgB1p3iT38L3AXcCBQBtwCFwMV9rWRm1wLXAkycOJGSkpJBFx2uxsbGqG6/txf2dfLjjW1MyjI+vzSVg9te5eC2mL18zNsbNLV3ZEu09sZSvF363LMn9ivn3Nf8n0vMLBn4ppnNd8692Xsl59xqYDVAUVGRKy4ujlqBJSUlRHP7oe57oZx7NpSyYlY+91xeRG5GakxeN1Qs2zscqL0jW6K1N5aCOIxWQ997MP3t8YQ67E//3mv+Gn+69ATqiis/fuYt/uPhUs6fP4GfX7U8kKAREQlXEHs2pXjnbXpbAGwOY93j6R5SRXHmrrXb+e6abbx38WTuvPRU0lJ07xkRGd6C+JZ6BFhhZrN6ZpjZDOAM/3fH81e8/jkX9Jp/oT9dH6EahyXnHN9bs5XvrtnGB5YW8P3LFDQiEh+C+Ka6FygHHjazVWa2EngYqADu6VnIzKabWaeZ9ZybwTl3GPgG8Ekz+7qZnW9mNwJfA34Rejn1SHTX2h38cO0OLltWyHc/tIQU3U1TROJEzA+jOeeazOxc4E7gPsCAJ4DPOecaQxY1IJljA/EWoAH4FPBFYD/wHeDWKJceqJ88W8b3/r6ND5xWwNcvXkxSkjprikj8CORqNOfcbuCSAZYpxwuc3vMdXsfP43X+HFEeeGk3tz66mYsWTeLbl5yioBGRuKPjMMPcoxv28ZWHNnL2SeP5/mVLdehMROKSvrmGsZfKqvn8b9+gaPoY7v7Y6boYQETilr69hqmdhxq59r71TM3P4N6PF2nkZhGJawqbYaiqsY2rfvYyyWb8/Mrl5GWmBV2SiMgJibfhaka81o4urv7Feg42tPKba1YwbWxm0CWJiJwwhc0w4pzjpoc28XpFLXd/7DSWThtoXFIRkfigw2jDyH3rdvGHV/dww3lzuXDR5KDLERGJmCGFjZk9bWaTQ57r7PUJeqmsmlv+vJnz50/ghvPmBl2OiEhEDXXPZoxzbn/I8zwzeyISBSWiyrpWPnX/K0zLz+SOS09Vp00RGXGGGjYdZnZkXX/MsrGRKSmxdHU7bnjgNZrbu1j98dPJSdetAkRk5Blq2DwGfNfMDMDMUoCEvS3zibhr7Q5eLKvm1lWLmDMhO+hyRESiYqhXo92Cd2vm183sWbyblq2NWFUJ4sW3DvP9J7zbBVxy+tSgyxERiZohhY1zrgVYaWbvBIqAp4HfR7Kwka6mqZ3P/fZ1puVncsv7FwVdjohIVIUVNma2DJgBHAI2+udocM49jRc0MgjOOW7600aqGtv44/VnMHqUujuJyMg24Lecmd0O3MjR4f6dmVUAr/iP9cCrzrmqqFU5wjy6YT9/2VjJly6cx+KpuUGXIyISdeH8l/pK4A/ATUA6cApwGnA68GUgGy+A9gCvOOc+EJ1SR4ZDDW187eFNLCnM49qzZg28gojICBBO2GQDdzvntvvPNwL39/zSzE7CC58ifyr9cM7x1T9tpKm9i+996BTdm0ZEEkY4YfMkMK2/XzrntgHbgAciVdRI9cgb+3is9ABfuehkXeYsIgklnP9a/z/gk2aWF+VaRrS65g5ufXQzSwrzuFqHz0QkwYQTNi8Dc4ENZnajmS0zs1FRrmvE+fZjW6huauf29y8iWcPRiEiCCSdsvo53KK0duB1YBzSY2UYz+6WZ/Vs0CxwJXttdw69f2s2V75jJogJdfSYiiWfAczbOudt6fjazbOBUvBEDeqYfBu6MTnnxr7Orm5se2sSE7FF8/t0nBV2OiEggBtWb0DnXADzjPwAwM40ceRy/WreLzfvr+d+PnqbOmyKSsE742lvnXEckChmJapvbufPx7ZwxZywXLZoUdDkiIoFRR48o+sETO2ho7eCr712AP0C2iEhCUthEyVuHGvnlC+VcuqyQ+ZNzgi5HRCRQCpso+eZftzAqJYl/e5cuChARUdhEwbq3DrNm8wE+dc4cJmSnB12OiEjgAgkbMys0swfNrM7M6s3sj2bW75A4x9nOjWbm/Bu4DQvOOb772FYm5aTziTNnBl2OiMiwEPOwMbNMvLt6ngxcAVyON0LBk2aWNYjtzAK+ChyMRp1DtbGqi/W7avjXc+eQnpocdDkiIsNCEB0/rgFmAfOcczsAzGwDsB24DrgjzO38CG/06XkE045jOOd4aHsHBXkZXFpUGHQ5IiLDRhCH0VYC63qCBsA5VwY8B6wKZwNm9hG82xl8JSoVDtHjbx6krL6bG86bS1qKToeJiPQI4htxIbCpj/mlwIKBVjazMXjD43zJOVcd4dqGrLvbccfftzEx0/jAaQVBlyMiMqwEcfgpH6jpY341MCaM9b+Dd/+cn4f7gmZ2LXAtwMSJEykpKQl31bCtr+zkzf1tfPwkx7PPPB3x7Q9XjY2NUXk/hyu1d2RLtPbG0rA41xEuMzsL+DhwmnPOhbuec241sBqgqKjIFRcXR7Qu5xx3/s9zTB+bTPFMiPT2h7OSkhK1dwRTeyVSgjiMVkPfezD97fGEugf4CbDHzPL8G7qlAMn+80Dus/PCzsO8saeO6945myQNSyMicowgwqYU77xNbwuAzQOsOx/4JF4o9TzOAFb4P18fuTLD96OndjJu9CidqxER6UcQh9EeAb5rZrOcc28BmNkMvNC4cYB1z+lj3n8DycBngB19/D6qNu6p45ntVXz5wpPVr0ZEpB9BhM29wKeBh83sq4ADbgUq8A6TAWBm04GdwC3OuVsAnHMlvTdmZrVASl+/i4W7n95J9qgUPrpi0AMgiIgkjJgfRnPONQHn4l1Rdh9ex8wy4FznXGPIooa3xzJsO6zsqWnmrxv385EV08hJ1z3kRET6E8jVaM653cAlAyxTjhc4A22rODJVDd5963YB8PG3zwiqBBGRuDBs9xqGu9aOLn77cgXvXjCJgryMoMsRERnWFDZD9PDre6lt7uDKM2YEXYqIyLCnsBkC5xw/e66ckydl87aZ+UGXIyIy7ClshuClsmq2VDZwxTtmYOrEKSIyIIXNENz/4m5y0lN4/6nqxCkiEg6FzSDVNrfzt9JKLl5aQEaaOnGKiIRDYTNIf3ptL+2d3Vy6TJ04RUTCpbAZBOccD7xcweKCXBZMyQm6HBGRuKGwGYSNe+vYUtnApct0y2cRkcFQ2AzCAy9XkJ6axMpTpwRdiohIXFHYhKmlvYs/v76P9yyerHHQREQGSWETpsffPEBDWycfPH1q0KWIiMQdhU2YHn59H5Ny0nnbzLFBlyIiEncUNmGobW7nqW0Hed+SySQnacQAEZHBUtiE4a+bKunocqzSiAEiIkOisAnDw6/vZdb4LBaqb42IyJAobAawv66FF8uqWbWkQINuiogMkcJmAI++sR/nUN8aEZEToLAZwKMb9nHK1FxmjssKuhQRkbilsDmOfbUtvLGnjgsXTQq6FBGRuKawOY41pZUAXLhQYSMiciIUNsfxt9JK5k4Yzazxo4MuRUQkrils+lHd1M5LZdU6hCYiEgEKm348vvkA3Q4u0CE0EZETprDpx2OllRTkZagjp4hIBChs+tDc3skzO6q4YOEkdeQUEYkAhU0fnt9xmPbObs6fPyHoUkRERoRAwsbMCs3sQTOrM7N6M/ujmU0LY70iM1ttZlvMrNnMdpvZ/WY2M5L1Pbn1IFlpyRTNyI/kZkVEElbMw8bMMoG1wMnAFcDlwFzgSTMbqJv+ZcBC4AfARcCNwGnAejMrjER9zjlKth7izLnjSEvRjp+ISCSkBPCa1wCzgHnOuR0AZrYB2A5cB9xxnHW/5Zw7FDrDzJ4Dyvztfu1Ei9t+sJG9tS185tw5J7opERHxBfFf95XAup6gAXDOlQHPAauOt2LvoPHn7QIOARG52cyTWw4CUDxP52tERCIliLBZCGzqY34psGCwGzOz+cAE4M0TrAvwztfMn5zDpNz0SGxOREQIJmzygZo+5lcDYwazITNLAe7G27P5yYkWVt/awfryGs6ZN/5ENyUiIiGCOGcTSXcB7wDe65zrK8AAMLNrgWsBJk6cSElJSZ/LvXKgk85uR27zXkpKKodUUGNjY7/bH4nU3pFN7ZVICSJsauh7D6a/PZ4+mdk38QLkCufcmuMt65xbDawGKCoqcsXFxX0u98SfNpGZtoerVp4z5CvRSkpK6G/7I5HaO7KpvRIpQYRNKd55m94WAJvD2YCZ3QR8GfiMc+6+SBX2/M4qls/M1yXPIiIRFsS36iPACjOb1TPDzGYAZ/i/Oy4z+yxwG3CTc+6uSBVVWdfKzkNNnDF7XKQ2KSIiviDC5l6gHHjYzFaZ2UrgYaACuKdnITObbmadZva1kHmXAf8N/A1Ya2YrQh6DvpIt1PM7qwB4++yxJ7IZERHpQ8wPoznnmszsXOBO4D7AgCeAzznnGkMWNSCZfwzEC/35F/qPUE8BxUOt67kdhxmTmcqCyRrlWUQk0gK5Gs05txu4ZIBlyvGCJXTelcCVUaiH53dW8fbZY0lK0ijPIiKRpjPhQFlVE/vrWnmHzteIiESFwgZ4budhAM6Yo7AREYkGhQ3wUlk1k3LSmTE2M+hSRERGpIQPG+ccL5dVs2xmvu7KKSISJQkfNntrW6isb2XZjEENyyYiIoOQ8GGzvtwbIadouu7KKSISLQkfNi+XV5M9KoV5k7KDLkVEZMRK+LBZX17DadPHkKz+NSIiUZPQYVPX3MHWAw06XyMiEmUJHTav7K4GoGiGzteIiERTQofNy+U1pCYbS6bmBV2KiMiIltBhs768mkUFuWSkJQddiojIiJawYdPR1c3GvXUsLdT5GhGRaEvYsNl2oIHWjm6WFOYGXYqIyIiXsGHzRkUdAKcW5gVbiIhIAkjgsKklLzOVafkafFNEJNoSN2z21LJkap4G3xQRiYGEDJvm9k62HWhgyVSdrxERiYWEDJtNe+vpdrBE52tERGIiIcNmw55aAE5RZ04RkZhIyLB5vaKWgrwMxmePCroUEZGEkJBh88aeWvWvERGJoYQLm85uR0V1iw6hiYjEUMKFTWtHFwCLC7RnIyISKwkXNi1+2MyfnBNwJSIiiSPhwqa1vYvJuenkZ6UFXYqISMJIuLBp6ehigfZqRERiKuHCpq2zmwVTFDYiIrEUSNiYWaGZPWhmdWZWb2Z/NLNpYa6bbmbfMbP9ZtZiZi+Y2TsH8/rasxERia2Yh42ZZQJrgZOBK4DLgbnAk2aWFcYmfgJcA3wN+CdgP/CYmZ0abg3asxERia2UAF7zGmAWMM85twPAzDYA24HrgDv6W9HMlgAfAf7FOfczf95TQClwC7ByoBdPMqNwjG4rICISS0EcRlsJrOsJGgDnXBnwHLAqjHU7gN+GrNsJPABcYGYDjj+TkZpEUpJuKyAiEktBhM1CYFMf80uBBWGsW+aca+5j3TRgzkAvnp6aHE6NIiISQUEcRssHavqYXw2MOYF1e35/DDO7FrjWf9poZlvDqHOoxgFVUdz+cKP2jmxq78g2L1YvFETYxJxzbjWwOhavZWbrnXNFsXit4UDtHdnU3pHNzNbH6rWCOIxWQ997MP3ttYS7LhzdwxERkWEkiLApxTv30tsCYHMY6870L5/uvW47sOPYVUREJGhBhM0jwAozm9Uzw8xmAGf4vzuePwOpwIdC1k0BLgXWOOfaIl7t4MXkcN0wovaObGrvyBaz9ppzLlav5b2g13HzDaAF+CrggFuBbOAU51yjv9x0YCdwi3PulpD1HwAuAP4dKAOux+vc+Q7n3KsxbIqIiIQp5ns2zrkm4FxgG3AfcD9eaJzbEzQ+A5L7qPEq4GfAbcD/AYXAhQoaEZHhK+Z7NiIikngSbtTncJlZsZm5Ph61vZYbY2Y/NrMqM2sys8fNbHEf2zvhAUQjxcymmtkP/Rqa/XbN6GO5sGo2syQz+4qZlZtZq5m9YWaX9PPa15jZFjNrM7OtZvbJKDSx92uG296+/t6u97h7cdDeD5rZH8xsl/9322pm3zCz7F7LRfSzO5j3JZLCaa+ZzTjO3zcvztp7gZmtNbNK/3O1x8x+Z2YLei0X1oDHkf4c9Ms5p0cfD6AY73zSZ4AVIY+ikGUMeBbYA/wzcCHwFF6nsKm9tnc/UIs3Ntx5wB/xzludGlDbDgB/AR7z2zmjj+XCqhm4HWgDvgicA9wDdAPv6bXcNf782/3lbvOfXz9M2uvwDtGu6PXIjLP2rgN+B3wUOBv4nP93XAckReuzG+77ElB7Z/h/36/38fdNjrP2/jPwHeCDfnsvx7tStx6Y7i+TiTfe5Cbg/XhDgW3EOw+eFbKtmH2HRe0NifcHR8Pm/OMss8pf5pyQebl4/X1+EDJvib/cVSHzUoCtwCMBtC0p5Oer6ePLN9yagQn+P7j/6rX+E8CGXuseBH7Ra7mf+h/s1CDb6//OAbcNsK14aO/4PuZ93G/fudH47Ib7vgTY3hn+86sH2Nawb28/dc/z6/6C//wGoAuYE7LMTKAT+HzIvJh9h+kw2olZCexzzj3ZM8M5V4d3ifaqXsud0ACikeSc6w5jsXBrvgBvXLpf9Vr/V8BiM5vpP387ML6P5e4DxgJnDqYNgxFme8MVD+091Mfsl/1pgT+N9Gc33Pcl4sJsb7iGfXv7cdifdvrTcAc8jtl3mMJmYPebWZeZHTazX/c65nm8QUWnmdnokOXK3AkMIBqAcGteiPc/vN4dakv96YKQ5eDY96v3ckG73j8O3uwfFz+r1+/jtb1n+9M3/WmkP7vhvi+x0ru9Pb5hZp3+eYxH+jg3ETftNbNkM0szs7l4h/Aqgd+E1BfOgMcx+w5T2PSvDvge3mGXc/H6Ap0PvGBmE/xlBhoYdEyYy/U5gGjAwq05H6h1/n71AMvRxzaH03vwK+BTeH/na/H2QNaaWXHIMnHXXjMrwLvf0+POuZ6xsCL92Q33fYm6ftrbhveFfB3e+ZUvAouB581sfsjq8dTeF/HatQ04Be+Q4cGQ1++vHaFDfsXsOywhBuIcCufca8BrIbOeMrOngZeAz+J1SJURxDl3ecjTZ8zsYbz/9d1GFA97RZP/P9OH8Q6vXBVwOVHXX3udc/uB0CsBnzGzv+H9z/wm4GOxrDNCLgdy8G5G+UXg72Z2pnOuPNCq+qE9m0FwXsfRbcAyf1YNxx8YtCbM5YbjAKLh1lwD5JlZ7zvS9bUcfWxz2L4HzrkGvI7Dy0Jmx017zSwD79j7LOAC59yekF9H+rMb7vsSNQO09xjOuQq8K7F6/33jor3OuTedcy86536Dd3XYaODGkPrCGfA4Zt9hCpuh6dl1Pt6gorvd0RER4nEA0XBrLgVGAbP7WA6ODq7acyy79/vVe7nhKPRQSVy018xSgQeBIrxLcTf2WiTSn91w35eoCKO9x9P77zvs29ubc64Wr7aecyfhDngcu++wWF+iF88PvA9yF954beBdv+6As0OWycG7MuSHIfOW+stdETIvBe/k5Z8DblN/lz6HVTPeJaDtwH/2Wv9xYGPI81TgEPCzXsv92H+/0oJsbz/L5gC7gafjqb14/4n8HV4fiPP6WSain91w35eg2tvPetPw+qb8Mp7a209bJgKNwD3+88/hHUqcFbLMDLwryr4Qrc/BcWuM9ZsSLw+8Dky3AR/Au0DgC3j9I3YD4/xlkoDngQrgMrzLIUvwdikLe23vAbxd0avxdnkfBFqB0wJq3wf9x4/8D9H1/vOzB1sz8E1//ufx+if9CK9z2z/1Wu6T/vzb/OVu8Z//a9DtxTvmfS/wEb+2K/A6wbUDZ8VTe0PaeBvHdmCcGq3PbrjvS0Dt/R5wJ/BhvAsEPgnswuukOC/O2vsQ8B94lyafg3fRwxa/LSf5y2Th7W1s9JdbiTcA8lvA6JBtxew7LKr/wOP5AXwF2IB3VVqH/8dYDUzutVw+Xke9aqAZr1PXkj62lwHcgXd5YivelSTFAbbP9fMoGWzNeAOmftX/x9vmv28f7Od1r8M779WG18P5U8OhvcD78PogVPl/78N4t7xYHm/tBcqP096bo/XZHcz7Euv2Av+C1/emxv/7VgK/plfQxEl7vwy8ghcuzXgdK+/h2KMT04A/4O29NQB/6r1MND4H/T00EKeIiESdLhAQEZGoU9iIiEjUKWxERCTqFDYiIhJ1ChsREYk6hY2IiESdwkZERKJOYSMiIlGnWwyIRIGZhdNbepdzbka0axEZDhQ2ItHx9l7PH8Ibm+rmkHltMatGJGAKG5EocM6tC31uZm1AVe/5IolC52xERCTqFDYiIhJ1ChsREYk6hY1IgMws28zqzSw56FpEoklhIxKs04FXnXNdQRciEk0KG5FgLQf2mtnvzWyzmb1kZvlBFyUSaQobkWAtw7st71XOuQV4t/A9P9iSRCJP/WxEgrUcuMg51+g/TwUOBViPSFRoz0YkIGY2Ach0zm32nycBpwKvBlmXSDRoz0YkBvoZA2058HLI85OBA865upgUJRJD2rMRCc4yYH3I86Jez0VGDHMunMFpRUREhk57NiIiEnUKGxERiTqFjYiIRJ3CRkREok5hIyIiUaewERGRqFPYiIhI1ClsREQk6v4/mxSjff4cHcMAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "Tc = 300\n",
    "Th = np.linspace(300,3000,1000)\n",
    "η = 1-Tc/Th\n",
    "\n",
    "plt.rc('font', size=16)\n",
    "plt.plot(Th, η)\n",
    "plt.grid()\n",
    "plt.xlabel(r'T$_h$')\n",
    "plt.ylabel(r'$\\eta_c$')\n",
    "plt.xlim([300,3000])\n",
    "plt.ylim([0,1.0]);\n",
    "\n",
    "plt.tight_layout()\n",
    "plt.savefig('carnot_eff.png')\n",
    "plt.savefig('carnot_eff.pdf')\n",
    "# use convert -density 300 -quality 100 carnot_eff.pdf carnot_eff.png"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "9e8f4abd-6e91-489b-be49-3846d41db731",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "1.5034722222222223"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "433/288\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "7842777e-aa29-4bb2-8c2b-d826df7e5c1b",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "666.6666666666666"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "1000/1.5\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "d29cc830-4dad-413f-8db0-5503487131d6",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.9.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
