The matrices in 17,27 should appear as: 
x(t)=(αγcos(ωtϕ)ωαγsin(ωtϕ)ωtmod2π)T
17
 
dxdτ=f(x)=(εv2ζεvεx+εV2(1cos(2θ))2(1x)21)T
27