  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. 
Reconstruction of the optical potential from scattering data