We construct complex integral step by step in this article.

Riemann Sums

We first review the classical Riemann integral and see how this motivates complex integral. Given an interval [a,b][a, b], the Riemann integral involves dividing an interval, denoted by [a,b][a, b], into several smaller sub-intervals.

Definition. A finite sequence P={x0,x1,…,xn}P = \{x_0, x_1, \ldots, x_n \} is a partition of [a,b][a, b], if

a=x0<x1<…<xn=b. a = x_0 < x_1 < \ldots < x_n = b.

The length of PP, denoted as ∥P∥\Vert P \Vert, is given by

∥P∥≐max⁡i∈[n]Δi \Vert P \Vert \doteq \max_{i \in [n] } \Delta_i

where Δi≐xi−xi−1\Delta_i \doteq x_i - x_{i - 1} for each i∈[n]i \in [n].

Given a function f:[a,b]→Rf : [a, b] \rightarrow \R, the Riemann integral computes the sum of the products of the length of each sub-interval and the value of ff at an arbitrary chosen point on that sub-interval. Such a sum is well-defined if it converges when the chosen partition is "refined" enough.

Definition. A number J∈RJ \in \R is the definite integral of ff over [a,b][a, b], if
∀ϵ>0\forall \epsilon > 0, ∃δ>0\exists \delta > 0 s.t., for every partition P={x0,x1,…,xn}P = \{x_0, x_1, \ldots, x_n \} of [a,b][a, b] with ∥P∥<δ\Vert P \Vert < \delta, and ∀ck∈[xk−1,xk]\forall c_k \in [x_{k - 1}, x_k], it holds that

∣∑k=1nf(ck)Δk−J∣<ϵ \left| \sum_{k = 1}^n f(c_k) \Delta_k - J \right| < \epsilon

There is a dedicated notation for the limit of Riemann sum

∫abf(x)dx=lim⁡∥P∥→0∑k=1nf(ck)Δk. \int_a^b f(x) dx = \lim_{\Vert P \Vert \rightarrow 0} \sum_{k = 1}^n f(c_k) \Delta_{k}.

Complex Integral

Complex integration is performed over curves, specifically over a set of curves that have favorable properties.

Smooth Curve

A point set γ\gamma in the C\mathbb{C} is a ``smooth curve'', if it is either a smooth arc or a smooth closed curve, defined as thus.

Definition. A point set γ\gamma in C\mathbb{C} is said to be a smooth arc if it is the range of some continuous complex-valued function z=z(t),t∈[a,b]z = z(t), t \in [a, b], s.t.,

  1. z(t)z(t) has a continuous derivative on [a,b][a, b],
  2. z(t)≠0,∀t∈[a,b]z(t) \neq 0, \forall t \in [a, b],
  3. z(t)z(t) is one-to-one on [a,b][a, b].

A point set γ\gamma is called a smooth closed curve if it is the range of some continuous function z=z(t),t∈[a,b]z = z(t), t \in [a, b], satisfying conditions (1) and (2) and the following:

  • z(t)z(t) is one-to-one on the half-open interval [a,b)[a, b), but z(b)=z(a)z(b) = z(a) and z′(b)=z′(a)z' (b) = z' (a).

Remark: in general, we refer to a function z=z(t)z = z(t), where t∈[a,b]t \in [a, b], as "admissible" if it satisfies the conditions (1), (2), and (3) / (⋅\cdot).

To explain the definition, let us consider the scenario where a curve γ\gamma is being drawn on a piece of paper using a pen. The position of the pen at each time t∈[a,b]t\in [a, b] is given by the function z(t)z(t), which parameterizes the curve.

To ensure that the curve is smooth, we must adhere to certain guidelines:

  1. We must not lift the pen from the paper during the drawing process, which implies that z(γt)z(\gamma_t) must be continuous.

  2. The velocity of the pen point must be well-defined and finite, and it must also vary continuously. Therefore, z′(t)z'(t) must exist and be continuous on the interval [a,b][a, b].

  3. To avoid sharp corners in the curve (which would require a sudden break in the drawing process), we require that z′(t)z'(t) is never zero.

  4. We also require that no point on the curve is drawn more than once, except for the initial and terminal points, which may coincide. In other words, z(t)z(t) must be one-to-one, except for the endpoints aa and bb.

