function F = paramfun(x, tau, w, omega, m, Ty_sp, Ty_sa, r_b, by_sa,dy, xi_1, xi_2, Ky, Ny, alpha, psi, nu, Lambda)

F(1) = x(1)-(1-tau)*(w*x(2)+(1-omega)*(1-m))-Ty_sp;
F(2) = x(3)-(1-tau)*(w*x(4)+(1-omega)*(1-m))-Ty_sa-r_b*by_sa+dy+xi_1*abs(dy)^xi_2*(Ky/Ny)^(-alpha*(1-xi_2))*x(5)^(xi_2-1);
F(3) = x(2)-(w*(1-tau)*x(1)^(-1)/psi)^(1/nu)*(Ky/Ny)^(-alpha*(1+nu)/nu)*x(5)^(-(1+nu)/nu);
F(4) = x(4)-(w*(1-tau)*x(3)^(-1)/psi)^(1/nu)*(Ky/Ny)^(-alpha*(1+nu)/nu)*x(5)^(-(1+nu)/nu);
F(5) = Ny-Lambda*x(2)-(1-Lambda)*x(4);


