Integral equation method is powerful to handle the PDE problem with nice local symmetry. Usually, it is coupled with some other numerical techniques to deal with more complicated PDE problems. In this talk, I will present a discrete version of integral equation method, formulating which through an operator point of view. Its application to PDE problems and lattice problems will be reported.