function outvalue=solve_calibration(par_vec,rhos,betta,s,kappa,nup,nu,zeta5,zeta6,zeta7,gammac,sigmm,sigmaa)

J=par_vec(1);
f=par_vec(2);
wp=par_vec(3);
U=par_vec(4);
V=par_vec(5);
A=par_vec(6);
Q=par_vec(7);

outvalue(1)=J-nup+wp; % (eq 15)
outvalue(2)=J-(s/Q+kappa); % (eq 18)
outvalue(3)=V-wp-A; % (eq 19) 
outvalue(4)=(1-(1-rhos)*betta)*A-rhos*betta*(f*V+(1-f)*U); % (eq 20)
outvalue(5)=0.4*(1-(1-rhos)*betta)*wp+betta*f*V+(betta*(1-f)-1)*U; % (eq 21) plugged with (p8) at steadystate.m
outvalue(6)=J-zeta5*(V-U)-zeta7*nu+(zeta6*gammac+zeta7*0.4*(1-(1-rhos)*betta)*wp); % (eq 25) plugged with (p8) at steadystate.m
outvalue(7)=Q-sigmm^(1/(1-sigmaa))*f^(-sigmaa/(1-sigmaa)); % (eq 14), (eq 22), and (eq 23)
end