We present a sample of 1148 ab-type RR Lyrae (RRLab) variables identified from Catalina Surveys Data Release 1, combined with SDSS DR8 and LAMOST DR4 spectral data. We firstly use a large sample of 860 Galactic halo RRLab stars and derive the circular velocity distributions for the stellar halo. With the precise distances and carefully determined radial velocities (the center-of-mass radial velocities) by considering the pulsation of the RRLab stars in our sample, we can obtain a reliable and comparable stellar halo circular velocity curve. We take two different prescriptions for the velocity anisotropy parameter {beta} in the Jeans equation to study the circular velocity curve and mass profile. We test two different solar peculiar motions in our calculation. Our best result with the adopted solar peculiar motion 1 of (U, V, W) = (11.1, 12, 7.2) km/s is that the enclosed mass of the Milky Way within 50 kpc is (3.75 +/- 1.33) *10^11Msun based on beta = 0 and the circular velocity 180 +/- 31.92 (km/s) at 50 kpc. This result is consistent with dynamical model results, and it is also comparable to the previous similar works.