Integracao por Monte Carlo — rnp_monte_carlo • rnp

Estima a integral de integrando em \[a, b\] por amostragem uniforme, com erro-padrao e intervalo de confianca.

Usage

rnp_monte_carlo(
  integrando,
  limites = c(0, 1),
  n = 10000L,
  conf = 0.95,
  seed = 42L,
  digits = 4L
)

Arguments

integrando: Funcao vetorizada a integrar.
limites: Vetor c(a, b) com os limites de integracao.
n: Inteiro. Numero de pontos de Monte Carlo.
conf: Nivel de confianca do IC.
seed: Inteiro. Semente.
digits: Inteiro. Casas decimais.

Value

tibble com estimativa, erro_padrao, ic_inf, ic_sup, n.

See also

Other probabilidade: rnp_bayes(), rnp_cadeia_markov(), rnp_distribuicao_conjunta(), rnp_esperanca_condicional(), rnp_lei_grandes_numeros(), rnp_passeio_aleatorio(), rnp_processo_poisson(), rnp_simula_aceitacao_rejeicao(), rnp_simula_inversao(), rnp_tcl_simulacao()

Examples

# integral de x^2 em [0, 1] = 1/3
rnp_monte_carlo(function(x) x^2, limites = c(0, 1), n = 1e5)
#> # A tibble: 1 × 5
#>   estimativa erro_padrao ic_inf ic_sup      n
#>        <dbl>       <dbl>  <dbl>  <dbl>  <dbl>
#> 1      0.334      0.0009  0.333  0.336 100000