% % this a script file that reads in dh_dzp, di_dzp, dq_dzp, dc_dzp, % d2h_dzp2, d2i_dzp2, d2q_dzp2, d2c_dzp2 % dh_dzp = zeros(nz+np,1); di_dzp = zeros(nz+np,1); dq_dzp = zeros(nz+np,1); dc_dzp = zeros(nz+np,1); d2h_dzp2 = zeros(nz+np,nz+np); d2i_dzp2 = zeros(nz+np,nz+np); d2q_dzp2 = zeros(nz+np,nz+np); d2c_dzp2 = zeros(nz+np,nz+np); dr_dz(1,1)=0; dr_dz(2,1)=(1-alpha)*exp(A)*exp(k)^(-alpha)*h1*(exp(n1)*exp(w))^alpha; dr_dz(3,1)=0; dr_dz(4,1)=-(1-alpha)*exp(A)*exp(k)^(-alpha)*alpha*h1*(exp(n1)*exp(w))^alpha; dr_dz(5,1)=(1-alpha)*exp(A)*exp(k)^(-alpha)*alpha*h1*(exp(n1)*exp(w))^alpha; dr_dz(6,1)=0; dr_dz(7,1)=(1-alpha)*exp(A)*exp(k)^(-alpha)*alpha*h1*(exp(n1)*exp(w))^alpha; dr_dz(8,1)=0; dh_dzp(16,1)=1/h1; dh_dzp(5,1)=1; di_dzp(19,1)=exp(kp1)/(gamma*exp(kp1)-(1-delta)*exp(k)); di_dzp(4,1)=-(1-delta)*exp(k)/(gamma*exp(kp1)-(1-delta)*exp(k)); di_dzp(8,1)=gamma*exp(kp1)/(gamma*exp(kp1)-(1-delta)*exp(k)); di_dzp(11,1)=exp(k)/(gamma*exp(kp1)-(1-delta)*exp(k)); dq_dzp(16,1)=1/h1; dq_dzp(2,1)=1; dq_dzp(4,1)=1-alpha; dq_dzp(5,1)=alpha; dq_dzp(7,1)=alpha; dq_dzp(9,1)=(-exp(A)*... exp(k)^(1-alpha)*log(exp(k))*h1*(exp(n1)*... exp(w))^alpha+exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*log(exp(n1)*... exp(w)))/exp(A)/(exp(k)^(1-alpha))/h1/((exp(n1)*... exp(w))^alpha); dc_dzp(16,1)=exp(A)*... exp(k)^(1-alpha)*(exp(n1)*exp(w))^alpha/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g)); dc_dzp(2,1)=exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g)); dc_dzp(3,1)=-exp(g)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g)); dc_dzp(19,1)=-exp(kp1)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g)); dc_dzp(4,1)=(exp(A)*... exp(k)^(1-alpha)*(1-alpha)*h1*(exp(n1)*... exp(w))^alpha+(1-delta)*exp(k))/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g)); dc_dzp(5,1)=exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... alpha/(exp(A)*exp(k)^(1-alpha)*h1*... (exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g)); dc_dzp(7,1)=exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*alpha/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g)); dc_dzp(8,1)=-gamma*... exp(kp1)/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g)); dc_dzp(9,1)=(-exp(A)*exp(k)^(1-alpha)*... log(exp(k))*h1*(exp(n1)*exp(w))^alpha+exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... log(exp(n1)*exp(w)))/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g)); dc_dzp(11,1)=-exp(k)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g)); d2h_dzp2(16, 16)=-1/h1^2; d2i_dzp2(19, 11)=-exp(kp1)/(gamma*... exp(kp1)-(1-delta)*exp(k))^2*exp(k); d2i_dzp2(4, 8)=gamma*... exp(kp1)/(gamma*exp(kp1)-(1-delta)*... exp(k))^2*(1-delta)*exp(k); d2i_dzp2(8, 8)=gamma*... exp(kp1)/(gamma*exp(kp1)-(1-delta)*... exp(k))-gamma^2*exp(kp1)^2/(gamma*... exp(kp1)-(1-delta)*exp(k))^2; d2i_dzp2(8, 19)=exp(kp1)/(gamma*... exp(kp1)-(1-delta)*exp(k))-exp(kp1)^2/(gamma*... exp(kp1)-(1-delta)*exp(k))^2*gamma; d2i_dzp2(19, 4)=exp(kp1)/(gamma*... exp(kp1)-(1-delta)*exp(k))^2*(1-delta)*... exp(k); d2i_dzp2(11, 8)=-exp(k)/(gamma*... exp(kp1)-(1-delta)*exp(k))^2*gamma*... exp(kp1); d2i_dzp2(11, 19)=-exp(kp1)/(gamma*... exp(kp1)-(1-delta)*exp(k))^2*exp(k); d2i_dzp2(19, 19)=-exp(kp1)^2/(gamma*... exp(kp1)-(1-delta)*exp(k))^2; d2i_dzp2(4, 11)=exp(k)/(gamma*... exp(kp1)-(1-delta)*exp(k))+exp(k)^2/(gamma*... exp(kp1)-(1-delta)*exp(k))^2*(1-delta); d2i_dzp2(8, 11)=-exp(k)/(gamma*... exp(kp1)-(1-delta)*exp(k))^2*gamma*... exp(kp1); d2i_dzp2(4, 19)=exp(kp1)/(gamma*... exp(kp1)-(1-delta)*exp(k))^2*(1-delta)*... exp(k); d2i_dzp2(11, 11)=-exp(k)^2/(gamma*... exp(kp1)-(1-delta)*exp(k))^2; d2i_dzp2(19, 8)=exp(kp1)/(gamma*... exp(kp1)-(1-delta)*exp(k))-exp(kp1)^2/(gamma*... exp(kp1)-(1-delta)*exp(k))^2*gamma; d2i_dzp2(4, 4)=-(1-delta)*... exp(k)/(gamma*exp(kp1)-(1-delta)*... exp(k))-(1-delta)^2*exp(k)^2/(gamma*... exp(kp1)-(1-delta)*exp(k))^2; d2i_dzp2(8, 4)=gamma*... exp(kp1)/(gamma*exp(kp1)-(1-delta)*... exp(k))^2*(1-delta)*exp(k); d2i_dzp2(11, 4)=exp(k)/(gamma*... exp(kp1)-(1-delta)*exp(k))+exp(k)^2/(gamma*... exp(kp1)-(1-delta)*exp(k))^2*(1-delta); d2q_dzp2(9, 7)=1; d2q_dzp2(9, 9)=(exp(A)*... exp(k)^(1-alpha)*log(exp(k))^2*h1*... (exp(n1)*exp(w))^alpha-2*exp(A)*exp(k)^(1-alpha)*... log(exp(k))*h1*(exp(n1)*exp(w))^alpha*... log(exp(n1)*exp(w))+exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*log(exp(n1)*... exp(w))^2)/exp(A)/(exp(k)^(1-alpha))/h1/((exp(n1)*... exp(w))^alpha)+(-exp(A)*exp(k)^(1-alpha)*... log(exp(k))*h1*(exp(n1)*exp(w))^alpha+exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... log(exp(n1)*exp(w)))/exp(A)/(exp(k)^(1-alpha))/h1/((exp(n1)*... exp(w))^alpha)*log(exp(k))-(-exp(A)*... exp(k)^(1-alpha)*log(exp(k))*h1*(exp(n1)*... exp(w))^alpha+exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*log(exp(n1)*... exp(w)))/exp(A)/(exp(k)^(1-alpha))/h1/((exp(n1)*... exp(w))^alpha)*log(exp(n1)*exp(w)); d2q_dzp2(7, 9)=(-exp(A)*... exp(k)^(1-alpha)*log(exp(k))*h1*(exp(n1)*... exp(w))^alpha*alpha+exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*alpha*log(exp(n1)*... exp(w))+exp(A)*exp(k)^(1-alpha)*h1*... (exp(n1)*exp(w))^alpha)/exp(A)/(exp(k)^(1-alpha))/h1/((exp(n1)*... exp(w))^alpha)-(-exp(A)*exp(k)^(1-alpha)*... log(exp(k))*h1*(exp(n1)*exp(w))^alpha+exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... log(exp(n1)*exp(w)))/exp(A)/(exp(k)^(1-alpha))/h1/((exp(n1)*... exp(w))^alpha)*alpha; d2q_dzp2(4, 9)=(-exp(A)*... exp(k)^(1-alpha)*(1-alpha)*log(exp(k))*... h1*(exp(n1)*exp(w))^alpha-exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha+exp(A)*... exp(k)^(1-alpha)*(1-alpha)*h1*(exp(n1)*... exp(w))^alpha*log(exp(n1)*... exp(w)))/exp(A)/(exp(k)^(1-alpha))/h1/((exp(n1)*... exp(w))^alpha)-(-exp(A)*exp(k)^(1-alpha)*... log(exp(k))*h1*(exp(n1)*exp(w))^alpha+exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... log(exp(n1)*exp(w)))/exp(A)/(exp(k)^(1-alpha))/h1/((exp(n1)*... exp(w))^alpha)*(1-alpha); d2q_dzp2(16, 16)=-1/h1^2; d2q_dzp2(16, 9)=(-exp(A)*... exp(k)^(1-alpha)*log(exp(k))*(exp(n1)*... exp(w))^alpha+exp(A)*exp(k)^(1-alpha)*... (exp(n1)*exp(w))^alpha*log(exp(n1)*... exp(w)))/exp(A)/(exp(k)^(1-alpha))/h1/((exp(n1)*... exp(w))^alpha)-(-exp(A)*exp(k)^(1-alpha)*... log(exp(k))*h1*(exp(n1)*exp(w))^alpha+exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... log(exp(n1)*exp(w)))/exp(A)/(exp(k)^(1-alpha))/h1^2/((exp(n1)*... exp(w))^alpha); d2q_dzp2(5, 9)=(-exp(A)*... exp(k)^(1-alpha)*log(exp(k))*h1*(exp(n1)*... exp(w))^alpha*alpha+exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*alpha*log(exp(n1)*... exp(w))+exp(A)*exp(k)^(1-alpha)*h1*... (exp(n1)*exp(w))^alpha)/exp(A)/(exp(k)^(1-alpha))/h1/((exp(n1)*... exp(w))^alpha)-(-exp(A)*exp(k)^(1-alpha)*... log(exp(k))*h1*(exp(n1)*exp(w))^alpha+exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... log(exp(n1)*exp(w)))/exp(A)/(exp(k)^(1-alpha))/h1/((exp(n1)*... exp(w))^alpha)*alpha; d2q_dzp2(9, 4)=-1; d2q_dzp2(9, 5)=1; d2c_dzp2(4, 5)=exp(A)*... exp(k)^(1-alpha)*(1-alpha)*h1*(exp(n1)*... exp(w))^alpha*alpha/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*alpha/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... (exp(A)*exp(k)^(1-alpha)*(1-alpha)*... h1*(exp(n1)*exp(w))^alpha+(1-delta)*... exp(k)); d2c_dzp2(8, 5)=gamma*... exp(kp1)/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*alpha; d2c_dzp2(16, 4)=exp(A)*... exp(k)^(1-alpha)*(1-alpha)*(exp(n1)*... exp(w))^alpha/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-exp(A)*exp(k)^(1-alpha)*... (exp(n1)*exp(w))^alpha/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*(exp(A)*exp(k)^(1-alpha)*... (1-alpha)*h1*(exp(n1)*exp(w))^alpha+(1-delta)*... exp(k)); d2c_dzp2(5, 4)=exp(A)*... exp(k)^(1-alpha)*(1-alpha)*h1*(exp(n1)*... exp(w))^alpha*alpha/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*alpha/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... (exp(A)*exp(k)^(1-alpha)*(1-alpha)*... h1*(exp(n1)*exp(w))^alpha+(1-delta)*... exp(k)); d2c_dzp2(2, 4)=exp(A)*... exp(k)^(1-alpha)*(1-alpha)*h1*(exp(n1)*... exp(w))^alpha/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-(exp(A)*exp(k)^(1-alpha)*... (1-alpha)*h1*(exp(n1)*exp(w))^alpha+(1-delta)*... exp(k))/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha; d2c_dzp2(7, 16)=exp(A)*... exp(k)^(1-alpha)*(exp(n1)*exp(w))^alpha*... alpha/(exp(A)*exp(k)^(1-alpha)*h1*... (exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-exp(A)^2*(exp(k)^(1-alpha))^2*... ((exp(n1)*exp(w))^alpha)^2/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... h1*alpha; d2c_dzp2(3, 5)=exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... alpha/(exp(A)*exp(k)^(1-alpha)*h1*... (exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*exp(g); d2c_dzp2(11, 5)=exp(k)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(A)*exp(k)^(1-alpha)*h1*(exp(n1)*... exp(w))^alpha*alpha; d2c_dzp2(9, 7)=-exp(A)*... exp(k)^(1-alpha)*log(exp(k))*h1*(exp(n1)*... exp(w))^alpha*alpha/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))+exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*log(exp(n1)*... exp(w))*alpha/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))+exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))-(-exp(A)*... exp(k)^(1-alpha)*log(exp(k))*h1*(exp(n1)*... exp(w))^alpha+exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*log(exp(n1)*... exp(w)))/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*alpha; d2c_dzp2(19, 11)=-exp(kp1)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(k); d2c_dzp2(16, 5)=exp(A)*... exp(k)^(1-alpha)*(exp(n1)*exp(w))^alpha*... alpha/(exp(A)*exp(k)^(1-alpha)*h1*... (exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-exp(A)^2*(exp(k)^(1-alpha))^2*... ((exp(n1)*exp(w))^alpha)^2/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... h1*alpha; d2c_dzp2(5, 5)=exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... alpha^2/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-exp(A)^2*(exp(k)^(1-alpha))^2*... h1^2*((exp(n1)*exp(w))^alpha)^2*alpha^2/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2; d2c_dzp2(2, 5)=exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... alpha/(exp(A)*exp(k)^(1-alpha)*h1*... (exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-exp(A)^2*(exp(k)^(1-alpha))^2*... h1^2*((exp(n1)*exp(w))^alpha)^2*alpha/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2; d2c_dzp2(3, 16)=exp(A)*... exp(k)^(1-alpha)*(exp(n1)*exp(w))^alpha/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(g); d2c_dzp2(11, 16)=exp(A)*exp(k)^(1-alpha)*... (exp(n1)*exp(w))^alpha/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*exp(k); d2c_dzp2(9, 8)=(-exp(A)*... exp(k)^(1-alpha)*log(exp(k))*h1*(exp(n1)*... exp(w))^alpha+exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*log(exp(n1)*... exp(w)))/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*gamma*exp(kp1); d2c_dzp2(7, 7)=exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... alpha^2/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-exp(A)^2*(exp(k)^(1-alpha))^2*... h1^2*((exp(n1)*exp(w))^alpha)^2*alpha^2/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2; d2c_dzp2(19, 2)=exp(kp1)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(A)*exp(k)^(1-alpha)*h1*(exp(n1)*... exp(w))^alpha; d2c_dzp2(4, 7)=exp(A)*... exp(k)^(1-alpha)*(1-alpha)*h1*(exp(n1)*... exp(w))^alpha*alpha/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*alpha/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... (exp(A)*exp(k)^(1-alpha)*(1-alpha)*... h1*(exp(n1)*exp(w))^alpha+(1-delta)*... exp(k)); d2c_dzp2(8, 7)=gamma*... exp(kp1)/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*alpha; d2c_dzp2(5, 16)=exp(A)*... exp(k)^(1-alpha)*(exp(n1)*exp(w))^alpha*... alpha/(exp(A)*exp(k)^(1-alpha)*h1*... (exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-exp(A)^2*(exp(k)^(1-alpha))^2*... ((exp(n1)*exp(w))^alpha)^2/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... h1*alpha; d2c_dzp2(2, 16)=exp(A)*exp(k)^(1-alpha)*... (exp(n1)*exp(w))^alpha/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-exp(A)^2*(exp(k)^(1-alpha))^2*... ((exp(n1)*exp(w))^alpha)^2/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... h1; d2c_dzp2(9, 9)=(exp(A)*... exp(k)^(1-alpha)*log(exp(k))^2*h1*... (exp(n1)*exp(w))^alpha-2*exp(A)*exp(k)^(1-alpha)*... log(exp(k))*h1*(exp(n1)*exp(w))^alpha*... log(exp(n1)*exp(w))+exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*log(exp(n1)*... exp(w))^2)/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-(-exp(A)*exp(k)^(1-alpha)*... log(exp(k))*h1*(exp(n1)*exp(w))^alpha+exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... log(exp(n1)*exp(w)))^2/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2; d2c_dzp2(9, 19)=exp(kp1)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... (-exp(A)*exp(k)^(1-alpha)*log(exp(k))*... h1*(exp(n1)*exp(w))^alpha+exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... log(exp(n1)*exp(w))); d2c_dzp2(7, 8)=gamma*... exp(kp1)/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*alpha; d2c_dzp2(19, 3)=-exp(kp1)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(g); d2c_dzp2(3, 7)=exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*alpha/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(g); d2c_dzp2(11, 7)=exp(k)/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*alpha; d2c_dzp2(4, 8)=gamma*... exp(kp1)/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*(exp(A)*exp(k)^(1-alpha)*... (1-alpha)*h1*(exp(n1)*exp(w))^alpha+(1-delta)*... exp(k)); d2c_dzp2(8, 8)=-gamma*... exp(kp1)/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-gamma^2*exp(kp1)^2/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2; d2c_dzp2(4, 16)=exp(A)*... exp(k)^(1-alpha)*(1-alpha)*(exp(n1)*... exp(w))^alpha/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-exp(A)*exp(k)^(1-alpha)*... (exp(n1)*exp(w))^alpha/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*(exp(A)*exp(k)^(1-alpha)*... (1-alpha)*h1*(exp(n1)*exp(w))^alpha+(1-delta)*... exp(k)); d2c_dzp2(16, 7)=exp(A)*... exp(k)^(1-alpha)*(exp(n1)*exp(w))^alpha*... alpha/(exp(A)*exp(k)^(1-alpha)*h1*... (exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-exp(A)^2*(exp(k)^(1-alpha))^2*... ((exp(n1)*exp(w))^alpha)^2/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... h1*alpha; d2c_dzp2(8, 19)=-exp(kp1)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))-exp(kp1)^2/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... gamma; d2c_dzp2(5, 7)=exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*alpha^2/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))-exp(A)^2*... (exp(k)^(1-alpha))^2*h1^2*((exp(n1)*... exp(w))^alpha)^2*alpha^2/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2; d2c_dzp2(8, 16)=exp(A)*... exp(k)^(1-alpha)*(exp(n1)*exp(w))^alpha/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... gamma*exp(kp1); d2c_dzp2(2, 7)=exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... alpha/(exp(A)*exp(k)^(1-alpha)*h1*... (exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-exp(A)^2*(exp(k)^(1-alpha))^2*... h1^2*((exp(n1)*exp(w))^alpha)^2*alpha/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2; d2c_dzp2(7, 9)=(-exp(A)*... exp(k)^(1-alpha)*log(exp(k))*h1*(exp(n1)*... exp(w))^alpha*alpha+exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*alpha*log(exp(n1)*... exp(w))+exp(A)*exp(k)^(1-alpha)*h1*... (exp(n1)*exp(w))^alpha)/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-(-exp(A)*exp(k)^(1-alpha)*... log(exp(k))*h1*(exp(n1)*exp(w))^alpha+exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... log(exp(n1)*exp(w)))/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*alpha; d2c_dzp2(7, 19)=exp(kp1)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(A)*exp(k)^(1-alpha)*h1*(exp(n1)*... exp(w))^alpha*alpha; d2c_dzp2(19, 4)=exp(kp1)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... (exp(A)*exp(k)^(1-alpha)*(1-alpha)*... h1*(exp(n1)*exp(w))^alpha+(1-delta)*... exp(k)); d2c_dzp2(3, 8)=-gamma*exp(kp1)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(g); d2c_dzp2(11, 8)=-exp(k)/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*gamma*exp(kp1); d2c_dzp2(4, 9)=(-exp(A)*... exp(k)^(1-alpha)*(1-alpha)*log(exp(k))*... h1*(exp(n1)*exp(w))^alpha-exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha+exp(A)*... exp(k)^(1-alpha)*(1-alpha)*h1*(exp(n1)*... exp(w))^alpha*log(exp(n1)*exp(w)))/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))-(-exp(A)*... exp(k)^(1-alpha)*log(exp(k))*h1*(exp(n1)*... exp(w))^alpha+exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*log(exp(n1)*... exp(w)))/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*(exp(A)*exp(k)^(1-alpha)*... (1-alpha)*h1*(exp(n1)*exp(w))^alpha+(1-delta)*... exp(k)); d2c_dzp2(8, 9)=(-exp(A)*... exp(k)^(1-alpha)*log(exp(k))*h1*(exp(n1)*... exp(w))^alpha+exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*log(exp(n1)*... exp(w)))/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*gamma*exp(kp1); d2c_dzp2(16, 8)=exp(A)*... exp(k)^(1-alpha)*(exp(n1)*exp(w))^alpha/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... gamma*exp(kp1); d2c_dzp2(16, 16)=-exp(A)^2*... (exp(k)^(1-alpha))^2*((exp(n1)*exp(w))^alpha)^2/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2; d2c_dzp2(5, 8)=gamma*... exp(kp1)/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*alpha; d2c_dzp2(9, 11)=exp(k)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... (-exp(A)*exp(k)^(1-alpha)*log(exp(k))*... h1*(exp(n1)*exp(w))^alpha+exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... log(exp(n1)*exp(w))); d2c_dzp2(2, 8)=gamma*... exp(kp1)/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha; d2c_dzp2(19, 5)=exp(kp1)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(A)*exp(k)^(1-alpha)*h1*(exp(n1)*... exp(w))^alpha*alpha; d2c_dzp2(3, 9)=(-exp(A)*... exp(k)^(1-alpha)*log(exp(k))*h1*(exp(n1)*... exp(w))^alpha+exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*log(exp(n1)*... exp(w)))/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*exp(g); d2c_dzp2(11, 9)=exp(k)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... (-exp(A)*exp(k)^(1-alpha)*log(exp(k))*... h1*(exp(n1)*exp(w))^alpha+exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... log(exp(n1)*exp(w))); d2c_dzp2(3, 19)=-exp(kp1)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(g); d2c_dzp2(11, 19)=-exp(kp1)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(k); d2c_dzp2(19, 16)=exp(kp1)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(A)*exp(k)^(1-alpha)*(exp(n1)*... exp(w))^alpha; d2c_dzp2(16, 9)=(-exp(A)*... exp(k)^(1-alpha)*log(exp(k))*(exp(n1)*... exp(w))^alpha+exp(A)*exp(k)^(1-alpha)*... (exp(n1)*exp(w))^alpha*log(exp(n1)*... exp(w)))/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-exp(A)*exp(k)^(1-alpha)*... (exp(n1)*exp(w))^alpha/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*(-exp(A)*exp(k)^(1-alpha)*... log(exp(k))*h1*(exp(n1)*exp(w))^alpha+exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... log(exp(n1)*exp(w))); d2c_dzp2(9, 2)=-exp(A)*... exp(k)^(1-alpha)*log(exp(k))*h1*(exp(n1)*... exp(w))^alpha/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))+exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*log(exp(n1)*... exp(w))/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-(-exp(A)*exp(k)^(1-alpha)*... log(exp(k))*h1*(exp(n1)*exp(w))^alpha+exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... log(exp(n1)*exp(w)))/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha; d2c_dzp2(19, 19)=-exp(kp1)^2/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2; d2c_dzp2(5, 9)=(-exp(A)*... exp(k)^(1-alpha)*log(exp(k))*h1*(exp(n1)*... exp(w))^alpha*alpha+exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*alpha*log(exp(n1)*... exp(w))+exp(A)*exp(k)^(1-alpha)*h1*... (exp(n1)*exp(w))^alpha)/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-(-exp(A)*exp(k)^(1-alpha)*... log(exp(k))*h1*(exp(n1)*exp(w))^alpha+exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... log(exp(n1)*exp(w)))/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*alpha; d2c_dzp2(5, 19)=exp(kp1)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(A)*exp(k)^(1-alpha)*h1*(exp(n1)*... exp(w))^alpha*alpha; d2c_dzp2(2, 9)=(-exp(A)*... exp(k)^(1-alpha)*log(exp(k))*h1*(exp(n1)*... exp(w))^alpha+exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*log(exp(n1)*... exp(w)))/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-(-exp(A)*exp(k)^(1-alpha)*... log(exp(k))*h1*(exp(n1)*exp(w))^alpha+exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... log(exp(n1)*exp(w)))/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha; d2c_dzp2(7, 11)=exp(k)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(A)*exp(k)^(1-alpha)*h1*(exp(n1)*... exp(w))^alpha*alpha; d2c_dzp2(2, 19)=exp(kp1)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(A)*exp(k)^(1-alpha)*h1*(exp(n1)*... exp(w))^alpha; d2c_dzp2(4, 11)=-exp(k)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))+exp(k)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... (exp(A)*exp(k)^(1-alpha)*(1-alpha)*... h1*(exp(n1)*exp(w))^alpha+(1-delta)*... exp(k)); d2c_dzp2(8, 11)=-exp(k)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... gamma*exp(kp1); d2c_dzp2(4, 19)=exp(kp1)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... (exp(A)*exp(k)^(1-alpha)*(1-alpha)*... h1*(exp(n1)*exp(w))^alpha+(1-delta)*... exp(k)); d2c_dzp2(9, 3)=(-exp(A)*exp(k)^(1-alpha)*... log(exp(k))*h1*(exp(n1)*exp(w))^alpha+exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... log(exp(n1)*exp(w)))/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*exp(g); d2c_dzp2(7, 2)=exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... alpha/(exp(A)*exp(k)^(1-alpha)*h1*... (exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-exp(A)^2*(exp(k)^(1-alpha))^2*... h1^2*((exp(n1)*exp(w))^alpha)^2*alpha/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2; d2c_dzp2(16, 2)=exp(A)*... exp(k)^(1-alpha)*(exp(n1)*exp(w))^alpha/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))-exp(A)^2*... (exp(k)^(1-alpha))^2*((exp(n1)*exp(w))^alpha)^2/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... h1; d2c_dzp2(3, 11)=-exp(k)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(g); d2c_dzp2(11, 11)=-exp(k)^2/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2; d2c_dzp2(4, 2)=exp(A)*... exp(k)^(1-alpha)*(1-alpha)*h1*(exp(n1)*... exp(w))^alpha/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-(exp(A)*exp(k)^(1-alpha)*... (1-alpha)*h1*(exp(n1)*exp(w))^alpha+(1-delta)*... exp(k))/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha; d2c_dzp2(8, 2)=gamma*... exp(kp1)/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha; d2c_dzp2(19, 7)=exp(kp1)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(A)*exp(k)^(1-alpha)*h1*(exp(n1)*... exp(w))^alpha*alpha; d2c_dzp2(16, 11)=exp(A)*... exp(k)^(1-alpha)*(exp(n1)*exp(w))^alpha/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(k); d2c_dzp2(9, 4)=(-exp(A)*exp(k)^(1-alpha)*... (1-alpha)*log(exp(k))*h1*(exp(n1)*... exp(w))^alpha-exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha+exp(A)*... exp(k)^(1-alpha)*(1-alpha)*h1*(exp(n1)*... exp(w))^alpha*log(exp(n1)*exp(w)))/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))-(-exp(A)*... exp(k)^(1-alpha)*log(exp(k))*h1*(exp(n1)*... exp(w))^alpha+exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*log(exp(n1)*... exp(w)))/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*(exp(A)*exp(k)^(1-alpha)*... (1-alpha)*h1*(exp(n1)*exp(w))^alpha+(1-delta)*... exp(k)); d2c_dzp2(5, 11)=exp(k)/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*alpha; d2c_dzp2(7, 3)=exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... alpha/(exp(A)*exp(k)^(1-alpha)*h1*... (exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*exp(g); d2c_dzp2(2, 11)=exp(k)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(A)*exp(k)^(1-alpha)*h1*(exp(n1)*... exp(w))^alpha; d2c_dzp2(3, 2)=exp(g)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(A)*exp(k)^(1-alpha)*h1*(exp(n1)*... exp(w))^alpha; d2c_dzp2(11, 2)=exp(k)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(A)*exp(k)^(1-alpha)*h1*(exp(n1)*... exp(w))^alpha; d2c_dzp2(4, 3)=(exp(A)*... exp(k)^(1-alpha)*(1-alpha)*h1*(exp(n1)*... exp(w))^alpha+(1-delta)*exp(k))/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(g); d2c_dzp2(8, 3)=-gamma*... exp(kp1)/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*exp(g); d2c_dzp2(19, 8)=-exp(kp1)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))-exp(kp1)^2/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... gamma; d2c_dzp2(5, 2)=exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... alpha/(exp(A)*exp(k)^(1-alpha)*h1*... (exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-exp(A)^2*(exp(k)^(1-alpha))^2*... h1^2*((exp(n1)*exp(w))^alpha)^2*alpha/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2; d2c_dzp2(9, 5)=-exp(A)*... exp(k)^(1-alpha)*log(exp(k))*h1*(exp(n1)*... exp(w))^alpha*alpha/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))+exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*log(exp(n1)*... exp(w))*alpha/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))+exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))-(-exp(A)*... exp(k)^(1-alpha)*log(exp(k))*h1*(exp(n1)*... exp(w))^alpha+exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*log(exp(n1)*... exp(w)))/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*alpha; d2c_dzp2(2, 2)=exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))-exp(A)^2*... (exp(k)^(1-alpha))^2*h1^2*((exp(n1)*... exp(w))^alpha)^2/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2; d2c_dzp2(7, 4)=exp(A)*exp(k)^(1-alpha)*... (1-alpha)*h1*(exp(n1)*exp(w))^alpha*... alpha/(exp(A)*exp(k)^(1-alpha)*h1*... (exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*alpha/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... (exp(A)*exp(k)^(1-alpha)*(1-alpha)*... h1*(exp(n1)*exp(w))^alpha+(1-delta)*... exp(k)); d2c_dzp2(3, 3)=-exp(g)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))-exp(g)^2/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2; d2c_dzp2(11, 3)=-exp(k)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(g); d2c_dzp2(4, 4)=(exp(A)*... exp(k)^(1-alpha)*(1-alpha)^2*h1*(exp(n1)*... exp(w))^alpha+(1-delta)*exp(k))/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))-(exp(A)*... exp(k)^(1-alpha)*(1-alpha)*h1*(exp(n1)*... exp(w))^alpha+(1-delta)*exp(k))^2/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2; d2c_dzp2(16, 19)=exp(kp1)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(A)*exp(k)^(1-alpha)*(exp(n1)*... exp(w))^alpha; d2c_dzp2(8, 4)=gamma*... exp(kp1)/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*(exp(A)*exp(k)^(1-alpha)*... (1-alpha)*h1*(exp(n1)*exp(w))^alpha+(1-delta)*... exp(k)); d2c_dzp2(19, 9)=exp(kp1)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... (-exp(A)*exp(k)^(1-alpha)*log(exp(k))*... h1*(exp(n1)*exp(w))^alpha+exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... log(exp(n1)*exp(w))); d2c_dzp2(16, 3)=exp(A)*... exp(k)^(1-alpha)*(exp(n1)*exp(w))^alpha/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(g); d2c_dzp2(5, 3)=exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... alpha/(exp(A)*exp(k)^(1-alpha)*h1*... (exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*exp(g); d2c_dzp2(9, 16)=-exp(A)*... exp(k)^(1-alpha)*log(exp(k))*(exp(n1)*... exp(w))^alpha/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))+exp(A)*exp(k)^(1-alpha)*... (exp(n1)*exp(w))^alpha*log(exp(n1)*... exp(w))/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))-exp(A)*exp(k)^(1-alpha)*... (exp(n1)*exp(w))^alpha/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*(-exp(A)*exp(k)^(1-alpha)*... log(exp(k))*h1*(exp(n1)*exp(w))^alpha+exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha*... log(exp(n1)*exp(w))); d2c_dzp2(2, 3)=exp(g)/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(A)*exp(k)^(1-alpha)*h1*(exp(n1)*... exp(w))^alpha; d2c_dzp2(7, 5)=exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha*alpha^2/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))-exp(A)^2*... (exp(k)^(1-alpha))^2*h1^2*((exp(n1)*... exp(w))^alpha)^2*alpha^2/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2; d2c_dzp2(3, 4)=(exp(A)*... exp(k)^(1-alpha)*(1-alpha)*h1*(exp(n1)*... exp(w))^alpha+(1-delta)*exp(k))/(exp(A)*... exp(k)^(1-alpha)*h1*(exp(n1)*exp(w))^alpha-gamma*... exp(kp1)+(1-delta)*exp(k)-exp(g))^2*... exp(g); d2c_dzp2(11, 4)=-exp(k)/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))+exp(k)/(exp(A)*exp(k)^(1-alpha)*... h1*(exp(n1)*exp(w))^alpha-gamma*exp(kp1)+(1-delta)*... exp(k)-exp(g))^2*(exp(A)*exp(k)^(1-alpha)*... (1-alpha)*h1*(exp(n1)*exp(w))^alpha+(1-delta)*... exp(k));