K. Boshkayev, T. Konysbayev, Ye. Kurmanov, M. Muccino, H. Quevedo

We numerically analyze quasiperiodic oscillations (QPOs) using a well-established spacetime model with neutron star sources. Within the framework of general relativity, we present expressions for the fundamental frequencies of test particles in the gravitational field of a slowly rotating and slightly deformed compact object defined by the Hartle-Thorne (HT) metric. Using the Relativistic Precession Model (RPM) formulated by Stella and Morinsk, we examine quasiperiodic oscillation data from eight neutron stars in low-mass X-ray binary systems. Employing Markov Chain Monte Carlo analyses with the Metropolis-Hastings algorithm, we estimate 1-$\sigma$ and 2-$\sigma$ error bars. Finally, we compare our results with predictions from the Schwarzschild, Lense-Thirring, and Kerr metrics, demonstrating that three of the eight sources can be well explained within the Hartle-Thorne model.