In order to demonstrate that a collection of points γ\gamma in the complex plane C\mathbb{C} forms a smooth curve, we must present an admissible parametrization function z(t)z(t) whose range is γ\gamma. There can be infinite number of such parametrizations if γ\gamma is smooth.

Directed Smooth Curve

Contour

Definition. A contour is either a single point z0z_0 or a finite sequence of directed smooth curves (γ1,γ2,…,γn)(\gamma_1, \gamma_2,\ldots,\gamma_n) such that the terminal point of γk\gamma_k coincides with the initial point of γk+1\gamma_{k+1} for each k=1,2,…,n−1k = 1, 2, \ldots, n − 1. In this case one can write Γ=γ1+γ2+⋯+γn.\Gamma = \gamma_1 + \gamma_2 + \cdots + \gamma_n.

The (undirected) point set underlying a contour is known as a piecewise smooth curve

is said to be a closed contour or a loop if its initial and terminal points coincide. A simple closed contour is a closed contour with no multiple points other than its initial-terminal point; in other words, if z=z(t),a≤t≤bz = z(t), a \le t \le b, is a parametrization of the closed contour, then z(t)z(t) is one-to-one on the half-open interval [a,b)[a, b)

Contour Integral

We are ready to define the definite integral of a complex-valued function ff along a contour Γ\Gamma.

Directed Smooth Curve

We start by defining the integral along a directed smooth curve γ\gamma (with initial point α\alpha and terminal point β\beta), and then defining integral along a contour in terms of the integrals along its smooth components.

As in the construction of the Riemann integrals, we repeat the partition and product idea: let P={z0,z1,…,zn}⊆γP = \{z_0, z_1, \ldots, z_n\} \subseteq \gamma be a partition of the curve such that z0=αz_0 = \alpha and zn=βz_n = \beta. For each k∈[n]k \in [n], define Δk≐zk−zk−1\Delta_k \doteq z_k - z_{k - 1}, and let mesh of partition PP be the maximum length of the Δk\Delta_k:

∥P∥≐max⁡k∈[n]∥Δk∥. \Vert P \Vert \doteq \max_{k \in [n]} \Vert \Delta_k \Vert.

A Riemann sum associated with PP is just what it is in the real case:

∑k=1nf(dk)⋅Δk, \sum_{k = 1}^n f(d_k) \cdot \Delta_{k},

where dkd_k is an arbitrary point on the arc zk−1,zk^\widehat{z_{k - 1}, z_k}.

Definition. A function f:γ→Cf: \gamma \rightarrow \mathbb{C} is integrable over a directed smooth curve γ\gamma, if ∃J∈C\exists J \in \mathbb{C}, s.t., ∀ϵ>0\forall \epsilon > 0, ∃δ>0\exists\delta > 0 s.t., for every partition PP with ∥P∥<δ\Vert P \Vert < \delta and every choice of dkd_k in zk−1,zk^\widehat{z_{k - 1}, z_k}, it holds that

∣∑k=1nf(dk)Δk−J∣<ϵ. \left| \sum_{k = 1}^n f(d_k) \Delta_{k} - J \right| < \epsilon.

We call JJ the integral of ff over γ\gamma, and write it as

∫γf(x)dx≐J=lim⁡∥P∥→0∑k=1nf(dk)Δk. \int_\gamma f(x) dx \doteq J = \lim_{\Vert P \Vert \rightarrow 0} \sum_{k = 1}^n f(d_k) \Delta_{k}.
Properties

Before establishing conditions under which the path integral exists, we study its properties directly, via its definition.

Linearity.

  1. Constant Multiple: ∫γcf(x)dx=c∫γf(x)dx\int_\gamma c f(x) dx = c \int_\gamma f(x) dx.

  2. Sum: ∫γf(x)+g(x)dx=∫γf(x)dx+∫γg(x)dx\int_\gamma f(x) + g(x) dx = \int_\gamma f(x) dx + \int_\gamma g(x) dx.

Proof. For every partition PP of γ\gamma and for every choice of the points dkd_k,

