  We present a brief review of general results about non-rotating neutron stars
in simple R-square gravity and its extension with a scalar axion field.
Modified Einstein equations are presented for metrics in isotropical
coordinates. The mass-radius relation, mass profile and dependence of mass from
central density on various equations of state are given in comparison to
general relativity.
