function F = consumption_trial(x, whs,  kgs, l_v, kfp_v, ss,  igs, Vps, Vws)
global mmunr kkappa tcbar eepsilon oomega cgbar mmur   twbar   ttheta ddelta aalpha_g aalpha_1 aalpha_2 nxybar; 

cnrhat = x(2).^((eepsilon-1)/eepsilon)+ oomega*cgbar^((eepsilon-1)/eepsilon);
lambdanrs= (x(2).^(-1/eepsilon))/((1+tcbar)*cnrhat);
lhnrs = 1-(mmunr/(kkappa*whs*lambdanrs));
F(2)=x(2)- (kkappa*whs*lhnrs)/(1+tcbar);
crhat = x(1).^((eepsilon-1)/eepsilon)+ oomega*cgbar^((eepsilon-1)/eepsilon);
lambdars = (1/((1+tcbar)*crhat))*x(1).^(-1/eepsilon);
lhrs = 1- (mmur/(whs*(1-twbar)*lambdars));
lhs = ttheta*lhrs + (1-ttheta)*kkappa*lhnrs;
ls= lhs/Vws;
vs= (1/l_v)*ls;
kpfs = kfp_v*vs; % Firms Demand for Capital
kps = kpfs/ttheta; % Capital Stock per Ricardian household
ips= ddelta*kps;% Investment per Ricardian household
yfs = kgs^aalpha_g*kpfs^aalpha_1*vs^aalpha_2*ls^(1-aalpha_1-aalpha_2);
y1= yfs/Vps;
nxs= nxybar*y1;
F(1) = x(1)-((y1-(1-ttheta)*x(2)-ttheta*ss*ips - cgbar-ss*igs-nxs)/ttheta); 
