function outvalue=solve_calibration(par_var,gammaa, xii, omegaa, rhoo, alphaa, nuu, deltaa, betaa, nn, xx, Gamm)

Z=par_var(1);
R=par_var(2);
epsilontheta=par_var(3);
theta=par_var(4);
om=par_var(5);
c=par_var(6);
cr=par_var(7);
cw=par_var(8);
l=par_var(9);
lw=par_var(10);
lr=par_var(11);
lambda=par_var(12);
a=par_var(13);
hr=par_var(14);
hw=par_var(15);
y=par_var(16);
ii=par_var(17);
w=par_var(18);
k=par_var(19);

outvalue(1)=Z-1;
outvalue(2)=R-1+deltaa-alphaa*y/k;
outvalue(3)=epsilontheta-1+R^rhoo/(1-rhoo)*gammaa*betaa^1/(1-rhoo)*(1/(1+xx))^rhoo*(1-nuu)/(1-rhoo); 
outvalue(4)=theta-1+(om)^(rhoo/(1-rhoo))*betaa*(1/(1+xx))^rhoo*(1-nuu)/(1-rhoo);
outvalue(5)=om-omegaa-(1-omegaa)*(epsilontheta/theta)^(rhoo-1)/rhoo*xii^(nuu-1); 
outvalue(6)=cr-epsilontheta*(R*lambda*k+hr);
outvalue(7)=cw-(epsilontheta/theta)*(R*(1-lambda)*a+hw);
outvalue(8)=c-cr-cw;
outvalue(9)=l-lw-xii*lr;
outvalue(10)=lr-Gamm+(1-nuu)*cr/nuu*xii*w;
outvalue(11)=lw-1+(1-nuu)*cw/nuu*w;
outvalue(12)=lambda-(omegaa*((1-epsilontheta)*hr+w*xii*lr-hr)+(1-omegaa)*a*(1+xx+nn))/a*(1+xx+nn-omegaa*(1-epsilontheta)*R);
outvalue(13)=a-k;
outvalue(14)=hr-w*xii*lr/(1-gammaa*(1+xx)/R);
outvalue(15)=hw-(w*lw+(1-omegaa)*(epsilontheta/theta)^(rhoo-1)/rhoo*(1/xii)^(1-nuu)*hr*(1+xx)/R*om)/(1-omegaa*(1+xx)/R*om);
outvalue(16)=y-l^(1-alphaa)*k^alphaa;
outvalue(17)=ii-y+c;
outvalue(18)=w-(1-alphaa)*y/l;
outvalue(19)=k-ii/(xx+nn+deltaa);
end