void setup(){
  size(800,600);
  rectMode(CENTER_DIAMETER);
  noStroke();
}
void loop(){
  background(0);
  int limite = 40;
  translate(400,300,-8000);
  rotateY(linea(mouseX,0,width,0,TWO_PI));
  rotateZ(HALF_PI/2);

  float ancho=200;

  Movim m=new Movim();
  for(int j=0;j<limite;j++){
    for(int i=0;i<limite;i++){
      m.reiniciar();
      cilindro(i*TWO_PI/limite,j*100, linea(abs(mouseY-j*15),0,40,ancho,ancho*2) ,m);
      Movim m2=new Movim();
      for(int k=0;k<6;k++){
        int paso=330;
        translate(0,0,paso);
        m2.totZ+=paso;
        fill(linea(i,0,limite,0,255),linea(k,0,6,0,255),linea(i,0,limite,255,0),250);
        rect(0,0,70,70);        
      }
      m2.deshacer();
      m.deshacer();
    }
    //ancho*=1.08;
  }
}
class Movim{
  float totX,totY,totZ;
  float angX,angY,angZ;
  Movim(){
    reiniciar();
  }
  void reiniciar(){
    totX=totY=totZ=angX=angY=angZ=0;
  }
  void deshacer(){
    translate(-totX,-totY,-totZ);
    rotateX(-angX);
    rotateY(-angY);
    rotateZ(-angZ);
    reiniciar();
  }
}
void cilindro(float angulo, float altura, float radio, Movim m){
  rotateY(angulo);
  m.angY+=angulo;
  translate(0,altura,radio);
  m.totY+=altura;
  m.totZ+=radio;
}
float linea(float x,float x1, float x2, float y1, float y2){
  return (x-x1)/(x2-x1)*(y2-y1)+y1;
}