%
function [hbin, h] = mkhisto(rgen, xmin, xmax, ncan, nensayos)
%
% Realiza histograma de un generador aleatorio pasado como argumento.
%
  sf = ncan / (xmax - xmin);
  
  %hbin(ncan + 1) = xmax;
  for i = ncan:-1:1
    hbin(i) = xmin + (i - 1) / sf;
    h(i)    = 0;
  end
%
  for n = 1:nensayos
    x = rgen();
    if ((x >= xmin) && (x < xmax))
      z   = sf * (x - xmin);
      ihi = floor(z) + 1;
      h(ihi) = h(ihi) + 1;
    end
  end
 %
 % Normalización.
 %
 h = h .* (sf / nensayos);
end

