Encontrar la raíz de una función con fsolve

Posted on jue 14 febrero 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:

17MtNybhdkA9GV3UNS6BTwPcuhjXoPrSzV