function out=monod(x,p);
out=(p(1)*x)./(p(2)+x);