  We investigate f-mode oscillations of static anisotropic stable neutron stars
within the framework of full general relativity. We present equations governing
unperturbed stellar structures and oscillations with an ansatz to account for
the anisotropy. We solve those equations for two different equations of states.
We see that, moderately anisotropic neutron stars with the tangential pressure
larger than the radial pressure can give more massive neutron stars than the
isotropic or very anisotropic ones. We find that the frequency of the f-mode
exhibits a linear relationship with the square root of the average density of
the stars and the slope of the fit depends on the anisotropic strength. For any
given value of the anisotropic strength, the frequency increases with the
increase of the mass of the neutron star, linearly for lower masses, and
rapidly at higher masses. However, this non-linear rise in the frequency with
the mass is not prominent when the radial pressure is larger than the
tangential pressure. For a fixed value of a small mass, higher anisotropy leads
to a larger value of the frequency, but when the fixed mass is above a
threshold value, higher anisotropy leads to a smaller value of the frequency.
The nature of the variation in the frequency with the change in the anisotropic
strength is similar for the two equations of state, but for a fixed mass and
the same amount of the anisotropy, the softer equations of state gives higher
frequency. We also find that the damping time of the f-mode oscillation
decreases as the mass of the neutron star increases for all values of the
anisotropic strength. For a fixed mass of the neutron star and for the same
amount of the anisotropy, the value of the damping time is lower for the softer
equation of state, but the nature of the variation in the damping time with the
change in the anisotropic strength is similar.
