function coupled_langevin(f, x, delta) {
var noise = [Gaussian(), Gaussian()];
var next = [];
for (var i = 0; i < x.length; ++i) {
var d = df(x[i]);
next.push([x[i][0] - delta * d[0] + Math.sqrt(2*delta) * noise[0],
x[i][1] - delta * d[1] + Math.sqrt(2*delta) * noise[1]]);
}
return next;
}