The nonlinear Wannier functions in square Kronig-Penney potential are investigated with the help of a set of exact nonlinear Bloch solutions. The nonlinear interaction makes the Wannier functions fall off as non-exponential law with distance and enhances the tunneling coupling between the neighbor sites.
