diff --git a/doc/source/models/cubics.ipynb b/doc/source/models/cubics.ipynb
index feb63f18ff32f48217477faaac057576721fa39f..3fa6f80057b33e1289ce1b1f7d300dabeff173e3 100644
--- a/doc/source/models/cubics.ipynb
+++ b/doc/source/models/cubics.ipynb
@@ -28,7 +28,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 2,
+   "execution_count": null,
    "id": "9a884dab",
    "metadata": {},
    "outputs": [],
@@ -59,21 +59,10 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 3,
+   "execution_count": null,
    "id": "a8c87890",
    "metadata": {},
-   "outputs": [
-    {
-     "data": {
-      "text/plain": [
-       "<teqp.teqp.GenericCubic at 0x18f604c3a70>"
-      ]
-     },
-     "execution_count": 3,
-     "metadata": {},
-     "output_type": "execute_result"
-    }
-   ],
+   "outputs": [],
    "source": [
     "k_12 = 0.01\n",
     "kmat = [[0,k_12,0],[k_12,0,0],[0,0,0]]\n",
@@ -95,24 +84,34 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 7,
+   "execution_count": null,
    "id": "1b02ef72",
    "metadata": {},
-   "outputs": [
-    {
-     "data": {
-      "text/plain": [
-       "(30846.392909514052, 42.480231719002326)"
-      ]
-     },
-     "execution_count": 7,
-     "metadata": {},
-     "output_type": "execute_result"
-    }
-   ],
+   "outputs": [],
    "source": [
     "teqp.canonical_PR([Tc_K[0]], [pc_Pa[0]], [acentric[0]]).superanc_rhoLV(100)"
    ]
+  },
+  {
+   "cell_type": "markdown",
+   "id": "da00ff20",
+   "metadata": {},
+   "source": [
+    "## a and b\n",
+    "\n",
+    "For the cubic EOS, it can be useful to obtain the a and b parameters directly. The b parameter is particularly useful because 1/b is the maximum allowed density in the EOS"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "b5f6aaea",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "import numpy as np\n",
+    "modelPR.get_a(140, np.array([0.3, 0.4, 0.3])), modelPR.get_b(140, np.array([0.3, 0.4, 0.3]))"
+   ]
   }
  ],
  "metadata": {