read MVPoisson:

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

quit: