Dérivée trop longue sur sympy.

Bonsoir.

Je cherche à plot une dérivée assez longue don voici le code:

Code:

1
2
3
4
5
6
7
8
9
from sympy import *
from sympy.abc import x
 
sa4 = (sqrt(2)*(sqrt(3)*(-sqrt(3)*cot(2*x + pi/3) + 1)*(-2*cot(2*x + pi/3)**2 - 2)/2 - sqrt(3)*(-2*cot(2*x + pi/3)**2 - 2)/2)*(1.5*sqrt(3)*(-sqrt(3)*cot(2*x + pi/3) + 1)*(-2*cot(2*x + pi/3)**2 - 2) - 1.5*sqrt(3)*(-2*cot(2*x + pi/3)**2 - 2))*((-sqrt(3)*cot(2*x + pi/3) + 1)**2/2 + sqrt(3)*cot(2*x + pi/3) + 1)**(-2.5) + sqrt(2)*(-sqrt(3)*(-sqrt(3)*cot(2*x + pi/3) + 1)*(-4*cot(2*x + pi/3)**2 - 4)*cot(2*x + pi/3) + sqrt(3)*(-4*cot(2*x + pi/3)**2 - 4)*cot(2*x + pi/3) + (-2*cot(2*x + pi/3)**2 - 2)*(3*cot(2*x + pi/3)**2 + 3))*((-sqrt(3)*cot(2*x + pi/3) + 1)**2/2 + sqrt(3)*cot(2*x + pi/3) + 1)**(-1.5))**2/(-sqrt(2)*(sqrt(3)*(-sqrt(3)*cot(2*x + pi/3) + 1)*(-2*cot(2*x + pi/3)**2 - 2)/2 - sqrt(3)*(-2*cot(2*x + pi/3)**2 - 2)/2)*(1.5*sqrt(3)*(-sqrt(3)*cot(2*x + pi/3) + 1)*(-2*cot(2*x + pi/3)**2 - 2) - 1.5*sqrt(3)*(-2*cot(2*x + pi/3)**2 - 2))*(2.5*sqrt(3)*(-sqrt(3)*cot(2*x + pi/3) + 1)*(-2*cot(2*x + pi/3)**2 - 2) - 2.5*sqrt(3)*(-2*cot(2*x + pi/3)**2 - 2))*((-sqrt(3)*cot(2*x + pi/3) + 1)**2/2 + sqrt(3)*cot(2*x + pi/3) + 1)**(-3.5) - sqrt(2)*(sqrt(3)*(-sqrt(3)*cot(2*x + pi/3) + 1)*(-2*cot(2*x + pi/3)**2 - 2)/2 - sqrt(3)*(-2*cot(2*x + pi/3)**2 - 2)/2)*(-3.0*sqrt(3)*(-sqrt(3)*cot(2*x + pi/3) + 1)*(-4*cot(2*x + pi/3)**2 - 4)*cot(2*x + pi/3) + 3.0*sqrt(3)*(-4*cot(2*x + pi/3)**2 - 4)*cot(2*x + pi/3) + (-2*cot(2*x + pi/3)**2 - 2)*(9.0*cot(2*x + pi/3)**2 + 9.0))*((-sqrt(3)*cot(2*x + pi/3) + 1)**2/2 + sqrt(3)*cot(2*x + pi/3) + 1)**(-2.5) - 2*sqrt(2)*(1.5*sqrt(3)*(-sqrt(3)*cot(2*x + pi/3) + 1)*(-2*cot(2*x + pi/3)**2 - 2) - 1.5*sqrt(3)*(-2*cot(2*x + pi/3)**2 - 2))*(-sqrt(3)*(-sqrt(3)*cot(2*x + pi/3) + 1)*(-4*cot(2*x + pi/3)**2 - 4)*cot(2*x + pi/3) + sqrt(3)*(-4*cot(2*x + pi/3)**2 - 4)*cot(2*x + pi/3) + (-2*cot(2*x + pi/3)**2 - 2)*(3*cot(2*x + pi/3)**2 + 3))*((-sqrt(3)*cot(2*x + pi/3) + 1)**2/2 + sqrt(3)*cot(2*x + pi/3) + 
1)**(-2.5) - sqrt(2)*((-sqrt(3)*cot(2*x + pi/3) + 1)**2/2 + sqrt(3)*cot(2*x + pi/3) + 1)**(-1.5)*(-sqrt(3)*(-sqrt(3)*cot(2*x + pi/3) + 1)*(-4*cot(2*x + pi/3)**2 - 4)*(-2*cot(2*x + pi/3)**2 - 2) + 4*sqrt(3)*(-sqrt(3)*cot(2*x + pi/3) + 1)*(-4*cot(2*x + pi/3)**2 - 4)*cot(2*x + pi/3)**2 + (-6*cot(2*x + pi/3)**2 - 6)*(-4*cot(2*x + pi/3)**2 - 4)*cot(2*x + pi/3) + 3*(-4*cot(2*x + pi/3)**2 - 4)*(-2*cot(2*x + pi/3)**2 - 2)*cot(2*x + pi/3) + sqrt(3)*(-4*cot(2*x + pi/3)**2 - 4)*(-2*cot(2*x + pi/3)**2 - 2) - 2*(-4*cot(2*x + pi/3)**2 - 4)*(3*cot(2*x + pi/3)**2 + 3)*cot(2*x + pi/3) - 4*sqrt(3)*(-4*cot(2*x + pi/3)**2 - 4)*cot(2*x + pi/3)**2)) + 1
 
sa5 = sa4.diff(x)
 
plot(sa5, (x, 0, pi/12), ylim=[0, 10])

Et j'obtien cette erreur:

Code:

RecursionError: maximum recursion depth exceeded while calling a Python object

Je me demande si j'ai atteint un plafond ou si je peux pousser la machine un peu plus loin?

Salut,

Entrer le message d'erreur dans votre moteur de recherche préféré donnera des réponses.

- W

Je ne comprends pas trop l'intéret d'écrire une fonction comme ca. Plein de chose se simplifie ...
La fonction peut s'écrire simplement (code donné en python classique ici, je vous laisse regarder en sympy) :
Code:

1 2 3 4 5 6 from math import sqrt,pi from numpy import arctan as cot def sa4(x) : cc = cot(2*x+pi/3) return 1 + sqrt(3)*4/3*(cc**2 + 1)**(-0.5)/cc