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 < jmcontre...@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@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es > -- Saludos, Carlos Ortega www.qualityexcellence.es [[alternative HTML version deleted]] _______________________________________________ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es