We evaluate the Majorana phases for a general 3×3 complex symmetric neutrino mass matrix on the basis of Mohapatra–Rodejohann's phase convention using the three rephasing invariant quantities I12, I13 and I23 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 (Σimi) and the neutrinoless double beta decay (ββ0ν) parameter |m11|. 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 cases 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.