# Random Walk Solutions to the Dirichlet Problem for the Laplace Equation

Requires a Wolfram Notebook System

Interact on desktop, mobile and cloud with the free Wolfram Player or other Wolfram Language products.

The Dirichlet problem seeks to find the solution to a partial differential equation inside a domain , with prescribed values on the boundary of . In 1944, Kakutani showed that the Dirichlet problem for the Laplace equation can be solved using random walks as follows. Given a point in the interior of , generate random walks that start at and end when they reach the boundary of . Then compute the average of the values of the given function at these boundary points. This average value is approximately equal to the value of the solution to the Dirichlet problem at the point .

[more]
Contributed by: Cameron Nachreiner (March 2019)

Based on an undergraduate research project at the Illinois Geometry Lab by Yuheng Chang, Baihe Duan, Yirui Luo, Yitao Meng, Cameron Nachreiner and Yiyin Shen and directed by A. J. Hildebrand.

Open content licensed under CC BY-NC-SA

## Snapshots

## Details

The random walks are generated using steps given by a two-dimensional normal distribution with standard deviation equal to the "scale" parameter. As the scale approaches 0, these random walks approach a Brownian motion starting at the point and ending at the boundary of , and the average value of the function at the endpoints of these random walks approaches the expected value of this stopped Brownian motion as the number of random walks approaches infinity. See [1].

Reference

[1] R. Hersh and R. J. Griego, "Brownian Motion and Potential Theory," *Scientific American*, 220(3), 1969 pp. 66–74. doi:10.1038/scientificamerican0369-66.

## Permanent Citation