0

Después de descubrir que sympy sí funciona en mi pc jajaja, se da a lugar a un nuevo problema, necesito derivar un polinomio para después integrarlo (Necesito sacar longitud de arco de una función) **Lo que ya tengo: Almaceno el polinomio desde el teclado en un arreglo Se deriva por diff (Instrucción de sympy) **Mi duda: Usando el metodo de Simpson 1/3, lo tengo que integrar, sin embargo me tira el error

'int' object has no attribute 'diff'

¿Cómo puedo tratar de integrar esta derivada? Muchas gracias a todos por ayudarme Adjunto las operaciones del metodo de Simpson

def simpson13(g,a,b):
    m=(a+b) / 2
    integral= (b - a) / 6 *(g(a) + 4 * g(m) + g(b))
    return integral 


n = 100
dx=(b - a) / n
suma = 0
for i in range(n):  
    b = a + dx
    area= simpson13(g, a, b)
    suma = suma + area 
    a = b

print(suma)
def simpson13(g,a,b):
    m=(a+b) / 2
    integral= (b - a) / 6 *(g(a) + 4 * g(m) + g(b))
    return integral 


n = 100
dx=(b - a) / n
suma = 0
for i in range(n):  
    b = a + dx
    area= simpson13(g, a, b)
    suma = suma + area 
    a = b

print(suma)
5
  • Hola Darío, antes de intentar reinventar la rueda XD ¿que versión de Python tienes? ¿Por que dices que no funciona con tu verión de Python? ¿Algún error al instalar? ¿Sistema operativo?
    – FJSevilla
    Commented el 4 abr. 2020 a las 20:16
  • @FJSevilla Jajaja, tengo la versión 3.7 de python, perdón por no ponerla, intenté instalarlo y todo bien, de hecho en el prompt de anaconda, me aparece como instalado, pero al quererlo usar me aparece que no existe, después de volver a intentar de todo, encontré en un blog de anaconda, que no funciona en las versiones nuevas de este mismo. Commented el 4 abr. 2020 a las 20:23
  • Dos cositas: la pregunta ahora es totalmente diferente que el original. Eso no es la función de editar preguntas, por que ahora la respuesta original ya no encaja para nada. Tendria que ser un post nuevo. Otra cosa: hace falta importar librería, dar valores iniciales a todo (a, b, g) para obtener un ejemplo completo y quitar el copy-paste doble. Si hay un error, al menos hay que explicar en que linea y con que variables pasa.
    – JohanC
    Commented el 5 abr. 2020 a las 1:15
  • No olvides votar para respuestas útiles. Y aceptar una respuesta si esta respuesta resuelve la pregunta.
    – JohanC
    Commented el 5 abr. 2020 a las 1:41
  • @DaríoValerio, si te ha sido de utilidad la respuesta proprocionada, puedes marcarla como válida. Saludos.
    – Adrian
    Commented el 5 sept. 2020 a las 11:31

1 respuesta 1

0

Primero de todo, la librería simbólica de Python se llama 'sympy'. La librería 'simpy' tiene que ver con simulaciones.

Una diferencia entre Python y muchos programas de matemáticas simbólicas, es que en Python el simbolo ^ solo funciona como "Disyunción exclusiva" (XOR). Para potencias hay que emplear **.

La función Derivative sirve para derivar. Sin nada más, solo represente el derivado sin evaluar. Para evaluarlo, hay que añadir la función .doit().

Alternativamente, la función diff calcula el derivado sin necesidad de .doit().

Un ejemplo:

from sympy import Symbol, Derivative, diff, sin, sympify

x = Symbol('x')
print(Derivative(x**2).doit()) # 2*x
print(diff(sin(x**2)) # 2*x*cos(x**2)
print(diff(1/sin(x**2))) # -2*x*cos(x**2)/sin(x**2)**2

No se puede emplear diff con una expresión que no sea de tipo sympy. Pero si con el resultado de otro diff. Y tambien con constantes convertido en tipo sympy:

print(diff(7*x)) # 7
print(diff(diff(7*x))) # 0
print(diff(sympify(23))) # 0

Para derivar un polinomio:

def derivar_polinomio(coef, x):
    p = sum([c * x**i  for i, c in enumerate(coef[::-1])])
    return diff(p)

print(derivar_polinomio([7, 3, 5], x)) # 14*x + 3

Para la pregunta nueva, supongo que la idea seria que g es derivado de una función x**2 y a y b son 2 constantes. Para evaluar g con x substituido por a, en sympy se calcula como g.subs(x, a). Sympy lo mire como expresiones más que como funciones.

from sympy import diff, Symbol

def simpson13(g, a, b):
    m = (a + b) / 2
    integral = (b - a) / 6 * (g.subs(x, a) + 4 * g.subs(x, m) + g.subs(x, b))
    return integral

x = Symbol('x')
g = diff(x ** 2)
a = 2
b = 3

n = 100
dx = (b - a) / n
suma = 0
for i in range(n):
    b = a + dx
    area = simpson13(g, a, b)
    suma += area
    a = b

print(suma)
7
  • muchas gracias por contestar, una pregunta más, se puede derivar con derivate un polinomio almacenado con la instrucción 'np.poly1d'? Commented el 4 abr. 2020 a las 22:29
  • El asunto es que el usuario tiene que insertar un polinomio, este se almacena en un arrelo y de ahí no sé como puedo derivarlo e integrarlo, ya intenté pasarlo a función, y ahora estoy intentando con 'np.poly1d' Commented el 4 abr. 2020 a las 22:33
  • Gracias nuevamente por contestar, @JohanC, disculpa, habrá alguna manera de usar esa derivada de sympy, con operaciones aritmeticas? Es decir, que me dé el resultado de una cierta derivada y esta pasarla a función para evaluarla y poder integrar? Especificamente hablando de la regla o metodo de simpson Gracias Commented el 4 abr. 2020 a las 22:38
  • Ya que al querer integrar con operaciones aritmeticas la derivada que me tira diff, me lanza este error: 'int' object has no attribute 'diff' Commented el 4 abr. 2020 a las 22:40
  • gracias nuevamente por seguir contestando Commented el 4 abr. 2020 a las 22:42

Tu Respuesta

By clicking “Publica tu respuesta”, you agree to our terms of service and acknowledge you have read our privacy policy.

¿No es la respuesta que buscas? Examina otras preguntas con la etiqueta o formula tu propia pregunta.