function [ youtput ] = FunctionToFit( xparams, xinput )
% produces the output of the analytic expression for Monod occupancy rate
% takes a 3-vector (KR, L, C) "xparams"
% returns a vector of y-values corresponding to the x-values in xinput
kr = xparams(1); L = xparams(2); kt = xparams(3);
n = 2;

for i=1:length(xinput)
  youtput(i) = (1+xinput(i)/kr)^n / (L*(1+xinput(i)/kt)^n + (1+xinput(i)/kr)^n);
end

end

