Ordinary differential equation algorithms for a frequency-domain water wave Green’s function

Our study investigated an algorithm for a second-order ordinary differential equation for the frequency-domain Green’s function of free-surface waves in water of infinite depth. A power-series method is introduced if the wave frequency ω<1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega <1$$\end{document}; otherwise, if ω≥1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega \ge 1$$\end{document}, then the trigonometrically fitted block Numerov-type method (TBNM) is employed. The calculation precision of the power-series method and the TBNM reached 10-7\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10^{-7}$$\end{document} and 10-6\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10^{-6}$$\end{document}, respectively. The two methods have a high calculation efficiency compared with calculating the Green’s function using the series expansion representation approach. The calculation speed for these two methods is 15 times faster using the same computing codes.