∑k∈[n]cf(dk)Δk=c∑k∈[n]f(dk)Δk,∑k∈[n](f(dk)+g(zk))Δk=∑k∈[n]f(dk)Δk+∑k∈[n]g(zk)Δk. \begin{aligned} \sum_{k \in [n]} c f( d_k ) \Delta_k &= c \sum_{k \in [n]} f( d_k ) \Delta_k, \\ \sum_{k \in [n]} (f( d_k ) + g( z_k ) )\Delta_k &= \sum_{k \in [n]} f( d_k ) \Delta_k + \sum_{k \in [n]} g( z_k ) \Delta_k. \end{aligned}

Taking limit of both side, we get the desired answer.

Remark: given two sequence {a1,a2,…,}\{ a_1, a_2, \ldots, \} and {b1,b2,…,}\{b_1, b_2, \ldots, \}, if ai=bia_i = b_i for arbitrary i≥0i \ge 0 and lim⁡ibi=S\lim_i b_i = S, then lim⁡iai\lim_i a_i exits and lim⁡iai=S\lim_i a_i = S.

□\square

Existence

We provide a sufficient condition under which the path integral exists.

Lemma. Let γ\gamma be a directed smooth curve, and r:[a,b]→Cr: [a, b] \rightarrow \mathbb{C} be a parameterization of γ\gamma. For each partition P={z0,…,zn}P = \{z_0, \ldots, z_n\} of γ\gamma, define the inverse partition of [a,b][a, b] as

r−1(P)≐{r−1(z0),…,r−1(zn)}. r^{-1} (P) \doteq \{ r^{-1} (z_0), \ldots, r^{-1} (z_n) \}.

Then

∥P∥→0implies∥r−1(P)∥→0. \Vert P \Vert \rightarrow 0 \qquad \text{implies} \qquad \Vert r^{-1}(P) \Vert \rightarrow 0.

Remark: When γ\gamma is a directed smooth arc, then rr is a bijection and r−1(z)r^{-1}(z) is well-defined for each z∈γz \in \gamma. But when γ\gamma is a directed smooth closed curve, then rr is no longer a bijection, since r(a)=r(b)r(a) = r(b). But there is still a unique partition (t0,…,tn)⊆[a,b](t_0, \ldots, t_n) \subseteq [a, b], s.t., a=t0≤⋯≤tn=ba = t_0 \le \cdots \le t_n = b and r(ti)=zi,r(t_i) = z_i, for each i∈[0 . . n]i \in [0 \, . \, . \, n]. In this case, r−1(zi)≐ti,r^{-1}(z_i) \doteq t_i, for each i∈[0 . . n]i \in [0 \, . \, . \, n]. Therefore, r−1(P)r^{-1}(P) is well-defined.

Proof. Since rr is continuous on [a,b][a, b], it is uniformly continuous.

□\square

Theorem. If f:γ→Cf: \gamma \rightarrow \mathbb{C} is continuous on the directed smooth curve γ\gamma, then ff is integrable over γ\gamma and

∫γf(x) dx=∫abf(r(t)) r′(t) dt, \int_\gamma f (x) \, d x = \int_a^b f \big( r(t) \big) \, r'(t) \, d t,

where r:[a,b]→Cr : [a, b] \rightarrow \mathbb{C} is an parameterization of γ\gamma.

Proof. We first prove that f(r(t))⋅r′(t)f (r (t)) \cdot r'(t) is integrable, so the RHS of the equation is well-defined. Then we show that ∑k=1nf(dk)Δk\sum_{k = 1}^n f(d_k) \Delta_{k} converges to the RHS as ∥P∥→0\Vert P \Vert \rightarrow 0. Hence the limit lim⁡∥P∥→0∑k=1nf(dk)Δk\lim_{\Vert P \Vert \rightarrow 0} \sum_{k = 1}^n f(d_k) \Delta_{k} exists.

