Python pour une edo

Bonjour,
j'essaye de programmer le système d'edo suivant avec Python:
\begin{align*}
U'(t)&=-a_1 U V_1 - a_2 U V_2;
\\
I_1'(t)&=a_1 U V_1 - \beta_1 I_1;
\\
I_2'(t)&=a_2 U V_2 - \beta_2 I_2;
\\
V_1'(t)&=b_1 I_{1_{\tau_1}}-\sigma_1 V_1;
\\
V_2'(t)&= \varepsilon I_{1_{\tau_1}}+b_2 I_{2_{\tau_2}}-\sigma_2 V_2,

\end{align*} L'objectif est de calculer $\log(V_1)$ et $\log(V_2)$ où $V_1$ et $V_2$ sont les solutions des deux dérnières équations.
Voici ce que j'ai écrit, mais j'ai l'erreur suivante.
RuntimeWarning: divide by zero encountered in log10 W2= np. log10(V2)
Je ne comprends pas cette erreur.
Voici le code.
Je vous remercie d'avance

import matplotlib.pyplot  as plt
import numpy as np
import os

path = PATH = str(os.path.abspath(''))

# Constantes
tau=2.
a1=10.
a2=10.
b1=5.
b2=5.
beta1=0.
beta2=0.
sigma1=0.4
sigma2=0.4
eps=0.1

# Pas de temps
h=0.01
ntau=int(tau/h)

#temps final
T=25

#M=nb de pas de temps
M=int(T/h)+1

# pre-dimensionnement du tableau & initialisation
U = np.zeros(M+1)
I1 = np.copy(U)
I2 = np.copy(U)
V1 = np.copy(U)
V2 = np.copy(U)
t = np.copy(U)

U[0]=1.
I1[0]=0.
I2[0]=0.
V1[0]=1.
V2[0]=0.

for n in range(1,M+1): # on commense au premier
    if n>ntau:
        V1[n] = (V1[n-1] + h*b1*I1[n-ntau])/(1+h*sigma1)
        V2[n]= (V2[n-1]+eps*h*I1[n-ntau]+b2*h*I2[n-ntau])/              (1+sigma2*h)
    else:
        V1[n] = V1[n-1]/(1+h*sigma1)
        V2[n]= V2[n-1]/(1+h*sigma2)

    U[n]=U[n-1]/(1+a1*h*V1[n]+a2*h*V2[n])
    I1[n] = (I1[n-1] + a1*h*U[n]*V1[n])/(1+beta1*h)
    I2[n] = (I2[n-1] + a2*h*U[n]*V2[n])/(1+beta2*h)
    t[n] = t[n-1] + h

#fichier = 'u.txt'
#tableau = np.transpose(np.vstack((t,U)))
#np.savetxt(path + '/' + fichier, tableau, delimiter=' ')

# passage de V1 à log10(V1)
W1 = np.log10(V1)

# passage de V2 à log10(V2)
W2 = np.log10(V2)

#fichier = 'i.txt'
#tableau = np.transpose(np.vstack((t,I)))
#np.savetxt(path + '/' + fichier, tableau, delimiter=' ')

#fichier = 'v.txt'
#tableau = np.transpose(np.vstack((t,V)))
#np.savetxt(path + '/' + fichier, tableau, delimiter=' ')

fichier = 'w1.txt'
tableau = np.transpose(np.vstack((t,W1)))
np.savetxt(path + '/' + fichier, tableau, delimiter=' ')

fichier = 'w2.txt'
tableau = np.transpose(np.vstack((t,W2)))
np.savetxt(path + '/' + fichier, tableau, delimiter=' ')