Skip to main content
\(\newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\) \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\) \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\) \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\) \( \newcommand{\Span}{\mathrm{span}}\)
Mathematics LibreTexts

2.3: Nonlinear Equations in Two Variables

  • Page ID
  • \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \) \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)\(\newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\) \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\) \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\) \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\) \( \newcommand{\Span}{\mathrm{span}}\) \(\newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\) \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\) \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\) \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\) \( \newcommand{\Span}{\mathrm{span}}\)

    Here we consider equation


    where \(z=u(x,y)\), \(p=u_x(x,y)\), \(q=u_y(x,y)\) and \(F\in C^2\) is given such that \(F_p^2+F_q^2\not=0\).

    In contrast to the quasilinear case, this general nonlinear equation is more complicated.

    Together with (\ref{nonlinear1}) we will consider the following system of ordinary equations which follow from considerations below as necessary conditions, in particular from the assumption that there is a solution of (\ref{nonlinear1}).


    Definition. Equations (\ref{char1})--(\ref{char5}) are said to be characteristic equations of equation (\ref{nonlinear1}) and a solution


    of the characteristic equations is called characteristic strip or Monge curve.

    Gaspard Monge (Panthèon, Paris)

    Figure 2.3.1: Gaspard Monge (Panthèon, Paris)

    We will see, as in the quasilinear case, that the strips defined by the characteristic equations build the solution surface of the Cauchy initial value problem.

    Let \(z=u(x,y)\) be a solution of the general nonlinear differential equation (\ref{nonlinear1}).

    Let \((x_0,y_0,z_0)\) be fixed, then equation (\ref{nonlinear1}) defines a set of planes given by
    \((x_0,y_0,z_0,p,q)\), i. e., planes given by \(z=v(x,y)\) which contain the point \((x_0,y_0,z_0)\)
    and for which \(v_x=p\), \(v_y=q\) at \((x_0,y_0)\). In the case of quasilinear equations these set of planes is a bundle of planes which all contain a fixed straight line, see Section 2.1. In the general case of this section the situation is more complicated.

    Consider the example


    where \(f\) is a given positive function. Let \(E\) be a plane defined by \(z=v(x,y)\) and which contains \((x_0,y_0,z_0)\). Then the normal on the plane \(E\) directed downward is

    $${\bf N}=\frac{1}{\sqrt{1+|\nabla v|^2}}(p,q,-1),$$

    where \(p=v_x(x_0,y_0)\), \(q=v_y(x_0,y_0)\). It follows from (\ref{monge1}) that the normal \({\bf N}\) makes a constant angle with the \(z\)-axis, and the \(z\)-coordinate of \({\bf N}\) is constant, see Figure 2.3,2.

    Monge cone in an example

    Figure 2.3.2: Monge cone in an example

    Thus the endpoints of the normals fixed at \((x_0,y_0,z_0)\) define a circle parallel to the \((x,y)\)-plane, i. e., there is a cone which is the envelope of all these planes.

    We assume that the general equation (\ref{nonlinear1}) defines such a Monge cone at each point in \(\mathbb{R}^3\). Then we seek a surface \(S\) which touches at each point its Monge cone, see Figure 2.3.3.

    Monge cones

    Figure 2.3.3: Monge cones

    More precisely, we assume there exists, as in the above example, a one parameter \(C^1\)-family

    $$p(\lambda)=p(\lambda;x,y,z),\ q(\lambda)=q(\lambda;x,y,z)$$

    of solutions of (\ref{nonlinear1}). These \((p(\lambda),q(\lambda))\) define a family \(\Pi(\lambda)\) of planes.


    $${\bf x}(\tau)=(x(\tau),y(\tau),z(\tau))$$

    be a curve on the surface \(S\) which touches at each point its Monge cone, see Figure 2.3.4. Thus we assume that at each point of the surface \(\mathcal S\) the associated tangent plane coincides with a plane from the family \(\Pi(\lambda)\) at this point.

    Figure 2.3.4: Monge cones along a curve on the surface

    Consider the tangential plane \(T_{{\bf x}_0}\) of the surface \(S\) at \({\bf x}_0=(x(\tau_0),y(\tau_0),z(\tau_0))\). The straight line

    $${\bf l}(\sigma)={\bf x}_0+\sigma {\bf x}'(\tau_0),\qquad -\infty<\sigma<\infty,$$

    is an apothem (in German: Mantellinie) of the cone by assumption and is contained in the tangential plane \(T_{{\bf x}_0}\) as the tangent of a curve on the surface \(S\). It is defined through

    {\bf x}'(\tau_0)={\bf l}'(\sigma).

    The straight line \({\bf l}(\sigma)\) satisfies


    since it is contained in the tangential plane \(T_{{\bf x}_0}\) defined by the slope \((p,q)\). It follows


    Together with (\ref{xl}) we obtain


    The above straight line \({\bf l}\) is the limit of the intersection line of two neighboring planes which envelopes the Monge cone:


    On the intersection one has


    Let \(h\to 0\), it follows


    Since \(x=l_1(\sigma)\), \(y=l_2(\sigma)\) in this limit position, we have


    and it follows from (\ref{xl}) that


    From the differential equation \(F(x_0,y_0,z_0,p(\lambda),q(\lambda))=0\) we see that

    Assume \(x'(\tau_0)\not=0\) and \(F_p\not=0\), then we obtain from (\ref{pxqy}), (\ref{Fpq})
    and from (\ref{zpq}) (\ref{Fpq}) that
    It follows, since \(\tau_0\) was an arbitrary fixed parameter,
    {\bf x}'(\tau)&=&(x'(\tau),y'(\tau),z'(\tau))\\
    i. e., the tangential vector \({\bf x}'(\tau)\) is proportional to \((F_p,F_q,pF_p+qF_q)\).
    where \(F=F(x(\tau),y(\tau),z(\tau),p(\lambda(\tau)),q(\lambda(\tau)))\). Introducing the new parameter \(t\) by the inverse of \(\tau=\tau(t)\), where
    t(\tau)=\int_{\tau_0}^\tau\ a(s)\ ds,
    we obtain the characteristic equations (\ref{char1})--(\ref{char3}).
    Here we denote \({\bf x}(\tau(t))\) by \({\bf x}(t)\) again.
    From the differential equation (\ref{nonlinear1}) and from (\ref{char1})--(\ref{char3}) we get equations (\ref{char4}) and (\ref{char5}). Assume the surface \(z=u(x,y)\) under consideration is in \(C^2\), then
    F_x+F_zp+F_pp_x+F_qp_y&=&0,\ \ (q_x=p_y)\\
    since \(p=p(x,y)=p(x(t),y(t))\) on the curve \({\bf x}(t)\). Thus equation (\ref{char4}) of the characteristic system is shown. Differentiating the differential equation (\ref{nonlinear1}) with respect to \(y\), we get finally equation (\ref{char5}).

    Remark. In the previous quasilinear case
    the first three characteristic equations are the same:
    x'(t)=a_1(x,y,z),\ y'(t)=a_2(x,y,z),\ z'(t)=a_3(x,y,z).
    The point is that the right hand sides are independent on \(p\) or \(q\). It follows from Theorem 2.1 that there exists a solution of the Cauchy initial value problem provided the initial data are non-characteristic. That is, we do not need the other remaining two characteristic equations.

    The other two equations (\ref{char4}) and (\ref{char5}) are satisfied in this quasilinear case automatically if there is a solution of the equation, see the above derivation of these equations.

    The geometric meaning of the first three characteristic differential equations (\ref{char1})--(\ref{char5}) is the following one. Each point of the curve
    {\mathcal A}:\quad (x(t),y(t),z(t))
    corresponds a tangential plane with the normal direction \((-p,-q,1)\) such that
    This equation is called strip condition. On the other hand, let \(z=u(x,y)\) defines a surface, then \(z(t):=u(x(t),y(t))\) satisfies the strip condition, where \(p=u_x\) and \(q=u_y\), that is, the ''scales'' defined by the normals fit together.

    Proposition 2.3. \(F(x,y,z,p,q)\) is an integral, i. e., it is constant along each characteristic curve.



    Corollary. Assume \(F(x_0,y_0,z_0,p_0,q_0)=0\), then \(F=0\) along characteristic curves with the initial data \((x_0,y_0,z_0,p_0,q_0)\).

    Proposition 2.4.
    Let \(z=u(x,y)\), \(u\in C^2\), be a solution of the nonlinear equation (\ref{nonlinear1}). Set
    z_0=u(x_0,y_0,)\ p_0=u_x(x_0,y_0),\ q_0=u_y(x_0,y_0).
    Then the associated characteristic strip is in the surface \(\mathcal{S}\), defined by \(z=u(x,y)\).
    where \((x(t),y(t),z(t),p(t),q(t))\) is the solution of the characteristic system
    (\ref{char1})--(\ref{char5}) with initial data \((x_0,y_0,z_0,p_0,q_0)\)

    Proof. Consider the initial value problem
    with the initial data \(x(0)=x_0\), \(y(0)=y_0\). We will show that
    is a solution of the characteristic system. We recall that the solution exists and is uniquely determined.

    Set \(z(t)=u(x(t),y(t))\), then \((x(t),y(t),z(t))\subset\mathcal{S}\), and
    Set \(p(t)=u_x(x(t),y(t)),\ q(t)=u_y(x(t),y(t))\), then
    Finally, from the differential equation \(F(x,y,u(x,y),u_x(x,y),u_y(x,y))=0\) it follows