// BEADS / Increase '99.12
// By Tsutomu HIGO  E-mail: nj2t-hg@asahi-net.or.jp

#version 3.1;
global_settings {
  max_trace_level 6
  assumed_gamma 2.2
  }

light_source {<40, 20, -10> color <1,1,1>}

camera {
  location <-1, -4, -100>
  angle 92
  look_at <1, -1, 0>
  }

#declare Imax=5000;
#declare Cp=array[Imax][2]
#declare I=0; #while (I<Imax)
  #declare J=0; #while (J<2)
    #declare Cp[I][J]=0;
  #declare J=J+1; #end
#declare I=I+1; #end

#macro Acosc (Rc,Dy,Ac)
  #declare Ac=acos(Rc);
  #if (Dy<0) #declare Ac=2*pi-Ac; #end
  #if (Ac=0) #declare Ac=2*pi; #end
#end

#declare I=0; #while (I<6)
  #declare Pe=1+I;
  #declare Cp[Pe][0]=Cp[0][0]+2*cos(pi/3*I);
  #declare Cp[Pe][1]=Cp[0][1]+2*sin(pi/3*I);
#declare I=I+1; #end

#declare Ap=0;
#declare As=0;
#declare Al=0;
#declare Ae=0;
#declare Ad=0;
#declare Rndm = seed(89);

#macro Spculc (N,P,Pe)
  #declare Rs=sqrt(pow(Cp[P][0]-Cp[N][0],2)+pow(Cp[P][1]-Cp[N][1],2))/4;
  Acosc((Cp[P][0]-Cp[N][0])/4/Rs,Cp[P][1]-Cp[N][1],Ap)
  #declare Re=sqrt(pow(Cp[N+1][0]-Cp[N][0],2)+pow(Cp[N+1][1]-Cp[N][1],2))/4;
  Acosc((Cp[N+1][0]-Cp[N][0])/4/Re,Cp[N+1][1]-Cp[N][1],Ae)
  Acosc(Rs,1,As)
  Acosc(Re,1,Al)
  #if (rand(Rndm)>0.1) #declare Ad=rand(Rndm); #end
  #if (Ae-Ap<0) #declare Ap=Ap-2*pi; #end
  #if (Ae-Ap-As-Al-Ad>0)
    #declare I=0; #while (I<int((Ae-Ap-As-Al-Ad)*3/pi)+1)
      #declare Pe=P+1+I;
      #declare Cp[Pe][0]=Cp[N][0]+2*cos(Ap+As+Ad+pi/3*I);
      #declare Cp[Pe][1]=Cp[N][1]+2*sin(Ap+As+Ad+pi/3*I);
    #declare I=I+1; #end
    #declare Ad=0;
  #else
    #if (Ae-Ap-As-Al>0)
      #declare I=0; #while (I<int((Ae-Ap-As-Al)*3/pi)+1)
        #declare Pe=P+1+I;
        #declare Cp[Pe][0]=Cp[N][0]+2*cos(Ap+As+pi/3*I);
        #declare Cp[Pe][1]=Cp[N][1]+2*sin(Ap+As+pi/3*I);
      #declare I=I+1; #end
    #end
  #end
#end

#declare Ps=6;
#declare Jd=1; #while (Jd<4000)
  Spculc (Jd,Ps,Pe)
  #declare Ps=Pe;
#declare Jd=Jd+1; #end

#declare I=0; #while (I<5000)
  sphere { <Cp[I][0],Cp[I][1],0> 1
    pigment { color rgbt <0,0.5,0.3,0.98> }
    finish { reflection 1.5 ambient 0.2 phong 0.3 brilliance 0.5 diffuse 2 }
    }
#declare I=I+1; #end

background {color rgbt <0,0.03,0>}