By applying finite difference formula to time discretization and the cubic B-SPLINEs for spatial variable, a numerical method for solving the inverse system of Burgers equations is presented. Also, the convergence analysis and stability for this problem are investigated and the order of convergence is obtained. By using two test problems, the accuracy of presented method is verified. Additionally, obtained numerical results of the cubic B-SPLINE method are compared to trigonometric cubic B-SPLINE method, exponential cubic Bspline method and radial basis function method. Implementation simplicity and less computational cost are the main advantages of proposed scheme compared to previous proposals.