gera_polino<-function(x,n)
  # x : vetor de dados da amostra
  # n : grau do ultimo polinomio
  #
  #Retorna um array com dimensões (N, n+1, n) onde N é o comprimento
  #do vetor x de entrada.
  #
  #Cada polinomio esta na posição P[,,i], onde i é o grau do polinomio
  {

  N<-1*length(x)
  Pols<-array(,c(N,n+1,n))
  M<-matrix(0,nrow=N,ncol=n+1)
  for(i in 1:n)
  {
    j<-1
    while(j<=i+1)
    {
      aux<-x^(j-1)
      M[,(n+2)-j]<-aux
      j<-j+1
    }
    Pols[,,i]<-M
    }
  Pols
}