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

background{ color rgb < 0.5, 0.5, 0.5 > }



#declare branco = 
  texture{
    pigment{ color rgb < 1 1 1 > }
    finish{ diffuse 0.9 ambient 0.1 }
  }
#declare preto = 
  texture{
    pigment{ color rgb < 0 0 0 > }
    finish{ diffuse 0.9 ambient 0.1 }
  }


// ======================================================================
// DESCRI��O DA CENA 




#declare cubo =  
  box{ 
    <-1,-1,-1>,<1,1,1> texture{branco}
  }
#declare esfera =
  sphere{ 
    <0,0,0>,1 texture{branco}
  }
#declare cone_invertido =
  cone{
    <0,0,0.0>,0,<0,0,3>,1 texture{branco}
  }
#declare cone_normal=
  cone{
    <0,0,1>,1,<0,0,3>,0 texture{branco}
  }
#declare base=
  cylinder{
    <0,0,0>,<0,0,0.5>,1.5 texture{preto}
  }
#declare esfera=
  sphere{
    <0,0,0>,1 texture{branco}
  }
#declare mini_esfera=
  sphere{
    <0,0,0>,0.7 texture{branco}
  }
#declare esferinha=
  sphere{
    <0,0,0>,0.1 texture{preto}
  }
#declare cano=
  cylinder{
    <0,0,0>,<0.5,0,0>,0.1 texture{preto}
  }
#macro tanque1()
  union{
    object{esfera} object{cano translate<1,0,0> }
  }
#end
#macro tanque2()
  union{
    object{esfera} object{cone_invertido} object{cano translate<1,0,0> } object{cano translate<-1.5,0,0>} 
  }
#end
#macro tanque4()
  union{
    object{cubo} object{cone_normal} object{cano translate<1,0,0>} object{cano translate<1,0,0> rotate 90*z} object{cano translate<1,0,0> rotate 180*z} object{cano translate<1,0,0> rotate 270*z}

  }
#end
#macro tanque3()
  union{
    object{esfera} object{mini_esfera translate<0,0,1.5>} object{cano translate<1,0,0> } object{cano translate<-1.5,0,0>} object{cano translate<0.7,0,1.5> rotate 90*z}
  }
#end
#macro interpola1(T,t0,v0,t1,v1)
  #local ss=(T-t0)/(t1-t0);
  #local vv=(1-ss)*v0+ss*v1;
  vv
#end
#macro interpola3(T,ta,tb,v0,v1,v2,v3)
  #local v01=interpola1(T,ta,v0,tb,v1);
  #local v12=interpola1(T,ta,v1,tb,v2);
  #local v23=interpola1(T,ta,v2,tb,v3);
  #local v012=interpola1(T,ta,v01,tb,v12);
  #local v123=interpola1(T,ta,v12,tb,v23);
  #local v0123=interpola1(T,ta,v012,tb,v123);
  v0123
#end
#macro testa_interpola1(obj,t0,v0,t1,v1,n)
  object{obj translate<t0,v0,0>}
  #local cont=1;
  #while(cont<n-1)
    #local position=t0+(t1-t0)*cont/(n-1);
    object{obj translate<position,interpola1(position,t0,v0,t1,v1),0>} 
    #local cont=cont+1;
  #end
  object{obj translate<t1,v1,0>}
#end

#macro testa_interpola3_multi(obj,N,Pini,P1,P2,Pfim,n)
  #local t0=0;
  #local t1=N;
  object{obj translate Pini}
  #local cont=1;
  #while(cont<n-1)
    #local position=t0+(t1-t0)*cont/(n-1);
    object{obj translate interpola3_multi(position,Pini,N,P1,P2,Pfim)} 
    #local cont=cont+1;
  #end
  object{obj translate Pfim}
#end

#macro interpola3_multi(T,Pini,n,P1,P2,Pfin)
  #local k=int(T);
  #if(k=0)
    #local Q0=Pini;
  #else
    #local Q0=(P2[k-1]+P1[k])/2;
  #end
  #local Q1=P1[k];
  #local Q2=P2[k];
  #if(k=n-1)
    #local Q3=Pfin;
  #else
    #local Q3=(P2[k]+P1[k+1])/2;
  #end
  #local Pf=interpola3(T,k,k+1,Q0,Q1,Q2,Q3);
  Pf
#end

#declare roleta=seed(34567);
#declare ponta11=<1.5,0,0>;
#declare ponta12=<-1.5,0,0>;
#declare ponta13=<0,0,3.5>;
#declare ponta21=<1.5,0,0>;
#declare ponta22=<0,1.5,0>;
#declare ponta23=<0,0,3.5>;
#declare ponta31=<1.5,0,0>;
#declare ponta32=<-1.5,0,0>;
#declare ponta33=<0,0,2.7>;



#include "eixos.inc"

// Aqui est� a cena, finalmente:


// testa_interpola3(esferinha,0+ponta11.x,20+ponta11.x,ponta11,ponta11+<5,0,0>,ponta12+<20,20,0>-<5,0,0>,ponta12+<20,20,0>,10000)
// object{tanque2()}
// object{tanque2() translate<20,20,0>}

#local Pini=<0,0,0>;
#local n=3;
#local P1=array[n];
#local P2=array[n];
#local P1[0]=<0,10,0>;
#local P2[0]=<10,10,0>;
#local P1[1]=<10,20,0>;
#local P2[1]=<20,20,0>;
#local P1[2]=<20,30,0>;
#local P2[2]=<30,30,0>;
#local Pfim=<30,40,0>;
#local Pini2=<-10,10,0>;
#local n2=1;
#local P12=array[n2];
#local P22=array[n2];
#local P12[0]=<-10,20,0>;
#local P22[0]=<0,20,0>;
#local Pfim2=<0,30,0>;
testa_interpola3_multi(esferinha,n,Pini,P1,P2,Pfim,100000)
testa_interpola3_multi(esferinha,n2,Pini2,P12,P22,Pfim2,100000) 
#include "camlight.inc"
#declare centro_cena = < 20 ,20, -3>;
#declare raio_cena = 30;
#declare dir_camera = < 2, 2, 2>;
#declare dist_camera = 5*raio_cena;
#declare intens_luz = 1.20;
camlight(centro_cena, raio_cena, dir_camera, dist_camera , z, intens_luz)