// Exemplo de arquivo de descricao de cena para POV-ray
// Last edited on 2010-03-04 15:44:01 by stolfi

// ======================================================================
// CORES E TEXTURAS

background{ color rgb < 0.75, 0.80, 0.85 > }

#declare tx_plastico = 
  texture{
    pigment{ color rgb < 0.10, 0.80, 1.00 > }
    finish{ diffuse 0.8 ambient 0.1 specular 0.5 roughness 0.005 }
  }
#declare tx_espelho = 
  texture{
    pigment{ color rgb < 1.00, 0.85, 0.30 > }
    finish{ diffuse 0.2 reflection 0.7*< 1.00, 0.85, 0.30 > ambient 0.1 }
  }

// ======================================================================
// DESCRIÇÃO DA CENA 

#declare raio = 1.000;
#include "eixos.inc"

// Aqui está a cena, finalmente:


#macro interpola(f0,v0,f1,v1,f)
  #local ss=(f-f0)/(f1-f0);
  #local rr=1-ss;
  
  (rr*v0 +ss*v1)
#end

#macro interpola4(t0,v0,v1,v2,v3,t3,tt)
#local v01 = interpola(t0,v0,t3,v1,tt);
#local v12 = interpola(t0,v1,t3,v2,tt);
#local v23 = interpola(t0,v2,t3,v3,tt);

#local v012 = interpola(t0,v01,t3,v12,tt);
#local v123 = interpola(t0,v12,t3,v23,tt);

interpola(t0,v012,t3,v123,tt);
#end

#macro grafico(x0,y0,y1,y2,y3,x3,N)
 #local delta = (x3-x0)/N;
 #local i = 0;
 #local xx = x0;
 #while(i <= N)
   #local yy = interpola4(x0,y0,y1,y2,y3,x3,xx);
   object{ sphere { yy, 0.10 } texture { tx_plastico } translate <-3,0,2.5>}
   #local xx = xx + delta;
   #local i = i + 1;
 #end
#end

#macro grafico4(x0,y1,y2,h1,h2,u1,u2,p1,p2,x1,n)
   #local h0 = (h1+y2)/2;
   #local y3=h0;
   #local u0 = (u1+h2)/2;
   #local h3=u0;
   #local p0 = (p1+u2)/2;
   #local u3=p0;

   #local y0 = (y1+p2)/2;
   #local p3 = y0;

   #local k = (x1-x0)/4;

	grafico(x0,y0,y1,y2,y3,x0+k,n)
	grafico(x0+k,h0,h1,h2,h3,x0+k*2,n)	
	grafico(x0+k*2,u0,u1,u2,u3,x0+k*3,n)
	grafico(x0+k*3,p0,p1,p2,p3,x0+k*4,n)

#end
#include "retalho.inc"

#macro ret(p1,p2,p3,p4,p5,p6,p7,p8,p9,p10,p11,p12,p13,p14,p15,p16,p17,p18,p19,p20,p21,p22,p23,p24,r,tx1,tx2,tx3,tx4)

	#local pp1 = (p9+p13)/2;
	#local pp2 = (p10+p14)/2;
	#local pp3 = (p11+p15)/2;
	#local pp4 = (p12+p16)/2;

   retalho(p1,p2,p3,p4,p5,p6,p7,p8,p9,p10,p11,p12,pp1,pp2,pp3,pp4,r,tx1,tx2)
   retalho(pp1,pp2,pp3,pp4,p13,p14,p15,p16,p17,p18,p19,p20,p21,p22,p23,p24,r,tx3,tx4)

#end
union {
object{ eixos(3.00) }
grafico4(0,<3,-2,5>,<5,3,1>,<1,-5,-4>,<-4,1,4>,<4,3,2>,<-4,-1,1>,<4,-3,2>,<-5,-4,1>,5,400)
union{ret(<0,0,0>,<0,1,0>,<0,2,0>,<0,3,0>,<1,0,3>,<1,1,3>,<1,2,3>,<1,3,3>,<2,0,-1>,<2,1,-1>,<2,2,-1>,<2,3,-1>,<4,0,2>,<4,1,2>,<4,2,2>,<4,3,2>,<5,0,2>,<5,1,2>,<5,2,2>,<5,3,2>,<6,0,-3>,<6,1,-3>,<6,2,-3>,<6,3,-3>,0.0000001,tx_espelho,tx_plastico,tx_plastico,tx_espelho) translate <2,0,0>}

}



#include "camlight.inc"
#declare centro_cena = < 1.00, 1.00, 1.00 >;
#declare raio_cena = 10.0;
#declare dir_camera = < 1.00, -2.00, -0.50 >;
#declare dist_camera = 80.0;
#declare intens_luz = 1.00;
camlight(centro_cena, raio_cena, dir_camera, dist_camera , z, intens_luz)