Solitary capillary-gravity wave

This function computes the steady irrotational surface solitary (classical and generalized, depending on the Bond number <> 1/3) capillary-gravity wave solutions of the full Euler equations with free surface (homogeneous, incompressible and perfect fluids). The wave is defined by its initial Froude and Bond numbers (Fr, Bo) and the result is about twelve digits accurate. The method works for all but the highest waves.
REFERENCE: D. Clamond, D. Dutykh & A. Duran. A plethora of generalised solitary gravity-capillary water waves. J. Fluid Mech., 2015 (https://hal.archives-ouvertes.fr/hal-01081798/)