agisga

library(dplyr)
library(tidyr)
library(ggplot2)
library(rTensor)

#--- generate data

# bell-shaped spacial component with different means
space_index <- seq(-1, 1, l = 100)
case1 <- matrix(rep(dnorm(space_index, mean = 0, sd = 0.3), 10), 100, 100)