Discussion :

[Singezare]

Bonsoir.

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

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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?

[wiztricks]

Salut,

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

- W

[lg_53]

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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