Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

I was solving Laplace equation in MATLAB numerically. However I have problems when the domain is not rectangular.

The equation is as follows:

$$ \frac{\partial^2 u}{\partial x^2}+\frac{\partial^2 u}{\partial y^2}=0 $$

domain is circular

$$ x^2 + y^2 < 16 $$

and boundary condition $$ u(x,y)= x^2y^2 $$

How should I start with solving this equation numerically ?

share|cite|improve this question
I have the same problem as you had. Could you help?. It's urgent. I need this program to solve some problems, but I'm bad using matlab. – user40196 Sep 13 '12 at 2:54
@BlancaHernándezGalván Please don't post comments as answers. – Bill Dubuque Sep 13 '12 at 3:16
up vote 2 down vote accepted

Perhaps begin by rewriting the problem in polar coordintates: $$\frac{\partial^{2}u}{\partial r^{2}}+\frac{1}{r}\frac{\partial u}{\partial r}+\frac{1}{r^{2}}\frac{\partial^{2}u}{\partial\theta^{2}}=0$$ $$r^2<16$$ $$\left. u(r,\theta)\right|_D=r^4\cos^2\theta\sin^2\theta$$

share|cite|improve this answer
But be careful around $r=0$ and that the function is periodic in $\theta$, $u(r,\theta)=u(r,\theta+2\pi)$. – Daryl Sep 1 '12 at 8:45

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.