  We propose a method for reconstruction of the optical potential from
scattering data. The algorithm is a two-step procedure. In the first step the
real part of the potential is determined analytically via solution of the
Marchenko equation. At this point we use a diagonal Pad\'{e} approximant of the
corresponding unitary $S$-matrix. In the second step the imaginary part of the
potential is determined via the phase equation of the variable phase approach.
We assume that the real and the imaginary parts of the optical potential are
proportional. We use the phase equation to calculate the proportionality
coefficient. A numerical algorithm is developed for a single and for coupled
partial waves. The developed procedure is applied to analysis of $^{1}S_{0}$
$NN$, $^{3}SD_{1}$ $NN$, $P31$ $\pi^{-} N$ and $S01$ $K^{+}N$ data.
