The Davydov's soliton propagation in the linear polymer chain is analyzed numerically by solving the discrete equations of motion for exciton and phonon amplitudes. The main difference with respect to the results of continual approximation is that two exciton–phonon coupling constants influence separately the soliton behavior. Their influence is studied both in the ideal chain and the chain with single impurity.