想做一个矩阵乘方的程序。有两个序列a<-c(1,2,3), b<-c(1,2),现在想做 c=a^b, 最后得到一个3*2矩阵c=(1,2,3, 1,4,9)，程序如下

a=c(1,2,3)
b=c(1,2)
d<- matrix(, nrow=3, ncol=2)
for (i in 1:3){
d=a^b
}
d后面是中括号i，不知为什么打不上


但是总是得到number of items to replace is not a multiple of replacement length的错误，不知道该怎么修改程序，请大家指点，谢谢啦。