%--------------------------------------
function yinit = Finitwithproduction_nop(x)
global rho sigma sigmasqr lambda gamma A phi r elbarnew wlbarnew qinf pinf plbarnew qlbarnew eubar

%no debt case
% beta1=(B*deltabar-Bfb*deltab)/(deltabar-deltab);
%  yinit = [  Bfb*(x-deltab)
%            beta1];
 qprime=(qinf-qlbarnew)/(eubar-elbarnew);
  yinit = [ qlbarnew+qprime*x
           qprime];
       
 
%  prime=0;
%   yinit = [ wlbarnew+prime*x
%       prime];