% Cobb-Douglas production function
% Y=F(K,N)=A*K^{alpha}*N^{1-alpha}
% In per capita form y=A*k^{alpha}
% Need to get the capital share parameter and scale prameter from main.m

function y=cobb(k)
global A alpha
if alpha>0 & alpha<1
   y=A*k.^alpha;
else
   y=-inf*abs(k);
   disp('invalid capital share')
end
