“…The equations of motion for the coupled system, Equation (33), and all the expressions for the forces and moments on the right‐hand side of Equation (33), constitute an initial value problem for a pair of nonlinear second order differential equations in
,
,
, and
. Similar initial value problems have previously been solved with Runge Kutta method,
23 Adams–Bashforth–Moulton (multi‐step) methods, Newton Raphson method,
3‐5 Quasi Newton methods (Broyden),
7,29 and modified extended backward differentiation method
20‐22,30 . In this article, ode15s is used to solve the stiff nonlinear initial value problem.…”