こちらのpythonをRに直す a <- 2.8 * 10^(-4) b <- 5 * 10^(-3) tau <- 0.1 k <- -0.005 size <- 100 dx <- 2/size T <- 10.0 dt <- 0.9 * dx^2/2 n <- floor(T/dt) U <- matrix(runif(size^2),size,size) V <- matrix(runif(size^2),size,size) ctr <- 2:(s…