1. include "colors.inc"
2. include "functions.inc"
3. include "rand.inc"

/*

1. include "textures.inc"
2. include "glass.inc"
3. include "metals.inc"
4. include "golds.inc"
5. include "stones.inc"
6. include "woods.inc"
7. include "shapes.inc"
8. include "shapes2.inc"
9. include "math.inc"
10. include "transforms.inc"

• /

1. include "stars.inc"

global_settings { assumed_gamma 1.0 }

camera {

 location  <0,0,120>

// location  <0,200,0> 
 
 
//  location  <0,100,50> 
  
 
  
// location  <30,40,60> 
 

//location <15,15,0>

 up x sky z
 
 look_at   <0.0, 0.0, 0.0> 
 angle 40

}

/*

// stars sphere { 0,1

    texture{ Starfield2 scale 0.25
              } // 

scale 10000

}

*/  
  
  

1. declare sun1= sphere {

   <0,0,0>,2
   texture {
   
   pigment {
       color White 
       transmit 0
   }

   finish {ambient 10}
   }   

}

light_source {

   <0, 0, 0>

   color Gray75 *10
   fade_distance 5
   fade_power 2  
   
 looks_like {sun1}

}

1. declare fn_Pigm=function {

   pigment {
     agate
     color_map {
       [0 color rgb 0]
       [1 color rgb 1]
     }
   }
 }

1. declare motikka= isosurface {

//function { f_sphere(x, y, z, 0.5) - f_noise3d(x, y, z) } // function { f_sphere(x, y, z, 0.5) - f_noise3d(x * 5, y * 5, z * 5) * 0.5 }

 function { f_sphere(x, y, z, 0.5)-fn_Pigm(x/2, y/2, z/2).gray*0.5 }
 threshold 0
 accuracy 0.005

texture { pigment {color Gray }

       finish{
       diffuse 1.0
       ambient 0.0
       specular 0.0
       reflection 0.1
     }

} }

1. declare rmax=20;
2. declare rmin=10;
3. declare rdelta=rmax-rmin;
4. declare rmid=(rmax+rmin)/2;

1. declare fiidelta=10;

1. local nn=0;

1. declare Seed1=seed(1);
2. declare Seed2=seed(2);
3. declare Seed3=seed(3);
4. declare Seed4=seed(4);
5. declare Seed5=seed(5);

1. declare theta=0.0;
2. declare fii=0.0;

1. while (nn<30000)

   object 
   {
       motikka    
       

       #declare theta=radians( rand(Seed2)*360 );
       
       
       #declare seed1=seed(nn*15);
       #declare seed2=seed(nn*15143);
 
       #declare a1=Rand_Normal( 1.0, 0.33, seed1 );

       #declare rc1=a1*rdelta+rmid;   
        
       #declare a2=Rand_Normal( 1.0, 0.5, seed2 );
        
       
       #declare zc1= rand(Seed4)*rc1/10-rc1/20; 
       
       #declare zc2=zc1*a2;
                                  
                                    
       #declare x1=cos(theta)*rc1;
       #declare y1=sin(theta)*rc1;
       #declare z1=zc2;

       #declare skale1=0.05+0.05*Seed2; 
       
       
       scale skale1*0.5 
       
       rotate <x1,y1,z1>  
 
       translate <x1,y1,z1>
    }
     
 #local nn=nn+1;

1. end