FTCS.java

for(int i=0;i<F.length;++i) {
  for(int j=0;j<F[0].length;++j) {
    F_next[i][j] = F[i][j]+deltaT*diffusionRate/deltaXSq*(F[(i+1)%F.length][j]+
                                                          F[(i-1+F.length)%F.length][j]+
                                                          F[i][(j+1)%F[0].length]+
                                                          F[i][(j-1+F[0].length)%F[0].length]-
                                                          4*F[i][j]);
  }
}