diff --git a/FloryHugg_DiffUnbleached.ipynb b/FloryHugg_DiffUnbleached.ipynb
index b8ed00758c9555df5cb71c02206c340cc9b54116..da0718cd718bef5f6b8ea6b4539900bb8e887ad5 100644
--- a/FloryHugg_DiffUnbleached.ipynb
+++ b/FloryHugg_DiffUnbleached.ipynb
@@ -88,12 +88,17 @@
     "c0_1 = create_func(F, 'x[0]<0.1 ? 0 :' + p_tot(0.9, 0.9), 1)\n",
     "Ga0_1 = create_func(F, p_tot(1, 1/9), 1)\n",
     "\n",
+    "c_tot_2 = create_func(F, p_tot(0.9, 0.45), 1)\n",
+    "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_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, 10, 2)\n",
-    "c0_9 = calc_sim(c0_9, c_tot_9, Ga0_9, 0.001, 10, 2)"
+    "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)"
    ]
   },
   {
@@ -104,8 +109,20 @@
    "source": [
     "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_9, x))\n",
-    "plt.xlim(0, 0.2)"
+    "plt.xlim(0.0, 0.125)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# Check partitioning\n",
+    "cp = eval_func(c0_9, x)\n",
+    "print('Partitioning: ' + str(np.max(cp[1:1050])/np.min(cp[999:1050])))"
    ]
   },
   {