Solution of the Laplace Equation for Temperature Distribution in a Square

Consider the Laplace equation in a square region, where . We wish to solve for the temperature distribution , subject to the following Dirichlet boundary conditions:
BC1: ,
BC2: ,
BC3: ,
BC4: .
This Demonstration solves the problem using a spectral method [1] for a user-set value of the number of the Chebyshev–Gauss–Lobatto collocation points. The colored dots reflect the intensity of the scalar field at the various nodes in the square domain. You can clearly recognize the form of a contour plot of in the third snapshot.


  • [Snapshot]
  • [Snapshot]
  • [Snapshot]


[1] L. N. Trefethen, Spectral Methods in MATLAB, Philadelphia: SIAM, 2000.
    • Share:

Embed Interactive Demonstration New!

Just copy and paste this snippet of JavaScript code into your website or blog to put the live Demonstration on your site. More details »

Files require Wolfram CDF Player or Mathematica.