function [x] = metodo_exclusao( N )
%  [x] = metodo_exclusao( N )
%
% Para gerar um vetor de N dados (N linhas, 1 coluna) 
% Obs: se omitido, N=1

if nargin==0, N=1, end

% -------
% exemplo do guia sobre metodo exclusao: f(x)=2*x com x entre 0 e 1;
xmin = 0;
xmax = 1;
ymax = 2;
f = @(x)( 2*x );  %funcao densidade de probabilidade
% -------

i = 0;
x = zeros( N, 1 );
while( i<N )
  x_cand = xmin + (xmax-xmin)*rand;
  y_test = ymax*rand;
  if y_test <= f(x_cand)
    i=i+1;
    x(i) = x_cand;
  end
end