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

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

#include "textures.inc"

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

#declare tx_chao =
  texture{
    pigment{ color rgb < 0.00, 0.70, 0.00 > }
    scale 2.0
  }

#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_fosca = 
  texture{
    pigment{ color rgb < 1.00, 0.80, 0.10 > }
    finish{ diffuse 0.9 ambient 0.1 }
  }

#declare tx_prateado = 
  texture{
    pigment{ color rgb < 0.30, 0.30, 0.30 > }
    finish{ diffuse 0.2 reflection 0.6 brilliance 0.1 }
  }

#declare tx_vidro = 
  texture{
    pigment{ color rgb < 0.85, 0.95, 1.00 > filter 0.70 }
    finish{ diffuse 0.03 reflection 0.25 ambient 0.02 specular 0.25 roughness 0.005 }
  }

#declare tx_xadrez =
  texture{
    pigment{ checker color rgb < 0.10, 0.32, 0.60 >, color rgb < 1.00, 0.97, 0.90 > }
    finish{ diffuse 0.9 ambient 0.1 }
    scale 2.0
  }

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

#declare raio = 2.000;

// Partes da cena:

#declare chao = 
  box{ <-20,-20,-1>, <+20,+20,0> }

#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 { <0.00, xx, yy>, 0.10 } texture { tx_plastico }}
    #local xx = xx + delta;
    #local i = i + 1;
  #end
#end

#macro grafico2(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 }}
    #local xx = xx + delta;
    #local i = i + 1;
  #end
#end

#declare graphs = array[4];
#declare graphs[0] = array[2];
#declare graphs[1] = array[2];
#declare graphs[2] = array[2];
#declare graphs[3] = array[2];

#declare graphs[0][0] = 2;
#declare graphs[0][1] = 1;
#declare graphs[1][0] = 3;
#declare graphs[1][1] = 2;
#declare graphs[2][0] = -2;
#declare graphs[2][1] = 1;
#declare graphs[3][0] = 0;
#declare graphs[3][1] = 2;

#declare graphs1 = array[4];
#declare graphs1[0] = array[2];
#declare graphs1[1] = array[2];
#declare graphs1[2] = array[2];
#declare graphs1[3] = array[2];

#declare graphs1[0][0] = <1,2,2>;
#declare graphs1[0][1] = <-3,2,2>;
#declare graphs1[1][0] = <2,-2,2>;
#declare graphs1[1][1] = <1,2,-2>;
#declare graphs1[2][0] = <3,2,2>;
#declare graphs1[2][1] = <1,2,5>;
#declare graphs1[3][0] = <1,7,2>;
#declare graphs1[3][1] = <2,-8,2>;

#macro graficoTotal(x0,y0,y3,x3,N,Npontos)
  #local delta = (x3-x0)/N;
  #local xx = x0;
  
  #local yCalculado = (graphs[0][1] + graphs[1][0])/2;
  grafico(x0,y0,graphs[0][0],graphs[0][1],yCalculado,xx+delta,Npontos)
  #local xx = xx + delta;
  #local yCalculado1 = (graphs[1][1] + graphs[2][0])/2;
  grafico(xx,yCalculado,graphs[1][0],graphs[1][1],yCalculado1,xx+delta,Npontos)
  #local xx = xx + delta;
  #local yCalculado = (graphs[2][1] + graphs[3][0])/2;
  grafico(xx,yCalculado1,graphs[2][0],graphs[2][1],yCalculado,xx+delta,Npontos)
  #local xx = xx + delta;
  grafico(xx,yCalculado,graphs[3][0],graphs[3][1],y3,x3,Npontos)
#end

#macro graficoTotal2(x0,y0,y3,x3,N,Npontos)
  #local delta = (x3-x0)/N;
  #local xx = x0;
  
  #local pCalculado = (graphs1[0][1] + graphs1[1][0])/2;
  grafico2(x0,y0,graphs[0][0],graphs1[0][1],pCalculado,xx+delta,Npontos)
  #local xx = xx + delta;
  #local pCalculado1 = (graphs1[1][1] + graphs1[2][0])/2;
  grafico2(xx,pCalculado,graphs1[1][0],graphs1[1][1],pCalculado1,xx+delta,Npontos)
  #local xx = xx + delta;
  #local pCalculado = (graphs1[2][1] + graphs1[3][0])/2;
  grafico2(xx,pCalculado1,graphs1[2][0],graphs1[2][1],pCalculado,xx+delta,Npontos)
  #local xx = xx + delta;
  grafico2(xx,pCalculado,graphs1[3][0],graphs1[3][1],y3,x3,Npontos)
#end

#include "eixos.inc"

// Aqui está a cena, finalmente:

union{
  //object{ chao  translate < 0,0,-5 > texture{ tx_xadrez } }
  object{ eixos(3.00) }
  graficoTotal(0,2,2,9,16,100)
  //graficoTotal2(0,2,2,9,16,100)
}

#include "camlight.inc"
#declare centro_cena = < 0.00, 0.00, 1.00 >;
#declare raio_cena = 9;
#declare dir_camera = < 10.00, 0.00, 5.00>;
#declare dist_camera = 16;
#declare intens_luz = 1.00;
camlight(centro_cena, raio_cena, dir_camera, dist_camera , z, intens_luz)