function y = expfit(beta,x)

y = beta(1) + beta(2)*exp(beta(3)*x(:,1));