fix nb

Author: dbrakenhoff

docs/steady/04benchmarks/test_polygon_areasink.ipynb

@@ -83,7 +83,7 @@
     "d2hdx2 = (ml.head(x + d, y) - 2 * ml.head(x, y) + ml.head(x - d, y)) / (d**2)\n",
     "d2hdy2 = (ml.head(x, y + d) - 2 * ml.head(x, y) + ml.head(x, y - d)) / (d**2)\n",
     "d2hdx2 + d2hdy2\n",
-    "aqin = ml.aq.inhomlist[0]\n",
+    "aqin = ml.aq.inhomdict[p1.name]\n",
     "print(\"recharge from numerical derivative: \", np.sum(aqin.T * (d2hdx2 + d2hdy2)))\n",
     "h = ml.head(x, y)\n",
     "print(\"leakage from aq0 to aq1 from head difference: \", (h[1] - h[0]) / aqin.c[1])\n",
@@ -92,6 +92,13 @@
     "    aqin.T[1] * (d2hdx2[1] + d2hdy2[1]),\n",
     ")"
    ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": []
   }
  ],
  "metadata": {