Skip to content
Prev 261947 / 398503 Next

kronecker sum

On 11-06-05 9:36 AM, Lara Poplarski wrote:
Why not do it in the obvious way?

A <- matrix1[1,2]
a <- matrix1[2,1]
B <- matrix2[1,2]
b <- matrix2[2,1]
matrix(c(-A-B,a,b,0,  A,-a-B,0,b, B,0,-A-b,a, 0,B,A,-a-b), 4,4)

Duncan Murdoch