Harmonic functions on graphs and manifolds describe stationary distribution of random processes such as random walk and Brownian motion respectively. There are discrete, given by difference equation and continuous, given by differential equations models of definitions of harmonic functions, on graphs and manifolds respectively.
A function u on the vertices of a graph w/boundary is harmonic if its value at every interior vertex is the average of its values at neighboring vertices. That is,
a harmonic function on a manifold is a twice continuously differentiable function u : U → R (where U is an open subset of R2) which satisfies Laplace's equation:
Definition of the Dirichlet-to-Neumann operator for a domain
The harmonic functions satisfy the following properties:
mean-value property
The value of a harmonic function is a weighted average of its values at the neighbor vertices.
maximum principle
Corollary: the maximum (and the minimum) of a harmonic functions occurs on the boundary of the graph or the manifold.