pH_acetic_acid_versus_ionicstrength.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 50,
   "metadata": {},
   "outputs": [],
   "source": [
    "% some housekeeping stuff\n",
    "register_graphics_toolkit(\"gnuplot\");\n",
    "available_graphics_toolkits();\n",
    "graphics_toolkit(\"gnuplot\")\n",
    "clear\n",
    "% end of housekeeping"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# ionic strength effect on pH of acetic acid\n",
    "\n",
    "see notes for description"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 51,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/svg+xml": [
       "<svg xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\" width=\"600px\" height=\"500px\" viewBox=\"0 0 560 420\">\n",
       "\n",
       "<title>Gnuplot</title>\n",
       "<desc>Produced by GNUPLOT 5.2 patchlevel 7 </desc>\n",
       "\n",
       "<g id=\"gnuplot_canvas\">\n",
       "\n",
       "<rect x=\"0\" y=\"0\" width=\"560\" height=\"420\" fill=\"none\"/>\n",
       "<defs>\n",
       "\n",
       "\t<circle id=\"gpDot\" r=\"0.5\" stroke-width=\"0.5\" stroke=\"currentColor\"/>\n",
       "\t<path id=\"gpPt0\" stroke-width=\"0.222\" stroke=\"currentColor\" d=\"M-1,0 h2 M0,-1 v2\"/>\n",
       "\t<path id=\"gpPt1\" stroke-width=\"0.222\" stroke=\"currentColor\" d=\"M-1,-1 L1,1 M1,-1 L-1,1\"/>\n",
       "\t<path id=\"gpPt2\" stroke-width=\"0.222\" stroke=\"currentColor\" d=\"M-1,0 L1,0 M0,-1 L0,1 M-1,-1 L1,1 M-1,1 L1,-1\"/>\n",
       "\t<rect id=\"gpPt3\" stroke-width=\"0.222\" stroke=\"currentColor\" x=\"-1\" y=\"-1\" width=\"2\" height=\"2\"/>\n",
       "\t<rect id=\"gpPt4\" stroke-width=\"0.222\" stroke=\"currentColor\" fill=\"currentColor\" x=\"-1\" y=\"-1\" width=\"2\" height=\"2\"/>\n",
       "\t<circle id=\"gpPt5\" stroke-width=\"0.222\" stroke=\"currentColor\" cx=\"0\" cy=\"0\" r=\"1\"/>\n",
       "\t<use xlink:href=\"#gpPt5\" id=\"gpPt6\" fill=\"currentColor\" stroke=\"none\"/>\n",
       "\t<path id=\"gpPt7\" stroke-width=\"0.222\" stroke=\"currentColor\" d=\"M0,-1.33 L-1.33,0.67 L1.33,0.67 z\"/>\n",
       "\t<use xlink:href=\"#gpPt7\" id=\"gpPt8\" fill=\"currentColor\" stroke=\"none\"/>\n",
       "\t<use xlink:href=\"#gpPt7\" id=\"gpPt9\" stroke=\"currentColor\" transform=\"rotate(180)\"/>\n",
       "\t<use xlink:href=\"#gpPt9\" id=\"gpPt10\" fill=\"currentColor\" stroke=\"none\"/>\n",
       "\t<use xlink:href=\"#gpPt3\" id=\"gpPt11\" stroke=\"currentColor\" transform=\"rotate(45)\"/>\n",
       "\t<use xlink:href=\"#gpPt11\" id=\"gpPt12\" fill=\"currentColor\" stroke=\"none\"/>\n",
       "\t<path id=\"gpPt13\" stroke-width=\"0.222\" stroke=\"currentColor\" d=\"M0,1.330 L1.265,0.411 L0.782,-1.067 L-0.782,-1.076 L-1.265,0.411 z\"/>\n",
       "\t<use xlink:href=\"#gpPt13\" id=\"gpPt14\" fill=\"currentColor\" stroke=\"none\"/>\n",
       "\t<filter id=\"textbox\" filterUnits=\"objectBoundingBox\" x=\"0\" y=\"0\" height=\"1\" width=\"1\">\n",
       "\t  <feFlood flood-color=\"white\" flood-opacity=\"1\" result=\"bgnd\"/>\n",
       "\t  <feComposite in=\"SourceGraphic\" in2=\"bgnd\" operator=\"atop\"/>\n",
       "\t</filter>\n",
       "\t<filter id=\"greybox\" filterUnits=\"objectBoundingBox\" x=\"0\" y=\"0\" height=\"1\" width=\"1\">\n",
       "\t  <feFlood flood-color=\"lightgrey\" flood-opacity=\"1\" result=\"grey\"/>\n",
       "\t  <feComposite in=\"SourceGraphic\" in2=\"grey\" operator=\"atop\"/>\n",
       "\t</filter>\n",
       "</defs>\n",
       "<g fill=\"none\" color=\"white\" stroke=\"currentColor\" stroke-width=\"1.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "</g>\n",
       "<g fill=\"none\" color=\"black\" stroke=\"currentColor\" stroke-width=\"1.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "\t<g stroke=\"none\" shape-rendering=\"crispEdges\">\n",
       "\t\t<polygon fill=\"rgb(255, 255, 255)\" points=\"80.2,362.4 534.9,362.4 534.9,18.2 80.2,18.2 \"/>\n",
       "\t</g>\n",
       "</g>\n",
       "<g fill=\"none\" color=\"black\" stroke=\"currentColor\" stroke-width=\"1.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "\t<path stroke=\"rgb(  0,   0,   0)\" d=\"M80.2,362.4 L92.7,362.4 M535.0,362.4 L522.5,362.4  \"/>\t<g transform=\"translate(71.9,365.7)\" stroke=\"none\" fill=\"rgb(38,38,38)\" font-family=\"Arial\" font-size=\"10.00\" text-anchor=\"end\">\n",
       "\t\t<text><tspan font-family=\"Arial\">3.372</tspan></text>\n",
       "\t</g>\n",
       "</g>\n",
       "<g fill=\"none\" color=\"black\" stroke=\"currentColor\" stroke-width=\"1.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "\t<path stroke=\"rgb(  0,   0,   0)\" d=\"M80.2,313.2 L92.7,313.2 M535.0,313.2 L522.5,313.2  \"/>\t<g transform=\"translate(71.9,316.5)\" stroke=\"none\" fill=\"rgb(38,38,38)\" font-family=\"Arial\" font-size=\"10.00\" text-anchor=\"end\">\n",
       "\t\t<text><tspan font-family=\"Arial\">3.374</tspan></text>\n",
       "\t</g>\n",
       "</g>\n",
       "<g fill=\"none\" color=\"black\" stroke=\"currentColor\" stroke-width=\"1.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "\t<path stroke=\"rgb(  0,   0,   0)\" d=\"M80.2,264.0 L92.7,264.0 M535.0,264.0 L522.5,264.0  \"/>\t<g transform=\"translate(71.9,267.3)\" stroke=\"none\" fill=\"rgb(38,38,38)\" font-family=\"Arial\" font-size=\"10.00\" text-anchor=\"end\">\n",
       "\t\t<text><tspan font-family=\"Arial\">3.376</tspan></text>\n",
       "\t</g>\n",
       "</g>\n",
       "<g fill=\"none\" color=\"black\" stroke=\"currentColor\" stroke-width=\"1.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "\t<path stroke=\"rgb(  0,   0,   0)\" d=\"M80.2,214.8 L92.7,214.8 M535.0,214.8 L522.5,214.8  \"/>\t<g transform=\"translate(71.9,218.1)\" stroke=\"none\" fill=\"rgb(38,38,38)\" font-family=\"Arial\" font-size=\"10.00\" text-anchor=\"end\">\n",
       "\t\t<text><tspan font-family=\"Arial\">3.378</tspan></text>\n",
       "\t</g>\n",
       "</g>\n",
       "<g fill=\"none\" color=\"black\" stroke=\"currentColor\" stroke-width=\"1.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "\t<path stroke=\"rgb(  0,   0,   0)\" d=\"M80.2,165.7 L92.7,165.7 M535.0,165.7 L522.5,165.7  \"/>\t<g transform=\"translate(71.9,169.0)\" stroke=\"none\" fill=\"rgb(38,38,38)\" font-family=\"Arial\" font-size=\"10.00\" text-anchor=\"end\">\n",
       "\t\t<text><tspan font-family=\"Arial\">3.38</tspan></text>\n",
       "\t</g>\n",
       "</g>\n",
       "<g fill=\"none\" color=\"black\" stroke=\"currentColor\" stroke-width=\"1.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "\t<path stroke=\"rgb(  0,   0,   0)\" d=\"M80.2,116.5 L92.7,116.5 M535.0,116.5 L522.5,116.5  \"/>\t<g transform=\"translate(71.9,119.8)\" stroke=\"none\" fill=\"rgb(38,38,38)\" font-family=\"Arial\" font-size=\"10.00\" text-anchor=\"end\">\n",
       "\t\t<text><tspan font-family=\"Arial\">3.382</tspan></text>\n",
       "\t</g>\n",
       "</g>\n",
       "<g fill=\"none\" color=\"black\" stroke=\"currentColor\" stroke-width=\"1.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "\t<path stroke=\"rgb(  0,   0,   0)\" d=\"M80.2,67.3 L92.7,67.3 M535.0,67.3 L522.5,67.3  \"/>\t<g transform=\"translate(71.9,70.6)\" stroke=\"none\" fill=\"rgb(38,38,38)\" font-family=\"Arial\" font-size=\"10.00\" text-anchor=\"end\">\n",
       "\t\t<text><tspan font-family=\"Arial\">3.384</tspan></text>\n",
       "\t</g>\n",
       "</g>\n",
       "<g fill=\"none\" color=\"black\" stroke=\"currentColor\" stroke-width=\"1.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "\t<path stroke=\"rgb(  0,   0,   0)\" d=\"M80.2,18.1 L92.7,18.1 M535.0,18.1 L522.5,18.1  \"/>\t<g transform=\"translate(71.9,21.4)\" stroke=\"none\" fill=\"rgb(38,38,38)\" font-family=\"Arial\" font-size=\"10.00\" text-anchor=\"end\">\n",
       "\t\t<text><tspan font-family=\"Arial\">3.386</tspan></text>\n",
       "\t</g>\n",
       "</g>\n",
       "<g fill=\"none\" color=\"black\" stroke=\"currentColor\" stroke-width=\"1.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "\t<path stroke=\"rgb(  0,   0,   0)\" d=\"M80.2,362.4 L80.2,349.9 M80.2,18.1 L80.2,30.6  \"/>\t<g transform=\"translate(80.2,383.7)\" stroke=\"none\" fill=\"rgb(38,38,38)\" font-family=\"Arial\" font-size=\"10.00\" text-anchor=\"middle\">\n",
       "\t\t<text><tspan font-family=\"Arial\">-6</tspan></text>\n",
       "\t</g>\n",
       "</g>\n",
       "<g fill=\"none\" color=\"black\" stroke=\"currentColor\" stroke-width=\"1.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "\t<path stroke=\"rgb(  0,   0,   0)\" d=\"M171.2,362.4 L171.2,349.9 M171.2,18.1 L171.2,30.6  \"/>\t<g transform=\"translate(171.2,383.7)\" stroke=\"none\" fill=\"rgb(38,38,38)\" font-family=\"Arial\" font-size=\"10.00\" text-anchor=\"middle\">\n",
       "\t\t<text><tspan font-family=\"Arial\">-5</tspan></text>\n",
       "\t</g>\n",
       "</g>\n",
       "<g fill=\"none\" color=\"black\" stroke=\"currentColor\" stroke-width=\"1.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "\t<path stroke=\"rgb(  0,   0,   0)\" d=\"M262.1,362.4 L262.1,349.9 M262.1,18.1 L262.1,30.6  \"/>\t<g transform=\"translate(262.1,383.7)\" stroke=\"none\" fill=\"rgb(38,38,38)\" font-family=\"Arial\" font-size=\"10.00\" text-anchor=\"middle\">\n",
       "\t\t<text><tspan font-family=\"Arial\">-4</tspan></text>\n",
       "\t</g>\n",
       "</g>\n",
       "<g fill=\"none\" color=\"black\" stroke=\"currentColor\" stroke-width=\"1.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "\t<path stroke=\"rgb(  0,   0,   0)\" d=\"M353.1,362.4 L353.1,349.9 M353.1,18.1 L353.1,30.6  \"/>\t<g transform=\"translate(353.1,383.7)\" stroke=\"none\" fill=\"rgb(38,38,38)\" font-family=\"Arial\" font-size=\"10.00\" text-anchor=\"middle\">\n",
       "\t\t<text><tspan font-family=\"Arial\">-3</tspan></text>\n",
       "\t</g>\n",
       "</g>\n",
       "<g fill=\"none\" color=\"black\" stroke=\"currentColor\" stroke-width=\"1.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "\t<path stroke=\"rgb(  0,   0,   0)\" d=\"M444.0,362.4 L444.0,349.9 M444.0,18.1 L444.0,30.6  \"/>\t<g transform=\"translate(444.0,383.7)\" stroke=\"none\" fill=\"rgb(38,38,38)\" font-family=\"Arial\" font-size=\"10.00\" text-anchor=\"middle\">\n",
       "\t\t<text><tspan font-family=\"Arial\">-2</tspan></text>\n",
       "\t</g>\n",
       "</g>\n",
       "<g fill=\"none\" color=\"black\" stroke=\"currentColor\" stroke-width=\"1.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "\t<path stroke=\"rgb(  0,   0,   0)\" d=\"M535.0,362.4 L535.0,349.9 M535.0,18.1 L535.0,30.6  \"/>\t<g transform=\"translate(535.0,383.7)\" stroke=\"none\" fill=\"rgb(38,38,38)\" font-family=\"Arial\" font-size=\"10.00\" text-anchor=\"middle\">\n",
       "\t\t<text><tspan font-family=\"Arial\">-1</tspan></text>\n",
       "\t</g>\n",
       "</g>\n",
       "<g fill=\"none\" color=\"black\" stroke=\"currentColor\" stroke-width=\"1.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "</g>\n",
       "<g fill=\"none\" color=\"black\" stroke=\"currentColor\" stroke-width=\"1.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "\t<g transform=\"translate(24.9,190.3) rotate(-90)\" stroke=\"none\" fill=\"rgb(38,38,38)\" font-family=\"Arial\" font-size=\"11.00\" text-anchor=\"middle\">\n",
       "\t\t<text><tspan font-family=\"Arial\">pH</tspan></text>\n",
       "\t</g>\n",
       "</g>\n",
       "<g fill=\"none\" color=\"black\" stroke=\"currentColor\" stroke-width=\"1.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "\t<g transform=\"translate(307.6,411.0)\" stroke=\"none\" fill=\"rgb(38,38,38)\" font-family=\"Arial\" font-size=\"11.00\" text-anchor=\"middle\">\n",
       "\t\t<text><tspan font-family=\"Arial\">log(</tspan><tspan font-family=\"Symbol\">m</tspan><tspan font-family=\"Arial\">)</tspan></text>\n",
       "\t</g>\n",
       "</g>\n",
       "<g fill=\"none\" color=\"black\" stroke=\"currentColor\" stroke-width=\"1.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "</g>\n",
       "\t<g id=\"gnuplot_plot_1a\"><title>gnuplot_plot_1a</title>\n",
       "<g fill=\"none\" color=\"black\" stroke=\"currentColor\" stroke-width=\"1.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "\t<use xlink:href=\"#gpPt6\" transform=\"translate(80.2,63.2) scale(6.00)\" color=\"rgb(  0,   0, 255)\"/>\n",
       "\t<use xlink:href=\"#gpPt6\" transform=\"translate(98.4,63.1) scale(6.00)\" color=\"rgb(  0,   0, 255)\"/>\n",
       "\t<use xlink:href=\"#gpPt6\" transform=\"translate(116.6,63.0) scale(6.00)\" color=\"rgb(  0,   0, 255)\"/>\n",
       "\t<use xlink:href=\"#gpPt6\" transform=\"translate(134.8,62.9) scale(6.00)\" color=\"rgb(  0,   0, 255)\"/>\n",
       "\t<use xlink:href=\"#gpPt6\" transform=\"translate(153.0,62.8) scale(6.00)\" color=\"rgb(  0,   0, 255)\"/>\n",
       "\t<use xlink:href=\"#gpPt6\" transform=\"translate(171.2,62.7) scale(6.00)\" color=\"rgb(  0,   0, 255)\"/>\n",
       "\t<use xlink:href=\"#gpPt6\" transform=\"translate(189.4,62.5) scale(6.00)\" color=\"rgb(  0,   0, 255)\"/>\n",
       "\t<use xlink:href=\"#gpPt6\" transform=\"translate(207.5,62.3) scale(6.00)\" color=\"rgb(  0,   0, 255)\"/>\n",
       "\t<use xlink:href=\"#gpPt6\" transform=\"translate(225.7,62.1) scale(6.00)\" color=\"rgb(  0,   0, 255)\"/>\n",
       "\t<use xlink:href=\"#gpPt6\" transform=\"translate(243.9,61.9) scale(6.00)\" color=\"rgb(  0,   0, 255)\"/>\n",
       "\t<use xlink:href=\"#gpPt6\" transform=\"translate(262.1,61.6) scale(6.00)\" color=\"rgb(  0,   0, 255)\"/>\n",
       "\t<use xlink:href=\"#gpPt6\" transform=\"translate(280.3,61.4) scale(6.00)\" color=\"rgb(  0,   0, 255)\"/>\n",
       "\t<use xlink:href=\"#gpPt6\" transform=\"translate(298.5,61.2) scale(6.00)\" color=\"rgb(  0,   0, 255)\"/>\n",
       "\t<use xlink:href=\"#gpPt6\" transform=\"translate(316.7,61.3) scale(6.00)\" color=\"rgb(  0,   0, 255)\"/>\n",
       "\t<use xlink:href=\"#gpPt6\" transform=\"translate(334.9,61.6) scale(6.00)\" color=\"rgb(  0,   0, 255)\"/>\n",
       "\t<use xlink:href=\"#gpPt6\" transform=\"translate(353.1,62.6) scale(6.00)\" color=\"rgb(  0,   0, 255)\"/>\n",
       "\t<use xlink:href=\"#gpPt6\" transform=\"translate(371.3,64.5) scale(6.00)\" color=\"rgb(  0,   0, 255)\"/>\n",
       "\t<use xlink:href=\"#gpPt6\" transform=\"translate(389.5,67.9) scale(6.00)\" color=\"rgb(  0,   0, 255)\"/>\n",
       "\t<use xlink:href=\"#gpPt6\" transform=\"translate(407.7,73.6) scale(6.00)\" color=\"rgb(  0,   0, 255)\"/>\n",
       "\t<use xlink:href=\"#gpPt6\" transform=\"translate(425.8,82.8) scale(6.00)\" color=\"rgb(  0,   0, 255)\"/>\n",
       "\t<use xlink:href=\"#gpPt6\" transform=\"translate(444.0,96.9) scale(6.00)\" color=\"rgb(  0,   0, 255)\"/>\n",
       "\t<use xlink:href=\"#gpPt6\" transform=\"translate(462.2,118.0) scale(6.00)\" color=\"rgb(  0,   0, 255)\"/>\n",
       "\t<use xlink:href=\"#gpPt6\" transform=\"translate(480.4,148.3) scale(6.00)\" color=\"rgb(  0,   0, 255)\"/>\n",
       "\t<use xlink:href=\"#gpPt6\" transform=\"translate(498.6,190.4) scale(6.00)\" color=\"rgb(  0,   0, 255)\"/>\n",
       "\t<use xlink:href=\"#gpPt6\" transform=\"translate(516.8,246.9) scale(6.00)\" color=\"rgb(  0,   0, 255)\"/>\n",
       "\t<use xlink:href=\"#gpPt6\" transform=\"translate(535.0,320.1) scale(6.00)\" color=\"rgb(  0,   0, 255)\"/>\n",
       "</g>\n",
       "\t</g>\n",
       "\t<g id=\"gnuplot_plot_2a\"><title>gnuplot_plot_2a</title>\n",
       "<g fill=\"none\" color=\"black\" stroke=\"currentColor\" stroke-width=\"1.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "\t<use xlink:href=\"#gpPt5\" transform=\"translate(80.2,63.2) scale(6.00)\" color=\"rgb(  0,   0,   0)\"/>\n",
       "\t<use xlink:href=\"#gpPt5\" transform=\"translate(98.4,63.1) scale(6.00)\" color=\"rgb(  0,   0,   0)\"/>\n",
       "\t<use xlink:href=\"#gpPt5\" transform=\"translate(116.6,63.0) scale(6.00)\" color=\"rgb(  0,   0,   0)\"/>\n",
       "\t<use xlink:href=\"#gpPt5\" transform=\"translate(134.8,62.9) scale(6.00)\" color=\"rgb(  0,   0,   0)\"/>\n",
       "\t<use xlink:href=\"#gpPt5\" transform=\"translate(153.0,62.8) scale(6.00)\" color=\"rgb(  0,   0,   0)\"/>\n",
       "\t<use xlink:href=\"#gpPt5\" transform=\"translate(171.2,62.7) scale(6.00)\" color=\"rgb(  0,   0,   0)\"/>\n",
       "\t<use xlink:href=\"#gpPt5\" transform=\"translate(189.4,62.5) scale(6.00)\" color=\"rgb(  0,   0,   0)\"/>\n",
       "\t<use xlink:href=\"#gpPt5\" transform=\"translate(207.5,62.3) scale(6.00)\" color=\"rgb(  0,   0,   0)\"/>\n",
       "\t<use xlink:href=\"#gpPt5\" transform=\"translate(225.7,62.1) scale(6.00)\" color=\"rgb(  0,   0,   0)\"/>\n",
       "\t<use xlink:href=\"#gpPt5\" transform=\"translate(243.9,61.9) scale(6.00)\" color=\"rgb(  0,   0,   0)\"/>\n",
       "\t<use xlink:href=\"#gpPt5\" transform=\"translate(262.1,61.6) scale(6.00)\" color=\"rgb(  0,   0,   0)\"/>\n",
       "\t<use xlink:href=\"#gpPt5\" transform=\"translate(280.3,61.4) scale(6.00)\" color=\"rgb(  0,   0,   0)\"/>\n",
       "\t<use xlink:href=\"#gpPt5\" transform=\"translate(298.5,61.2) scale(6.00)\" color=\"rgb(  0,   0,   0)\"/>\n",
       "\t<use xlink:href=\"#gpPt5\" transform=\"translate(316.7,61.3) scale(6.00)\" color=\"rgb(  0,   0,   0)\"/>\n",
       "\t<use xlink:href=\"#gpPt5\" transform=\"translate(334.9,61.6) scale(6.00)\" color=\"rgb(  0,   0,   0)\"/>\n",
       "\t<use xlink:href=\"#gpPt5\" transform=\"translate(353.1,62.6) scale(6.00)\" color=\"rgb(  0,   0,   0)\"/>\n",
       "\t<use xlink:href=\"#gpPt5\" transform=\"translate(371.3,64.5) scale(6.00)\" color=\"rgb(  0,   0,   0)\"/>\n",
       "\t<use xlink:href=\"#gpPt5\" transform=\"translate(389.5,67.9) scale(6.00)\" color=\"rgb(  0,   0,   0)\"/>\n",
       "\t<use xlink:href=\"#gpPt5\" transform=\"translate(407.7,73.6) scale(6.00)\" color=\"rgb(  0,   0,   0)\"/>\n",
       "\t<use xlink:href=\"#gpPt5\" transform=\"translate(425.8,82.8) scale(6.00)\" color=\"rgb(  0,   0,   0)\"/>\n",
       "\t<use xlink:href=\"#gpPt5\" transform=\"translate(444.0,96.9) scale(6.00)\" color=\"rgb(  0,   0,   0)\"/>\n",
       "\t<use xlink:href=\"#gpPt5\" transform=\"translate(462.2,118.0) scale(6.00)\" color=\"rgb(  0,   0,   0)\"/>\n",
       "\t<use xlink:href=\"#gpPt5\" transform=\"translate(480.4,148.3) scale(6.00)\" color=\"rgb(  0,   0,   0)\"/>\n",
       "\t<use xlink:href=\"#gpPt5\" transform=\"translate(498.6,190.4) scale(6.00)\" color=\"rgb(  0,   0,   0)\"/>\n",
       "\t<use xlink:href=\"#gpPt5\" transform=\"translate(516.8,246.9) scale(6.00)\" color=\"rgb(  0,   0,   0)\"/>\n",
       "\t<use xlink:href=\"#gpPt5\" transform=\"translate(535.0,320.1) scale(6.00)\" color=\"rgb(  0,   0,   0)\"/>\n",
       "</g>\n",
       "\t</g>\n",
       "\t<g id=\"gnuplot_plot_3a\"><title>gnuplot_plot_3a</title>\n",
       "<g fill=\"none\" color=\"white\" stroke=\"rgb(  0,   0,   0)\" stroke-width=\"2.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "</g>\n",
       "<g fill=\"none\" color=\"black\" stroke=\"currentColor\" stroke-width=\"2.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "\t<path stroke=\"rgb(  0,   0,   0)\" stroke-dasharray=\"20.0,20.0\" d=\"M80.2,63.4 L98.4,63.4 L116.6,63.4 L134.8,63.4 L153.0,63.4 L171.2,63.4 L189.4,63.4 L207.5,63.4   L225.7,63.4 L243.9,63.4 L262.1,63.4 L280.3,63.4 L298.5,63.4 L316.7,63.4 L334.9,63.4 L353.1,63.4   L371.3,63.4 L389.5,63.4 L407.7,63.4 L425.8,63.4 L444.0,63.4 L462.2,63.4 L480.4,63.4 L498.6,63.4   L516.8,63.4 L535.0,63.4  \"/></g>\n",
       "\t</g>\n",
       "<g fill=\"none\" color=\"black\" stroke=\"currentColor\" stroke-width=\"2.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "</g>\n",
       "<g fill=\"none\" color=\"black\" stroke=\"black\" stroke-width=\"1.50\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "\t<path stroke=\"rgb( 38,  38,  38)\" d=\"M80.2,362.4 L535.0,362.4 M80.2,18.1 L535.0,18.1 M80.2,362.4 L80.2,18.1 M535.0,362.4 L535.0,18.1    \"/></g>\n",
       "<g fill=\"none\" color=\"black\" stroke=\"rgb( 38,  38,  38)\" stroke-width=\"1.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "</g>\n",
       "<g fill=\"none\" color=\"black\" stroke=\"currentColor\" stroke-width=\"1.00\" stroke-linecap=\"butt\" stroke-linejoin=\"miter\">\n",
       "</g>\n",
       "</g>\n",
       "</svg>"
      ],
      "text/plain": [
       "<IPython.core.display.SVG object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%plot -s 600,500 -f 'svg'\n",
    "logmu=-6:0.2:-1; mu=10.^logmu; AT=0.01; Ka=10^-4.75; sizeH=900; sizeAc=450;\n",
    "\n",
    "%zero ionic strength\n",
    "a=1; b=Ka; c=-Ka*AT; t=roots([a b c]); t(imag(t)==0); t=t(t>0); pHzero=-log10(t);\n",
    "\n",
    "for i=1:length(logmu)\n",
    "    loggammaH=-0.51*sqrt(mu(i))/(1+((sizeH/350)*sqrt(mu(i)))); gammaH=10^loggammaH;\n",
    "    loggammaAc=-0.51*sqrt(mu(i))/(1+((sizeAc/350)*sqrt(mu(i)))); gammaAc=10^loggammaAc;\n",
    "    Kaprime=Ka/(gammaH*gammaAc);\n",
    "    a=1; b=Kaprime; c=-Kaprime*AT;\n",
    "    t=roots([a b c]); \n",
    "    t(imag(t)==0); %sets any imaginary roots to zero\n",
    "    t=t(t>0); %t=positive real roots\n",
    "    pH(i)=-log10(t*gammaH);\n",
    "end\n",
    "\n",
    "plot(logmu,pH,'ko','markersize',4,'markerfacecolor','b')\n",
    "set(gca,'linewidth',1.5)\n",
    "xlabel('log(\\mu)'); ylabel('pH')\n",
    "hold on\n",
    "plot(logmu,pHzero*ones(size(logmu)),'k--','linewidth',2)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Octave",
   "language": "octave",
   "name": "octave"
  },
  "language_info": {
   "file_extension": ".m",
   "help_links": [
    {
     "text": "GNU Octave",
     "url": "https://www.gnu.org/software/octave/support.html"
    },
    {
     "text": "Octave Kernel",
     "url": "https://github.com/Calysto/octave_kernel"
    },
    {
     "text": "MetaKernel Magics",
     "url": "https://metakernel.readthedocs.io/en/latest/source/README.html"
    }
   ],
   "mimetype": "text/x-octave",
   "name": "octave",
   "version": "6.2.0"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}