I've never studied infinite matrix products, but would have expected that for one to converge, its general term would need to approach the identity. No? (Or perhaps some entries "diverge to 0" ?) --Dan << . . . prod(matrix([-(k+1/2)^3/(8*(k+1)^3),-(k+1/2)^2/(16*(k+1)^2),30*(k+7/30)],[-63*k*(k+1/2)^2/(16*(k+1)^3),0,0],[0,0,1]),k,0,inf) = matrix([0,0,24/%pi],[0,0,0],[0,0,1]) [ 1 3 1 2 ] [ (k + -) (k + -) ] [ 2 2 7 ] [ - ---------- - ----------- 30 (k + --) ] inf [ 3 2 30 ] [ 24 ] /===\ [ 8 (k + 1) 16 (k + 1) ] [ 0 0 --- ] | | [ ] [ %pi ] | | [ 1 2 ] = [ ] | | [ 63 k (k + -) ] [ 0 0 0 ] k = 0 [ 2 ] [ ] [ - ------------- 0 0 ] [ 0 0 1 ] [ 3 ] [ 16 (k + 1) ] [ ] [ 0 0 1 ] . . .

...

Those who sleep faster get more rest.