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