\documentclass[a4paper,12pt]{article}
\newcommand{\ds}{\displaystyle}
\newcommand{\un}{\underline}
\newcommand{\undb}{\underbrace}
\newcommand{\pl}{\partial}
\parindent=0pt
\begin{document}

{\bf Question}

Consider the equation

$$u_{xx}+8u_{xy}+7u_{yy}=0$$

\begin{description}
\item[(i)]
Classify and calculate the characteristics of this equation.

\item[(ii)]
Let the boundary data

$$u(\tau,\tau)=f(\tau),\ u_x(\tau,\tau)=g(\tau)$$

be given on the curve $C:(x,y)=(X(\tau),Y(\tau))$ where
$X(\tau)=Y(\tau)=\tau$ where $f$ and $g$ are sufficiently
differentiable functions.

Write down the general form of the Taylor expansion of the
solution $u(x,y)$ about a general given point $(x_0,y_0)$ to
second order in the derivatives wrt $x$ and $y$.

By differentiating the boundary data and using the equation
itself, show that it is not possible to calculate the second
derivatives anywhere, and hence it is not possible to solve the
equation with this boundary data. Comment on this in relation to
the answers of part (i).

\item[(iii)]
If now the boundary data is

$$u(X(\tau),Y(\tau))=f(\tau),\ u_y(X(\tau),Y(\tau))=g(\tau)$$

and is given on the curve $C:(x,y)=(X(\tau),Y(\tau))=(\tau,0)$,
calculate the expansion of $u(x,y)$ about $(0,0)$ up to second
order.
\end{description}

\medskip

{\bf Answer}

\begin{description}
\item[(i)]
Second order, linear, homogeneous, constant coefficients, PDE

$a=1,\ b=4,\ c=7$

$\Rightarrow b^2-ac=16-7=9>0 \Rightarrow$ hyperbolic everywhere

Characteristics given by

$$\ds\frac{dy}{dx}=\ds\frac{4 \pm \sqrt{16-7}}{1}=\ds\frac{4 \pm
3}{1}=7\ \rm{or}\ 1$$

$$\Rightarrow un{y=7x+c\ \rm{or}\ y=x+d}$$

where $c,\ d$ are arbitrary constants

(Note that general solution from bookwork is

$u(x,y)=P(y-7x)+Q(y-x)$ for $P,\ Q$ arbitrary functions)

\item[(ii)]
\begin{eqnarray*} u(x,y) & = & \left.u(x_0,y_0)+(x-x_0)\ds\frac{\pl u}{\pl
x}\right|_{(x_0,y_0)}+\left.(y-y_0)\ds\frac{\pl u}{\pl y}
\right|_{(x_0,y_0)}\\ & + &
\ds\frac{1}{2}\left.(x-x_0)^2\ds\frac{\pl^2 u}{\pl
x^2}\right|_{(x_0,y_0)}+\left.(x-x_0)(y-y_0)\ds\frac{\pl^2 u}{\pl
x\pl y}\right|_{(x_0,y_0)}\\ & + &
\ds\frac{1}{2}\left.(y-y_0)^2\ds\frac{\pl^2 u}{\pl
y^2}\right|_{(x_0,y_0)}+\cdots \end{eqnarray*}

Taylor expansion to 2nd order:

From lectures we use the following

$$u(x,y)=u(X(\tau)=Y(\tau))\ \rm{on}\ C$$

But when $X(\tau)=Y(\tau)=\tau$

\begin{eqnarray*} u(X,Y) & = & f(\tau)\ \ \ (1)\\ u_x(X,Y) & = &
g(\tau)\ \ \ (2)
\end{eqnarray*}

Differentiate both wrt $\tau$:

$\begin{array} {rlclcl} \ds\frac{d(1)}{d\tau}:= & u_x\dot{X} & + &
u_y\dot{Y} & = & \dot{f} \ \ \ (3)\\ \ds\frac{d(2)}{d\tau}:= &
u_{xx}\dot{X} & + & u_{xy}\dot{Y} & = & \dot{g}\ \ \ (4)
\end{array}$

($\dot{} \equiv \ds\frac{d}{d\tau}$)

Therefore on $C$ we have $\dot{X}=1,\ \dot{Y}=1$ and so

$(3): \un{\left.u_x\right|_C+\left.u_y\right|_C=\dot{f}}$

$(4): \un{\left.u_{xx}\right|_C+\left.u_{xy}\right|_C=\dot{g}}$

From $(1)$ and $(2)$ we know that $\un{\left.u\right|_C=f}$ and
$\un{\left.u_x\right|_C=g}$

Thus from $(3)$ we have

$$\un{\left.u_y\right|_C=\dot{f}-g}$$

Now seek $\left.u_{xx}\right|_C,\ \left.u_{xy}\right|_C,\
\left.u_{yy}\right|_C$:

Need 3 equations.

PDE itself is one equation, $(4)$ is another

The third comes from

$\ds\frac{d^2(1)}{d\tau^2}=\ds\frac{d(3)}{d\tau}:=$

$u_{xx}\dot{X}^2+u_{xy}\dot{X}\dot{Y}+u_x\ddot{X}
+u_{yx}\dot{Y}\dot{X}+u_{yy}\dot{Y}^2+u_y\ddot{Y}=\ddot{f}$

