Exact numerical methods and stochastic simulation methods are developed to study the force stretching single RNA issue on the secondary structure level in equilibrium. By computing the force-extension curves on the constant force and the constant extension ensembles, we find the two independent methods agree with each other quite well. To show the precision of our methods in predicting unfolding experiments, the unfolding forces of different RNA molecules under different experimental conditions are calculated. We find that the ionic corrections on the RNA free energies alone might not account for the apparent differences between the theoretical calculations and the experimental data; an ionic correction to the persistent length of single-stranded RNA should be necessary.