/*
 processing でRayleigh's Theoremを図示する
(Very) short proof of Rayleigh's Theorem (and extensions)
http://arxiv.org/abs/1007.4870
<Time-stamp: <2019-01-28 17:01:59 nakata>

アニメーションで何回か繰り返す
*/

 //初期パラメータの指定
int sq_size = 600; // 正方形のサイズ
int r = 40; // ランダムウォークのサイズ
int step = 100; // ランダムウォークのステップ数
int num_rw = 2;  // ランダムウォークの個数 
int center = sq_size/2; // 中央の座標
int c_size = r*2; // 円の大きさ
int stroke_wei = 9; // 線の太さ

void setup() {
  size(600,600); // 図の大きさを指定
  frameRate(20); //1秒間に20回リフレッシュする
}



void draw() {
  // draw が呼びだされるたび色をリセットする	0 黒、255 白
  background(255); 
  ellipse(center,center,c_size,c_size); // はじめに単位円を書く

  float theta, ix, iy, tx, ty; // 角度、始点、終点
  theta =random(-PI,PI); // 0から-πからπ=3.14までの一様乱数

  for(int i=0; i<num_rw; i++) {
    ix = center; iy = center;
    tx = ix + r*cos(theta);
    ty = iy + r*sin(theta);
    int fr = (int)random(256); //乱数で色を決める
    int fg = (int)random(256);
    int fb = (int)random(256);
    
    for(int j=0; j<step; j++) {  
      stroke(fr, fg, fb);	//輪郭線の色を指定
      strokeWeight(stroke_wei); //線の太さを指定
      line(ix,iy,tx,ty);
      ix = tx; iy = ty;
      theta =random(-PI,PI);
      tx = ix + r*cos(theta);
      ty = iy + r*sin(theta);
    }

  }
  delay(1000);  // プログラムを1000 ミリ秒停止させる
}