squares best fit $\lambda \sim 0.04452$;

-   the one term MAE result based on $ C = \Delta_1(R_e) = (\log(4/R_e) - \gamma + 1/2)^{-1}$;

-   the two term MAE result based on $ C = \Delta_1(R_e)- 0.8669 \Delta_1^3(R_e)$.

The choice of constants in the MAE approach is itself somewhat arbitrary. That made by Kaplun (1957), giving Eq. (\[nearcyl\]), is little more than convenience. These results give support to the new model and the agreement with experiment we get by taking $\lambda$ about one half that implied by the MAE approach is rather tantalizing.

Iteration: a Stokes-Whitehead anomaly?
======================================

Another use for the type of low Reynolds numbers solution developed here is to provide a basis for [*iteration*]{} of the solution in powers of $R_e$. The presence of the paradox obstructing a global representation in the standard approach makes this impossible. Furthermore, we know from a corresponding analysis of the 3D spherical problem that even when the base solution makes sense, even the first iteration may fail: the [*Stokes-Whitehead*]{} paradox emerges. So it is of considerable interest to investigate the iteration of the system. In doing so we will keep $\epsilon$ general initially.

We summarize the solution developed thus far. With 