  The supersymmetric solutions of PT-symmetric and Hermitian/non-Hermitian
forms of quantum systems are obtained by solving the Schrodinger equation for
the Exponential-Cosine Screened Coulomb potential. The Hamiltonian hierarchy
inspired variational method is used to obtain the approximate energy
eigenvalues and corresponding wave functions.
