Collocation method with transition matrix

We modify the program to aggregate model and others in general as well.

All we have to do is change the functional setting of fc and the related part to call.

Originally, it is fc(k,k1,theta). We change it to fc(K,K1,A,Q) or fc(K,K1,para).

Here, para consists of A and Q horizontally merged. If we consider non-aggregate

model, para=theta.