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

 //初期パラメータの指定
int sq_size = 600; // 正方形のサイズ
int r = 40; // ランダムウォークのサイズ
int step = 100; // ランダムウォークのステップ数
int num_rw = 3;  // ランダムウォークの個数


size(600, 600); // 図の大きさを指定

int center = sq_size/2; // 中央の座標
int c_size = r*2; // 円の大きさ

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);	//輪郭線の色を指定
    line(ix,iy,tx,ty);
    ix = tx; iy = ty;
    theta =random(-PI,PI);
    tx = ix + r*cos(theta);
    ty = iy + r*sin(theta);
  }
}

/*
Rayleigh.tifという名前で保存。
Sketch->Show Sketch Folder で保存場所を確認
saveFrame("Rayleigh.tif");
*/