[R-es] resolución de ecuación

Carlos Ortega cof en qualityexcellence.es
Vie Oct 23 13:52:17 CEST 2015


Hola,

Otra alternativia que tienes, por la complejidad de la ecuación que la
tienes expresada como un sumatorio (sin desarrollar) es utilizar "R" unido
al paquete de matemática simbólica "Yacas", mediante "RYacas":

https://cran.r-project.org/web/packages/Ryacas/index.html

Saludos,
Carlos Ortega
www.qualityexcellence.es

El 23 de octubre de 2015, 10:58, José Miguel Contreras García <
jmcontreras en ugr.es> escribió:

> Hola a todos
>
> Tengo una duda existencial, como siempre, jejeje
>
> Quiero resolver una ecuación (F=0) que depende de un parámetro z, pero
> este viene en función de una suma que depende de un vector de tiempos ti
> (1:25) y unos valores de un vector x
>
> Lo he intentado resolver mediante una función para luego utilizar solve,
> pero me he atrancado.
>
> La función es esta:
>
> ti<-1:25
> x<-sample(1:10, size=25, replace=T)
> beta<-0.1
>
> F<-function(z) {
> suma<-0
> for(i in 2:n){
>
> sum<-sum+exp(1)^(-beta*ti[i])/(z-beta)^2*(log(x[i])-exp(1)^(-beta)*log(x[i-1]))*(exp(1)^(ti[i]*(beta-z))*(1-ti[i]*(beta-z))-exp(1)^(-ti[i-1]*(beta-z))*(1-ti[i-1]*(beta-z)))
> }
> }
>
> Quiero calcular el valor de z para que la suma valga cero. ¿He metido
> mucho la pata?
>
> Gracias de antemano
>
> Saludos
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>



-- 
Saludos,
Carlos Ortega
www.qualityexcellence.es

	[[alternative HTML version deleted]]



Más información sobre la lista de distribución R-help-es