We present a method for retrieving of single-active electron potential in an atom or molecule from a given momentum distribution of photoelectrons ionized by a strong laser field. In this method the potential varying within certain limits is found as the result of the optimization procedure aimed at reproducing the given momentum distribution. The optimization using numerical solution of the time-dependent Schrödinger equation for ionization of a model one-dimensional atom shows the good accuracy of the potential reconstruction method. This applies to different ways used for representing of the potential under reconstruction, including a parametrization and determination of the potential by specifying its values on a spatial grid.