Integrability of f(r(t))⋅r′(t)f (r (t)) \cdot r'(t):

  • Observe that f∘r(t)=f(r(t)):[a,b]→Cf \circ r (t) = f(r(t)): [a, b] \rightarrow \mathbb{C} is continuous.

  • By assumption r′(t):[a,b]→Cr'(t) : [a, b] \rightarrow \mathbb{C} is continuous.

  • Therefore, f(r(t))⋅r′(t):[a,b]→Cf (r (t)) \cdot r'(t) : [a, b] \rightarrow \mathbb{C} is continuous.

  • Both the real and imaginary parts of f(r(t))⋅r′(t)f (r (t)) \cdot r'(t) are continuous, and hence integrable.

  • By definition of line integral,

    ∫abf(r(t)) r′(t) dt≐∫abRe(f(r(t)) r′(t)) dt+i⋅∫abIm(f(r(t)) r′(t)) dt \int_a^b f \big( r(t) \big) \, r'(t) \, d t \doteq \int_a^b \mathbb{Re} \Big( f \big( r(t) \big) \, r'(t) \Big) \, d t + i \cdot \int_a^b \mathbb{Im} \Big( f \big( r(t) \big) \, r'(t) \Big) \, d t

    is well-defined.

Proving The Equation:

Claim One. For each z∈γz \in \gamma, there is a unique point t∈[a,b]t \in [a, b] for which r(t)=zr(t) = z.

Proof of the Claim One.

□\square

We continue to prove the equation. Let P={z0,z1,…,zn}P = \{z_0, z_1, \ldots, z_n\} be a partition of γ\gamma. By Claim One, there exists unique points Q={a=t0<t1<…<tn=b}Q = \{ a = t_0 < t_1< \ldots < t_n = b \}, s.t., r(k)=zkr(k) = z_k, for k∈[0. . n]k \in [0 .\, .\, n].

□\square

To compute such integral, we can choose a description $\gamma(t) = x(t) + i y(t) : [a,b]→C[a, b] \rightarrow \mathbb{C} of curve γ\gamma. We partition γ\gamma by partitioning the interval [a,b][a, b] in usual way Pt={a=t0<t1<…<tn=b}P_t = \{ a = t_0 < t_1< \ldots < t_n = b \}. Then γ(Pt)≐{γ(t0),γ(t1),…,γ(tn)}\gamma(P_t) \doteq \{ \gamma(t_0), \gamma(t_1), \ldots, \gamma(t_n) \} is a partition of γ\gamma.

If γ\gamma is continuous, it is uniformly continuous on [a,b][a, b]. Therefore, as ∥Pt∥→0\Vert P_t\Vert \rightarrow 0, ∥γ(Pt)∥→0\Vert \gamma(P_t)\Vert \rightarrow 0. If ∫γf(z)dz\int_\gamma f(z) dz does exist, then

∫γf(z)dz=lim⁡∥Pt∥→0∑j=1nf(γ(ej))(γ(tj)−γ(tj−1)) \int_\gamma f(z) dz = \lim_{\Vert P_t\Vert \rightarrow 0} \sum_{j = 1}^n f \big( \gamma( e_j ) \big) \big(\gamma(t_j) - \gamma(t_{j - 1} ) \big)

for arbitrary choice ej∈[tj−1,tj]e_j \in [t_{j - 1}, t_j].

When the derivative of γ\gamma exists and is continuous, it is uniformly continuous on the closed bounded interval [a,b][a, b]. We now can rewrite the Riemann sum as

∑j=1nf(γ(ej))γ(tj)−γ(tj−1)tj−tj−1(tj−tj−1)=∑j=1nf(γ(ej))x(tj)−x(tj−1)+i(y(tj)−y(tj−1))tj−tj−1(tj−tj−1) \begin{aligned} &\sum_{j = 1}^n f \big( \gamma( e_j ) \big) \frac{\gamma(t_j) - \gamma(t_{j - 1} ) }{t_j - t_{j - 1} }(t_j - t_{j - 1} ) \\ &= \sum_{j = 1}^n f \big( \gamma( e_j ) \big) \frac{x(t_j) - x(t_{j - 1} )+ i (y(t_j) - y_(t_{j - 1} ) ) }{t_j - t_{j - 1} }(t_j - t_{j - 1} ) \end{aligned}

By mean value theorem, as x′(t)x'(t) and y′(t)y'(t) exist, there are some ej′∈[tj−1,tj]e_j' \in [t_{j - 1}, t_j] and ej′′∈[tj−1,tj]e_j'' \in [t_{j - 1}, t_j], such that

x(tj)−x(tj−1)tj−tj−1=x′(ej′),y(tj)−y(tj−1)tj−tj−1=x′(ej′′). \frac{x(t_j) - x(t_{j - 1} ) }{t_j - t_{j - 1} } = x'(e_j'), \qquad \frac{y(t_j) - y(t_{j - 1} ) }{t_j - t_{j - 1} } = x'(e_j'').

