A high‐order harmonic polynomial method for solving the Laplace equation with complex boundaries and its application to free‐surface flows. Part I: Two‐dimensional cases. (30th June 2020)