library(AIUQ)
M = 10
len_t = 50
sigma_fbm = 2
H = 0.3
sigma_ou = 2
rho = 0.95
pos0 = matrix(100/8+0.75*100*runif(M*2),nrow=M,ncol=2)
pos = fbm_ou_particle_intensity(pos0=pos0, M=M, len_t=len_t,
sigma_fbm=sigma_fbm, sigma_ou=sigma_ou, H=H, rho=rho)
Run the code above in your browser using DataLab