An efficient B-spline difference method for solving system of nonlinear parabolic PDEs