Encontrar la raíz de una función con fsolve
Posted on Thu 14 February 2013 in Tutorial Python • 1 min read
Este artículo explica como calcular la raíz de una función utilizando la función fsolve.
El artículo se basa en un artículo en Inglés "How to find the rooots of a function with fsolve".
La función fsolve retorna la raíces de una ecuación no lineal definida por f(x) = 0.
Para este caso se calculará la raíz de la función f(x) = x^3.
A continuación se muestra el código:
#Import fsolve para calcular la raiz de la funcion x^3
from scipy.optimize import fsolve
#Importar pylab
import pylab
#importar numpy
import numpy
#se calcula la potencia 3 de x con la funcion lambda
potencia3 = lambda x : x**3
#Se calcula la raiz de x^3 iniciando con x = 10
resultado = fsolve(potencia3,10) # starting from x = 10
print resultado
#Se define 400 valores de x entre -4 a 4
x = numpy.linspace(-4,4,400)
#Se genera la grafica, pasando el valor de x
#la potencia 3era de x, el valor de resultado, la potencia 3era de resultado
pylab.plot(x,potencia3(x),resultado,potencia3(resultado),'ro')
#Se define el grid
pylab.grid(b=1)
#Se muestra la grafica
pylab.show()
La gráfica muestra el punto donde se encuentra la raíz de la función:
¡Haz tu donativo! Si te gustó el artículo puedes realizar un donativo con Bitcoin (BTC) usando la billetera digital de tu preferencia a la siguiente dirección: 17MtNybhdkA9GV3UNS6BTwPcuhjXoPrSzV
O Escaneando el código QR desde la billetera:
