read MVPoisson:

A:=[[1,  0,  0,  1,  1], [0,  1,  1,  1,  1], [0, 1, 1, 0, 1]]:
RecsV(A,[3,2,5,3,8],b):

quit: