直接乘就好了，本来就是矢量化的。

> x <- 1:10
> y <- 2:11
> x
[1]  1  2  3  4  5  6  7  8  9 10
> y
[1]  2  3  4  5  6  7  8  9 10 11
> x * y
[1]   2   6  12  20  30  42  56  72  90 110

啊啊！！~~谢谢~虽然我用对角阵做了…看来我绕路了