omment">% x (mg / L COD)</span>
data.ydata = [0.053 0.060 0.112 0.105 0.099 0.122 0.125]'; <span class="comment">% y (1 / h)</span>
</pre><p>Here is a plot of the data set.</p><pre class="codeinput">figure(1); clf
plot(data.xdata,data.ydata,<span class="string">'s'</span>);
xlim([0 400]); xlabel(<span class="string">'x [mg/L COD]'</span>); ylabel(<span class="string">'y [1/h]'</span>);
</pre><img vspace="5" hspace="5" src="monodex_01.png" alt=""> <p>For the MCMC run we need the sum of squares function. For 