Consider the
Problem. Let $\gamma$ be a parameterized curve in $\mathbb{R}^{2}$ by $\gamma: I \to \Omega$, where $I$ is an interval of $\mathbb{R}$ and $\Omega$ an open in $\mathbb{R}^{2}$. Let $a,b,c: \Omega \to \mathbb{R}$ be given functions. Determine a function $\varphi(x,y)$ solution of the equation $$a(x,y)\frac{\partial \varphi}{\partial x} + b(x,y)\frac{\partial \varphi}{\partial y} = c(x,y) \tag{Eq 1}$$ where $\varphi(\gamma(t)) = \varphi_{0}(t)$ with $\varphi_{0}: I \to \mathbb{R}$ is a given function.
Proof. (just a loosely idea) Fixed a point $\gamma_{0} = \gamma_{0}(s_{0}) = \gamma_{0}(x_{0},y_{0})$ of $\gamma$, consider the curve $\Gamma(t) = (x(t),y(t))$ passing through $\gamma_{0}$, that is, $\Gamma(0) = \gamma_{0}$. Define $z(t) = \varphi(x(t),y(t))$ where $\varphi$ is a solution of Eq 1. If $\Gamma$ is differentiable, by the Chain Rule, $$\frac{dz}{dt} = \langle \Gamma'(t), \nabla\varphi(\Gamma(t)) \rangle = \frac{dx}{dt}\frac{d\varphi}{dx} + \frac{dy}{dt}\frac{d\varphi}{dy}.$$ Therefore, if $\Gamma$ satisfies the system of ODE $$\begin{cases} \frac{dx}{dt} = a(x,y),&x(0) = x_{0}\\ \frac{dy}{dt} = b(x,y),&y(0) = y_{0}, \end{cases}\tag{Sy 1}$$ we can the solution $\varphi$ solving $$\frac{dz}{dt} = c(x,y),\quad z(0)=\varphi(s_{0}).$$ If we repeat the previous argument for all points $\gamma(s)$, $s \in I$, we obtain a family of curves on which the solution $\varphi$ can be determined.
The solutions of the Sy 1 define a change of variables, that is, a function $$f: \mathbb{R}^{2} \to \mathbb{R}^{2}$$ $$(t,s) \mapsto (x,y)$$ and the Inverse Function Theorem ensures a solution if $(a,b)$ is transversal to curve $\gamma$ where $(a,b)$ represents $(x,y) \mapsto (a(x,y),b(x,y))$.
My objective is to apply the Method of Characteristics (the Problem) to solve the one-dimensional linear Wave Equation that is given by
$$\frac{\partial^{2} u}{\partial^{2} x} - c_{0}^{2}\frac{\partial^{2} u}{\partial^{2} y} = 0$$
and we write $$\frac{\partial^{2} u}{\partial^{2} x} - c_{0}^{2}\frac{\partial^{2} u}{\partial^{2} y} = \left(\frac{\partial }{\partial x} + c_{0}\frac{\partial }{\partial y}\right)\left(\frac{\partial }{\partial x} - c_{0}\frac{\partial }{\partial y}\right)u = 0.$$ Let $$\left(\frac{\partial }{\partial x} - c_{0}\frac{\partial }{\partial y}\right)u = v(x,y) = v.$$ Then it's enough to solve
$$\underbrace{\frac{\partial u}{\partial x} - c_{0}\frac{\partial u}{\partial y} = v}_{(1)}\quad\text{and}\quad\underbrace{\frac{\partial v}{\partial x} + c_{0}\frac{\partial v}{\partial y} = 0}_{(2)}$$
I broke the Wave Equation in (1) and (2) because the Problem is only for order $1$, since the partial derivatives are of order $1$. Then I think that solving the two equations, I'll find a solution for Wave Equation. This is the way that I found to use the Problem.
Now, I want to use the Problem for solve (2) and so, use the solution for (2) to solve (1).
I couldn't apply the above problem to get a solution to the equations (1) and (2). Hans Lundmark gave me a good reference, but in it, the author does the construction with some different details and my professor wants me to use the problem, exactly as it is, to get the solution. I wish someone could help me with this.