Bonjour à tous,

Débutante en Python, je dois modéliser un problème biologique, à savoir la concentration de mélatonine. Le graph se présente sous forme de concentration en fonction du temps sur une durée de 48h.
Mon code se lance, mais la fenêtre supposé afficher le graph se ferme immédiatement.. Je ne comprends donc pas... Serait ce un soucis d'index out of range ?

Voici le code :

Merci beaucoup d'avance!

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
import numpy as np
import matplotlib.pyplot as plt
 
def func(N,n_steps):
    alpha=1.5
    k=0
    B=0
    for i in range(1,int(n_steps)):
        N[i]= N[i-1]+(alpha+k*N[i-1]-B*N[i-1])*dt
        if 0 < i < 4000 or 10000 < i <= 15000:
            k=0
            B+= 1/4000
 
        elif 4000 < i < 6000:
            k+=1/2000
            B=0
 
        else :
            k-=1/4000
            B=0
 
 
simulated_time= 48
dt= 0.001
n_steps = round(simulated_time / dt)
N=np.zeros(n_steps)
N[0]=8
time_range = dt*np.arange(n_steps)
func(N,n_steps)
plt.plot(time_range,N)