This question is rather simple but I have some difficulties to find code. Let us suppose that I wrote a routine, in a given language, that computes the evolution of a particle doing Brownian motion in one dimension. Given all the steps, how should I write an algorithm to derive the heat kernel from them?
Some good references are fine as an answer.