用Julia、Python、R语言计算方型矩阵的n次方

用Julia、Python、R语言计算矩阵的n次方分别如下
Julia语言，直接用a^n即可：

a = Int[[1,1] [1,0]]
r = a^50
r

Python语言，只有dot函数，需要自定义一个矩阵的n次方函数来计算

import numpy  as np
#自定义函数dotx
def dotx(a,n):
    for i in range(1,n+1):
        if i == 1:
            b = a
        else :
            b = np.dot(a,b)
    return b

R语言也一样，只有"%*%"运算符，需要自定义一个矩阵的n次方函数来计算

#自定义"%^%"
"%^%" = function(a, n) {
  for (i in 1:n) {
    if (i == 1){
        b = a}
    else{
        b = a %*% b}}
    return = b
}

# 计算矩阵a的N次方
a = matrix(c(1,1,1,0),nrow=2,ncol=2) 
b = a %^% 50
b