function J = beetle_bh(p0,pndata,pn1data)
% Least Squares fit to Logistic Growth
N = length(pndata);
err = [pn1data - (p0(1)*pndata./(1 + p0(2)*pndata))];
J = err*err'; % Sum of square errors
end

