dv= dpz − drx + pqρ sin θ + qrρ sin θ,
dt dtdt

(16)

dw= dqx − q2ρ cos θ.
dt dt

(17)

The dynamical equations (4)...(9) can now be reduced to

X= − dqz − p2z cotθ + q2 (x − ρ sin θ) + prx cot θ − g sin θ,
M dt

(18)

Y= dpz − drx − pq (x + z cot θ − ρ sin θ) + qrp cos θ,
M dtdt

(19)

Z= dqx + q2 (z − ρ cos θ) + p2z − prx − g cos θ,
M dt