From example 3.5.4, page 212 nonlinear pde’s by Lokenath Debnath, 3rd edition.
First order PDE of three unknowns. Solve for \(u(x,y,z)\)
Mathematica ✓
ClearAll["Global`*"]; pde = (y - z)*D[u[x, y, z], x] + (z - x)*D[u[x, y, z], y] + (x - y)*D[u[x, y, z], z] == 0; sol = AbsoluteTiming[TimeConstrained[DSolve[pde, u[x, y, z], {x, y, z}], 60*10]];
Maple ✓
restart; pde :=(y-z)*diff(u(x,y,z),x)+(z-x)*diff(u(x,y,z),y)+(x-y)*diff(u(x,y,z),z)=0; cpu_time := timelimit(60*10,CodeTools[Usage](assign('sol',pdsolve(pde,u(x,y,z),'build')),output='realtime'));
Hand solution
Solve
The following method works only when the sum of the coefficients of \(u_{x},u_{y},u_{z}\) is zero. This is the case here. Hence we write
We need one more equation which is
And since we know that \(C_{1}=F\left ( C_{2},C_{3}\right ) \) where \(F\) is arbitrary function, then the solution is
Again, the above method only worked because of the special value of the coefficients. If the PDE was \(\left ( 2y-z\right ) u_{x}+\left ( z-x\right ) u_{y}+\left ( x-y\right ) u_{z}=0\) for example, then the above method will not work and we have to use the Lagrange-Charpit equations. But since the equations will be coupled, this would make the solution harder.
___________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________