On the other hand, as x′(t)x'(t) and y′(t)y'(t) are uniformly continuous, for arbitrary ϵ>0\epsilon > 0, there is a δ\delta such that if ∥Pt∥<δ\Vert P_t\Vert < \delta, ∣x′(ej′)−x′(ej)∣<ϵ|x'(e_j') - x'(e_j)| < \epsilon and ∣y′(ej′′)−y′(ej)∣<ϵ|y'(e_j'') - y'(e_j)| < \epsilon. It follows that

x(tj)−x(tj−1)+i(y(tj)−y(tj−1))tj−tj−1=x′(ej)+iy′(ej)±ϵ(1+i)∑j=1nf(γ(ej))γ(tj)−γ(tj−1)tj−tj−1(tj−tj−1)=∑j=1nf(γ(ej))γ′(ej)(tj−tj−1)±ϵ(1+i)⋅∑j=1nf(γ(ej))(tj−tj−1), \begin{aligned} \frac{x(t_j) - x(t_{j - 1} )+ i (y(t_j) - y_(t_{j - 1} ) ) }{t_j - t_{j - 1} } &= x'(e_j) + i y'(e_j) \pm \epsilon(1 + i)\\ \sum_{j = 1}^n f \big( \gamma( e_j ) \big) \frac{\gamma(t_j) - \gamma(t_{j - 1} ) }{t_j - t_{j - 1} }(t_j - t_{j - 1} ) &= \sum_{j = 1}^n f \big( \gamma( e_j ) \big) \gamma'(e_j) (t_j - t_{j - 1} ) \\ &\quad \pm \epsilon (1 + i) \cdot \sum_{j = 1}^n f \big( \gamma( e_j ) \big) (t_j - t_{j - 1} ), \end{aligned}

Taking the limit of ∥Pt∥→0\Vert P_t\Vert \rightarrow 0, we get

∫γf(z)dz=∫abf(γ(t))γ′(t)dt, \int_\gamma f(z) dz = \int_a^b f(\gamma(t) )\gamma'(t) dt,

Remark 1: Can we relax the condition that γ′(t)\gamma'(t) is continuous and that [a,b][a, b] is a closed interval?

Remark 2: We have also showed that the integration is independent of the choice of γ\gamma, i.e., any representation of γ\gamma that is continuously differentiable gives the same integration value.

Anti-Derivative

Now suppose that F(z)=u(x,y)+iv(x,y)F(z) = u(x, y) + i v(x, y) and F′(z)=f(z)F'(z) = f(z).

Theorem. ∫γf(z)dz=F(γ(b))−F(γ(a))\int_\gamma f(z) dz = F(\gamma(b) ) - F(\gamma(a) ).

First, F′(γ(t))=∂u∂x∂x∂t+∂u∂y∂y∂t+i(∂v∂x∂x∂t+∂v∂y∂y∂t)F'(\gamma(t) ) = \frac{\partial u}{\partial x} \frac{\partial x}{\partial t} + \frac{\partial u}{\partial y} \frac{\partial y}{\partial t} + i ( \frac{\partial v}{\partial x} \frac{\partial x}{\partial t} + \frac{\partial v}{\partial y} \frac{\partial y}{\partial t} ). Therefore,

∫γf(z)dz=∫ab(∂u∂x∂x∂t+∂u∂y∂y∂t) dt+i∫ab(∂v∂x∂x∂t+∂v∂y∂y∂t) dt \int_\gamma f(z) dz = \int_a^b \left(\frac{\partial u}{\partial x} \frac{\partial x}{\partial t} + \frac{\partial u}{\partial y} \frac{\partial y}{\partial t} \right) \, dt + i \int_a^b \left( \frac{\partial v}{\partial x} \frac{\partial x}{\partial t} + \frac{\partial v}{\partial y} \frac{\partial y}{\partial t} \right) \, dt

