File:Star that have asteroid zone 2.png

      //

// // star and asteroid zone // Povray 3.7 code // // 11.7.2023 0000.0001

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  <60,80,120> 
 

//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 3}
   }   

}

light_source {

   <0, 0, 0>

   color Gray75 *50
   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 numasteroids=1000;

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 theta=0.0;
2. declare fii=0.0;

1. while (nn<numasteroids)

   object 
   {
       motikka    
       

       #declare seed1=seed(nn*1);
       #declare seed2=seed(nn*15143);
       #declare seed3=seed(nn*76359);
       #declare seed4=seed(nn*8987547);
       #declare seed5=seed(nn*1543); 
       
       #declare seed6=seed(nn*76923439);
       #declare seed7=seed(nn*87639);
       #declare seed8=seed(nn*76312); 
       
       #declare theta=radians( rand(seed1)*360 );
 
       #declare a1=Rand_Normal( 1.0, 0.33, seed2 );

       #declare rc1=a1*rdelta+rmid;   
        
       #declare a2=Rand_Normal( 1.0, 0.5, seed3 );
        
       
       #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 de1=rand(seed5);
       
       // fictional size distribution
       #declare skale1=pow(de1,-0.111)*0.2;  
         
       #declare rot1=rand(seed6)*360;
       #declare rot2=rand(seed7)*360;     
       #declare rot3=rand(seed8)*360;
         
     
       scale skale1 
       
       rotate <rot1,rot2, rot3>  
 
       translate <x1,y1,z1>
    }
     
 #local nn=nn+1;

1. end