#include "colors.inc"
global_settings {
  assumed_gamma 1.0 
//  max_trace_level  20

}


#declare icount = 0;
#declare R_rad = seed(0);
#declare R_tilt = seed(12345);
#declare R_twist = seed(7777);

#declare R_phase = seed(6331);
#declare omeger = function(rr) { 8/pow(rr,1.5) } 
#declare transer = function(tt){.5*(1-cos(pi*max(min(tt,1),0)))}
#declare lclock =clock+0.0 ;

#while(icount<2000)
        #declare rad=rand(R_rad)+.02;
        
        sphere { <0, 0, 0>, 0.0035 
                pigment {color Grey}
                finish {ambient .52 diffuse 0.0 }
        translate rad*x
        rotate y*(-omeger(rad)*lclock+360*rand(R_phase)) //orbit and phase
        rotate x*(90*(rand(R_phase)-1))*max(.05,exp(-.5*lclock)) //orbit tilt
        rotate y*(360*rand(R_twist)) //twist of orbit tilt
        }
    
    
        #declare icount=icount+1;    
#end
/*  */
sphere {0, 1.2
    pigment{color rgbf <1,1,1,1>}
    texture {pigment {color rgbf 1}}
    hollow
    interior {
        media {
            emission color rgbf < 0.85, 0.0, 0,.95>
            scattering {5}
            intervals 2 samples 16
            method 3
            aa_threshold 0.3 aa_level 2
            density {cylindrical scale 1 turbulence .25}
        }
    }
 scale <1,max(.05,exp(-.5*lclock)),1>  
        rotate y*(omeger(1)*lclock+360*rand(R_phase)) //orbit and phase
 
}                                             
 
/*
#declare r0=.2;
#declare dr=.2;
#while(r0<1.4)
        intersection{
                difference{
                        sphere{0,1.2     pigment{color rgbf <1,1,1,1>}}
                        cylinder{-5*y,5*y,r0-dr+.001     pigment{color rgbf <1,1,1,1>}}
                }
                cylinder{-5.1*y,5.1*y,r0-.001     pigment{color rgbf <1,1,1,1>}}
                pigment{color rgbf <1,1,1,1>}
                texture {pigment {color rgbf 1}}
                hollow
                interior {
                        media {
                            emission color rgbf < 0.85, 0.0, 0,1>
                            scattering {5}
                            intervals 2 samples 16
                            method 3
                            aa_threshold 0.3 aa_level 2
                            density {cylindrical scale exp(.25*lclock) turbulence .2+.1*lclock}
                        }
                }
        scale <1,exp(-.5*lclock),1>
        rotate y*(omeger(rad)*lclock+360*rand(R_phase)) //orbit and phase
        }
        #declare r0=r0+dr; 
#end 
*/   
// perspective (default) camera
#declare cscaler = .04;
camera {
//  location  vrotate(-40*z,0-0*transer((lclock-5)/5)*x)*cscaler*1.
    location  vrotate(-40*z,5+(1-transer((lclock-8)/5))*35*x)*cscaler*1.
  look_at   <0.0, 0.0,  0.0>
  right     x*image_width/image_height
  angle 3/cscaler
}

#declare icount = 1;
#declare rrx=0;
#declare rry=0;
#declare rrz=0;
#while(icount<10000)
        #declare rrx=rrx+abs(sin(93*icount));
        #declare rry=rry+abs(sin(57*icount));
        #declare rrz=rrz+abs(sin(27*icount));      
        sphere { <0, 0, 0>, 10 
                pigment { color rgb 1*<.5+.5*sin(icount),.5+.2*sin(2*icount),.7+.3*cos(icount)>}           
                finish {ambient 1.0 diffuse 0.0 }
        translate 4000
        rotate 3570*<rrx,rry,rrz>
        }
    
    
        #declare icount=icount+1;    
#end