read MVPoisson: g:=[1$nops(A2[1])]: b:=[3,4]: print(`The matrix is`, A2, ` and the vector b is `, b): print(`The correlation matrix with parameters`, g, `is : `): print(evalf(CorMf(A2,g,b))): quit: