R语言建立矩阵计算双重积分求和

你好同学，代码供参考，如有帮助还望题主给个宝贵的采纳支持答主答题哟，谢谢啦（づ￣3￣）づ╭❤～

constructMat<-function(n) {
  m<-matrix(0, ncol=n, nrow=n)#生成n行n列的0矩阵
  for (i in 1:n){
    for (j in 1:i){
        m[i,j]<-j*j/(10+4*i*i*i)#得到矩阵的每个元素
    }
  }
  return (m);#返回矩阵
}
a <- constructMat(5)
print(a)
fn <- sum(a)
print(fn)#计算最终的fn

输出结果

> a <- constructMat(5)
> print(a)
            [,1]        [,2]       [,3]       [,4]       [,5]
[1,] 0.071428571 0.000000000 0.00000000 0.00000000 0.00000000
[2,] 0.023809524 0.095238095 0.00000000 0.00000000 0.00000000
[3,] 0.008474576 0.033898305 0.07627119 0.00000000 0.00000000
[4,] 0.003759398 0.015037594 0.03383459 0.06015038 0.00000000
[5,] 0.001960784 0.007843137 0.01764706 0.03137255 0.04901961
> fn <- sum(a)
> print(fn)#计算最终的fn
[1] 0.5297454