We investigate the upper mass limit predictions of baryonic for static neutron stars in context $f(R)$ gravity. use most popular gravity model, namely $R^2$ gravity, and calculate maximum baryon adopting several realistic equations state one ideal equation state, that causal limit. Our motivation is based on fact with masses larger than star configurations inevitably collapse to black holes. Th...