Finding of a numerical solution tothe Cauchy - Dirichlet problem for Boussinesq - Lo`ve equation using finitedifferences method

The article is devoted to the numerical investigation of Boussinesq - L`ove mathematical model. Algorithm for finding numerical solution to the Cauchy - Dirichlet problem for Boussinesq - Lo`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.

Boussinesq – L`ove equation, Cauchy – Dirichlet problem, finite differences method, Sobolev type equation, phase space, conditions of data consistency, system of difference equations, Thomas algorithm