[R-br] matematica com o R

[Samuel Carvalho]

Caros, boa tarde!
Alguém poderia por gentileza auxiliar com a duvida relativa ao codigo abaixo?

# RCMR
y <- expression(exp(a+b/x))
D(y,'x') #primeira derivada
dy.dx <- expression(-(exp(a + b/x) * (b/x^2))) #resultado da primeira derivada

D(dy.dx, 'x') #segunda derivada
fx <- function(a,b,x) exp(a + b/x) * (b * (2 * x)/(x^2)^2) + exp(a + b/x) * (b/x^2) * (b/x^2) #resultado da segunda derivada
optim(c(1), fx, a=3.3,b=-1.2,method="BFGS")


# A ideia é encontrar o ponto de inflexao para o modelo representado pelo objeto y. Por definição o ponto de inflexao é o valor de x que faz a segunda derivada da 
# função igual a zero. Montei este rcmr porém não estou certo se a otimização de fx está igualando a zero, até porque os resultados não estão condizentes.
# Gostaria se possível também de alguma dica de algum livro que tenha dicas de cálculos matemáticos com o R pois gostaria de adaptar algumas rotinas do maple 
# para o R, até porque recentemente precisei isolar um termo de um modelo não linear que manualmente não é tão simples e gostaria de fazer estes calculos com o 
# R e nao mais com o maple. Segue um exemplo simples de um calculo que gostaria de reproduzir com o R

y = b0 + b1x
b0=?
b0 = y - b1x # esta é a saída que gostaria de encontrar com o R. Ate cheguei a ver um pouco sobre o pacote Ryacas mas sem muito sucesso.

Ja deixo aqui meus agradecimentos


====================================
Samuel P. C. Carvalho
Engº Florestal [UFLA]
Mestre em Ciências Florestais [UFLA]
Doutor em Recursos Florestais [ESALQ/USP]
=============================================
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130912/2ef927cc/attachment.html>