We study a radial basis function collocation method (RBFCM) to discretize a coupled nonlinear Schr$ddot{o}$dinger equation (CNLSE) that governs a two dimensional rotating Bose-Einstein condensate (BEC) with an angular momentum rotation term. We exploit a RBFCM-continuation method (RBFCM-CM) to trace the solution curve of the CNLSE. We compare the performance of the RBFCM-CM with the FEM-CM. We observe that the RBFCM-CM is very robust in a coarse grid for resolving the ground state solution with many vortices when the angular momentum rotation is close to the limit. Numerical results demonstrate the efficiency and accuracy of the RBFCM-CM for computing the superfluid density of the ground level of the BEC.