“…> Y := proc(y, z) synchronization: average over a period of first equation according to (22) > local a, b, c, a1, a2, a3; > global V e , v e , 0 , e 0 , e n , n , A; > a := H (v)/v; > b := diff (H (v), v); > c := unapply(a − b, v); > a1 := λ( 0 , n , e 0 , e n , A); > a2 := α( n , e 0 , e n ); > a3 := (1/a2 + 1 + v e )/y; > −(a1/( pi * c(v e ) * y 2 * a2 2 )) * I J (z, a3); > end proc; > Z := proc(y, z) synchronization: average over a period of second equation according to (23) > local a, b, c, a1, a2, a3; > global η, V e , v e , 0 , e 0 , e n , n , A; (H (v), v); > c := unapply(a − b, v); > a1 := λ( 0 , n , e 0 , e n , A); > a2 := α( n , e 0 , e n ); > a3 := (1/a2 + 1 + v e )/y; > (η/2) − (a1/(2 * pi * c(v e ) * y 3 * a2 2 )) * I I (z, a3); > end proc; > S := proc() procedure giving the solution y 0 = (y 0 , z 0 ) of equations of synchronization and matrix S = S(y 0 ) = S(y 0 , z 0 ) of synchronization defined by (12) > local sol, a, b, g, g1, g2, g3, g4, g5, g6, c; > global Y, Z , h, G, H, v e , η; > sol := fsolve(Y (y, z), Z (y, z), y, z); > a := eval(y, sol); First component y 0 of y 0 > b := eval(z, sol); Second component z 0 of y 0 > c := G()(v e ) > g := unapply(η * y * cos(z + φ) − (2 * H (v e + y * cos(z + φ), φ)/c), y, z, φ); > g5 := unapply(sin(z + φ) * g(y, z, φ), y, z, φ); > g6 := unapply((1/y) * cos(z + φ) * g(y, z, φ), y, z, φ); > g1 := unapply(D [1](g5)(a, b, φ), φ); > g1 := value(Int(g1(φ), φ = 0..Pi)/(2 * Pi)); > g1 := evalf(%); > g2 := unapply(D [2](g5)(a, b, φ), φ); > g2 := value(Int(g2(φ), φ = 0..Pi)/(2 * Pi)); > g2 := evalf(%); > g3 := unapply(D [1](g6)(a, b, φ), φ); > g3 := value(Int(g3(φ), φ = 0..Pi)/(2 * Pi)); > g3 := evalf(%); > g4 := unapply(D [2](g6)(a, b, φ), φ); > g4 := value(Int(g4(φ), φ = 0..Pi)/(2 * Pi)); > g4 := evalf(%); > Matrix([[g1, g2], [g3, g4]]); Matrix S according to (12) > end proc;…”