I want to solve the partial differential equation $uu_{xy} = u_xu_y$. It is known that a solution is $u(x,y) = f(x)g(y)$ for all pairs of (differentiable) functions $f$ and $g$ of one variable (Strauss Section 1.1 Problem 11).
This is what Mathematica 12 gives, which I'm unsure is correct since DSolve
, etc. are seemingly buggy:
DSolveValue[{u[x, y] D[u[x, y], x, y] == D[u[x, y], x] D[u[x, y], y]}, u[x, y], {x, y}]
NDSolveValue[{u[x, y] D[u[x, y], x, y] == D[u[x, y], x] D[u[x, y], y]}, u[x, y], {x, y}]
Simple PDE
This is not simple PDE. It is non-linear It requires finding first integral, then use that to convert the pde to set of ODE's. It is not always easy to find first integral of a PDE, in particular, a non-linear one. Of course the question is asking to verify, and not to solve it. Big difference,. it is very easy to verify. Just plugin the assumed solution into the PDE. Mathematica can do this very easily. btw, what you wrote in Latex and what you wrote in the Mathematica code is not the same thing. $\endgroup$