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.