Beyond curse of dimensionality with transition matrix

We now set the aggregate model with A and Q changing according to transition matrix.

We set

fc(K,K1,A,Q)=A*K^alpha-(K1-(1-delta)*K)/Q

rather than

fc(k,k1,theta)=theta*A*k^alpha+(1-delta)*k-k1

for non-aggregate model.

We also change the value of kstar by Kstar.

All we have to modify is to change the calling part of fc in terms of A[j] and Q[j] rather

than theta[j].