  We analyze the spatial and energy dependence of the local density of states
in a SNS junction. We model our system as a one-dimensional tight-binding chain
which we solve exactly by numerical diagonalization. We calculate the
dependence of the Andreev bound states on position, phase difference, gate
voltage, and coupling with the superconducting leads. Our results confirm the
physics predicted by certain analytical approximations, but reveal a much
richer set of phenomena beyond the grasp of these approximations, such as the
metamorphosis of the discrete states of the normal link (the normal bound
states) into Andreev bound states as the leads become superconducting.
