Harmonic Functions

Work in progress. This note is still being written and incomplete.

Let DD be a domain in the xyxy-plane. A real-valued function u(x,y)u(x,y) is harmonic in DD iff:

  • uu has continuous second-order partial derivatives in DD.
  • uu satisfies Laplace’s equation throughout DD:
∇2u=uxx+uyy=0\nabla^2 u = u_{xx} + u_{yy} = 0

Component Functions

If f(z)=u(x,y)+iv(x,y)f(z) = u(x,y) + iv(x,y) is analytic in a domain DD, then uu and vv are harmonic in DD.

Harmonic Conjugate

Let uu and vv be real-valued functions on a region D⊆R2D \subseteq \mathbb{R}^2. vv is a harmonic conjugate of uu iff f(z)=u(x,y)+iv(x,y)f(z) = u(x,y) + iv(x,y) is analytic on DD.

A harmonic conjugate of uu exists only if uu is harmonic on DD. A harmonic uu does not always have one.

If DD is simply connected, a harmonic conjugate always exists and is unique up to an additive constant. Otherwise it may fail to exist.

For example: u=ln⁡∣z∣u = \ln|z| is harmonic on C∖{0}\mathbb{C} \setminus \{0\}, but the only candidate v=arg⁡zv = \arg z is multivalued, so no harmonic conjugate exists there.

Finding a Conjugate

Given a harmonic u(x,y)u(x,y) on a simply connected domain:

  1. Check that uxx+uyy=0u_{xx} + u_{yy} = 0.
  2. Compute uxu_x and uyu_y.
  3. Use the Cauchy-Riemann equation vy=uxv_y = u_x. Integrate with respect to yy to get v=∫ux dy+ϕ(x)v = \int u_x \, dy + \phi(x), where ϕ(x)\phi(x) is an unknown function of xx alone.
  4. Differentiate this vv with respect to xx and equate it to −uy-u_y, from the Cauchy-Riemann equation vx=−uyv_x = -u_y. All terms containing yy cancel because uu is harmonic, which leaves an expression for ϕ′(x)\phi'(x).
  5. Integrate ϕ′(x)\phi'(x) to get ϕ(x)\phi(x), including an arbitrary constant CC.
  6. Substitute ϕ(x)\phi(x) into vv. Then f=u+ivf = u + iv is analytic, and can be written in terms of zz.

Integrating vx=−uyv_x = -u_y with respect to xx first and adding an unknown function of yy works the same way.

Worked Example

Find a harmonic conjugate vv of u(x,y)=y3−3x2yu(x,y) = y^3 - 3x^2y, and the corresponding analytic function ff.

Solution Reveal Hide
ux=−6xy,uy=3y2−3x2u_x = -6xy, \qquad u_y = 3y^2 - 3x^2

Cauchy-Riemann vy=uxv_y = u_x:

vy=−6xy  ⟹  v=−3xy2+ϕ(x)v_y = -6xy \implies v = -3xy^2 + \phi(x)

Cauchy-Riemann vx=−uyv_x = -u_y:

−3y2+ϕ′(x)=−3y2+3x2  ⟹  ϕ′(x)=3x2-3y^2 + \phi'(x) = -3y^2 + 3x^2 \implies \phi'(x) = 3x^2ϕ(x)=x3+C\phi(x) = x^3 + Cv(x,y)=x3−3xy2+Cv(x,y) = x^3 - 3xy^2 + Cf(z)=(y3−3x2y)+i(x3−3xy2+C)=iz3+iCf(z) = (y^3 - 3x^2y) + i(x^3 - 3xy^2 + C) = iz^3 + iC
Written by October 17, 2026 2 min read
Was this helpful?