In order to solve the nonlinear inverse Camassa-Holm equations, a numerical method is developed by applying finite difference formula to time discrimination and collocation of polynomial scaling functions for spatial ...