or $\un{\left.u_x\right|_C+2\left.u_{xy}\right|_C
+\left.u_{yy}\right|_C=\ddot{f}}$\ \ \ $(5)$

Thus we have (dropping $|_C$ but remembering it's still there) on
$C$:

$$\left\{\begin{array} {rcrcrcl} u_{xx} & + & 2u_{xy} & + & u_{yy}
& = & \ddot{f}\\ u_{xx} & + & u_{xy} & & & = & \dot{g}\\ u_{xx} &
+ & 8u_{xy} & + & 7u_{yy} & = & 0 \end{array} \right.$$

or

$$\left(\begin{array}{ccc} 1 & 2 & 1\\ 1 & 1 & 0\\ 1 & 8 &
7\end{array}\right)\left(\begin{array}{c}u_{xx}\\u_{xy}\\u_{yy}\end{array}\right)
=\left(\begin{array}{c}\ddot{f}\\ \ddot{g}\\ 0
\end{array}\right)$$

To solve this we must check first that

$\bigtriangleup=\left|\begin{array}{ccc} 1 & 2 & 1\\ 1 & 1 & 0\\ 1
& 8 & 7\end{array}\right| \ne0$. However

$\bigtriangleup=1\left|\begin{array}{cc} 1&0\\ 8 & 7
\end{array}\right|-2\left|\begin{array}{cc} 1&0\\ 1 & 7
\end{array}\right|+1\left|\begin{array}{cc} 1&1\\ 1 & 8
\end{array}\right|=1\times 7-2 \times 7 +1\times(8-1)=\un{0}$

Therefore no solution \un{exists} on $C$.

Hence cannot find Taylor expansion about $C$. Hence cannot find
any solution with this boundary data. The reason is that
$(x,y)=(\tau,\tau)$ is a characteristic curve, cf (i), where
$y=x+d$ are the characteristics. Consequently $C$ is more than
just tangential to a characteristic!

\newpage
\item[(iii)]
Now given data on $C: (\tau,0)$

$\begin{array}{rclrcl} X(\tau) & = & \tau & Y(\tau) & = & 0\\
\dot{X} & = & 1 & \dot{Y} & = & 0\\ \ddot{X} & = & 0 & \ddot{Y} &
= & 0 \end{array}\ \ \ (\star)$

Given boundary data:

\begin{eqnarray*} u(X,Y) & = & f(\tau)\ \ \ (6)\\ u_y(X,Y) & = &
g(\tau)\ \ \ (7) \end{eqnarray*}

$\begin{array} {rlclcl} \ds\frac{d(6)}{d\tau}:= & u_x\dot{X} & + &
u_y\dot{Y} & = & \dot{f} \Rightarrow
\un{\left.u_x\right|_C=\dot{f}}\ \ \ (8)\\ \ds\frac{d(7)}{d\tau}:=
& u_{yx}\dot{X} & + & u_{yy}\dot{Y} & = & \dot{g} \Rightarrow
\un{\left.u_{xy}\right|_C=\dot{g}}\ \ \ (9) \end{array}$

Thus \un{$\left.u\right|_C=f,\ \left.u_x\right|_C=\dot{f},\
\left.u_y\right|_C=g$}

Need to find $\left.u_{xx}\right|_C,\ \left.u_{xy}\right|_C,\
\left.u_{yy}\right|_C$:

Have 2 equations from PDE itself and $(9)$. Get 3rd from

$\ds\frac{d(8)}{d\tau}:=$

$u_{xx}\dot{X}^2+u_{xy}\dot{X}\dot{Y}+u_x\ddot{X}
+u_{yx}\dot{Y}\dot{X}+u_{yy}\dot{Y}^2+u_y\ddot{Y}=\ddot{f}$

$\Rightarrow \un{u_{xx}=\ddot{f}}$

Thus we have

$\left(\begin{array}{ccc} 1 & 0 & 0\\ 0 & 1 & 0\\ 1 & 8 & 7
\end{array}\right)\left(\begin{array}{c}u_{xx}\\u_{yx}\\u_{yy}\end{array}\right)
=\left(\begin{array}{c}\ddot{f}\\ \dot{g}\\ 0 \end{array}\right)$

$\Rightarrow \left\{\begin{array} {rcl} u_{xx} & = & \ddot{f}\\
u_{xy} & = & \dot{g}\\ u_{yy} & = &
-\ds\frac{(\ddot{f}+8\dot{g}}{7} \end{array} \right.$

Thus expansion of solution about $(0,0)$ is given by expansion in
part (ii) for the appropriate value of $\tau$.

What is it? $(x_0,y_0)=(0,0) \equiv \un{\tau =0}$ from $(\star)$

Thus

\begin{eqnarray*} u(x,y) & = & f(0)+x\dot{f}(0)+yg(0)\\ & &
+\ds\frac{1}{2}x^2\ddot{f}(0)+xy\dot{g}(0)-\left(\ds\frac{\ddot{f}(0)
+8\ddot{g}(0)}{14}\right)y^2+\cdots \end{eqnarray*}
\end{description}

\end{document}
