Cooperativity-driven singularities in asymmetric exclusion (probabilistic version)

asep.js

var lambda = .5;
var p = new Array(10);

for(var i=0; i<p.length;i++) p[i]=0;
p[2]=1;


var niter = 3;

for (var k=0;k<niter;k++){
var p2 = [0,0];
for(var i=2; i<p.length-1;i++) {
  // The current particle stays + the previous one spontaneously bounces + the previous one is pushed
  p2[i] = (p[i]*p[i+1] + (1-p[i+1])*p[i]*(1-p[i-1])*(1-lambda))
       + (1-p[i-2])*p[i-1]*(1-p[i])*lambda
       + p[i-2]*p[i-1]*(1-p[i]);
}
p=p2;
}