[R-es] potencia fracional de un número negativo

[José Trujillo Carmona]

Olivier te dio las soluciones. Dos soluciones.

Una utilizar la función para solucionar polinomios de orden n:

a_1 x^n + a_2 x^(n-1) + ... + a_(n-1) x^2 + a_n x + a_(n+1) = 0

Se resuelve mediante:

polyroot(a_(n+1), a_n, a(n-1), ..., a_2, a_1)

La raiz enésima de un número b es un caso particular:

b^1/n = x --------> x^n - b = 0 ----------> x^n + 0 x^(n-1) + ... + 0 
x^2 + 0 x + b = 0

Que se resuelve mediante:

polyroot(b, 0, ..., 0, 0, 1)

Y la otra te fabricó una función para que aún sea más fácil:

sq<-function(a,q){
	x=as.complex(a)
	theta=(2*pi*(1:q)+Arg(x))/q
	complex(mod=Mod(x)^(1/q),arg=theta)
	}
	
>/  sq(-8,3)
/[1] -2+0.000000i  1-1.732051i  1+1.732051i
>/  sq(-.5,5)
/[1] -0.2690149+0.8279428i -0.8705506+0.0000000i -0.2690149-0.8279428i
[4]  0.7042902-0.5116968i  0.7042902+0.5116968i



Saludos

El 16/10/15 a las 01:44, Alex J. Zambrano escribió:
> Mirando los comentarios, realmente lo que deseo es encontrar la raíz real
> de (-0.5)^(1/5) la cual debería ser -0.87055056329. José me hace caer en
> cuenta que además de no encontrar la raiz real, tampoco da todas las raiz
> complejas. Habría alguna manera de que tuviera en cuenta?
>
>


	[[alternative HTML version deleted]]