function G=optpath2(k,a,b,gamma,beta,teta,kc1,kk1,kcT,kkT,T); kc=k(:,1); kk=k(:,2); g=(a*kc+kk).^(1-teta); A=4; for t=1:T-1; ic=gamma*kc(t+1)-b*kc(t); G(A*(t-1)+1)=ic*(ic<0)*100; ik=gamma*kk(t+1)-b*kk(t); G(A*(t-1)+2)=ik*(ik<0)*100; G(A*(t-1)+3)=gamma*g(t)/kc(t)-beta*(teta+b*g(t+1))/kc(t+1); if G(A*(t-1)+1)==0 & G(A*(t-1)+2)==0 G(A*(t-1)+4)=kk(t)/g(t)+b*(kc(t)+kk(t))-gamma*(kc(t+1)+kk(t+1)); elseif G(A*(t-1)+1)<0 G(A*(t-1)+4)=kk(t)/g(t)+b*kk(t)-gamma*kk(t+1); elseif G(A*(t-1)+2)<0 G(A*(t-1)+4)=kk(t)/g(t)+b*kc(t)-gamma*kc(t+1); end end G(A*(T-1)+1)=100000*(kk(1)-kk1); G(A*(T-1)+2)=100000*(kc(1)-kc1); %G(A*(T-1)+3)=100000*(kk(T)-kkT); %G(A*(T-1)+4)=100000*(kc(T)-kcT);