> a <- matrix(1:9, 3,3)
> b <- matrix(1:3, nrow=3)
> a
[,1] [,2] [,3]
[1,] 1 4 7
[2,] 2 5 8
[3,] 3 6 9
> b
[,1]
[1,] 1
[2,] 2
[3,] 3
> dim(a)
[1] 3 3
> dim(b)
[1] 3 1
> a %*% b
[,1]
[1,] 30
[2,] 36
[3,] 42
行列のかけ算は、%*%で無いと駄目らしいです。
道理で、先日、赤道座標から銀河座標への変換の計算がうまくいかなかったわけです。