We compare the results obtained from analytical perturbation theory and the AKM numerical code for an axistationary spacetime built from matching a rotating perfect fluid interior with the equation of state $\epsilon-3p=4B$ of the simple MIT bag model and an asymptotically flat exterior. We discuss the behaviour of the error in the metric components of the analytical approximation going to higher orders. Additionally, we check and comment the errors in multipole moments, central pressure and some other physical properties of the spacetime.