We develop the discrete derivatives representation method (DDR) to find the physical structures of the Schrödinger equation in which the interpolation polynomial of Bernstein has been used. In this paper the particle swarm optimization (PSO for short) has been suggested as a means to improve qualitatively the solutions. This approach is carefully handled and tested with a numerical example.