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)