Found!
BHJ = Table[ T[r[[p + 1]], r[[p]], x[[p + 1]] - x[[p]], 10], {p, 1, 2}] Apply[Dot, BHJ] // MatrixForm
Thanks