This research describes an efficient numerical method based on Wendland’s compactly supported functions to simulate the time-space fractional coupled nonlinear Schrödinger (TSFCNLS) equations. Here, the time and space fractional derivatives are considered in terms of Caputo and Conformable derivatives, respectively. The present numerical discussion is based on the following ways: we first approximate the Caputo fractional derivative of the proposed equation by a scheme order O(∆t2−α), 0 < α < 1 and then the Crank-Nicolson scheme is employed in the mentioned equation to discretize the equations. Second, applying a linear difference scheme to avoid solving nonlinear systems. In this way, we have a linear, suitable calculation scheme. Then, the conformable fractional derivatives of the Wendland’s compactly supported functions are established for the scheme. The stability analysis of the suggested scheme is also examined in a similar way to the classic Von-Neumann technique for the governing equations. The efficiency and accuracy of the present method are verified by solving two examples.