#declare Min_factor = 0.6;

#declare Star =
isosurface {
	function { x*x*y*y + x*x*z*z + y*y*z*z - 0.25 }
//	evaluate 356*Min_factor,  sqrt(356/(356*Min_factor)),  0.7
	max_gradient 130
	open
	all_intersections
	contained_by { box{-3.01, +3.01} }
	clipped_by { box{-3, +3} }
	bounded_by { clipped_by }

	pigment { boxed color_map{
				[0.02 color rgb<1.0, 0.2, 0.5>]
				[0.23 color rgb<1.0, 0.7, 0.2>]
				[1.00 color rgb<0.1, 0.8, 0.5>]
				}
				scale 3
			}
	finish { phong 0.6 reflection 0.2 }
	scale 1/3
}


#declare Z=-5;
#while(Z<=15)
	#declare X=-10;
	#while (X<=+10)
		object { Star translate<X,0,Z> }
	#declare X=X+2;
	#end

#declare Z=Z+2;
#end

light_source { <-33,66,-99> color 1
	area_light <0,5,0>, <0,0,5>, 4, 4 jitter orient circular adaptive 1
}

plane { y,-1.125
	pigment { checker color rgb 0.2 color rgb 0.9 scale 1.5}
	}

camera {
	location <1.5, 6, -10>
	look_at <0, .25, >
	angle 40
}