Collocation method with transition matrices

We generalize the program for any model in terms of x:

with

for all groups of parameters.

We continue to consider para1, para2, . . . and so on.

There are two kinds in the programming. One is to offset the state between the breaks.

The other is keep the generation of states across breaks.