在R中处理热数据

Question

我试图用R进行热模拟;我应该尝试在系列内部制作流程。例如，我应该多次运行这个过程：

 T[i,j]=.25*(T[i+1,j]+T[i,j+1]+T[i-1,j]+T[i,j-1])

我知道如何在C或FORTRAN中制作这样的东西，但在R中我仍然不清楚这一点。

Answer 1

设T为10 * 10阵列。迭代是：

T[2:9,2:9] <- (T[1:8,2:9] + T[3:10,2:9] + T[2:9,1:8] + T[2:9,3:10]) / 4
- and similar for edges