The program solves nonlinear Navier–Stokes equations by using the Galerkin method for analytical expressions as proposed in [1]. The method was first implemented in FORTRAN in [2].
References
[1] M. Kawaguti, "The Critical Reynolds Number for the Flow past a Sphere," Journal of the Physical Society of Japan, 10(8), 1955 pp. 694–699. doi:10.1143/JPSJ.10.694.
[2] C. Y. Chow, An Introduction to Computational Fluid Mechanics, New York: Wiley, 1979.