Computation Q matrix
q.matrix.diffregr4(b.mat, act.a, act.b, ss)
b.mat
no descr
act.a
act.b
ss
n.stadler