% Input: z0, punto iniziale del ramo
%        r0, direzione del ramo (punto finale = z0+r0)
%        c, vettore di lunghezza N_rami dei fattori moltiplicativi
%           rho e^( j alpha ) per ogni sottoramo
%        livello, intero del livello del ramo considerato
%        livello_max, massimo numero di livelli
function ramo( z0 , r0 , c , livello , livello_max )
  z1 = z0 + r0 ;
  plot( [z0 z1] ) ;
  if livello < livello_max
    for sottoramo = 1 : length(c)
      r1 = r0 * c(sottoramo) ;
      ramo( z1 , r1 , c , livello+1 , livello_max ) ;
    end
  end
end