Value function iterations with transition matrix

We generalize the program for any model and any parameter. That is, we consider

where

together with

by replacing both k and k’ by x and x’ and in terms of matrix para. Here,

we replace A and Q by para[1] and para[2].

We can set to para alone in the non-aggregate model.