J'ai programmé la méthode d'euler implicite en python et maintenant j'aimerais l'appliquer sur des vecteurs pour le systèmes d'équations différentielle
Sauf que ça marche pas j'obtiens une liste de termes tous constant

def euler_vectorielle(F,t0,tmax,y0,n) :
    h = (tmax-t0)/(n-1) 
    t = t0
    y1,y2 = y0
    T = [t0 + k*h for k in range(0, n)]
    Y = [y0 for k in range(0, n)]
    for k in range(1, len(T)) :
        y1 = y1 + h*F([y1,y2],t)[0]
        y2 = y2 + h*F([y1,y2],t)[1]
        t = t+h
        Y[k][0] = y1 
        Y[k][1] = y2
        print(Y)
    return (Y,T)

regarde si t'as pas pris des conditions initiales qui font que c'est constant?
(j'ai vaguement check ton code je crois pas voir d'erreur, si ce n'est possiblement F du coup mais c'est extérieur à ta définition)

Bah F j'avais déjà utilisé cette fonction et les mêmes conditions initiale pour résoudre l'équation avec le module odeint et ça marchait bien

je te conseille de taper ce qu'il y a dans ta définition ligne par ligne histoire de voir à quel moment ça bloque

C'est Euler explicite ça

Oui je voulais dire explicite my bad