read MVPoisson:

A:=[[1,  0,  0,  1,  1], [0,  1,  1,  1,  1], [0, 1, 1, 0, 1]]:
RecsV(A,g,b):

quit: