radial-shear, azimuthal-shear and axial-inflow)axial-inflow row stores the positive constant $-H_\infty$, so that the dimensional far-field velocity is $-\sqrt{\nu\Omega}\,(-H_\infty)$.from mpmath import mp
mp.dps = 50
# The attached generate.py solves the two shooting equations for F'(0)
# and G'(0), then reads H at the same finite endpoint.The values are computed by shooting on $F'(0)$ and $G'(0)$ in the finite-interval boundary-value problem.
The generator solves the same shooting equations at two working precisions, requires the 30 stored digits to agree, checks the first derivative of the integrated $H$ against $H'=-2F$, and checks the rounded values against the constants printed by White, Mohr and Newman [3] and the no-slip case of Miklavčič and Wang [4].