Methods of finding quasi-normal modes of non-rotating relativistic stars have been well established, however, none of the existing treatments which take spacetime and fluid oscillations fully into account can determine modes of long decay time, e.g., the p and g mode series, or the f modes for stars with low compactness ratio (M/R). In this paper we show how the quasi-normal modes of long lifetime can be determined through refinements of a treatment originally due to Detweiler and Lindblom. The determination of the p mode series has been argued in the literature to have implication on the life time of gravitational wave sources and stellar stability. In this paper we 1) provide detailed steps in our treatment to facilitate future effort in this direction; 2) correct mistakes in the literature on the formulation; and 3) analyse the accuracy of the quasi-normal mode frequencies obtained and the limitations of the treatment.