Réponse :
Voici mes matrices (mu et nu) si vous voulez les voir.
#Paramétrisation des espaces de départ et d'arrivée
s1 = np.linspace(-5,5,N, False)
s2 = np.linspace(-5,5,N, False)
S1,S2 = np.meshgrid(s1,s2) #points de discrétisation de l'espace de départ
t1 = np.linspace(-5,5,N, False)
t2 = np.linspace(-5,5,N, False)
T1,T2 = np.meshgrid(t1,t2) #points de discrétisation de l'espace de d'arrivé
#Définition des mesures
def n1(x1,x2):
if np.exp(-(x1**2+x2**2)/2) ==0:
return epsilon
else:
return np.exp(-(x1**2+x2**2)/2)/2*np.pi
mu = fmu(S1,S2)
mu = mu/np.sum(mu)
def n1(x1,x2):
if np.exp(-(x1**2+x2**2)/2) ==0:
return epsilon
else:
return np.exp(-(x1**2+x2**2)/4)/4*np.pi
nu = fnu(T1,T2)
#print(np.sum(nu))
nu = nu/np.sum(nu)
