  We evaluate the Majorana phases for a general $3\times3$ complex symmetric
neutrino mass matrix on the basis of Mohapatra-Rodejohann's phase convention
using the three rephasing invariant quantities $I_{12}$, $I_{13}$ and $I_{23}$
proposed by Sarkar and Singh. We find them interesting as they allow us to
evaluate each Majorana phase in a model independent way even if one eigenvalue
is zero. Utilizing the solution of a general complex symmetric mass matrix for
eigenvalues and mixing angles we determine the Majorana phases for both the
hierarchies, normal and inverted, taking into account the constraints from
neutrino oscillation global fit data as well as bound on the sum of the three
light neutrino masses ($\Sigma_im_i$) and the neutrinoless double beta decay
($\beta\beta_{0\nu}$) parameter $|m_{11}|$. This methodology of finding the
Majorana phases is applied thereafter in some predictive models for both the
hierarchical cases (normal and inverted) to evaluate the corresponding Majorana
phases and it is shown that all the sub cases presented in inverted hierarchy
section can be realized in a model with texture zeros and scaling ansatz within
the framework of inverse seesaw although one of the sub case following the
normal hierarchy is yet to be established. Except the case of quasi degenerate
neutrinos, the methodology obtained in this work is able to evaluate the
corresponding Majorana phases, given any model of neutrino masses.
