Dear all,
I am trying to maximize an utility function subject to asset accumulation. I am not able to obtain the asset profile using conopt.
Thanks,
Fabio
utility(ii,t).. sum(j,beta**(ord(j)-1)*pi(j)* (((c(ii,j,t+(ord(j)-1))**(1/1-rho)+alpha*(1-le(ii,j,t+(ord(j)-1)))**(1-1/rho))**(1/(1-1/rho)))**(1-1/gamma))/(1-(1/gamma)))=e=u;
assets(ii,j,t).. (1+r)*as(ii,j-1,t-1)+w(t)*e(ii,j)*(1-le(ii,j,t))+p(ii,j,t)+s(ii,j,t)+b(ii,j,t)-c(ii,j,t)-tax(ii,j,t)=e=as(ii,j,t);
tax_household(ii,j,t).. 0.2*(w(t)*e(ii,j)*(1-le(ii,j,t))+r*as(ii,j-1,t-1))+tau_c(t)*c(ii,j,t)=e=tax(ii,j,t) ;
pension(ii,j,t)$(ord(j) ge 45)..max(min(pmax,pmax-theta*(r*as(ii,j-1,t-1)+0.5*w(t)*e(ii,j)*(1-le(ii,j,t))-ybar)),0)=e=p(ii,j,t) ;
*private pension
super_asset(ii,j,t).. [1+(tau_r)*r]*super_hat(ii,j-1,t-1)+(1-tau_s)*nu*w(t)*e(ii,j)*(1-le(ii,j,t))=e=super_hat(ii,j,t);
super_payout(ii,j,t)$(ord(j) eq 40).. s(ii,j,t)=e=super_hat(ii,j,t);
model house /utility assets tax_household pension super_asset super_payout /;
option dnlp=conopt;
solve house using dnlp maximizing u;
comparing 4.gms (9.73 KB)