Gaussian Distribution

Translated by Will Secor

Based on: Gaussian Distribution by Kerry Jones

Category: direct

        
/*
Part of the ReCode Project (http://recodeproject.com)
Based on "Gaussian Distribution" by Kerry Jones
Originally published in "Computer Graphics and Art" v2n2, 1977
Copyright (c) 2013 Will Secor - OSI/MIT license (http://recodeproject/license).
*/

let w = 760;
let h = 330;

function setup() {
  createCanvas(w, h);
  noLoop();
}

function draw() {
  background(255, 255, 255);
  fill(0, 0, 0);
  push();
  translate(0, h / 2);
  boxmuller();
  translate(w / 4, 0);
  boxmuller();
  translate(w / 4, 0);
  boxmuller();
  translate(w / 4, 0);
  boxmuller();
  translate(w / 4, 0);
  boxmuller();
  pop();
}

// Box-Muller transform!  This still seems like magic to me.

function boxmuller() {
  for (let j = 0; j < 5500; j++) {
    let u1 = random(0, 1);
    let u2 = random(0, 1);
    let z0 = sqrt(-2 * log(u1)) * cos(2 * PI * u2);
    let z1 = sqrt(-2 * log(u1)) * sin(2 * PI * u2);
    ellipse(z0 * w / (1.8 * PI), z1 * h / (1.8 * PI), 2, 2);
  }
}

function mousePressed() {
  redraw();
}