flotox - Bibliothèque Ecole Centrale Lyon - École Centrale de Lyon

De l'ensemble des relations (1), (2), (3), on déduitdQ K2L2+Q 2 B2L2 K1L1Q2= 2 S2 (L2expK2 L2Q) - L.,8iL1exp C- ) (4)K1 L1 + Qlors de l'effacement, l'équation différentielle s'obtientde la même manièredQK1L1-i-Q2= a'1 S1 CL1expL$1 iK1 L1Q) - '2 S2L22expB' L22K2 L2Q)les équations différentielles (4) et (5) sont résolues de manièrenumérique par une méthode de Runge Kutta à quatre approximations.(4) et (5) s'écriventdQ = f (t, Q)dtsoient Q0 to les conditions initiales et t la raison de laprogression arithmétique en t pour laquelle on secalculer yon pose k1 = t f(t0, Q)propose dek2 = t f(t0 + 1/2 t, Q0 + 1/2 k1)k3 = t f(t + 1/2 t, Q0 + 1/2 k2)k = ,t f(t + t4 0 Q0 + k3)le résultat approché est alorsQ(t + t) = Q + (k1 + 2 k2 + 2 k3 + k4)