In class we sketched the essentials of the method, which is like the method of characteristics for quasilinear equations, except with two more equations. These two new equations describe how p(=ux) and q(=uy) evolve along the characteristics.
JOHN has a particularly clear discussion of why this method works, and does indeed generate a smooth surface.
I will not ask about this fully nonlinear case on the tests or exam.