Kelvin Ship-Wave Pattern
The displacement produced by a ship moving at constant speed produces a Kelvin wave pattern behind it.
The Kelvin ship-wave pattern can be calculated from an integral of the form , where and are the surface coordinates of the wave pattern. The integrand is oscillatory with stationary points.
Although Mathematica is able to perform the integration, here it is performed numerically using a midpoint formula to avoid the singularities at the limits. The stationary points are automatically handled with an asymptotic formula derived from analytical solutions at the stationary points.
The result is very robust and calculations can be done in seconds rather than minutes on a core i3 laptop.
 National Instutute of Standards and Technology, "Digital Library of Functions." (Jul 1, 2010) http://dlmf.nist.gov/.
"Kelvin Ship-Wave Pattern"
Wolfram Demonstrations Project
Published: September 8 2011