diff --git a/FloryHugg_DiffUnbleached.ipynb b/FloryHugg_DiffUnbleached.ipynb
index 04e3b89ec7e788783d046da12ecce802356298df..3d430d78c2bc66e7f688afdecb4bf59ffdc3238b 100644
--- a/FloryHugg_DiffUnbleached.ipynb
+++ b/FloryHugg_DiffUnbleached.ipynb
@@ -92,13 +92,18 @@
     "c0_2 = create_func(F, 'x[0]<0.1 ? 0 :'+p_tot(0.9, 0.45), 1)\n",
     "Ga0_2 = create_func(F,p_tot(1, 0.08), 1)\n",
     "\n",
+    "c_tot_3 = create_func(F, p_tot(0.9, 0.3), 1)\n",
+    "c0_3 = create_func(F, 'x[0]<0.1 ? 0 :'+p_tot(0.9, 0.3), 1)\n",
+    "Ga0_3 = create_func(F,p_tot(1, 0.145), 1)\n",
+    "\n",
     "c_tot_9 = create_func(F, p_tot(0.9, 0.1), 1)\n",
     "c0_9 = create_func(F, 'x[0]<0.1 ? 0 :'+p_tot(0.9, 0.1), 1)\n",
     "Ga0_9 = create_func(F,p_tot(1, 1), 1)\n",
     "\n",
-    "c0_1 = calc_sim(c0_1, c_tot_1, Ga0_1, 0.001, 50, 0)\n",
-    "c0_2 = calc_sim(c0_2, c_tot_2, Ga0_2, 0.001, 50, 0)\n",
-    "c0_9 = calc_sim(c0_9, c_tot_9, Ga0_9, 0.001, 50, 0)"
+    "c0_1 = calc_sim(c0_1, c_tot_1, Ga0_1, 0.001, 10, 0)\n",
+    "c0_2 = calc_sim(c0_2, c_tot_2, Ga0_2, 0.001, 10, 0)\n",
+    "c0_9 = calc_sim(c0_9, c_tot_9, Ga0_9, 0.001, 10, 0)\n",
+    "c0_3 = calc_sim(c0_3, c_tot_3, Ga0_3, 0.001, 10, 0)"
    ]
   },
   {
@@ -110,8 +115,15 @@
     "x = np.linspace(0, 1, 10000)\n",
     "plt.plot(x, eval_func(c0_1, x))\n",
     "plt.plot(x, eval_func(c0_2, x))\n",
+    "plt.plot(x, eval_func(c0_3, x))\n",
     "plt.plot(x, eval_func(c0_9, x))\n",
-    "plt.xlim(0.0, 0.125)"
+    "plt.xlim(0.0, 0.125)\n",
+    "plt.show()\n",
+    "# Diffusion versus partitioning\n",
+    "plt.plot([Ga0_1(1)*(1-c_tot_1(1)), Ga0_2(1)*(1-c_tot_2(1)), Ga0_3(1)*(1-c_tot_3(1)), Ga0_9(1)*(1-c_tot_9(1))],\n",
+    "         [c_tot_1(0)/c_tot_1(1), c_tot_2(0)/c_tot_2(1), c_tot_3(0)/c_tot_3(1), c_tot_9(0)/c_tot_9(1)])\n",
+    "plt.ylim(0, 9); plt.xlim(0, 1)\n",
+    "plt.plot(np.linspace(0, 1, 100), 9.5*np.sqrt(np.linspace(0, 1, 100)))"
    ]
   },
   {