But

u(x(b),y(b))−u(x(a),y(a))=∑j=1n[u(x(tj),y(tj))−u(x(tj−1),y(tj−1))]=∑j=1n[u(x(tj),y(tj))−u(x(tj),y(tj−1))+u(x(tj),y(tj−1))−u(x(tj−1),y(tj−1))]=∑j=1nu(x(tj),y(tj))−u(x(tj),y(tj−1))y(tj)−y(tj−1)y(tj)−y(tj−1)tj−tj−1(tj−tj−1)+∑j=1nu(x(tj),y(tj−1))−u(x(tj−1),y(tj−1))x(tj)−x(tj−1)x(tj)−x(tj−1)tj−tj−1(tj−tj−1)=∑j=1n[uy′(x(tj),yj∗)y′(ej′)+ux′(xj∗,y(tj−1))x′(ej′′)](tj−tj−1) \begin{aligned} &u(x(b), y(b) ) - u(x(a), y(a) ) \\ &= \sum_{j = 1}^n [u(x(t_j), y(t_j) ) - u(x(t_{j - 1} ), y(t_{j - 1} ) ) ] \\ &= \sum_{j = 1}^n [u(x(t_j), y(t_j) ) - u(x(t_j), y(t_{j - 1} ) ) + u(x(t_j), y(t_{j - 1} ) ) - u(x(t_{j - 1} ), y(t_{j - 1} ) ) ] \\ &= \sum_{j = 1}^n \frac{u(x(t_j), y(t_j) ) - u(x(t_j), y(t_{j - 1} ) ) }{y(t_j) - y(t_{j - 1} ) } \frac{y(t_j) - y(t_{j - 1} ) }{ t_j - t_{j - 1} }(t_j - t_{j - 1} )\\ &\quad +\sum_{j = 1}^n \frac{u(x(t_j), y(t_{j - 1} ) ) - u(x(t_{j - 1} ), y(t_{j - 1} ) ) }{x(t_j) - x(t_{j - 1} ) } \frac{x(t_j) - x(t_{j - 1} ) }{t_j - t_{j - 1} } (t_j - t_{j - 1} ) \\ &= \sum_{j = 1}^n [u_y' (x(t_j), y^*_j) y'(e_j') + u_x' (x^*_j, y(t_{j - 1} ) ) x'(e_j'') ] (t_j - t_{j - 1} )\\ \end{aligned}

where yj∗∈[y(tj−1),y(tj)]y^*_j \in [ y( t_{j - 1} ), y( t_j ) ], xj∗∈[x(tj−1),x(tj)]x^*_j \in [ x( t_{j - 1} ), x( t_j ) ], ej′,ej′′∈[tj−1,tj]e_j', e_j'' \in [ t_{j - 1}, t_j ]. If ux′,uy′,x′(t),y′(t)u'_x, u'_y, x'(t), y'(t) are continuous, then as ∥Pt∥→0\Vert P_t\Vert \rightarrow 0,

∑j=1n[uy′(x(tj),y(tj))y′(tj)+ux′(x(tj),y(tj))x′(tj)](tj−tj−1)=∑j=1n∂u∂t(tj−tj−1)→u(x(b),y(b))−u(x(a),y(a)). \begin{aligned} &\sum_{j = 1}^n [u_y' ( x( t_j ), y( t_j ) ) y'( t_j ) + u_x' ( x( t_j ), y( t_j ) ) x'( t_j ) ] ( t_j - t_{j - 1} ) \\ &= \sum_{j = 1}^n \frac{\partial u}{\partial t} ( t_j - t_{j - 1} ) \rightarrow u(x(b), y(b) ) - u(x(a), y(a) ). \end{aligned}

An immediate consequence is that

F(γ(b))−F(γ(a))=∫γF′(γ(t))dt F(\gamma(b) ) - F(\gamma(a) ) = \int_\gamma F'(\gamma(t) ) dt

On the other hand, since F(z)F(z) is differential,

lim⁡h→0F(z+h)−F(z)h=∂u∂x+i∂v∂x=lim⁡ih→0F(z+ih)−F(z)ih=−i∂u∂y+∂v∂y. \lim_{h \rightarrow 0} \frac{F(z + h) - F(z) }{h} = \frac{\partial u}{\partial x} + i \frac{\partial v}{\partial x} = \lim_{ih \rightarrow 0} \frac{F(z + ih) - F(z) }{ih} = -i \frac{\partial u}{\partial y} + \frac{\partial v}{\partial y}.

It holds that ∂u∂x=∂v∂y\frac{\partial u}{\partial x} = \frac{\partial v}{\partial y} and ∂v∂x=−∂u∂y\frac{\partial v}{\partial x} = -\frac{\partial u}{\partial y}. Therefore,

F′(γ(t))=∂u∂x∂x∂t−∂v∂x∂y∂t+i(∂v∂x∂x∂t+∂u∂x∂y∂t)=(∂u∂x+i∂v∂x)(∂x∂t+i∂y∂t)=F′(γ(t))γ′(t)=f(γ(t))γ′(t). \begin{aligned} F'(\gamma(t) ) &= \frac{\partial u}{\partial x} \frac{\partial x}{\partial t} - \frac{\partial v}{\partial x} \frac{\partial y}{\partial t} + i \left( \frac{\partial v}{\partial x} \frac{\partial x}{\partial t} + \frac{\partial u}{\partial x} \frac{\partial y}{\partial t} \right) \\ &= \left( \frac{\partial u}{\partial x} + i \frac{\partial v}{\partial x} \right) \left( \frac{\partial x}{\partial t} + i \frac{\partial y}{\partial t} \right) \\ &= F'(\gamma(t) ) \gamma'(t) \\ &= f(\gamma(t) ) \gamma'(t). \end{aligned}

and

F(γ(b))−F(γ(a))=∫γf(γ(t))γ′(t)dt F(\gamma(b) ) - F(\gamma(a) ) = \int_\gamma f(\gamma(t) ) \gamma'(t) dt

Remark: That ux′,uy′u_x', u_y' are continuous is a strong condition. What we need in the proof is just that *

u(x(tj),y(tj))=u(x(tj−1),y(tj−1))+ux′(x(tj−1),y(tj−1))(x(tj)−x(tj−1))+uy′(x(tj−1),y(tj−1))(y(tj)−y(tj−1))+o(x(tj)−x(tj−1),y(tj)−y(tj−1)) \begin{aligned} u(x(t_j), y(t_j) ) &= u(x(t_{j - 1} ), y(t_{j - 1} ) ) \\ &+ u_x'(x(t_{j - 1} ), y(t_{j - 1} ) ) (x(t_j) - x(t_{j - 1} ) ) \\ &+ u_y'(x(t_{j - 1} ), y(t_{j - 1} ) ) (y(t_j) - y(t_{j - 1} ) ) \\ &+ o(x(t_j) - x(t_{j - 1} ), y(t_j) - y(t_{j - 1} ) ) \end{aligned}

or

u(x(tj−1)+δx,y(tj−1)+δy)=u(x(tj−1),y(tj−1))+ux′(x(tj−1),y(tj−1))δx+uy′(x(tj−1),y(tj−1))δy+o(δx,δy) \begin{aligned} u(x(t_{j - 1} ) + \delta_x, y(t_{j - 1} ) + \delta_y) &= u(x(t_{j - 1} ), y(t_{j - 1} ) ) \\ &+ u_x'(x(t_{j - 1} ), y(t_{j - 1} ) ) \delta_x \\ &+ u_y'(x(t_{j - 1} ), y(t_{j - 1} ) ) \delta_y \\ &+ o(\delta_x, \delta_y) \end{aligned}

i.e., uu is differentiable at every point in the domain.

Question: Is is possible that uu is differentiable but its partial derivatives are not continuous? The answer is NO. First we consider an one dimension counter example:

f(x)={x2sin⁡1x∀x≠0,0,ifx=0 f(x) = \begin{cases} x^2 \sin \frac{1}{x} \quad \forall x \neq 0, \\ 0, \quad if \quad x = 0 \end{cases}

To show f(x)f(x) is differential at 00, we use the definition of derivative:

f′(0)=lim⁡h→f(h)−f(0)h−0=lim⁡h→0hsin⁡1h=0 f'(0) = \lim_{h \rightarrow } \frac{f(h) - f(0) }{h - 0} = \lim_{h \rightarrow 0} h \sin \frac{1}{h} = 0

On the other hand, f(x)f(x) is always differentiable from 00:

f′(x)=2xsin⁡1x+x2(−1/x2)cos1x=2xsin⁡1x−cos1x f'(x) = 2x \sin \frac{1}{x} + x^2 ( -1 / x^2) cos \frac{1}{x} = 2x \sin \frac{1}{x} -cos \frac{1}{x}

Hence, f(x)f(x) is differentiable everywhere, but f′(x)f'(x) is not continuous at 00.

The above example can be easily extended to two dimension case:

f(x,y)={(x2+y2)sin⁡1x2+y2∀(x,y)≠(0,0)0if(x,y)=(0,0) f(x, y) = \begin{cases} (x^2 + y^2) \sin \frac{1}{ \sqrt{x^2 + y^2} } \quad \forall (x, y) \neq (0, 0) \\ 0 \qquad \qquad \qquad \qquad \quad if \quad (x, y) = (0, 0) \end{cases}

By definition, f(x,y)f(x, y) is differential at (0,0)(0, 0),

lim⁡∥(x,y)∥→0∥f(x,y)−f(0)−(0,0)⋅(x,y)∥∥(x,y)∥≤lim⁡∥(x,y)∥→0∥x2+y2∥∥(x,y)∥=lim⁡∥(x,y)∥→0∥(x,y)∥=0 \lim_{\Vert (x, y)\Vert \rightarrow 0} \frac{\Vert f(x, y) - f(0) - (0,0) \cdot (x, y)\Vert }{\Vert (x, y)\Vert } \le \lim_{\Vert (x, y)\Vert \rightarrow 0} \frac{\Vert x^2 + y^2\Vert }{\Vert (x, y)\Vert } = \lim_{\Vert (x, y)\Vert \rightarrow 0} \Vert (x, y)\Vert = 0

However, it partial derivative from (0,0)(0, 0)

∂f∂x=2xsin⁡1x2+y2+(x2+y2)(−2xx2+y23)cos⁡1x+y=2xsin⁡1x2+y2−2xx2+y2cos⁡1x2+y2 \begin{aligned} \frac{\partial f}{\partial x} &= 2x \sin \frac{1}{ \sqrt{x^2 + y^2} } + (x^2 + y^2) \left( -\frac{ 2x }{ \sqrt{ x^2 + y^2 }^{3} } \right) \cos \frac{1}{ \sqrt{x + y} } \\ &= 2x \sin \frac{1}{ \sqrt{x^2 + y^2} } - \frac{ 2x }{ \sqrt{ x^2 + y^2 } } \cos \frac{1}{ \sqrt{x^2 + y^2 } } \end{aligned}

and

∂f∂x(0,0)=lim⁡h→h2sin⁡1∣h∣h=0 \frac{\partial f}{\partial x}(0, 0) = \lim_{h \rightarrow} \frac{h^2 \sin \frac{1}{|h|} }{h} = 0

Therefore, ∂f/∂x\partial f / \partial x is not continuous.

Line Integral

How do we extend integral in real line to C\mathbb{C}. The most simplest extension is of a real to complex function

γ(t)=x(t)+iy(t), \gamma(t) = x(t) + i y(t),

whose integral on an interval is defined by integrating its real part and imaginary parts separately.

Definition. Let γ(t)=x(t)+iy(t)\gamma(t) = x(t) + i y(t) be a function defined on [a,b][a, b]. If both x(x)x(x) and y(t)y(t) are Riemann integrable, then the integral of γ\gamma over [a,b][a, b] is defined as

∫abγ(t) dt≐∫abx(t)dt+i∫aby(t)dt, \int_a^b \gamma(t) \, dt \doteq \int_a^b x(t) dt + i\int_a^b y(t) dt,

Reference

[1] 常庚哲 and 史济怀, 数学分析教程 (第三版).

[2] E. M. Stein and R. Shakarchi, Complex analysis. Princeton University Press, 2003.