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);    