Finding of a numerical solution to the Cauchy–Dirichlet problem for Boussinesq–Lòve equation using finite differences method

The article is devoted to the numerical investigation of Boussinesq–Lòve mathematical model. Algorithm for finding numerical solution to the Cauchy–Dirichlet problem for Boussinesq–Lòve equation modeling longitudinal oscillations in a thin elastic rod with regard to transverse inertia was obtained on the basis of phase space method and by using finite differences method. This problem can be reduced to the Cauchy problem for Sobolev type equation of the second order, which is not solvable for arbitrary initial values. The constructed algorithm includes additional check if initial data belongs to the phase space. The algorithm is implemented as a program in Matlab. The results of numerical experiments are obtained both in regular and degenerate cases. The graphs of obtained solutions are presented in each case.