#!/usr/bin/perl -w

for (1..47_500_000){
    push @renta, rand() * 1_000_000_000;        # No cambia nada por el nominal
}
print "El índice Gini es ", gini(\@renta), "\n";
exit 1;



sub gini {   # Este se supone que es el código bueno para listas completas de rentas
    my $list = shift;
    my @l = @$list;
    @l = sort { $a <=> $b } @l;
    my $N = scalar @l;
    
    my (@pi, $sum_pi);
    for my $i (1..$N){
	$pi[$i] = $i/$N;
	$sum_pi += $pi[$i];
    }                                    # el denominador
       
    my $sum_RENT;
    for my $j (@l){
	$sum_RENT += $j;
    }
    my $counter = 0;
    my $accum;
    my @qi;
    for my $j (@l){
	$counter++;
	$accum += $j;
	$qi[$counter] = $accum / $sum_RENT;
    }
    
    my $resta;
    for my $i (1..$N){
	$resta += $pi[$i]-$qi[$i];
    }                                    # el numerador
			
    return ($resta/$sum_pi);
}

__END__
La definición del coeficiente o índice es curiosa, sobre
  la curva de Lorenz, pero no hace falta hacer la integral. 

