In the present meshfree boundary-domain integral equation
method, where the matrices are non-symmetric, convergence
cannot be guaranteed, but the inequalities shown in Eq. (26)
should equally be observed in order not to introduce and additional
source of instability. Thus seven kinds of combinations of
the parameters with a fixed time step of Δt¼0.005 [34] are investigated
in details to obtain the proper ones.