import matplotlib.pyplot as plt
import numpy as np
from py_expression_eval import *
def f(exp, var, x0):
p = Parser()
result = p.parse(exp).evaluate({var:x0})
return result
a = 0
b = 4
error = 10
i = 0
while(error>1e-8 and i!=100):
c = (a + b) / 2
fa = f('x^3 -2 * x^2 - 1', 'x', a)
fc = f('x^3 -2 * x^2 - 1', 'x', c)
if(fc == 0):
raiz = c
break
elif(fa * fc < 0):
b = c
else:
a = c
raiz = c
i += 1
error = abs(fc)
print("Iteracion",i,". Raiz aproximada:",raiz)
print(raiz)
print(i)
print(f('x^3 -2 * x^2 - 1', 'x',raiz))
x = np.linspace(0, 4, 101)
print(x)
plt.plot(x, f('x^3 -2 * x^2 - 1','x',x))
plt.plot(a, f('x^3 -2 * x^2 - 1', 'x', a), 'or')
plt.grid()
plt.show()
Me dice el error del título al tratar de graficar después de implementar el evaluador de expresiones, antes no me daba problema con el siguiente código:
import matplotlib.pyplot as plt
import numpy as np
import py_expression_eval
def f(x):
return x**3 -2 * x ** 2 - 1
a = 0
b = 4
error = 10
i = 0
while(error>1e-8 and i!=100):
c = (a + b) / 2
fa = f(a)
fc = f(c)
if(fc == 0):
raiz = c
break
elif(fa * fc < 0):
b = c
else:
a = c
raiz = c
i += 1
error = abs(fc)
print("Intervalo ("+str(a) + "," + str(b) +")" )
print("Iteracion",i,". Raiz aproximada:",raiz)
print(raiz)
print(i)
print(f(2.2055694311857224))
x = np.linspace(0, 4, 101)
plt.plot(x, f(x))
plt.grid()
plt.show()
Si alguien puede ayudarme a que funcione con el expression eval
se lo agradezco de antemano.
El mensage de error completo es este:
Traceback (most recent call last):
File "C:\Users\Daniel\Desktop\bisection3.py", line 34, in
plt.plot(x, f('x^3 -2 * x^2 - 1','x',x))
File "C:\Users\Daniel\Desktop\bisection3.py", line 7, in f
result = p.parse(exp).evaluate({var:x0})
File "C:\Users\Daniel\AppData\Local\Programs\Python\Python36-32\lib\site-packages\py_expression_eval__init__.py", line 122, in evaluate nstack.append(f(n1, n2))
TypeError: only length-1 arrays can be converted to Python scalars
Si trato de convertir x
a int
antes de pasarlo a la función que va graficar, me da el mismo error del titulo, es algo que tiene que ver con el tipo de dato pero nose que es :(