diff --git a/docs/usage/basics.ipynb b/docs/usage/basics.ipynb
index af3f37c8..edf7cc0f 100644
--- a/docs/usage/basics.ipynb
+++ b/docs/usage/basics.ipynb
@@ -1286,10 +1286,10 @@
     "fig.suptitle(\"Evolution of the parameter values during the fit\")\n",
     "ax1.set_title(\"Minuit2\")\n",
     "ax2.set_title(\"Scipy (Nelder-Mead)\")\n",
-    "ax1.get_shared_x_axes().join(ax1, ax3)\n",
-    "ax2.get_shared_x_axes().join(ax2, ax4)\n",
-    "ax1.get_shared_y_axes().join(ax1, ax2)\n",
-    "ax3.get_shared_y_axes().join(ax3, ax4)\n",
+    "ax1.sharex(ax3)\n",
+    "ax2.sharex(ax4)\n",
+    "ax1.sharey(ax2)\n",
+    "ax3.sharey(ax4)\n",
     "ax2.set_ylim(\n",
     "    1.02 * scipy_traceback[\"estimator_value\"].min(),\n",
     "    0.98 * scipy_traceback[\"estimator_value\"].max(),\n",
diff --git a/docs/usage/unbinned-fit.ipynb b/docs/usage/unbinned-fit.ipynb
index 31e4e245..35eb9c4a 100644
--- a/docs/usage/unbinned-fit.ipynb
+++ b/docs/usage/unbinned-fit.ipynb
@@ -128,8 +128,8 @@
     "empty.remove()\n",
     "ax2d.set_xlabel(\"$x$\")\n",
     "ax2d.set_ylabel(\"$y$\")\n",
-    "ax2d.get_shared_x_axes().join(ax2d, ax_x)\n",
-    "ax2d.get_shared_y_axes().join(ax2d, ax_y)\n",
+    "ax2d.sharex(ax_x)\n",
+    "ax2d.sharey(ax_y)\n",
     "for ax in [ax_x, ax_y]:\n",
     "    ax.set_xticks([])\n",
     "    ax.set_yticks([])\n",