[PYTHON] Système différentiel d'ordre 2

[Razorr]

Bonjour, j'ai un petit problème sur python, j'essaye de résoudre le système suivant

J'ai alors tapé le code suivant :

from scipy.integrate import odeint
from pylab import *

def epidemie2(syst, t):
    r=3 #r=beta
    a=1 #a=1/lambda
    D=1
    c=300
    syst=[S,I,y]
    dI=y
    dS=I*S*(r/c)
    dy=(-c*y+I*(a-r*S))/D
    return [dS,dI,dy]
    
t=linspace(0,30,0.001)

#conditions initales et résolution
S0=1
I0=10
y0=30
syst_CI=array([S0,I0,y0])
solution=odeint(epidemie2,syst_CI,t)


plot(solution,t)
show()

mais ça me renvoie aucune courbe, seulement un graphique blanc.
Si quelqu'un pouvait me venir en aide svp.
-----

[fred1599]

Que vaut la variable S dans la fonction epidemie2 ?

[Razorr]

Ah ouais y a un problème.
J'ai modifié

from scipy.integrate import odeint
from pylab import *

def epidemie2(syst, t):
    r=3 #r=beta
    a=1 #a=1/lambda
    D=1
    c=300
    S=syst[0]
    I=syst[1]
    dI=y
    dS=I*S*(r/c)
    dy=(-c*y+I*(a-r*S))/D
    return [dS,dI,dy]
    
t=linspace(0,30,0.001)

#conditions initales et résolution
S0=1
I0=10
y0=30
syst_CI=array([S0,I0,y0])
syst=odeint(epidemie2,syst_CI,t)


plot(syst,t)
show()

Toujours rien...

[fred1599]

dI = y, mais que vaut y ?

[A voir en vidéo sur Futura]

[Razorr]

En fait je veux mettre que y = la dérivée de I, mais je vois pas comment le rentrer en tant que variable. :/

[fred1599]

On fait des maths ou de l'informatique ? Parce-que si c'est des maths, je peux guère t'aider... Par contre si tu sais le faire en maths, je peux t'aider à le transcrire, mais dans ce cas il faut au moins que tu exprimes ce que tu veux algorithmiquement parlant, sinon ça rendra difficile la résolution de ton problème, si toi même tu n'arrives pas à résoudre cette équation sur papier.