Universidade Federal do Paraná
Curso de Estatística
CE 089 -
Estatística Computacional II - 2014/2
Prof. Dr. Walmes Marques Zeviani
Tabela de conteúdo
##-----------------------------------------------------------------------------
## Funções do slide:
## * f(x) = x^2+ln(x).
## * f(x) = x^2-8*x+15.
## Implementar o método de Newton-Raphson.
##-----------------------------------------------------------------------------
## Função 1.
f0 <- function(x, A, B, C){
x*(A+B*x)^(-C)
}
curve(f0(x, 1, 1, 10), 0, 2)
D(expression(x*(A+B*x)^(-C)), "x")
## (A + B * x)^(-C) + x * ((A + B * x)^((-C) - 1) * ((-C) * B))
f1 <- function(x, A, B, C){
(A + B * x)^(-C) + x * ((A + B * x)^((-C) - 1) * ((-C) * B))
}
curve(f1(x, 1, 1, 10), 0, 2)
abline(h=0, lty=2)
## Encontrar as raíz de f1 (que corresponde ao máximo de f0).
##-----------------------------------------------------------------------------
## Função 2.
## Gamma incompleta.
D(expression(A*x^B*exp(-C*x)), "x")
## A * (x^(B - 1) * B) * exp(-C * x) - A * x^B * (exp(-C * x) *
## C)