function [x] = gera_dados_polinomiais( G, L, N )
%  [x] = gera_dados_polinomiais( G, L, N )
%
% Para gerar um vetor de N dados (N linhas, 1 coluna) 
%  com funcao densidade de probabilidade:
%    f(x) = A*( 1 - abs(x/L).^G )
%
% valores padroes: G=2, L=10 e N=1000

if nargin<1  %se nao for passado nenhum argumento de entrada
  G = 2;
end
if nargin<2  %se nao for passado o segundo argumento de entrada
  L = 10;
end
if nargin<3  %se nao for passado o terceiro argumento de entrada
  N = 1000;
end
A = 1; % porque o que importa eh o valor relativo de f(x)
i = 0;
x = zeros( N, 1 );
while( i<N )
  x_cand = -L + 2*L*rand;
  y_test = A*rand;
  if y_test <= A*( 1 - abs(x_cand/L)^G )
    i=i+1;
    x(i) = x_cand;
  end
end